The bovine spongiform encephalopathy (BSE) epidemic and the emergence of a new human variant of Creutzfeldt-Jakob Disease (vCJD) have led to profound changes in the production and trade of agricultural goods. The rapid tests currently approved for BSE monitoring in slaughtered cattle are all based on the detection of the disease related isoform of the prion protein, PrPd, in brain tissue and consequently are only suitable for post-mortem diagnosis. Objectives: In instances such as assessing the health of breeding stock for export purposes where post-mortem testing is not an option, there is a demand for an ante-mortem test based on a matrix or body fluid that would permit easy access and repeated sampling. Urine and urine based analyses would meet these requirements.
Two dimensional differential gel eletrophoresis (2D-DIGE) and mass spectrometry analyses were used to identify proteins exhibiting differential abundance in the urine of BSE infected cattle and age matched controls over the course of the disease. Multivariate analyses of protein expression data identified a single protein able to discriminate, with 100% accuracy, control from infected samples. In addition, a subset of proteins were able to predict with 85% ± 13.2 accuracy the time post infection that the samples were collected.
These results suggest that in principle it is possible to identify biomarkers in urine useful in the diagnosis, prognosis and monitoring of disease progression of transmissible spongiform encephalopathy diseases (TSEs).
Bovine Spongiform Encephalopathy (BSE) was first described in the UK in 1985 and became an epidemic that peaked with 37, 280 cases reported in 1992. BSE and other Transmissible Spongiform Encephalopathy (TSE) diseases are untreatable, uniformly fatal degenerative syndromes of the central nervous system (CNS). A post-mortem examination revealing characteristic deposits of an insoluble host encoded protein, astrocytosis and spongiosis is required for definitive diagnosis. The characteristic protein deposits are formed by the accumulation of misfolded isoforms of a host-encoded protein, PrPc, or prion protein. The disease associated isoforms are derived from the host protein, PrPc, by a posttranslational process and are often distinguished by their partial resistance to proteinase K digestion. The term PrPd is used to denote the presence of abnormal accumulations or isoforms of PrP detected by any method without prejudice as to its biochemical properties, its infectivity, source or host range.
The new human variant of Creutzfeldt-Jakob Disease (vCJD) identified in 1996 is thought to have been caused by dietary exposure to BSE infected products. In contrast to typical cases of classical CJD, vCJD seems to affect predominantly young adults. Risk reduction measures implemented in response to BSE and the emergence of this new disease led to profound changes in the production and trade of agricultural goods. To minimize the risk of disease transmission to consumers specified risk materials (SRM), constituting tissues known to harbour high levels of infectivity such as the brain and spinal cord, have been removed from the food chain. In addition, the testing of risk animals and of all slaughtered animals above the age of thirty months for BSE is a requisite for access to many potential markets for beef products. The rapid tests currently approved for BSE monitoring in slaughtered cattle are all based on the detection of the disease related isoform of the prion protein, PrPd, in brain tissue and consequently are only suitable for post-mortem diagnosis.
A reliable ante-mortem test would provide an alternative to the routine culling of herds when a confirmed case of BSE is detected. More importantly, in instances such as assessing the health of breeding stock where post-mortem testing is not an option, there is a demand for ante-mortem tests based on a matrix or body fluid that would permit repeated sampling. The development of such assays, based on the detection of PrPd, have been complicated by the extremely low amounts of PrPd present in accessible tissues, or in body fluids such as cerebrospinal fluid (CSF), blood and urine [5-8]. Furthermore, the demonstration that most infectivity is associated with protease sensitive forms of PrPd also calls into question the reliability of tests reliant on the association of prion infectivity with the presence of a proteinase K resistant fragment that is measured by Western blotting, enzyme-linked immunosorbent assay, or immunohistochemistry.
The advent of the protein misfolding cyclic amplification assay (PMCA), offered a possible solution to this problem. Indeed, the use of modified PMCA assays using relatively defined components has resulted in the detection of PrPd in the CSF, serum and urine of terminal stage hamsters. However, PMCA based methods require further investigation and validation before they are ready for routine use [11-15]. Thus, the identification of alternative biomarkers in accessible tissues or body fluids applicable to the development of diagnostic tests remains a relevant approach.
Urine, due to its ease of collection and comparatively less complex protein profile, is perhaps the ideal matrix for surveillance provided a sufficiently sensitive and specific alternative biomarker for disease can be identified. Previously, it was demonstrated that the presence of protease resistant light chain immunoglobulin in urine may constitute a surrogate marker for prion diseases[6,16]. Building on these results we have used two dimensional differential gel electrophoresis (2D-DIGE) and mass spectrometry analyses to demonstrate that the relative abundance of other proteins in the urine of BSE infected cattle and age matched controls change over the course of the disease. These analyses, performed on biological replicates, identified a single protein able to discriminate between control and infected cattle throughout the course of the disease as well as a subset of proteins able to accurately identify the collection date of the samples. The results indicate that biomarkers in urine may be useful in the diagnosis, prognosis and monitoring of disease progression of transmissible spongiform encephalopathy diseases (TSEs).
Urine was collected from 4 Simmental cross-breed calves that were orally infected with BSE and 4 age matched controls at 8 month intervals throughout the course of the disease. The calves were infected at 4 months of age. All cattle were scored every second month for clinical signs. BSE was confirmed by immunohistochemistry of the obex. Urine samples were frozen immediately after collection and stored at -80°C until processing. This generated 4 infected and 4 control biological replicates at each of the 6 time points. One sample from this set was not obtained resulting in 47 samples in total.
Urine samples (80–90 mL) were thawed overnight at 4°C. Insoluble particles were removed by a pre-spin, 4°C at 415 g for 5 minutes. The soluble fraction was concentrated with a 5 K MWCO Centricon Plus-70 centrifugal filter unit (Millipore) in a swinging bucket rotor at 4°C and 3400 g for approximately 20 minutes or until volumes reached less than 4 mL. The urine was further concentrated with an Amicon 4 ml 5 K MWCO centrifugal filter unit (Millipore) at 4°C and 7000 g until volumes reached 200 μl.
The concentrated urine samples were purified using a 2D Clean Up kit (GE Healthcare) according to the manufacturer's recommendations. The resulting protein pellets were resuspended in 100 μL of Rehydration Buffer (0.03 M Tris, 8 M Urea, 2 M Thiourea, 2% Chaps, pH 8.5). Samples were adjusted to pH8.5 with the addition of 1–5 μl of 0.05 M NaOH. The concentration of each sample was determined using a 2D Quant kit (GE Healthcare) according to the manufacturer's recommendations. The pooled internal standard was created by combining 100 μg of each sample.
2D Gel Electrophoresis
CyDye™ (GE Healthcare) minimal labeling was performed as per the manufacturer's recommendations (400 pmol: 50 μg) with the Cy2 label reserved for the pooled sample. The control and infected samples were labeled in a randomized manner with either Cy3 or Cy5. An equal volume of 2× Rehydration Buffer (0.03 M Tris, 8 M Urea, 2 M Thiourea, 2% Chaps, pH 8.5, 4 mg/ml Dithiothreitol (DTT), 1% IPG buffer pH 4–7) was added to a mixture comprised of 30 μg each of the labeled Pooled, Control and Infected samples. After a 10 minute incubation on ice approximately 400 μl of 1× Rehydration Buffer (0.03 M Tris, 8 M Urea, 2 M Thiourea, 2% Chaps, pH 8.5, 2 mg/ml DTT, 0.5% IPG buffer pH 4–7) and 5 μl of 1% Bromophenol blue (10 mM TrisCl pH 8.5) were added to bring the volume up to a total of 450 μl. The labelled samples were loaded onto a reswelling tray and overlaid with a 24 cm Immobiline DryStrip pH 4–7 (GE Healthcare) and DryStrip Cover Fluid and allowed to rehydrate overnight at room temperature. The strip was transferred to a Manifold filled with 108 mL of DryStrip Cover Fluid and placed on an Ettan IPGphor3 focusing apparatus that was programmed as follows: Step 30 V 8 hrs, Step 500 V 1 hr, Step 1000 V 1 hr, Grad 10000 V 3 hrs and Step 10000 V 3 hrs. A final focusing program: Grad 10000 V 0:20 hr and Step 10000 V until the volt hours reached a total of 55000 was added as required.
Completed isoelectric focusing (IEF) runs were stored at -80°C until the second dimension was run. IEF strips were prepared for second dimension gels by incubating in two different Equilibration Buffer solutions (50 mM Tris-Cl pH8.8, 6 M Urea, 30% Glycerol, 2% SDS, 0.2% Bromophenol Blue, supplemented with either 65 mM DTT – 1st incubation or 135 mM Iodoacetamide – 2nd incubation) for 15 minutes each with gentle rocking.
The prepared IEF strips were placed on 15–20% gradient gels between low fluorescent glass plates (NextGen Sciences). After sealing in place with a 1% agarose solution, the gels were placed in the Ettan DALT6 unit (GE Healthcare) and run at 2 W overnight and then at 100 W until a total of 3100 Vhr was reached.
Data Acquisition and Analysis
Two randomly selected samples labeled with Cy5 and Cy3 were co-resolved with a Cy2 labeled pooled internal standard on each of 24 gels. Gels were scanned within 24 hours of being run on a Molecular Dynamics Typhoon 9400. Gel images were cropped using Molecular Dynamics Image Quant 5.2 software.
Upon visual inspection of the 24 gels, proteins on three gels were observed to not be well resolved. These three gels were immediately rerun using the same Cy2 labelled standard to obtain gel images suitable for analysis. Acquired gel images were first analyzed in the DeCyder™ Differential In-gel Analysis (DIA) module of the GE HealthCare DeCyder™ 2D Software version 6.5. The DIA generated identical spot feature detection patterns on all images derived from the same gel. This ensured that the internal standard and the sample spot features had identical spot boundaries. Quantification of spot features was achieved by normalizing spot feature volumes against the internal standard.
The DIA files were imported into the Biological Variation Analysis (BVA) module to match spot feature migration patterns and normalize abundance values using the unique signal of each spot feature from the pooled internal standard. The standardized abundance was derived from the normalized spot volume, standardized against the intra gel standard. To obtain a normal distribution around zero the spot feature standardized log abundance values were used for inter-gel spot comparisons.
Multivariate analysis was performed in the DeCyder™ Extended Data Analysis (EDA) module version 1.0. Marker selection and classifier creation were performed using partial least squares for the searching and ranking of spot features and K-nearest neighbor (KNN) to evaluate the spot feature set found.
Spot features of interest were manually excised using a Gilson P1000 Pipetman from SYPRO Ruby stained preparative gels and stored in 1% acetic acid. The ART pipet tips were cut with a razor blade to increase the pore size. The gel slices were washed a total of five times; first with sterile water, secondly with 25 mM ammonium bicarbonate and finally three consecutive washes with 25 mM ammonium bicarbonate/50% acetonitrile solution with the last wash being an overnight incubation at 10°C. Gel slices were dehydrated with 100% acetonitrile before adding trypsin (Trypsin Gold, Promega) at 20 μg/ml (in 40 mM ammonium biocarbonate/10% acetonitrile solution) and incubated at 37°C overnight. Tryptic peptides were extracted from the gel slices by washing with 0.1% formic acid and then 0.1% formic acid/50% acetonitrile solutions. The tryptic peptide extracts were vacuum-dried and reconstituted with 10 μl of 5% acetonitrile and 0.1% formic acid.
Nanoflow LC of tryptic peptide samples was performed with an Agilent 1100 nanoflow LC system equipped with a C18 pre-column (Zorbax 300SB-C18, 5 μm, 5 mm × 0.3 mm, Agilent) and a C18 analytical column (Zorbax 300SB-C18, 3.5 μm, 15 cm × 75 μm, Agilent). The aqueous mobile phase (solution A) contained 5% acetonitrile and 0.1% formic acid, and the organic mobile phase (solution B) contained 95% acetonitrile and 0.1% formic acid. Samples (5-μl injected) were loaded and washed on the pre-column for 5 minutes with solution A at 50 μl/min. Peptides were then eluted off the pre-column and through the analytical column with a 50 minute profile at 300 nl/min: 1 to 30% solution B over 30 min, 40% to 95% B over 5 minutes, 95% B over 5 minutes, and re-equilibrated for 10 minutes at initial conditions. Peptides were eluted directly into a QStar XL via a nanospray ion source (Applied Biosystems). The ion source was equipped with a 50-μm inner-diameter, fuse-silica needle with a 15-μm tip (PicoTip Emitter, New Objective). Data dependant acquisition was performed with a 10 second cycle: 1-second interval for acquiring intact peptide signal (MS), and three 3-second intervals for collision induced dissociation of the 3 most intense peptide signals in the initial 1-second interval (MS/MS). The MS m/z range was 350 to 1500, and the MS/MS m/z range was 70 to 2000. Collision energy was automatically determined by the data acquisition software (Analyst QS 1.1). MS/MS data were acquired for the entire LC run.
Mascot (version 2.2, Matrix Science) search engine was used to search the NCBI database with the MS/MS data. The search parameters were as follows: taxonomy was unrestricted, protein molecular weight was unrestricted, fixed modification was Carbamidomethyl (C), variable modification was Oxidation (M), peptide and fragment mass tolerance was ± 0.3 Da, and up to one missed cleavage was allowed. Individual ion scores > 54 indicate identity or extensive homology (p < 0.05).
Gel image acquisition
Urine was collected from 4 Simmental cattle orally infected with BSE and 4 age matched controls at 8 month intervals over the first 40 months of the disease. This generated 4 infected and 4 control biological replicates at each of the 6 time points. An internal gel standard was created by pooling equivalent amounts of protein from each sample. In total 27 gels, each comprised of the internal standard and two biological samples, were run to obtain gel images suitable for analysis. Ultimately, 46 gel images representing 46 biological samples and 24 gel images of the internal standard were used for analyses (Table 1).
Table 1. Sample loading and labelling matrix.
The acquired gel images were analyzed using the DeCyder™ DIA and BVA software modules. Manual landmarking of 12 spot features across 24 gel images of the internal standard was performed in order to improve the accuracy of the gel-to-gel matching process. This resulted in the detection, quantification, and matching of 1329 master spot features across the 24 gels, all containing the same internal standard (Figure 1).
Figure 1. Representative Cy2-labelled internal standard proteome gel image illustrating proteins resolved in the pH4-7 range. The gel image as loaded into the DIA module prior to spot detection (A). In panel B the 1329 spot features, including spots at the edges of the gel that were outside the pH range of the 1st dimension separation, are each denoted by a green dot. The position of the 16 spot features used in the class prediction classifier have been marked with yellow and the associated master gel spot feature number assigned to the same spot features on all gels by the DeCyder BVA module are shown.
Multivariate analyses of protein expression data derived from the BVA were performed using the DeCyder™ Extended Data Analysis Software (EDA). The gel images were first grouped such that the 6 samples from each individual animal formed a group. This resulted in 8 groups each representing one of 8 biological replicates. The data were filtered so that only the 36 spot features exhibiting statistically significant (ANOVA p < 0.01) changes in abundance and present on all 46 gel images were considered in the following analyses.
Principle component analysis (PCA) on this filtered data set was performed to identify the relative contributions of the inherent differences between individuals and disease state on the variance exhibited by the 8 biological replicates. The PCA analysis demonstrated that the cows generally segregated into infected and control groups indicating that the disease status of the animals was the primary factor affecting the differential abundance of urinary proteins (Figure 2). The exception was cow #54. At no time point did this animal cluster with the other infected animals and it further diverged from all other animals as the disease progressed.
Figure 2. Principle component analysis of the 8 biological replicates. The samples obtained from individual infected and control cows clustered together indicating that disease is the factor that most influences the differential abundance observed in the urine samples. The time of the 6 sample collections from cow #54 are given to illustrate that the urine proteome of this animal diverged further from all the other animals as the disease progressed. This analysis was based on the 36 spot features exhibiting statistically significant (ANOVA p < 0.01) changes in abundance and present on all 46 gel images. (PC1 = 36.3, PC2 = 15.2).
The reason for the atypical pattern demonstrated by this animal is not known, but it is interesting to note that it had an atypical phenotype as well. Cow #54 developed clinical signs, as determined by regular scoring for clinical signs (20 minutes/animal), between 40–46 months post infection (mpi), but then phenotypically recovered before reaching the terminal stage of the disease at 56 mpi. In contrast, cows #38 and #53 reached the terminal stage of the disease at 44 mpi and cow #40 at 48 mpi. The behaviour of cow #54 was atypical, not only with respect to the other three infected cattle considered in this study, but from the other 10 infected cows in the herd that were allowed to reach the terminal stage of the disease. For the purposes of this experiment all the data representing this animal were excluded from subsequent analyses. BSE infection was confirmed in all four cases by immunohistochemical analyses of the obex.
The 40 remaining gel images were then grouped according to disease state and months post infection. Control and infected cows at each time point formed a group. This resulted in 10 groups each representing either the 3 remaining infected cows or the 4 control cows at a particular point in the experiment. Samples collected from these 7 cows prior to the start of the experiment formed an 11th group (normal). The data were again filtered revealing 56 spot features that exhibited statistically significant (ANOVA p < 0.01) changes in abundance and were present on all 40 gel images for consideration in the following analyses.
PCA analysis on this filtered data set was performed to identify the relative contributions of time and disease state on the variance observed (Figure 3). The cows again segregated into infected and control groups further indicating that the disease status of the animals was the primary cause of the differential abundance in urinary proteins observed. Nonetheless, within the infected group it was also observed that the individual time points clustered together and generally moved down and to the right as disease progressed. A somewhat similar, but less pronounced pattern was observed in the control samples. This indicates that time did factor into the differential abundance of urinary proteins observed. That time played a role in the urinary protein profile was reinforced by the observation that the animals at 0 mpi formed a distinct group (normal).
Figure 3. Principle component analysis of different disease states followed throughout the disease progression. Ellipses have been drawn to illustrate the clustering of the 3 groups (BSE infected, control and normal). Within the infected group it can also be seen that the individual time points cluster together. A somewhat similar but less pronounced pattern is observed in the control samples. This analysis is based on the 56 spot features exhibiting statistically significant (ANOVA p < 0.01) changes in abundance and present on all 40 gel images. (PC1 = 38.6, PC2 = 23.0).
PCA of the 56 spot features that exhibited statistically significant (ANOVA p < 0.01) changes in abundance and were present on all 40 gels revealed 3 outliers representing either mismatched features or strongly differentially abundant features (Figure 4A). Visual inspections of all potential outliers were made in the corresponding BVA file containing the gel images. The 3 spot features clearly observed in the infected image were not visible in the control image of the gel (Figure 4B). Visualization of the shape of the peak representing spot feature 405 and those in the surrounding area demonstrated that the matches were legitimate (Figure 4C). The increase in abundance of 405 ranged from 17 to 77 fold over the course of the experiment. Graphical representation of the standardized log abundance obtained for the spot feature 405 demonstrated strong differential abundance of the spot feature throughout the experiment. Note the segregation of the control and infected biological replicates at each time point (Figure 4D).
Figure 4. Differentially abundant spot features. Principle component analysis of 56 spot features that exhibited statistically significant (ANOVA p < 0.01) changes in abundance and were present on all 40 gels (A). The red rectangle on the gel image shows the region on the gels where the potential outliers were situated. Two magnified views of this region showing infected and control images (B). 3D images of spot feature 405 showing the 23.68 fold increase in abundance observed at 8 mpi (C). Graphical representation of the standardized log abundance data obtained for spot feature 405 (D).
In order to determine whether the differences suggested by the PCA analyses were sufficient to discriminate control from infected animals the marker selection and classifier creation functions of the EDA module were used. Initially this was performed on the pooled data obtained from all control, infected and normal samples. The data, as before, were filtered so that only spot features exhibiting statistically significant (ANOVA p < 0.01) changes in abundance and present on all 40 gel images were considered.
The partial least squares method was used for the searching and ranking of spot features and K-nearest neighbor (KNN) was used to evaluate the set of biomarkers found. The use of 16 biomarkers and KNN classification demonstrated that the training data set could be classified with 83.3% ± 18.3 accuracy (Table 2). Removal of the confounding normal samples led to the identification of a single spot feature (405) that could discriminate between the remaining control and infected sample sets with 100% accuracy.
Table 2. Classification Matrix
In order to evaluate the information concerning disease progression contained in the urine samples the 14 gel images of the infected samples were considered separately. The data were filtered so that only the 25 spot features exhibiting statistically significant (ANOVA p < 0.01) changes in abundance and present on all 14 gel images were considered. The classifier created using the 16 spot features identified by marker selection was able to classify the time post infection that the infected samples were collected with 85% ± 13.2 accuracy. The two misclassified samples were placed into the immediately proceeding time point (Table 3).
Table 3. Disease Progression Matrix
In order to parse out those changes due to disease progression from those associated with aging a similar analysis was performed on the control samples. The classifier created using 16 of the spot features identified by marker selection was able to classify with 85% ± 19.1 accuracy the time post infection that the control samples were collected (Table 4). Two of the three misclassified samples were placed into the immediately proceeding time point.
Table 4. Aging Matrix
The 16 spot features used by the 3 classifiers are provided in Table 5. The rank assigned to the spot features denotes the relative contribution of each protein to the classification. The standardized abundance ratio of the top three ranked proteins in the disease progression classifier were plotted with respect to time post infection (Figure 5). The steady increase or decrease in abundance observed over the course of the experiment illustrated the utility of the relative abundance of the spot features in classifying the urine samples with respect to date post infection. Three of the spot features utilized to predict disease progression (597, 1022, 1041) were also used in disease class prediction indicating that these markers are disease specific. Another disease progression spot feature (710) also appeared among the control progression markers indicating that changes in the abundance of this spot feature were age related.
Figure 5. Proteins that exhibited a steady increase or decrease in abundance throughout disease progression. The average standardized abundance ratios of the top three ranked proteins used in the disease progression classier (437, 1041, 1022) are shown. The consistent increase or decrease in abundance over the course of the experiment illustrates the utility of the relative abundance of these spot features in classifying the urine samples with respect to date post infection.
Table 5. Biomarker sets used to create classifiers.
The 16 biomarkers making up the class prediction classifier (Figure 1B), designed to discriminate between control, infected and normal samples, were excised from a preparative gel and subjected to protein identification using mass spectrometry and database interrogation as described in Materials and Methods. MS analysis identified 5 unique proteins, not including redundancies likely due to post-translational modifications or proteolysis, and enabled proteins to be assigned to 13 of the 16 spot features. A summary of the results of this analysis are shown in Table 6. See Additional File 1 for additional protein statistics.
Additional File 1. Results of LC/MS/MS analyses and protein statistics. Thirteen of the 16 spot features included in the class prediction classifier were identified. For each identified protein the mascot score, the number of peptides, the % coverage and the corresponding NCBI identifier are provided. Individual ions scores > 54 indicate identity or extensive homology (p < 0.05). The average ratios at each time point are also given.
Format: DOC Size: 58KB Download file
This file can be viewed with: Microsoft Word Viewer
Table 6. Thirteen of the 16 spot features included in the class prediction classifier were identified. The average ratios at each time point are given.
The 5 different proteins representing the 13 spot features were: clusterin, Ig Gamma-2 chain C region, similar to GCAP-11/uroguanylin, cystatin E/M, and cathelicidin1. Some of the redundancies appeared to be due to post-translational modifications that created charge-related isoforms that had different iso-electric points, but indistinguishable molecular masses. For example, 2 spot features identified as clusterin (393, 405), appeared to meet this criteria (Figure 4B). It is also interesting to note that the differential abundance of one of the isoforms is much more robust than the other (Figure 4A).
Recent advances in 2-dimensional gel electrophoresis technologies, namely the introduction of fluorescent dyes, which allow multiple samples to be co-separated and visualized on one 2-D gel, have increased the utility of this methodology for the discovery of robust protein biomarkers for disease [18-20]. For example, power analysis for 2D-DIGE has demonstrated statistical power of >0.8 for detecting 2-fold changes at p ≤ 0.01 with 4 biological replicates. In our experiment urine samples representing 4 control and 4 age matched infected cattle were run with an internal standard. The internal standard, created by pooling aliquots of all biological samples in the experiment, was run on all gels within the experiment thereby creating an intrinsic link across all gels. Normalization of the internal standard across gels allowed the ratio of relative abundance of the same protein to be compared directly, separating gel-to-gel variation from biological variation. Differences in protein abundance were then determined by comparing the ratio obtained from one fluorescent labelled sample directly with another.
Changes in the amount of a given soluble protein in urine can result from a change in its concentration in the blood plasma, a change in the function of the glomerular filter, an alteration in the proximal tubule scavenging system or a change in local production and excretion. No proteinuria indicating a change in the glomerular filter was observed in response to BSE infection (Additional file 2). In addition, no significant difference (p < 0.01) in the standardized log of abundance for cystatin was observed at any time point when compared to the corresponding control (Additional file 3). Cystatin is a low molecular weight cysteine proteinase inhibitor that is freely filtered by the renal glomeruli and reabsorbed in the proximal tubules. Consequently cystatin serves as an indicator of the health of the proximal tubule scavenging system[23,24]. The absence of change in cystatin abundance indicated that the proximal tubule scavenging system was unaffected by BSE infection. Together these results suggest that the differential abundances observed in response to BSE infection were likely due to changes in the protein concentration in the blood plasma or a disease associated alteration in local production.
Additional file 2. Statistical Analysis of Protein Concentration in the Urine Samples. Protein concentrations of urine are evaluated to determine if there was any difference amongst cows and whether or not the concentration changed throughout the course of the disease.
Format: DOC Size: 41KB Download file
This file can be viewed with: Microsoft Word Viewer
Additional file 3. Statistical Analysis of Cystatin Abundance in Control and Infected Urine samples collected at each of the 5 time points. The relative cystatin abundances found in control and infected urine were compared.
Format: DOC Size: 42KB Download file
This file can be viewed with: Microsoft Word Viewer
One of the identified classifier proteins that exhibited increased abundance in response to infection was immunoglobulin Gamma-2 chain C region. Increased immunoglobulin has previously been observed in the urine of scrapie infected hamsters and in one report it was speculated that this was the result of a nephritic syndrome [5,6,16]. In contrast, the data presented suggests that the increased abundance of immunoglobulin in urine associated with BSE was probably due to a change in immunoglobulin concentration in the blood plasma. Nonetheless, the detection of differential abundance of another immunoglobulin protein in urine by an unbiased screen lends support to the suggestion that immunoglobulin light chain may constitute a surrogate marker for TSE diseases[6,16].
In addition to immunoglobulin Gamma-2 chain C region and cystatin, the 3 other identified classifier proteins were well known urinary proteins. One was the antimicrobial peptide cathelicidin that is produced in the kidney by the epithelial cells that line the urinary tract. When exposed to bacteria the levels of cathelcidin mRNA are known to rapidly increase, however, even in the absence of microbes the epithelial surface of the urinary tract are continuously bathed with cathelicidin[25,26]. A second protein was the natriuretic peptide uroguanylin that is produced in the small intestine and kidney. In response to salt loading no increase in circulating uroguanylin is observed indicating that the natriuretic effect of uroguanylin is in part mediated by increased renal production that inhibits tubular resorption of ions from the glomerular filtrate. These characteristics suggested that the decreased abundance of both cathelicidin, observed after 24 mpi, and uroguanylin throughout the experiment were probably the result of decreased renal production and excretion. The precise cause of the decreased production and the possible effects of the altered abundance on prion pathobiology are not known.
The third protein, identified as clusterin, was able to distinguish between infected animals and age matched controls with 100% accuracy throughout the experiment. Clusterin is a multifunctional glycoprotein with nearly ubiquitous tissue distribution. Increased abundance of clusterin in association with TSE diseases has been reported previously and has included increased expression in astrocytes as well as a significant accumulation in cerebrospinal fluid and blood plasma. Thus, increased amounts of circulating clusterin may have caused the increased clusterin abundance observed in the urine of the BSE infected animals. Despite the power of clusterin as a biomarker of BSE in this experiment, the increased clusterin abundance in CSF observed in models of other neurodegenerative diseases, such as Alzheimer's disease, and in response to a variety of renal insults raises doubts as to the specificity of clusterin per se as a biomarker of BSE in cattle [30-34]. However, the specificity of the particular isoform of clusterin observed to best discriminate between BSE infected and control cattle remains to be seen.
EDA analyses also demonstrated that the differential abundance of different subsets of proteins provided accurate measures of disease progression and aging. This was an unexpected result, but the ability to follow disease progression by monitoring the differential abundance of a subset of proteins has potential applications as a prognostic indicator or in the assessment of the therapeutic benefit of potential treatments. Furthermore, while markers of disease progression must be sensitive to changes in disease state and present in easily accessible tissues that permit repeated sampling, they do not require the same high disease specificity as diagnostic markers. The ultimate utility of these markers of disease progression will be determined by their identification and validation as well as their applicability to clinically relevant disease models. Significantly, the markers of disease progression demonstrated very little overlap with those identified as able to track age. This indicates that they were a measure of disease specific progression and that their identification may also provide insight into the pathology of these diseases.
The results demonstrate that in principle it is possible to identify biomarkers of TSE disease by analyzing changes in the urine protein profile provoked by the disease. Extending the present study to larger numbers of cattle and to those of other strains will test the value of the biomarkers identified. Even more promising markers may have been missed due to the bias of 2D-DIGE to the identification of abundant proteins. This shortfall will be addressed in future studies by utilizing a variety of pre-fractionation methods.
JDK, MG, SC, CG have made substantial contributions to conception and design of the experiment. SLRS, LL, MP, JL, UZ, MS performed the acquisition of data. JDK carried out the analysis and interpretation of data. JDK, SLRS, MP, LL, SC have been involved in drafting the manuscript or revising it critically for important intellectual content. All authors read and approved the final manuscript.
We would like to thank Michael Carpenter and Julie Boutillier for their advice and many helpful discussions. We would also like to thank Garrett Westmacott and Keding Cheng of the Proteomics Core Facility, Public Health Agency of Canada, for the Mass Spectrometry Analyses. This work was funded by the Alberta Prion Research Institute and the Public Health Agency of Canada.
Folia Neuropathol 2004, 42 Suppl A:85-91. PubMed Abstract
Furukawa H, Doh-Ura K, Okuwaki R, Shirabe S, Yamamoto K, Udono H, Ito T, Katamine S, Niwa M: A pitfall in diagnosis of human prion diseases using detection of protease-resistant prion protein in urine. Contamination with bacterial outer membrane proteins.
Hoffmann C, Ziegler U, Buschmann A, Weber A, Kupfer L, Oelschlegel A, Hammerschmidt B, Groschup MH: Prions spread via the autonomic nervous system from the gut to the central nervous system in cattle incubating bovine spongiform encephalopathy.
Alban A, David SO, Bjorkesten L, Andersson C, Sloge E, Lewis S, Currie I: A novel experimental design for comparative two-dimensional gel analysis: two-dimensional difference gel electrophoresis incorporating a pooled internal standard.
Methods Mol Biol 2007, 367:219-239. PubMed Abstract
Chromek M, Slamova Z, Bergman P, Kovacs L, Podracka L, Ehren I, Hokfelt T, Gudmundsson GH, Gallo RL, Agerberth B, Brauner A: The antimicrobial peptide cathelicidin protects the urinary tract against invasive bacterial infection.
Ghiggeri GM, Bruschi M, Candiano G, Rastaldi MP, Scolari F, Passerini P, Musante L, Pertica N, Caridi G, Ferrario F, Perfumo F, Ponticelli C: Depletion of clusterin in renal diseases causing nephrotic syndrome.