Protein profiling of plasma proteins in dairy cows with subclinical hypocalcaemia
Irish Veterinary Journal volume 70, Article number: 3 (2017)
Subclinical hypocalcaemia (SH) is an important metabolic disease in dairy cows that has a serious impact on production performance. The objective of this study was to investigate novel aspects of pathogenesis using proteomics technology to identify proteins that are differentially expressed in diseased and healthy animals. Dairy cows were divided into an SH group (T, n = 10) and a control group (C, n = 10) based on plasma calcium concentration. A total of 398 differentially expressed proteins were identified, of which 265 proteins were overlapped in the two parallel experiments. Of these, 24 differentially expressed proteins were statistically significant. Gene Ontology analysis yielded 74 annotations, including 7 cellular component, 55 biological process and 12 molecular function categories. Bioinformatics analysis indicated that calcium regulation, immune and inflammatory response, blood coagulation and complement pathway were all related to SH. Our iTRAQ/LC-MS/MS (isobaric tags for relative and absolute quantification/liquid chromatography-mass spectrometry/mass spectrometry) approach proved highly effective for plasma protein profiling of dairy cows with SH, and the results pave the way for further studies in this area.
Animal husbandry, and dairy farming in particular, has grown considerably in China in recent years, and is now a major industry. However, the increase in dairy farming has been accompanied by the proliferation of metabolic diseases that have a significant deleterious effect on the production performance of dairy cows. Subclinical hypocalcaemia (SH) is one such metabolic disease that has no obvious clinical signs and can therefore be easily ignored. Biochemical tests show that the blood calcium (Ca) concentration of dairy cows with SH is between 1.38 and 2.00 mmol/L. This disease increases the risk of other perinatal diseases due to the reduced blood Ca concentration and muscle contractive capacity . SH has been studied using proteomics approaches based on SELDI-TOF-MS (surface-enhanced laser desorption/ionization-time of flight-Mass spectrometry) [2, 3], 2D-DIGE (fluorescence two-dimensional differential gel electrophoresis) combined with MALDI-TOF-MS (matrix-assisted laser desorption/ionization-time of flight-mass spectrometry)  and metabolomics based on 1H-NMR (1hydrogen-nuclear magnetic resonance) . However, these techniques have their limitation on studying differentially expressed proteins in dairy cows with SH. And our iTRAQ/LC-MS/MS technique can obtain more comprehensive information on differentially expressed proteins, which can further explore the pathogenesis and preventive methods of SH in dairy cows. In the present study, we used iTRAQ (isobaric tags for relative and absolute quantitation) to isotopically label and simultaneously analyze protein expression in diseased and healthy animals . The objective of this study was to identify differentially expressed proteins in plasma from dairy cows with SH using iTRAQ/LC-MS/MS combined with bioinformatics analysis to provide novel insight into the pathogenesis of SH.
Materials and methods
This experiment was conducted in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health located in United States. All experimental animals were treated according to the International Guiding Principles for Biomedical Research Involving Animals.
All cows were obtained from an intensive dairy farm located in Heilongjiang Province. All diets were in accordance with the Chinese standards for cattle breeding. All cows were fed a total mixed ration during the dry period, which consisted of 2.02 kg of concentrated feed, 17.37 kg of silage maize, and 4.2 kg of leymus chinensis, with a dry matter content of 42.9%. The net energy for dry cow of this diet was 1.3 Mcal/kg DM. The nutritional content was 12.89% crude protein, 2.83% fat, 49.21% neutral detergent fibre, 27.47% acid detergent fibre, and 71.71 g of Ca and 42.75 g of P per cow per day. The major mineral content was 91.81 mg of zinc, 59.29 mg of manganese, 16.60 mg copper, 0.44 mg of iodine, 0.25 mg of cobalt and 0.79 mg of selenium per cow per day. The dietary cation anion difference (DCAD) was 91 meq/kg DM. All cows were fed the same diet during the dry period.
Experimental Holstein dairy cows of a similar age, parity and body condition were selected. In total, 10 animals with SH were placed in T group, and 10 healthy cows were placed in C group. Cows were considered to be suffering from SH if within 1–3 days postpartum they displayed no obvious clinical signs but exhibited a plasma Ca concentration between 1.38 and 2.00 mmol/L. In addition, both groups contained only clinically normal cows in every other respect. In order to minimize individual differences between samples, 5 plasma samples from each group were mixed to generate one hybrid sample, and this was labelled using an iTRAQ kit and used in subsequent experiments. Basic information on dairy cows and experimental design are listed in Table 1.
Collection of blood samples
All blood samples were obtained within 6 h of calving from tail veins of cows in the early morning, anticoagulated with heparin and immediately centrifuged at 3000 × g for 10 min at room temperature. The plasma were transferred to eppendorf tubes (1.5 mL) and stored at −80 °C until needed.
Measurement of plasma Ca concentration
The plasma Ca concentration was measured using an automatic biochemical analyser (modular DPP, Roche Diagnostics GmbH) and a Ca plasma kit (651564-01, Roche Diagnostics GmbH).
Elimination of high abundance proteins from plasma samples
In order to enrich low abundance proteins, high abundance proteins were eliminated according to the operating instructions supplied with the ProteoMiner kit (Bio-Rad, USA, catalog#s163-3007) following thawing and mixing of samples. Protein concentration was measured by the Bradford method to determine the loading quantity, and eliminated high abundance proteins were detected by SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis).
Plasma samples (100 μL) digested with protease were mixed with trypsase (dissolved in precooled deionized water) at a ratio of 20:1, digested again at 27 °C overnight, and centrifuged at 1000 × g for 30 s at 4 °C. iTRAQ reagent (in each tube) was dissolved completely in 70 μL of ethanol, and centrifuged at 1000 × g for 30 s at 4 °C. Treated plasma samples were labelled with treated iTRAQ reagent following standing for 2 h, and labelled plasma samples (1 μL) were vortexed and evaporated to dryness in vacuum. Information on iTRAQ labelling is shown in Table 1.
Strong cation exchange (SCX) and high pH reverse fractionation
Peptides in dried plasma samples were dissolved in loading buffer (10 mmol/L K3PO4 pH 3.0, 25% ACN) (acetonitrile) and the pH was adjusted to 3. A PolyLC SCX chromatographic column (10 mm id × 14 mm, 12 μm, 300 A; PolyLC Inc, Columbia, MD, USA) was equilibrated in loading buffer, and 400 μg of sample was loaded onto the column and eluted with 2 ml loading buffer. Eluted peptides were evaporated to dryness in vacuum and separated by high pH reverse fractionation. Following drying in vacuum, peptides were dissolved in 150 μl of 2% ACN/0.5% FA (formaldehyde) and stored at −80 °C until needed for MS.
Nano LC-MS/MS separation and MS identification
Samples were separated using an UltiMate 3000 RSLCnano nano-litre flow rate LC system (Thermo-Dionex, Sunnyvale, CA, USA), and MS was performed with a LTQ Orbitrap Velos (Thermo-Fisher Scientific, San Jose, CA, USA). Raw data were collected by Xcalibur 2.1 operation software (Thermo-Fisher Scientific), corrected using Proteome Discoverer 1.3 (Thermo) and exported as MGF formatted files. Mascot 2.3.02 (Matrix Science, Boston, MA, USA) was used for database searching, and processed data were used to search the NCBInr Bovine RefSeq database that contained 71,248 protein sequences when downloaded on April 11, 2011. The error rate of peptide and protein identification was <5%, and at least two unique peptides were identified for each protein. All potential SH-related proteins were subjected to quantitative analysis by Proteome Discoverer l.3 according to the pseudo-molecular- ion intensity obtained in the peptides report. The ratio of the quasi-molecular-ion intensity (114/116 and 115/117) was determined by each protein present in both T and C group.
The age, parity, body condition score  and plasma Ca concentration of cows between groups were analyzed by one-way ANOVA using SPSS version 19.0 (IBM, Armonk, NY, USA). Differentially expressed proteins were analyzed by Venn analysis (VENN 2.1, http://bioinfogp.cnb.csic.es/tools/venny/) and Gene Ontology (GO) analysis (GOTERM_BP_ALL,GOTERM_CC_ALL, GOTERM _MF_ALL) using DAVID bioinformatics analysis software (v6.7, http://david.abcc.ncifcrf.gov/home.jsp).
Information on experimental animals
The age, parity, body condition score and plasma Ca concentration of SH and control cows are listed in Table 1. The observed decrease in the plasma Ca concentration between T and C groups was the only significant difference.
Venn analysis of experimental results
Data (114/116 and 115/117) were collected in two parallel experiments, and 365 and 298 proteins were identified in each parallel experiment, which resulted in 398 differentially expressed proteins. Venn analysis of all 398 differentially expressed proteins confirmed 265 overlapping proteins, and 100 and 33 unique proteins in each parallel experiment (Fig. 1). Overlapping proteins accounted for 66.58% of all identified proteins.
Identification of differentially expressed proteins
Mean values of the fold change in expression were calculated for proteins overlapping (265) in the two parallel experiments. A mean value greater than 1.5 was considered up-regulated, and a mean value lower than 0.67 was considered down-regulated. As a result, 24 differentially expressed proteins were identified, comprising 21 that were up-regulated and 3 that were down-regulated (Tables 2 and 3).
Bioinformatics analysis of differentially expressed proteins
All 24 differentially expressed proteins were subjected to function enrichment analysis of cell component, molecular function and biological process categories (Fig. 2). In total, 74 annotations were obtained, comprising 7 cell components, 12 molecular functions and 55 biological processes.
GO analysis revealed that the most highly represented biological process terms were ‘inflammatory reaction, immunoreaction and resistance function’, while the top molecular function terms were ‘cell’ and ‘molecular’, and the dominant cell component terms were ‘extracellular domain’ and ‘lipoprotein’.
In this study, iTRAQ and LC-MS/MS were used for plasma protein profiling of SH in dairy cows. A total of 398 differentially expressed proteins were identified, of which 265 proteins overlapped in the two parallel experiments. Of these, 24 differentially expressed proteins were statistically significant, of which 21 were up-regulated and 3 were down-regulated. Based on bioinformatics analysis, most differentially expressed proteins were associated with inflammatory action and immunological processes.
Differentially expressed proteins related to blood Ca concentration
Two of the up-regulated proteins were identified as potentially related to blood Ca concentration, namely cadherin and periostin. Cadherin is a single chain transmembrane glycoprotein involved in cell adhesion, forming adherens junctions and binding cells together within tissues, and its activity is Ca dependent. Cadherin plays an important role in the formation of osteoblasts and bone by controlling adhesion and interaction between cells [6, 8, 9], and it has also been linked to cell recognition and sorting, boundary formation and maintenance, coordination of cell motility, and maintenance of cell structure . Meanwhile, the 811 amino acid periostin induces cell adhesion and separation, and members of the periostin family are characterized by a cysteine-rich N-terminal region . Studies have indicated that periostin is a novel member of the vitamin K-dependent gamma-carboxylate protein family . Periostin is involved in both normal physiological and pathological processes, including vascular disease , trauma repair , bone formation  and formation of tumours . Our GO analysis results indicated that cadherin and periosteal cells function in cell adhesion and protein binding.
In general, clinical hypocalcaemia (CH) is diagnosed based on clinical signs and blood Ca concentration. However, accurate diagnosis of SH in dairy cows, which has been linked to various diseases affecting productivity in transition dairy cows, depends on the measurement of blood calcium. Up-regulation of the Ca-associated cadherin and periostin in SH suggests these proteins may initiate a mechanism that stimulates blood Ca when levels are decreased. The exact mechanism through which these proteins operate in dairy cows with SH requires further study.
Differentially expressed proteins related to immune response and inflammatory reaction
SH is a well-characterised metabolic disease which has the high prevalence in the first week of lactation in dairy cows  and the high incidence of peripartum problems in dairy cows. Furthermore, immune cells in cows with SH display impaired function . Our bioinformatics analysis indicated that most proteins differentially expressed in diseased animals were related to immune response and inflammatory reaction. These included serum amyloid A protein and haptoglobin. Haptoglobin is an acidoglycoprotein that is highly expressed in blood and is believed to act as a cytokine following release from activated macrophages during the initial stages of pathogenic infection and tissue damage, which can stimulate the liver to produce acute phase reactive protein. Expression of haptoglobin in the blood is increased during stress or following pathogen-induced damage .
Serum amyloid A is mainly synthesized by the liver and can exist in acute phase amyloid A and structural amyloid protein forms. Expression of serum amyloid A can be up-regulated 100-fold within 1–2 days of an organism suffering from physical damage, stress or microbial infection . As mentioned above, dairy cows with SH show no obvious clinical signs. However, our results indicated that diseased animals do show hypoimmunity. Dairy cows with metabolic diseases such as SH are therefore more likely to contract other more debilitating infectious and non-infectious diseases that could greatly influence production performance.
Differentially expressed proteins related to blood coagulation and complement pathways
Previous studies on CH and SH in dairy cows reported that the occurrence of CH is related to blood coagulation and complement pathways [3, 21–23], and our results are consistent with this observation. A number of differentially expressed proteins identified in the present study were related to blood coagulation and complement pathways, specifically platelet factor, serpin and complement factor. Complement and coagulation pathway proteins such as serpins and complement factors are known to be highly dependent on the plasma Ca concentration [24, 25]. Serpins play an important role in coagulation, fibrinolysis and the protection of neurons, while platelet factors assist adhesion and aggregation of platelets, and accelerate blood coagulation. Many complement factors also stimulate the production of specific antibodies and therefore boost the immune system. Complement activation ultimately results in the formation of membrane attack complexes that attack and kill pathogenic microorganisms . Ca is an important factor in blood coagulation, and malfunctions can occur if the blood Ca concentration is reduced. In this study, we identified a number of differentially expressed proteins involved in coagulation and complement pathways. We speculate that reduced blood Ca in dairy cows with SH may alter key metabolic and immunological pathways, and lead to diminished production performance and hypoimmunity. These results therefore indicate that CH and SH could serve as recognizable risk factors for diseases affecting productivity in dairy cows such as mastitis, ketosis, retained placenta, displaced abomasum and uterine prolapse .
In this study, we identified proteins that were differentially expressed in dairy cows with SH through a combined proteomics and bioinformatics approach. A total of 398 differentially expressed proteins were obtained, of which 24 were statistically significant and confirmed by MS. These proteins were mainly associated with blood Ca concentration, coagulation and complement pathways. Our results provide novel insight into SH in dairy cows.
Dietary cation anion difference
Strong cation exchange
Vansaun RJ, Sniffen CJ. Nutritional management of the pregnant dairy cow to optimize health, lactation and reproductive performance. Anim Feed Sci Tech. 1996;59(1):13–26.
Wang PX, Shu S, Xia C, Wang Z, Wu L, Wang B, et al. Protein expression in dairy cows with and without subclinical hypocalcaemia. New Zeal Vet J. 2015;64(2):1–15.
Shu S, Xia C, Zhang HY, Sun ZL, Liu JN, Wang B. Plasma proteomics analysis of dairy cows with milk fever using SELDI-TOF-MS. Asian J Anim Vet Adv. 2014;9(1):1–12.
Xia C, Zhang HY, Wu L, Xu CC, Zheng JS, Yan Y, et al. Proteomic analysis of plasma from cows affected with milk fever using two-dimensional differential in-gel electrophoresis and mass spectrometry. Res Vet Sci. 2011;93(2):857–61.
Paweletz CP, Trock B, Pennanen M, Tsangaris T, Magnant C, Liotta LA, et al. Proteomic patterns of nipple aspirate fluids obtained by SELDI-TOF: potential for new biomarkers to aid in the diagnosis of breast cancer. Dis Markers. 2000;17(4):301–7.
Schambony A, Kunz M, Gradl D. Cross-regulation of Wnt signaling and cell adhesion. Differentiation. 2004;72(7):307–18.
Edmondson IJ, Lean IJ, Weaver LD, Farver T, Webster G. A body condition scoring chart for Holstein Dairy Cows. J Dairy Sci. 1989;72(1):68–78.
Bienz M. β-Catenin: a pivot between cell adhesion and Wnt signalling. Curr Biol. 2005;15(2):R64–7.
Brembeck FH, Rosário M, Birchmeier W. Balancing cell adhesion and Wnt signaling, the key role of β-catenin. Curr Opin Genet Dev. 2006;16(1):51–9.
Halbleib JM, Nelson WJ. Cadherins in development: cell adhesion, sorting, and tissue morphogenesis. Gene Dev. 2006;20(23):3199–214.
Liu GX, Xi HQ, Sun XY, Wei B. Role of periostin and its antagonist PNDA-3 in gastric cancer metastasis. World J Gastroenterol. 2015;21(9):2605–13.
Coutu DL, Wu JH, Monette A, Rivard G-É, Blostein MD, Galipeau J. Periostin, a member of a novel family of vitamin K-dependent proteins, is expressed by mesenchymal stromal cells. J Biol Chem. 2008;283(26):17991–8001.
Stanton LW, Garrard LJ, Damm D, Garrick BL, Lam A, Kapoun AM, et al. Altered patterns of gene expression in response to myocardial infarction. Circ Res. 2000;86(9):939–45.
Roy S, Patel D, Khanna S, Gordillo GM, Biswas S, Friedman A, et al. Transcriptome-wide analysis of blood vessels laser captured from human skin and chronic wound-edge tissue. Proc Natl Acad Sci U S A. 2007;104(36):14472–7.
Kruzynska FA, Wang JM, Maeda M, Rogers R, Krug E, Hoffman S, et al. Periostin is expressed within the developing teeth at the sites of epithelial–mesenchymal interaction. Dev Dynam. 2004;229(4):857–68.
Kudo Y, Siriwardena B, Hatano H, Ogawa I, Takata T. Periostin: novel diagnostic and therapeutic target for cancer. Histol Histopathol. 2007;22(10):1167–74.
Reinhardt TA, Lippolis JD, McCluskey BJ, Goff JP, Horst RL. Prevalence of subclinical hypocalcemia in dairy herds. Vet J. 2011;188:122–4.
Martinez N, Sinedino LDP, Bisinotto RS, Ribeiro ES, Gomes GC, Lima FS, et al. Effect of induced subclinical hypocalcemia on physiological responses and neutrophil function in dairy cows. J Dairy Sci. 2014;97(2):874–87.
Jung SM, Lee K, Lee JW, Namkoong H, Kim HK, Kim S, et al. Both plasma retinol-binding protein and haptoglobin precursor allele 1 in CSF: candidate biomarkers for the progression of normal to mild cognitive impairment to Alzheimer's disease. Neurosci Lett. 2008;436(2):153–7.
Bahk YY, Na BK, Cho S, Kim JY, Lim KJ, Kim TS. Proteomic analysis of haptoglobin and amyloid A protein levels in patients with vivax malaria. Korean J Parasitol. 2010;48(3):203–11.
Goff JP. The monitoring, prevention, and treatment of milk fever and subclinical hypocalcemia in dairy cows. Vet J. 2008;176(1):50–7.
Horst RL, Jorgensen NA. Elevated plasma cortisol during induced and spontaneous hypocalcemia in ruminants. J Dairy Sci. 1982;65(12):2332–7.
Ducusin RJT, Uzuka Y, Satoh E, Otani M, Nishimura M, Tanabe S, et al. Effects of extracellular Ca 2+ on phagocytosis and intracellular Ca 2+ concentrations in polymorphonuclear leukocytes of postpartum dairy cows. Res Vet Sci. 2003;75(1):27–32.
Beinrohr L, Murrayrust TA, Dyksterhuis L, Zavodszky P, Gal P, Pike RN, et al. Serpins and the complement system. Method Enzymol. 2011;499:55–75.
Amara U, Rittirsch D, Flierl M, Bruckner U, Klos A, Gebhard F, et al. Interaction between the coagulation and complement system. Adv Exp Med Biol. 2008;632(2):71–9.
Onat A, Can G, Rezvani R, Cianflone K. Complement C3 and cleavage products in cardiometabolic risk. Clin Chim Acta. 2011;412(13–14):1171–9.
DeGaris PJ, Lean IJ. Milk fever in dairy cows: a review of pathophysiology and control principles. Vet J. 2008;176(1):58–69.
The authors thank the participating farmers.
This study was supported by the National Science and Technology Foundation of China (Grant number 30972235), the Heilongjiang Province Nature Science and Technology Foundation of China (Grant number C200916), the Chinese National Science and Technology Project of China (Grant number 2013BAD21B01-2) and the Natural Science Foundation of Heilongjiang Province of China (Grant number C2016041).
Availability of data and supporting materials
ZLF carried out the final analysis and drafted the manuscript. CX assisted in the analysis and helped to draft the manuscript. SS, CCX, LYB, XHX and GW collected data at dairy farms and carried out initial analysis. HYZ, CX, LW and WY participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.
None of the authors of this paper has a financial or personal relationship with other people or organisations that could inappropriately influence or bias the content of the paper. The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the farm owner and all animal experiments were conducted according to the International Guiding Principles for Biomedical Research. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Heilongjiang Bayi Agricultural University.
Ziling Fan and Shi Shu are co-first authors.
About this article
Cite this article
Fan, Z., Shu, S., Xu, C. et al. Protein profiling of plasma proteins in dairy cows with subclinical hypocalcaemia. Ir Vet J 70, 3 (2017). https://doi.org/10.1186/s13620-017-0082-0