- Research article
- Open access
- Published:
Purification of protein C from canine plasma
BMC Veterinary Research volume 10, Article number: 251 (2014)
Abstract
Background
In order to characterize the functional properties of canine protein C (CnPC), the zymogen needs to be purified from plasma. The goals of this study were (1) to purify protein C from fresh frozen canine plasma by barium chloride and ammonium sulphate precipitation, followed by immunoaffinity chromatography using a monoclonal mouse antibody against human protein C (HPC4) and (2) to characterize this protein’s structure.
Results
The purified protein contained three glycosylated forms of a heavy chain (~49 kDa) and a glycosylated light chain (~25 kDa). Tandem mass spectra of the peptides obtained following trypsin digestion and liquid chromatography identified this protein to be protein C (vitamin K-dependent protein C precursor, gi|62078422) with 100% probability. Three glycosylation sites (Asn139, Asn202, and Asn350) were identified by detection of peptides containing an N-linked glycosylation consensus sequon with a 3-dalton increase in mass following incubation of the protein with PNGase F in 18O-labeled water. Following incubation with Protac (a specific activator of protein C), the heavy chain showed a slight decrease in molecular size and amidolytic activity measured by a synthetic chromogenic substrate containing an amide bond [H-D-(γ-carbobenzoxyl)-lysyl-prolyl-arginine-paranitroanilide diacetate salt]. The amidolytic activity was increased by ~303-fold in the final protein preparation compared to that in plasma. The purified protein showed concentration-dependent anti-factor V and anti-factor VIII activities in canine plasma in coagulometric factor assays.
Conclusions
These studies showed that CnPC could be purified from plasma using HPC4 and that this protein showed amidolytic and anti-coagulant properties upon activation with Protac.
Background
Activated protein C (APC) is a natural anticoagulant in plasma. It is derived from protein C zymogen, which is synthesized by hepatocytes and circulates in plasma at low concentration. Upon activation by thrombin, APC inhibits factor V (FV) and factor VIII (FVIII) activities by irreversible proteolysis [1],[2]. Recombinant human activated protein C (rhAPC) was previously available as drotrecogin alpha (activated) (Xigris, Eli Lilly & Company, Indianapolis, IN) and approved for treatment of sepsis in patients with high risk of death by the United States Food and Drug Administration (FDA) and European Medicines Agency (EMEA) based on results obtained from the Recombinant Human Activated Protein C Worldwide Evaluation in Severe Sepsis (PROWESS) study [3]. A subsequent trial, PROWESS-SHOCK, showed no difference in the 28-day or 90-day mortality rates between the Xigris and placebo groups [4]. Furthermore, a systematic review of six randomized clinical trials, including the two aforementioned studies, showed no efficacy of treatment with Xigris in sepsis [5]. The manufacturer of Xigris eventually announced the withdrawal of Xigris from the market worldwide on October 25, 2011 [6].
Despite the lack of clinical efficacy of Xigris for treating sepsis, the anti-inflammatory and cytoprotective effects of protein C are well established. It has also been shown that the anticoagulant effects of APC are not required for protecting septic mice against mortality [7]. An advantage of using non-coagulant APC variants, as opposed to the anticoagulant wild-type APC, in a clinical setting is reduced risk of bleeding, as the use of Xigris in human sepsis was associated with increased risk of serious bleeding [8]. Furthermore, some APC variants conferred a higher degree of cytoprotection than wild-type APC [9],[10].
The overall structure of the gene encoding for canine protein C (CnPC) zymogen is very similar to that for humans, consisting of nine exons and conserved exon/intron boundaries [11]. Computational analysis of CnPC cDNA predicted that CnPC had 456 amino acids and was 72.1% similar to human protein C [11]. Furthermore, in naturally occurring sepsis for both species, reduced circulating protein C activities were noted and early increase in protein C levels during hospitalization were of positive prognostic value [12]-[16]. The biological properties of CnPC are not as well established as those for human protein C. Given the therapeutic potential of rhAPC and human APC variants, the therapeutic potential of CnPC or canine activated protein C (CnAPC) warrants investigation. To facilitate such investigation, CnPC needs to be purified so that its precise structure and function may be determined. Knowledge of the CnPC pathway could potentially benefit canine patients and enhance our understanding of the protein C pathway in other species.
In this manuscript, we describe a method for purifying CnPC zymogen from fresh frozen plasma by salt precipitation followed by immunoaffinity chromatography using an antibody directed against a peptide sequence on human protein C (clone HPC4). This method was modified from a method described by Stenflo et al. for isolating protein C from fresh bovine plasma [17].
Results
A plasma protein was markedly enriched over the course of the purification process, which consisted of two major parts: salt precipitation (barium sulphate precipitation followed by ammonium sulphate precipitation) and immunoaffinity chromatography using a monoclonal antibody that targets the amino acid sequence EDQVDPRLIDGK of human protein C. On a silver-stained non-reducing sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE) of the different plasma protein fractions obtained at various stages of the purification process, this protein was not readily discernable until after immunoaffinity chromatography (Figure 1A). However, Western immunoblotting with a polyclonal antibody detected this protein in small quantities in the plasma protein fraction obtained after salt precipitation and in much greater (37.6-fold increase) amounts in the eluate obtained from immunoaffinity chromatography (Figure 1B-D). This protein showed strong immunoreactivity for a polyclonal antibody against human protein C, while no other plasma proteins in human protein C-deficient plasma or canine plasma showed such immunoreactivity (Figure 2). Furthermore, using domain-enhanced lookup time-accelerated BLAST (DELTA-BLAST) [18] to search the Swissprot database for canis lupus familiaris, only CnPC (accession number Q28278.2) was found to show significant alignment with the sequence EDQVDPRLIDGK.
The peptide mass spectra obtained following in-gel trypsin digestion and liquid chromatography followed by tandem mass spectral analysis (LC-MS/MS) identified protein C (vitamin K-dependent protein C precursor, gi|62078422) in all silver-stained protein bands with 100% probability. Table 1 shows a summary of the spectra assigned to identification of protein C in each protein band. Table 2 shows a summary of the mass spectral data for each excised band. Only protein C was identified in the three closely spaced bands seen at ~49 kDa (MW of the three bands were estimated to be: 49 kDa, 44 kDa, and 40 kDa) and a band at ~25 kDa (Figure 3). Protein C was also identified in other bands (~129 kDa, ~107 kDa, and ~75 kDa) (Figure 3). The 75-kDa band may represent the single chain form of protein C, which is also found in human plasma. However, Western immunoblotting using antibodies against human protein C did not detect immunoreactivity in the 75-kDa region in reduced CnPC (Figure 2B). No further studies were performed to assess whether or not single-chain form of protein C was in canine plasma. Apolipoprotein A-1 was present in a band at ~25 kDa. Complement C3 was present in two bands at ~129 kDa and ~75 kDa.
The structure of the purified protein showed many similarities with human protein C. The protein migrated as two tightly spaced bands of ~75 kDa under non-reducing conditions (Figures 1 and 4); and as three tightly spaced bands of ~49 kDa and a much smaller band of ~25 kDa (Figure 2). The results suggest that the protein consisted of two chains linked by at least one disulphide bond. To show the tightly spaced bands were different glycosylated forms, as opposed to different peptides or different degraded forms of the same peptide, the protein was reduced and then deglycosylated in a mixture of deglycosylation enzymes. The protein was first reduced to allow better access of the enzymes to their substrates. Following glycosylation, the three tightly spaced bands seen at ~49 kDa were no longer detected on SDS-PAGE (12%, silver-stained) (Figure 3). A smaller, single band was seen at ~27 kDa instead, suggesting that the three tightly spaced bands were indeed different glycosyated forms of the same peptide. The 25-kDa band that was seen with the parent form of the protein was also no longer detectable after deglycosylation. Instead, a new, smaller band was seen at ~23 kDa, suggesting that the parent form of this peptide was glycosylated as well (Figure 3). To determine the exact N-linked glycosylated sites of this protein, we searched the amino acid sequence of the protein for the consensus sequences (also known as sequons) N-X-T and N-X-S. The following four sites were identified: Asn139, Asn202, Asn289, and Asn350. It is recognized that these sequons are not the only requisite for N-linked glycosylation to occur [19],[20]. To confirm N-linked glycosylation at these four sites, isotope-coded glycosylation site-specific tagging (IGOT) [21] was performed. This method is based on the principle that while PNGase F removes a glycan from the Asn residue, it deamidates Asn, converting it into an Asp residue. During this process, an oxygen ion from a water molecule is incorporated into the new Asp residue, which makes it 1 dalton heavier than the predicted Asn residue. Given that natural deamidation of Asn residues is known to occur, a more definitive way to show that deamidation is due to the action of PNGase F is to incorporate a labeled oxygen ion, such as 18O-labelled water. Deamidation of Asn 18O-labelled water will result in incoporation of an oxygen ion that is 2 daltons heavier than the H216O, which is the more abundant form found in nature. This reaction will result in a 3-dalton mass shift, which is very unlikely to be a natural occurrence. By use of IGOT, three of the potential N-linked glycosylated sites were found to show a 3-dalton increase in mass: Asn139, Asn202, and Asn289 (Figures 5 and 6). The fourth N-linked glycosylated site (Asn350) was not captured on mass spectrometry on multiple attempts. Therefore, the glycosylation status of this site remains unknown.
Both human [22],[23] and CnPC [12],[13],[24] are known to show amidolytic activity following incubation with Protac. For human protein C, Protac removes a peptide from the heavy chain, thereby exposing the catalytic site, which cleaves FV and FVIII, as well as the synthetic chromogenic substrate Spectrozyme. To show that the protein purified in this study shared these properties, we investigated its amidolytic activity in a plasma-free system and effects on FV and FVIII activities in canine plasma. Following incubation with Protac, the protein showed amidolytic activities that were directly proportional to the protein concentration (Figure 7A and B). Such activities appeared to be enriched during the purification process as well. Initial attempts to detect protein C amidolytic activities in plasma and the salt precipitate using final protein concentrations of 10 μg/mL did not yield consistently detectable amidolytic activities. Activity was only detected when the final protein concentration was increased to 100 μg/mL (Figure 7C). Compared to plasma, the amidolytic activities of CnPC was increased by ~8.19-fold after salt precipitation and (after adjustment for the final protein concentration used in the assay) ~303-fold following immunoaffinity chromtagraphy (i.e. CnPC activity was enriched ~37.6 times in the final purified form compared to that in the salt precipitate). These results corroborated with results obtained by western immunoblotting using anti-protein C tag antibody.
The preservation of amidolytic activities of the purified protein suggested to us that the protein would most likely retain its anticoagulant properties as well. To confirm this, the effects of this protein (following activation in vitro by Protac) on the activities of FV and FVIII in canine plasma were determined. As predicted, activated form of the protein inhibited the activities of these coagulation factors in a concentration-dependent manner (Figure 8).
Other researchers have concluded that human and CnPC share structural similarities when comparing the expression sequences between the two species [11]. More specifically, they predicted that CnPC was synthesized as a single-chain protein, which was processed into a double-chain protein by removal of a dipeptide adjoining the two chains, accompanied by removal of a signal peptide at the N-terminus. We performed N-terminal sequencing of the light and heavy chains of the purified protein, as well as the Protac-activated heavy chain to determine if findings at the protein level would corroborate with predictions based on alignment of the canine and human expression sequences. Edman sequencing showed that the N-terminus of the light chain began as: ANSFL, and that of the heavy chain began as: DTXQT, where X could be a modified amino acid (a glycosylated asparagine in this case). These results corroborated with the predictions. Following incubation with Protac, CnPC showed a decrease in molecular weight (Figure 9), and the N-terminus of the heavy chain was determined to be GVNGK by Edman sequencing. However, CnPC cDNA obtained by us (data not shown) and others [11],[25] showed that the amino acid residue prior to VNGK should be a leucine, and LC-MS/MS results showed that this amino acid at this position in the peptide had a mass consistent with either leucine or isoleucine (the only two natural amino acids that share the exact same mass). We concluded that the N-terminal sequence of Protac-digested heavy chain was likely LVNGK, and the glycine detected was probably from a second protein or peptide present in the sample.
Discussion
In this study, we showed that a method involving salt precipitation followed by immunoaffinity chromatography using a monoclonal antibody directed against human protein C zymogen (clone HPC4) was useful for purifying protein C zymogen from canine fresh frozen plasma. This method is straightforward and does not require expensive equipment. The matrix and antibody are coupled as a commercially available product, which can be used up to 7 to 10 times, according to the manufacturer’s instructions, without the need for regeneration. The proteins purified from fresh frozen plasma using this method were identified to be CnPC based on tandem mass spectrometry and western immunoblotting. Using a method that quantified the amidolytic activity of the activated form of the protein, we showed that such activity was enriched 8.19-fold by salt precipitation and a further 37.59-fold by immunoaffinity chromatography. A specific antigenic assay for CnPC was not available at the time of this study, rendering absolute quantification of the yield and purification level of the final protein preparation impossible.
Clone HPC4 targets the peptide sequence EDQVDPRLIDGK in human protein C zymogen [26], and the corresponding sequence in CnPC, TDQIDPRLVNGK, differs at four of twelve amino acid residues. Despite these differences, HPC4 showed sufficient cross-reactivity with CnPC for protein purification and western immunoblotting purposes. It has been suggested that calcium ions expose the HPC4-binding site by inducing a conformational change in the protein [26]. Calcium ions may induce a conformational change in CnPC such that differences in primary amino acid sequence do not affect antibody binding.
In humans, protein C is produced by hepatocytes as a single-chain protein [27]. During post-translational modification (PTM), the pre-pro-peptide at the N-terminus is removed [28]. Additionally, an internal arginine-lysine dipeptide is removed, resulting in a heavy chain and light chain linked by a disulphide bond [29]. By performing densitometry on Western blots, Greffe et al. showed that in human adults, 90% of the plasma protein C was double-chained, while 10% was single-chained [30]. They also showed that single-chain human protein C had three isoforms: alpha (MW 62 kDa), beta (MW 58 kDa), and gamma (MW 55 kDa); that the heavy chain had three isoforms alpha (MW 41 kDa), beta (MW 37 kDa), and gamma (MW 34 kDa); and that the light chain had two isoforms (MW 21 kDa and MW 18 kDa).
Our findings show that, similar to human protein C, CnPC exists in different glycoforms as well. On reducing SDS-PAGE, CnPC showed three distinct bands of MW 49 kDa, 44 kDa, and 40 kDa. These bands likely represent three different glycosylated forms of the heavy chain in CnPC, as treatment with a mixture of deglycosylation enzymes reduced these three bands into one single band and resulted in a decreased in MW. A faint band was visible at ~75-kDa; this may represent the single-chain CnPC. Results from LC-MS/MS analysis did confirm the presence of protein C in this band, and no immunoreactivity with HPC4 was detected in this region on western immnoblotting of the reduced form of the protein. We concluded that further studies were necessary to confirm the presence of single-chain CnPC in plasma.
We further determined three N-linked glycosylation sites on the CnPC. However, no further studies were performed to determine the biological significance of the glycans present at these sites. In humans, glycans on protein C appear to affect the functional properties of activated protein C, and the effect of glycosylation at Asn329 is best described among all the N-linked glycosylated sites. This site is conserved across many mammalian species [10]. Compared to wild-type APC, a recombinant APC variant with the glycan eliminated at Asn329 (ACP-N329Q) showed enhanced endothelial barrier protective function by 6-fold and improved anti-apoptotic function by 30-fold in endothelial cells treated with staurosporine. The effects of Asn329 on anticoagulation are not as well defined. APC-N239Q increased the anticoagulant effects [10],[31], but naturally occurring protein C with a mutation at this site (N329T) showed decreased rate in deactivation of FVa following purification from plasma [32].
Canine protein C isolated by the method described here not only preserved the immunogenic properties of the protein, as demonstrated by Western immunoblotting, but also the amidolytic and anticoagulant properties. A commercially available protein C assay has been validated for measuring protein C activity in canine plasma [24]. This assay is based on the amidolytic activity of activated protein C on a chromogenic substrate, which yields a quantifiable change in absorbance upon cleavage by activated protein C. In this study, we showed that CnPC purified by the described method was activatable by Protac and retained its amidolytic properties. Similarly, following ex vivo activation with Protac, CnPC inhibited FV and FVIII activities in a dose-dependent manner, as expected. It is known that human APC inhibits FV and FVIII activities by proteolysis of these factors [1],[2]. The same mechanisms are presumably present in dogs. Compared to human protein C, relatively little is known about CnPC.
Currently, the main clinical application of protein C in veterinary medicine is the differentiation of microvascular dysplasia and portosystemic vascular anomaly in dogs [33]. Combined with other clinical and laboratory data, plasma protein C activity could be used to distinguish portosystemic vascular anomaly from microvascular dysplasia [33]. Although not routinely performed, plasma protein C levels were found to have prognostic value in canine sepsis [13]. Clinical application of purified CnPC is not shown in this manuscript or described elsewhere and warrants further investigation.
Conclusions
We showed that a functional protein C zymogen could be purified by salt precipitation followed by immunoaffinity chromatography using a commercially available matrix coupled to an antibody against human protein C. Upon activation with Protac, the purified protein retained its amidolyic, anti-FV, and anti-FVIII properties. We concluded that this method was useful for purifying CnPC for future functional studies.
Methods
Sources of CnPC
All procedures were carried out with the approval of the Animal Care Committee of the University of Guelph and in compliance with the guidelines provided by the Canadian Council on Animal Care, and informed client consent was obtained. Healthy adult dogs were selected for blood donation based on clinical history, physical examination, routine blood testing that included complete blood count and biochemistry, urinalysis, and testing for heartworm, erhrlichiosis, Lyme disease, anaplasmosis, and hemoplasmosis. Approximately 450 mL of whole blood was collected aseptically by a trained veterinary technician from the jugular vein into a bag containing citrate-phosphate-dextrose and Optisol (sodium, dextrose, adenine, and mannitol) (Teruflex Blood Bag, Terumo, Somerset, NJ) with constant mixing during collection. The bag was immediately centrifuged at 5,000 g at 4°C for 8 min, followed by immediate separation of the plasma from other blood components into the attached satellite bag. The plasma was stored at -20°C for up to one year, after which the fresh frozen plasma was moved to storage at -70°C for up to two years.
Salt precipitation
Unless specified, all procedures described below, including centrifugation, were performed at 4°C. All reagents and containers were pre-chilled. Where protease inhibitors were used, they were prepared by dissolving tablets containing a mixture of protease inhibitors (Complete, EDTA-free, Protease Inhibitor Cocktail Tablets, Roche, Laval, QC) in the solution specified in each of the following steps described below so that the final concentration during incubation, centrifugation, or washing was one tablet per 50 mL of solution, as suggested by the manufacturer.
A bag of fresh frozen plasma was thawed in a cool water bath and divided into 40-mL aliquots in 50-mL BD Falcon polypropylene conical tubes (BD Biosciences, Mississauga, ON). 3.2 mL of 1 M barium chloride (Fisher Scientific, Ottawa, ON) containing protease inhibitors was added one drop at a time to each aliquot with periodic mixing. The plasma mixture was allowed to equilibrate for an additional one hour with constant mixing, followed by centrifugation at 1,500 g for 15 min. The supernatant was discarded, and the precipitate was resuspended in 10 mL of 0.9% NaCl containing protease inhibitors. This plasma protein fraction was centrifuged at 1,500 g for 15 min. After two such washes, the precipitate from two of the original 40-mL plasma aliquots was resuspended in 10 mL of 0.2% EDTA containing protease inhibitors. This suspension was allowed to equilibrate for 30 min with constant mixing. Ammonium sulphate pellets (Fisher Scientific, Ottawa, ON) were added to the suspension in small quantities with periodic mixing to achieve a final saturation level of 40% (243 Mg per mL solution). This was followed by 10 min of constant mixing to ensure complete solubilization of the ammonium sulphate pellets and equilibration. The suspension was centrifuged at 1,500 g for 30 min. The supernatant was separated. Ammonium sulphate pellets were added in the same manner described above to achieve a final saturation level of 65% (168 Mg per mL solution). This was followed by 10 min of constant mixing and then centrifugation at 1,500 g for 30 min. The supernatant was removed. The precipitate was dissolved in 2 mL of water. A buffer exchange into 20 MM Tris, pH 7.5, containing 0.1 M sodium chloride, 1 MM calcium chloride, and protease inhibitors was performed immediately using 15-mL Amicon Ultra filter units with molecular weight cut-off of 30 kDa (Millipore, Billerica, MA). After the final spin, the retentate was resuspended in the same buffer containing protease inhibitors so that proteins isolated from two of the original 40-mL plasma aliquots was resuspended in 1 mL of solution.
Immunoaffinity chromatography
The following procedure was carried out at 4°C with chilled reagents and containers. The same protease inhibitor tablets were used in the same manner as described above. The matrix described in the following paragraphs refers to a commercially available product (Anti-Protein C Affinity Matrix; Roche, Laval, QC), which consists of agarose beads coupled to a monoclonal antibody (clone HPC4; mouse IgG1κ) against a 12-amino acid sequence (EDQVDPRLIDGK) of human protein C.
One mL of matrix was packed in the column by gravitation. The matrix was equilibrated with 10 bed-volumes of equilibration buffer (20 MM Tris, pH 7.5, containing 0.1 M sodium chloride, 1 MM calcium chloride). The buffer was drained. The bottom end of the column was capped and 1 mL of equilibration buffer containing protease inhibitors was applied to the matrix. One mL of the plasma protein suspension obtained after barium chloride and ammonium sulphate precipitation was applied to the matrix. The final concentration of calcium chloride was adjusted to 1 MM. The top of the column was capped. The column, now containing the proteins and the matrix, was inverted a few times to loosen the matrix, followed by gentle constant mixing on an end-to-end rocker for two hours. The matrix was then allowed to settle for 10 min. The bottom was uncapped and the unbound protein fraction was allowed to drain. A 23-g needle was applied to the bottom of the column. This set-up resulted in a flow rate of approximately 0.5 mL per minute. The matrix was washed with 30 bed-volumes of wash buffer (20 MM Tris, pH 7.5, containing 1.0 M sodium chloride, 1 MM calcium chloride) containing protease inhibitors. A clean cap was applied to the bottom of the column. Five mLs of elution buffer (20 MM Tris, pH 7.5, containing 0.1 M sodium chloride and 5 MM EDTA) was applied to the matrix, followed by 30 minutes of incubation without mixing. The bottom of the column was uncapped and the 5 mL of elution buffer, now containing the eluted proteins, was collected through a fresh 23-g needle. A buffer exchange into the equilibrium buffer was performed using 4-mL Amicon Ultra filter units with molecular weight cut-off of 30 kDa (Millipore, Billerica, MA). The retentate was concentrated to 250–500 μL in the last spin.
SDS-PAGE
The eluted proteins were subjected to one-dimensional SDS-PAGE under both denaturing and non-denaturing conditions. Buffers were reconstituted from stock solutions (Bio-Rad, Mississauga, ON) for a modular Mini-Protean II system (Bio-Rad, Mississauga, ON). Laemmli sample buffer (Bio-Rad, Mississauga, ON) and β-mercaptoehthanol (Sigma-Aldrich, Oakville, ON) were purchased and used according to the manufacturer’s instructions. Briefly, 1 part β-mercaptoethanol was added to 9 parts of Laemmli sample buffer. For denaturing SDS-PAGE, 5 μg of proteins was combined with this sample buffer at a 1:1 (v:v) ratio, followed by incubation at 65°C for 10 minutes. For non-denaturing SDS-PAGE, 5 μg of eluted proteins was directly combined with the Laemmli sample buffer at a 1:1 (v:v) ratio and kept at room temperature prior to electrophoresis without addition of β-mercaptoethanol. The protein samples were loaded into sample wells of a 1.0-mm thick, 15% Tris-HCl polyacrylamide gel and then subjected to SDS-PAGE at 220 V on ice for 75 min. The gels were then stained with silver for visualization of protein bands. Briefly, the gels were fixed in a solution of 40% ethanol and 10% glacial acetic acid, followed by rinsing with 20% ethanol and then water. The gels were then incubated with 0.2 g/L of sodium thiosulfate for one min, followed by rinsing with water, and then incubation with 2 g/L of silver nitrate for 30 minutes. Following a brief rinse with water, the gels were developed in a solution containing 30 g/L of sodium carbonate, 1.4 mL/L of 37% formaldehyde, and 10 Mg/L of sodium thiosulfate for up to 15 min until the desired intensity of the protein bands was achieved. Molecular weights of the visible bands were compared against markers of known molecular weight (Precision Plus Protein WesternC Standards, Bio-Rad, Mississauga, ON).
Mass spectrometry and protein identification
The silver-stained protein bands visible on SDS-PAGE were excised, placed in siliconized microcentrifuge tubes containing 1% v/v acetic acid, and submitted on dry ice for LC-MS/MS (The Advanced Protein Technology Centre, Hospital for Sick Children, Toronto, ON). In-gel tryptic digestion was performed as previously described [34]. Briefly, the silver-stained gel incubated with a solution containing 30 MM of potassium ferricyanide and 100 MM of sodium thiosulfate at a 1:1 (v/v) ratio for 15 minutes, followed by rinsing with water and then with 50 MM of ammonium bicarbonate. The gel was shrunk with 50% acetonitrile containing 25 MM of ammonium bicarbonate for 10 min. The proteins were then reduced in 10 MM of dithiothreitol (DTT) at 56°C for 30 min, followed by alkylation in 100 MM of iodoacetamide for 15 min. The gel was shrunk for the second time for 15 minutes. The gel was incubated overnight with 50 MM ammonium bicarbonate containing 13 ng/μL of porcine trypsin (Sequencing Grade Modified Trypsin, Promega, Madison, WI) at 37°C. The digested peptides were loaded onto a 150 μm ID pre-column (Magic C18, Michrom Biosciences, Auburn, CA) at 4 μL per min and separated over a 75 μm ID analytical column packed into an emitter tip containing the same packing material. The peptides were eluted over 60 minutes at 400 nL per minute using a 0–40% acetonitrile gradient in 0.1% formic acid into a linear ion trap mass spectrometer (LTQ-Orbitrap) (Thermo-Fisher, San Jose, CA) operated in a data-dependent mode. One MS followed by 6 MS/MS scans were obtained per cycle. Tandem mass spectra were extracted by Extract_msn.exe 4.0 and Bioworks 3.3 (Thermo-Fisher, San Jose, CA). Charge state deconvolution and deisotoping were performed. The generated data files were analyzed using Mascot 2.3 (Matrix Science, London, UK; version Mascot). Mascot was set up to search the NCBInr_20090719 database (selected for Canis lupus familiaris, 34574 entries) assuming that digestion enzyme was non-specific. Mascot was searched with a fragment ion mass tolerance of 0.50 Da and a parent ion tolerance of 20 ppm. Iodoacetamide derivative of cysteine was specified in Mascot as a fixed modification. Deamidation of asparagine and glutamine and oxidation of methionine were specified in Mascot as variable modifications. Scaffold (version Scaffold_3_00_08, Proteome Software Inc., Portland, OR) was used to validate MS/MS based peptide and protein identifications. Peptide identifications were accepted if they could be established at greater than 95.0% probability as specified by the Peptide Prophet algorithm [35]. Protein identifications were accepted if they could be established at greater than 99.9% probability and contained at least 5 identified peptides. Protein probabilities were assigned by the Protein Prophet algorithm [36]. Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony.
Western immunoblotting
Following SDS-PAGE, the separated proteins were electrotransfered at 100 V for 1 h onto Immobilon-P polyvinylidene fluoride (PVDF) membrane (Millipore, Billerica, MA). The membrane was blocked for 2 h with phosphate-buffered saline, pH7.4 containing 0.1% Tween-20 (PBS-T) and 5% non-fat skim milk. The membrane was then incubated with polyclonal rabbit anti-protein C (HPC4) tag antibody (GenScript, Piscataway, NJ) at 1:1,000 dilution (5 × 10-4 μg/μL) in PBS-T containing 5% non-fat skim milk for 60 min, followed by three 10-min washes in 20 mL of PBS-T. The membrane was then incubated with horseradish peroxidase (HRP)-conjugated polyclonal goat anti-rabbit IgG antibody (Dako, Denmark) at 1:2,000 dilution (1.25 × 10-4 μg/μL) in PBS-T containing 5% non-fat skim milk for 2 h and washed three times as described above. The membrane was then incubated with HRP conjugate for the MW standards (Precision Plus Protein StrepTactin-HRP Conjugate, Bio-Rad, Mississauga, ON) at 1:5,000 dilution (2 × 10-4 μg/μL) in PBS-T containing 5% non-fat skim milk for 20 min, followed by three washes as described above. The membrane was then incubated with the Amersham ECL Western Blotting Detection Reagents (GE Healthcare, Buckinghamshire, UK) according to the manufacturer’s instructions. Chemiluminesence was detected using the ChemiDoc XRS + System with Image Lab Software (BioRad, Mississauga, ON). Relative intensity of detected bands was determined by densitometry using Image Lab Software.
Additionally, the purified proteins were subjected to Western immunoblotting using a mouse monoclonal anti-human protein C heavy chain antibody, clone GMA-093 (Upstate, Lake Placid, NY), as primary antibody at 1:1,000 dilution (1 × 10-3 μg/μL) and a polyclonal rabbit anti-mouse IgG antibody coupled to horseradish peroxidase (Abcam, Cambridge, MA) as secondary antibody at 1: 1,000 dilution (5 × 10-4 μg/μL). Procedures for Western immunoblotting and detection of immunoreactive proteins were carried out in the same way as described above.
N-linked glycosylation mapping
For visualization of deglycosylated proteins on SDS-PAGE, the isolated proteins were incubated with a Protein Deglycosylation Mix (New England Biolabs, Mississauga, ON) according to the manufacturer’s instructions. Briefly, 9 μg of purified proteins were incubated with 40 MM DTT and 0.5% SDS at 100°C for 10 min and cooled on ice for 10’s. The denatured proteins were then incubated with 50 MM sodium phosphate, pH7.5, containing 1% NP-40, O-glycosidase (4 × 103 units/μL), beta-N-acetylglucosaminidase (4 × 10-1 units/μL), neuraminidase (5 units/μL), beta-1-4-galactosidase (8 × 10-1 units/μL), and PNGase F (5 × 10 units/μL) for 4 h at 37°C in a 50-μL reaction. Purified proteins subjected to the same procedures in the absence of the deglycosylation enzymes served as a negative control. Bovine fetuin (a glycosylated protein) was used as a positive control. SDS-PAGE and silver-staining were performed as described above. Proteins isolated from three different dogs were studied.
To identify potential N-linked glycosylated sites, the predicted amino acid sequence of CnPC was searched for the consensus sequences N-X-T and N-X-S [19]. For IGOT [21], 10 μg of purified protein was first concentrated and diafiltered into 10 μL of 18O-labelled water (97% atom labelled) (Sigma-Aldrich, Oakville, ON) using a Microcon Centrifugal Filter Device with molecular weight cut-off at 10 kDa (Millipore, Billerica, MA). Following addition of 2 μL of 400 MM DTT and 5% SDS, the proteins were denatured (total volume: 12 μL) at 100°C for 10 min and cooled on ice for 10’s. The following were then added to the mixture: 20 μL of 18O-labelled water, 4 μL of 500 MM sodium phosphate, pH7.5, 2 μL of 10% NP-40, and 2 μL of PNGase F (5 × 102 units/μL). This mixture was incubated for 4 h at 37°C with periodic mixing. At the end of the reaction period, the mixture was diafiltered into TBS constituted with Milli-Q (mostly H216O) water using a Microcon Centrifugal Filter Device with molecular weight cut-off at 10 kDa (Millipore, Billerica, MA). The proteins were subjected to trypsin digestion and LC-MS as described above. The software PEAKS (Bioinformatics Solutions, Waterloo, ON) was used to search for +3 mass unit shifts at the predicted N-linked glycosylated sites.
Digestion of CnPC with Protac
Five μg of purified proteins in TBS were incubated with 0.1 unit of Protac (Sekisui Diagnostics, Stamford, CT), an activator of human [22] and canine [24] protein C and a substance derived from the venom of Agkistrodon c. contortrix, or same volume of TBS at 37°C for 15 min. All samples were then subjected to non-reducing 12% SDS-PAGE, followed by silver-staining or Western immunoblotting using a polyclonal rabbit anti-protein C (HPC4) tag antibody as described above.
N-terminal sequence determination
The N-termini of the heavy chain, light chain, and Protac-digested heavy chain were analyzed by Edman degradation [37]. Purified proteins were subjected to reducing 12% SDS-PAGE and electrotransfered overnight at 30 V in CAPS buffer (Sigma-Aldrich, Oakville, ON) onto Immobilon-PSQ PVDF membrane (Millipore, Billerica, MA). The membrane was incubated for 10 mins in 0.2% (w/v) Ponceau S (Sigma-Aldrich, Oakville, ON) in 1% (v/v) acetic acid and washed three times in Milli-Q water. It was then stored in a moist chamber and submitted for sequencing (The Advanced Protein Technology Centre, Hospital for Sick Children, Toronto, ON). Stained protein bands were excised and analyzed by the Procise Protein Sequencing System (Applied Biosystems, Foster City, CA) as previously described [38].
Detection of amidolytic activity
The amidolytic activity of the purified protein was measured by a method modified from previously described studies [22]-[24],[39] using Protac and Spectrozyme [H-D-(γ-carbobenzoxyl)-lysyl-prolyl-arginine-paranitroanilide diacetate salt], a chromogenic substrate for activated protein C. Protac (Sekisui Diagnostics, Stamford, CT) and Sepctrozyme (Sekisui Diagnostics, Stamford, CT) were reconstituted with sterile water to 0.5 unit/mL and 10 MM, respectively. Spectrozyme was further diluted by combining 7 parts of 50 MM Tri-HCl, pH 8.4 and 1 part Spectrozyme. The purified protein was adjusted to 0, 20, 40, 60, 80, and 100 μg/mL by dilution with TBS. Fifty μL of purified protein was incubated with 50 μL of Protac (0.5 unit/mL) for 5 min at 37°C. Forty μL of this mixture was placed in each well in a 96-well plate. One hundred and eighty μL of Spectrozyme was added to each of all samples simultaneously using a multi-channel pipette. Final protein concentrations in the samples were 0, 2, 4, 6, 8, and 10 μg/mL. Absorbance at λ = 405 nm was measured at 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10 min following addition of Spectrozyme using a BioTek Synergy HT plate reader (BioTek, Winooski, VT) and Gen5 version 1.11 (BioTek, Winooski, VT). All samples were analyzed in duplicate and the average of the two measurements was used as one true replicate for statistical analysis using GraphPad Prism (GraphPad Software, La Jolla, CA). For each concentration of protein (final concentrations: 0, 2, 4, 6, 8, 10 μg/mL), three replicates were obtained using proteins isolated from three different dogs. The mean absorbance ± standard error of mean (SEM) for each time point was plotted against time. For each time point, the mean absorbance ± SEM was plotted against concentration of purified protein. Additionally, amidolytic activities in canine plasma and the salt precipitate obtained following barium chloride and ammonium sulphate precipitation were determined in a similar manner using final protein concentration of 100 μg/mL.
Canine factor V & VIII assays
Canine coagulation factor assays were performed based on published methods [40],[41] using the KC 4 Delta coagulation analyzer (Stago, Parsipanny, NJ). Canine plasma pool was generated by pooling 0.5-2 mL of citrated plasma from each of ten clinically healthy adult dogs. The plasma pool was stored in aliquots at -70°C. For calibration, the activity level of each coagulation factor present in the undiluted plasma pool was designated to be 100% in the calibration process; and the plasma pool was serially diluted five times by a factor of two to obtain samples with activity levels of 50%, 25%, 12.5%, 6.25%, and 3.125%. The isolated canine protein was activated with Protac as described above and prepared by dilution with TriniCLOT imidazole buffer (Diagnostica Stago, Toronto, ON) to obtain these concentrations of 1, 0.8, 0.6, 0.4, 0.2 and 0 μg/μL. Fifty μL of factor V-deficient plasma (Trinity Biotech, Jamestown, NY) or factor VIII-deficient plasma (Trinity Biotech, Jamestown, NY) was incubated with 50 μL of pooled canine plasma and 3 μL of purified protein for 60’s at 37°C with constant mixing, followed by incubation with 10 μL of TriniCLOT aPTT reagent S (Diagnostica Stago, Toronto, ON) at 37°C for 300 sec. One hundred μL of 0.025 M CaCl2 (Diagnostica Stago, Toronto, ON) was added to the mixture and the time from then to clot formation was recorded [referred to as modified activated partial thromboplastin time (APTT)]. All samples were analyzed in duplicate and the average of the two measurements was used as one true replicate for statistical analysis using GraphPad Prism (GraphPad Software, La Jolla, CA). The clotting times were compared against a standard curve generated within the same assay using canine pooled plasma diluted with TriniCLOT imidazole buffer. Percent activity of FV or FVIII was plotted against concentration of the purified canine protein.
Authors’ contributions
VMW designed the study, carried out the experimental work (except for mass spectrometry), analyzed and interpreted the data, and prepared the manuscript. DB, MAH, and RDW participated in the design of the study and critically revised the manuscript for important intellectual content. PT carried out the mass spectral analysis and interpreted the data. All authors read and approved the final manuscript.
Abbreviations
- FV:
-
Factor V
- FVIII:
-
Factor VIII
- FDA:
-
Food and Drug Administration
- EMEA:
-
European Medicines Agency
- PROWESS:
-
Recombinant human activated protein C worldwide evaluation in severe sepsis
- CnPC:
-
Canine protein C
- CnAPC:
-
Canine activated protein C
- SDS-PAGE:
-
Sodium dodecyl sulphate polyacrylamide gel electrophoresis
- DELTA-BLAST:
-
Domain-enhanced lookup time-accelerated BLAST
- LC-MS/MS:
-
Liquid chromatography followed by tandem mass spectral analysis
- IGOT:
-
Isotope-coded glycosylation site-specific tagging
- PTM:
-
Post-translational modification
- PVDF:
-
Polyvinylidene fluoride
- PBS-T:
-
Phosphate-buffered saline, pH 7.4 containing 0.1% Tween-20
- SEM:
-
Standard error of mean
- APTT:
-
Activated partial thromboplastin time
References
Van de Waart P, Bruls H, Hemker HC, Lindhout T: Functional properties of factor Va subunits after proteolytic alterations by activated protein C. Biochim Biophys Acta. 1984, 799: 38-44. 10.1016/0304-4165(84)90324-6.
Fulcher CA, Gardiner JE, Griffin JH, Zimmerman TS: Proteolytic inactivation of human factor VIII procoagulant protein by activated human protein C and its analogy with factor V. Blood. 1984, 63: 486-489.
Bernard GR, Vincent JL, Laterre PF, LaRosa SP, Dhainaut JF, Lopez-Rodriguez A, Steingrub JS, Garber GE, Helterbrand JD, Ely EW, Fisher CJ: Efficacy and safety of recombinant human activated protein C for severe sepsis. N Engl J Med. 2001, 344: 699-709. 10.1056/NEJM200103083441001.
Ranieri VM, Thompson BT, Barie PS, Dhainaut J-F, Douglas IS, Finfer S, Gårdlund B, Marshall JC, Rhodes A, Artigas A, Payen D, Tenhunen J, Al-Khalidi HR, Thompson V, Janes J, Macias WL, Vangerow B, Williams MD: Drotrecogin alfa (activated) in adults with septic shock. N Engl J Med. 2012, 366: 2055-2064. 10.1056/NEJMoa1202290.
Martí-Carvajal AJ, Solà I, Gluud C, Lathyris D, Cardona AF: Human recombinant protein C for severe sepsis and septic shock in adult and paediatric patients.Cochrane Database Syst Rev 2012, 12:CD004388.
Lilly Announces Withdrawal of Xigris Following Recent Clinical Trial Results [], [http://newsroom.lilly.com/releasedetail.cfm?releaseid=617602]
Kerschen EJ, Fernández JA, Cooley BC, Yang XV, Sood R, Mosnier LO, Castellino FJ, Mackman N, Griffin JH, Weiler H: Endotoxemia and sepsis mortality reduction by non-anticoagulant activated protein C. J Exp Med. 2007, 204: 2439-2448. 10.1084/jem.20070404.
Martí-Carvajal AJ, Solà I, Lathyris D, Cardona AF: Human recombinant activated protein C for severe sepsis.Cochrane Database Syst Rev 2011:CD004388.
Guo H, Singh I, Wang Y, Deane R, Barrett T, Fernández JA, Chow N, Griffin JH, Zlokovic BV: Neuroprotective activities of activated protein C mutant with reduced anticoagulant activity. Eur J Neurosci. 2009, 29: 1119-1130. 10.1111/j.1460-9568.2009.06664.x.
Ní Ainle F, O’Donnell JS, Johnson JA, Brown L, Gleeson EM, Smith OP, Preston RJS: Activated protein C N-linked glycans modulate cytoprotective signaling function on endothelial cells. J Biol Chem. 2011, 286: 1323-1330. 10.1074/jbc.M110.159475.
Leeb T, Kopp T, Deppe A, Breen M, Matis U, Brunnberg L, Brenig B: Molecular characterization and chromosomal assignment of the canine protein C gene. Mamm Genome. 1999, 10: 134-139. 10.1007/s003359900958.
de Laforcade AM, Freeman LM, Shaw SP, Brooks MB, Rozanski EA, Rush JE: Hemostatic changes in dogs with naturally occurring sepsis. J Vet Intern Med. 2003, 17: 674-679. 10.1892/0891-6640(2003)017<0674:HCIDWN>2.3.CO;2.
de Laforcade AM, Rozanski EA, Freeman LM, Li W: Serial evaluation of protein C and antithrombin in dogs with sepsis. J Vet Intern Med. 2008, 22: 26-30. 10.1111/j.1939-1676.2007.0021.x.
Fourrier F, Chopin C, Goudemand J, Hendrycx S, Caron C, Rime A, Marey A, Lestavel P: Septic shock, multiple organ failure, and disseminated intravascular coagulation. Compared patterns of antithrombin III, protein C, and protein S deficiencies. Chest. 1992, 101: 816-823. 10.1378/chest.101.3.816.
Lorente JA, García-Frade LJ, Landín L, de Pablo R, Torrado C, Renes E, García-Avello A: Time course of hemostatic abnormalities in sepsis and its relation to outcome. Chest. 1993, 103: 1536-1542. 10.1378/chest.103.5.1536.
Boldt J, Papsdorf M, Rothe A, Kumle B, Piper S: Changes of the hemostatic network in critically ill patients-is there a difference between sepsis, trauma, and neurosurgery patients?. Crit Care Med. 2000, 28: 445-450. 10.1097/00003246-200002000-00026.
Stenflo J: A new vitamin K-dependent protein. Purification from bovine plasma and preliminary characterization. J Biol Chem. 1976, 251: 355-363.
Boratyn GM, Schäffer AA, Agarwala R, Altschul SF, Lipman DJ, Madden TL: Domain enhanced lookup time accelerated BLAST.Biol Direct 2012, 7:12.
Gavel Y, Heijne von G: Sequence differences between glycosylated and non-glycosylated Asn-X-Thr/Ser acceptor sites: implications for protein engineering. Protein Eng. 1990, 3: 433-442. 10.1093/protein/3.5.433.
Apweiler R, Hermjakob H, Sharon N: On the frequency of protein glycosylation, as deduced from analysis of the SWISS-PROT database. Biochim Biophys Acta. 1999, 1473: 4-8. 10.1016/S0304-4165(99)00165-8.
Kaji H, Isobe T: Stable isotope labeling of N-glycosylated peptides by enzymatic deglycosylation for mass spectrometry-based glycoproteomics. Methods Mol Biol. 2013, 951: 217-227. 10.1007/978-1-62703-146-2_14.
Martinoli JL, Stocker K: Fast functional protein C assay using Protac, a novel protein C activator. Thromb Res. 1986, 43: 253-264. 10.1016/0049-3848(86)90145-3.
Madden RM, Oppenheimer C, Wydro R, Marlar RA: Recombinant human protein C: comparative functional studies with human plasma protein C. Thromb Res. 1990, 57: 425-435. 10.1016/0049-3848(90)90258-E.
Fry MM, Snyder KR, Tobias KM, Williamson BG, Reed GA: Protein C activity in dogs: adaptation of a commercial human colorimetric assay and evaluation of effects of storage time and temperature.Vet Med Int 2011, 2011:751849.
Leeb T, Pfeiffer I, Kopp T, Deppe A, Brenig B: Analysis of canine protein C gene polymorphisms. Anim Genet. 1999, 30: 237-238. 10.1046/j.1365-2052.1999.00404-17.x.
Stearns DJ, Kurosawa S, Sims PJ, Esmon NL, Esmon CT: The interaction of a Ca2 + -dependent monoclonal antibody with the protein C activation peptide region. Evidence for obligatory Ca2+ binding to both antigen and antibody. J Biol Chem. 1988, 263: 826-832.
Fair DS, Marlar RA: Biosynthesis and secretion of factor VII, protein C, protein S, and the protein C inhibitor from a human hepatoma cell line. Blood. 1986, 67: 64-70.
Foster DC, Yoshitake S, Davie EW: The nucleotide sequence of the gene for human protein C. Proc Natl Acad Sci U S A. 1985, 82: 4673-4677. 10.1073/pnas.82.14.4673.
Foster DC, Rudinski MS, Schach BG, Berkner KL, Kumar AA, Hagen FS, Sprecher CA, Insley MY, Davie EW: Propeptide of human protein C is necessary for gamma-carboxylation. Biochemistry. 1987, 26: 7003-7011. 10.1021/bi00396a022.
Greffe BS, Manco-Johnson MJ, Marlar RA: Molecular forms of human protein C: comparison and distribution in human adult plasma. Thromb Haemost. 1989, 62: 902-905.
Grinnell BW, Walls JD, Gerlitz B: Glycosylation of human protein C affects its secretion, processing, functional activities, and activation by thrombin. J Biol Chem. 1991, 266: 9778-9785.
Simioni P, Kalafatis M, Millar DS, Henderson SC, Luni S, Cooper DN, Girolami A: Compound heterozygous protein C deficiency resulting in the presence of only the beta-form of protein C in plasma. Blood. 1996, 88: 2101-2108.
Toulza O, Center SA, Brooks M, Erb HN, Warner KL, Deal W: Evaluation of plasma protein C activity for detection of hepatobiliary disease and portosystemic shunting in dogs. J Am Vet Med Assoc. 2006, 229: 1761-1771. 10.2460/javma.229.11.1761.
Shevchenko A, Wilm M, Vorm O, Mann M: Mass spectrometric sequencing of proteins silver-stained polyacrylamide gels. Anal Chem. 1996, 68: 850-858. 10.1021/ac950914h.
Keller A, Nesvizhskii AI, Kolker E, Aebersold R: Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002, 74: 5383-5392. 10.1021/ac025747h.
Nesvizhskii AI, Keller A, Kolker E, Aebersold R: A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003, 75: 4646-4658. 10.1021/ac0341261.
Edman P: A method for the determination of amino acid sequence in peptides.Arch Biochem 1949, 22:475.
Fan X, Tkachyova I, Sinha A, Rigat B, Mahuran D: Characterization of the biosynthesis, processing and kinetic mechanism of action of the enzyme deficient in mucopolysaccharidosis IIIC.PLoS One 2011, 6:e24951.
Takahashi H, Hanano M, Tatewaki W, Shibata A: Fast functional assay of protein C in whole plasma using a snake venom activator: evaluation in patients with congenital and acquired protein C deficiencies. Clin Chim Acta. 1988, 175: 217-225. 10.1016/0009-8981(88)90098-8.
Bauer N, Moritz A: Characterisation of changes in the haemostasis system in dogs with thrombosis. J Small Anim Pract. 2013, 54: 129-136. 10.1111/jsap.12037.
Mischke R: Optimization of coagulometric tests that incorporate human plasma for determination of coagulation factor activities in canine plasma. Am J Vet Res. 2001, 62: 625-629. 10.2460/ajvr.2001.62.625.
Acknowledgements
This study was generously funded by Pet Trust at the Ontario Veterinary College (OVC). Dr. Wong’s position was supported by the OVC PhD Fellowship. The authors are grateful for Miss Bette Anne Quinn and Dr. GenevièVe Labbè for their expert advice.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wong, V.M., Bienzle, D., Hayes, M.A. et al. Purification of protein C from canine plasma. BMC Vet Res 10, 251 (2014). https://doi.org/10.1186/s12917-014-0251-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12917-014-0251-2