Previous Article | Next Article ![]()
Journal of Virology, September 2009, p. 9464-9473, Vol. 83, No. 18
0022-538X/09/$08.00+0 doi:10.1128/JVI.00352-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Molecular Pathogenesis and Genetics,1 Department of Pathology, Veterinary Laboratories Agency, Woodham Lane, New Haw, Surrey KT15 3NB, United Kingdom,2 Institute of Neuropathology, Ludwig-Maximilians-University Munich, Munich, Germany3
Received 17 February 2009/ Accepted 29 June 2009
|
|
|---|
|
|
|---|
BSE was first reported in the United Kingdom in 1986, and its most likely cause was contaminated meat and bone meal, a dietary supplement for cattle (33, 34). The resulting BSE epidemic reached its maximum in 1992, and to date, more than 185,000 cattle succumbed to BSE in the United Kingdom. BSE is also considered to be the origin of the human disease variant CJD (15, 31). Worldwide, there have been almost 200 cases of variant CJD, and most have been found in the United Kingdom.
In recent years, few cases of unusual forms of BSE (bovine amyloidotic spongiform encephalopathy, lower-molecular-weight BSE, and higher-molecular-weight BSE [known as BASE, L-type BSE, and H-type BSE, respectively]) were identified mainly in Europe (4, 6, 12). The molecular signature of PrPSc of these cases differs from that of classical BSE, and the distribution of PrPSc in the brains of the affected animals also differ where data are available. These unusual (or atypical) forms of BSE mainly affect older cattle, and the etiology is currently unclear.
Many aspects of the pathogenesis of classical BSE have been studied extensively following oral infection of cattle (32). Calves orally exposed to BSE by ingestion of 100-g portions of homogenized brain stem samples from cattle with BSE developed BSE with a mean incubation time of 45 months (between 33 and 55 months). The earliest detection of infectivity in the central nervous system was 32 months postinoculation (34). The correlation between BSE pathogenesis and the transcriptional activities of genes is however not explored in any detail. Gene expression profiling studies of mice with scrapie have shown that these studies can provide valuable insights in the possible pathomechanisms of the disease (5, 25, 27, 29, 36, 37).
To determine the correlation between BSE pathogenesis and gene expression, we analyzed brain stem samples from cattle experimentally infected with the BSE agent over the time course of the infection by using microarray-based gene expression profiling.
|
|
|---|
Microarray analysis. The preparation of samples and reagents used were as specified by the Affymetrix GeneChip Expression Analysis manual (1). Briefly, total RNA was extracted from 0.2-g brain stem tissue samples by the Trizol method (Invitrogen), cleaned by RNeasy columns (Qiagen), and checked for integrity by gel electrophoresis. Total RNA (5 µg) was used as a template for T7-oligo(dT)-primed reverse transcription to cDNA according to the Affymetrix manual (1). The double-stranded cDNA was cleaned using the GeneChip sample cleanup module (Affymetrix). Synthesis of biotin-labeled cRNA from the cDNA was performed by Affymetrix one-cycle target labeling assays. The cRNA was cleaned with the GeneChip sample cleanup module and quantified according to the Affymetrix manual (1). The biotin-labeled cRNA (20 µg) was fragmented in fragmentation buffer at 94°C for 35 min (Affymetrix). The quality of the fragmented cRNA was examined by electrophoresis using 1% agarose gel (Sigma-Aldrich). The fragmented cRNA (15 µg) was used to prepare the hybridization cocktail (Affymetrix 49 format). To ensure that the RNA samples were of high quality, Test 3 probe arrays (Affymetrix) were used to evaluate seven RNA samples, and all of them showed satisfactory quality (3'/5' ratios of the probe sets for glyceraldehyde-3-phosphate dehydrogenase [GAPDH] < 3). The cRNA cocktails were hybridized to Affymetrix bovine genome GeneChips, which contain 24,128 probe sets, at 60 rpm for 16 h at 45°C in Affymetrix hybridization oven 640. The microarrays were washed and stained in Affymetrix fluidics station 450DX, and the data from the microarrays were collected by an Affymetrix scanner.
The raw data were first imported to Affymetrix GeneChip operating software version 1.4. After initial analysis of the data and inspection of the images, the pivot formatted data were analyzed further with GeneSpring version 7 software (Silicon Genetics). The data were normalized by three steps: (i) transform data transformation set measurements less than 0.01 to 0.01, (ii) normalize each chip to the 50th percentile, and (iii) normalize each gene to the median. These steps are the default settings for the GeneSpring package.
One-way analysis of variance (ANOVA) statistical analysis was used as a filter tool that compares mean expression levels between two or more groups of samples. The comparison is performed for each gene, and the genes with the most significant differential expression (smallest P value) are returned. The parameters for the analysis were as follows. (i) The P value cutoff was 0.05. (ii) The parametric test did not assume that variances were equal. (iii) Multiple testing correction was used. (iv) The post hoc tests were Student-Newman-Keuls. The post hoc tests were carried out in conjunction with ANOVA to determine which specific group pair(s) was statistically different.
Quantitative reverse transcription (RT)-PCR. The RNA samples were treated with DNA-free DNase treatment and removal reagents (Ambion) for 1 h at 37°C to remove any trace DNA. The treated RNA was then used as a template for cDNA synthesis with a TaqMan reverse transcription kit (Applied Biosystems). The real-time PCR was carried out by denaturing at 95°C for 15 s, annealing at 50°C for 2 min, and extension at 60°C for 1 min for 40 cycles using ABI Prism 7700 sequencing detector. The GAPDH gene was used as an internal control to normalize the expression levels of target mRNA.
The primer sets were chosen by using Primer Express 1.5 for TaqMan software. The sequences of the primer sets and probes follow. For Arg2 (arginase, type II), the primers were 5'-GGC AGT GGA CGT CAT TGC T and 5'-GGT CGT AGA CAA TAT GCC CTC C, and the probe was 5'-FAM-CGA GTT TCG GGC AGA CGA GGG A (FAM is 6-carboxyfluorescein). For ITGB5 (integrin, beta 5), the primers were 5'-TGC TCG TCA CCA TCC ACG and 5'-GGCCCT GGA TCG CTC ACT, and the probe was 5'-FAM-CCG GAG AGA GTT CGC CAA GTT CCA. For ITGA4 (integrin, alpha 4), the primers were 5'-TGC TTT CCT TCA TTT CTT ATG TTA TGT G and 5'-TCC AAC TGT CTC TTC TGT TTT CTT TTT, and the probe was 5'-FAM-AAG GCT GGC TTC TTC AAA AGA CAG TAC CAA TC. For BOLA-DMA (major histocompatibility complex, class II, DM alpha-chain), the primers were 5'-TGT GCG GCG TGG CC and 5'-AGA GGA CCA AGC CAA CAA TGA, and the probe was 5'-FAM-TTG GCC TGG GTG TGC TGG GC. For COL9A1 (collagen, type IX, alpha 1), the primers were 5'-GCC ACT GGG AAT CGA ACA AG and 5'-GGG CTG GAT GGA AGT CTC C, and the probe was 5'-FAM-AGG GCA AAT TCG AAA AGG CTG CAG TT. For GSTA2 (glutathione S-transferase, alpha 2), the primers were 5'-GAC CCT AGC CTT TTG GCC A and 5'-CCG GGA GAC TGC TGA CTC TG, and the probe was 5'-FAM-TTC CCT CTG CTG AAG GGC CTG AAA G. For GAPDH, the primers were 5'-TCA GCA ATG CCT CCT GCA C and 5'-CAG TCT TCT GGG TGG CAG TGA, and the probe was 5'-VIC-CCC CTG GCC AAG GTC ATC CAT.
|
|
|---|
Clustering analysis of the time course samples. To identify any relationship between the changes in gene expression of the various samples, clustering analysis was performed using the GeneSpring package to group individual samples according to the similarities in expression patterns using the condition tree with the Spearman correlation.
The results showed that the samples were divided into two major groups (Fig. 1). One group combined all the negative samples and all 6 mpi samples plus one 36 mpi sample. The other group was divided into two subgroups: one subgroup combined all 21 mpi samples with one 27 mpi sample and one 39 mpi sample; the other subgroup combined all positive samples together with two 27 mpi samples, two 39 mpi samples, and one 36 mpi sample. This clustering analysis indicated that there was a broad correlation between disease progression and expression patterns. It is also noteworthy that the highest degree of difference was between 21 mpi samples and the negative controls, a finding that was confirmed in further analyses (see below).
![]() View larger version (9K): [in a new window] |
FIG. 1. Condition tree of clustering analysis for BSE time course samples. The analysis was performed by GeneSpring using all 24,128 genes (probe sets) on Affymetrix bovine microarray GeneChips. Similarity was measured using the Spearman correlation with a value of 1 for the separation ratio and a value if 0.001 for the minimum distance to merge similar branches. Three positive-control animals (Positive-1, -2, and -3) and three negative-control animals (Negative-1, -2, and -3) were used in the study. 21m-2, sample from 21 mpi from animal 2.
|
|
View this table: [in a new window] |
TABLE 1. Relative levels of differentially expressed genes during the progression of BSE
|
![]() View larger version (10K): [in a new window] |
FIG. 2. Gene expression profiles of T-cell receptor gamma variable 3-1 (a), proteasome 26S subunit (b), a gene similar to 14-3-3 protein theta (c), a gene similar to metalloprotease 1 (d), nuclear receptor (NR1H3) (e), T-cell receptor delta chain variable region (f), acetylcholine receptor (nicotinic, beta 4) (g), and a gene similar to glycine dehydrogenase during the progression of BSE (h). neg, negative control; pos, positive control; 6m, 6 mpi.
|
![]() View larger version (22K): [in a new window] |
FIG. 3. Relative levels of expression of Arg2, BOLA-DMA, IGTA4, IGTB5, COL9A1, and GSTA2 during the progression of BSE. The graphs of quantitative RT-PCR and microarrays are shown side by side for comparison. The time points are 6 mpi (6m), 21 mpi (21m), 27 mpi (27m), 36 mpi (36m), and 39 mpi (39m). neg, negative control; pos, positive control.
|
|
View this table: [in a new window] |
TABLE 2. Relative levels of expression for genes previously associated with prion diseases
|
|
View this table: [in a new window] |
TABLE 3. Number of genes that pass ANOVA analysis and change (fold) tests at each time pointa
|
Differentially regulated genes correlated well with the time course of BSE. We next carried out clustering analysis to investigate whether those 205 probe sets from the ANOVA analysis reflect the progression of BSE. For this purpose, individual samples of each time point postincubation were grouped together. Figure 4a showed that according to gene expression profiles of the 205 probe sets, there was a good reflection of the BSE time course except for the samples at 36 mpi.
![]() View larger version (95K): [in a new window] |
FIG. 4. Condition trees of clustering analysis. The analysis was performed by GeneSpring using 205 genes found by ANOVA analysis. The individual samples were grouped by time point. Similarity was measured using the standard correlation with a value of 1 for the separation ratio and a value of 0.001 for the minimum distance to merge similar branches. Each colored bar represents a gene, and the color represents the level of expression. The relative levels of expression are displayed in different colors: red, 5; orange, 2; yellow, 1; dark yellow, 0.7; dark blue, 0.4; and blue, 0.1. (a) Samples taken from negative controls, samples taken from animals at 6 mpi, 21 mpi, 27 mpi, 36 mpi, and 39 mpi, and samples taken from positive controls. (b) Samples taken from negative controls and positive controls and samples taken from animals at 6 mpi and 45 mpi.
|
|
|
|---|
ANOVA analysis defined 205 significantly changed probe sets, of which 114 genes with a known function could be identified. The association of seven of those genes with prion diseases had been described previously. The expression of three T-cell receptors was found to have significantly changed. This may be due to T cells infiltrating the brain of an animal infected with the BSE agent in the same way as occurs in murine and human transmissible spongiform encephalopathies (19).
Also, in the immune response group, the expression of the MHC-II gene was upregulated with an expression maximum at 39 mpi. This may indicate an immune response in the brain either through microglia activation or T-cell infiltration (14, 19). In the brains of CJD patients, the increased expression of MHC-II correlates well with neuronal apoptosis (14).
The microarray data also showed that there was a 12-fold reduction in expression of acetylcholine receptor (nicotinic, beta 4) at 6 mpi and that expression remained low throughout the time course of BSE. The regulation of acetylcholine receptor is linked to many neurodegenerative diseases. In Parkinson's disease and Alzheimer's disease, one of the pathogenic features is the loss of several subunits of the nicotinic acetylcholine receptor (13).
Furthermore, the expression of NR1H3 was reduced from the 6 mpi level. At 36 mpi, the expression of NR1H3 reached its lowest level (–4.7-fold of the negative controls). This might be linked to the downregulation of cholesterol synthesis similar to the findings in mice with scrapie (36). The nuclear oxysterol receptors liver X receptor-alpha [LXRalpha (NR1H3)] and LXRbeta (NR1H2) regulate genes involved in cholesterol homeostasis in a coordinate manner. Perturbations of cholesterol metabolism are associated with the development of other neurodegenerative diseases (21).
There were not many genes used both in this study and in the study by Sawiris and colleagues using VM mice infected with the BSE agent (27). This may be due to three fundamental differences between these two studies. First, Sawiris and colleagues used mice inoculated intraperitoneally with the mouse-adapted BSE strain 301V, while in this study cattle were orally infected with BSE. Second, in the mouse system, the whole brain had to be used for the expression analysis, while in this study the samples were from the brain stem. The third major difference is that in the mouse model only the negative and positive samples at the time point of terminal disease were analyzed, while in this study the progression of BSE was analyzed by using time course samples.
In this study we found that at 21 mpi, more genes passed both the ANOVA and change (fold) filters than at any other time point. Many of them are known to be associated with prion diseases. This may explain why at 24 to 25 mpi some cattle have shown some early clinical signs of BSE even though there is no detectable PrPSc in the brain at this stage of the disease, and the earliest detection of infectivity of BSE in the central nervous system starts at 32 mpi (3, 32, 34). Our results indicate that the global changes of gene expression activities are prior to the PrPSc accumulation and the appearance of clinical signs. Some of the genes identified at 21 mpi may therefore be responsible for the subsequent pathological and clinical events. While there is no detectable pathology present at this time point, microglia activation might already occur as it coincides with the earliest changes in neuronal morphology (11, 35). Furthermore, in a recent systems biology study based on gene expression changes in mice with scrapie, mRNA levels of C1qa/b and C3ar1 reflecting activation of microglia and astrocytes are among the first detectable expression changes (10 weeks postinfection) and precede, e.g., clinical signs by many weeks (16). Not only are our results consistent with early changes of gene expression in mice infected with the TSE agent (16, 36), they are also consistent with the behavioral changes in the early stages of prion infection in mice (9, 10, 20). Prion-infected mice have shown significant inability to discriminate a novel object from 7 week postinoculation compared with healthy mice (20). In this study mice succumb to RML prion at around 13 weeks postinoculation. Early cognitive deficits are therefore preceding the clinically recognizable prion disease in this model. The authors reason that an unidentified transient neurotoxic species is generated within neurons when PrPC is converted to PrPSc, which rapidly impairs neuronal function and synaptic responses. Recently, Collinge and Clarke (8) have proposed that it is PrPL that is responsible for the neurotoxicity instead of PrPSc. PrPL is a templated intermediate or side product during the accumulation of PrPSc. It is possible that the early changes in gene expression in this study are the responses to neurotoxicity produced by PrPL. Other studies have shown that soluble, low-molecular-weight oligomers of the full-length prion protein (PrP) are neurotoxic both in vitro and in vivo (28). The pronounced changes in gene expression at 21 mpi might therefore be a consequence of the emergence of low-molecular-weight oligomers. These oligomers are probably very difficult to detect with current assay technology, but sensitive methods, such as protein misfolding cyclic amplification might be able to detect the presence of these oligomers in the future.
The differentially regulated genes identified by the ANOVA analysis correlated well with the development of BSE. Therefore, the profiles of those genes might be able to tell the infectivity status of any given sample as shown for the samples at 45 mpi that were still clinically and pathologically negative. Their expression profiles were very similar to those of the positive controls. It will be interesting to know whether the results from this study can be applied to field cases of BSE, which are far more complex than the samples from the time course experiments regarding genotype, breed, feed, age, and environment.
In summary, there are two major findings of this study. First, at 21 mpi, there were more changes in gene expression compared to the negative controls than at any other period during the time course of BSE. Second, the 205 probe sets found by ANOVA might be used to predict the incubation period using clustering analysis.
This work was supported by a development grant from the Veterinary Laboratories Agency and a grant from the European Network of Excellence Neuroprion (FOOD-CT-2004-506579).
Published ahead of print on 8 July 2009. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»