3.1 Assessment of gene expression normality for Mock and Utah-inoculated clones
Seven genes with a potential role in prion permissibility were assessed with RT-qPCR,
and normalized to 18s rRNA and hPRT1b. Gene expression data from Mock-inoculated clones showed a normal distribution (Goodness-of-Fit; p > 0.05). Utah-inoculated clones were not normally distributed for SDC4 (Goodness-of-Fit; p < 0.001) and SPARC (Goodness-of-Fit; p < 0.01) but were normally distributed for the remaining genes. For SDC4 and SPARC, raw expression values from Utah-inoculated clones were transformed with Box-Cox Y transformation test, and these values were deemed as normally distributed after assessment with histograms and Goodness of Fit test (p > 0.05) and subsequently utilized in parametric statistical analyses as such.
3.2 Transcript levels in Mock-inoculated clones
To determine if the RNA levels of the targeted genes significantly differed between clones, fold changes in normalized gene expression were statistically compared amongst all of the clones from at least 3 independent replicates. Statistically significant differences in expression were observed for BIRC5 (p < 0.001), SPARC (p < 0.001) and FSTL1 (p < 0.01) (Figure 1). For BIRC5, differential gene expression was significantly decreased for both 439Late
and 439Early compared to all other clones (p < 0.05), and this decreased expression was a four-to-eight-fold change compared to clone 438. Expression of FSTL1 in clones 439Late (p < 0.05) and 439Early (p < 0.05) was significantly decreased compared to clone 434, and clone 439Late was also significantly decreased compared to clone 438 (p < 0.05), but these changes were less than two-fold compared to clone 438. Additionally, clone 440 had significantly decreased expression of FSTL1 compared to clones 438 (p < 0.05) and 434 (p < 0.05), but these changes were also less than two-fold compared to clone 438. With SPARC, clones 439Late and 439Early had significantly diminished expression compared to clones 438 (p < 0.0001 and p < 0.01) and 434 (p < 0.001 and p < 0.05), and these changes were two-to-three-fold compared to clone 438. Additionally, 439Late SPARC expression was significantly decreased compared to clone 441 (p = 0.0378). Furthermore, clone 438 had significantly increased SPARC expression compared to clones 441 (p < 0.05) and 440 (p < 0.01), and clone 434 was also significantly increased compared to 440 (p < 0.05), but these changes were less than two-fold compared to clone 438. Significant differential gene expression was lacking for AXL, FN1, PRNP and SDC4 in all pairwise comparisons.
Since clone 439 exhibited variation in permissibility over time, differences in gene expression were also assessed for 439Late relative to the most permissive clone, 439Early in all evaluated genes. Only in BIRC5 was there a two-fold change expression in 439Late relative to 439Early, but this change was not statistically significant (p > 0.05).
There were no statistically significant correlations between mean normalized gene expression of Mock-inoculated clones and prion permissibility (p > 0.05).
3.3 Transcript levels in Utah-inoculated clones
As described in section 3.2, analysis of fold changes in gene expression and correlations with prion permissibility were also performed for Utah-inoculated clones for at least 3 independent replicates, with the exception of clone 439Early for BIRC5 in which data from only 2 out of 3 independent replicates was available to assess (one replicate had undetectable transcript levels, consistent with low expression levels detected in the other replicates). Individual, significant variations in gene expression were present for PRNP (p < 0.05) and SPARC (p < 0.01) (Figure 2). A single, significant increase in PRNP expression was
observed for clone 439Late as compared to clone 438 (p < 0.05); however, this change was less than two-fold compared to clone 438. With SPARC, expression levels were significantly increased for clone 439Late as opposed to clone 438 (p < 0.01), and clone 434 (p < 0.01), but this variation was also less than two-fold compared to clone 438. Although a two-fold change in BIRC5 expression was present for 439Late and 439Early relative to clone 438, these differences were not statistically significant (p > 0.05) due to the low expression levels in these clones. BIRC5 tended to amplify late in RT-qPCR with Utah-inoculated clones. Significant differential gene expression was also lacking for AXL, FN1, FSTL1, and SDC4 in all pairwise comparisons. No significant differential expression existed for clones 439Late versus 439Early in any of the examined genes.
As also described for the Mock-inoculated clones, significant correlations were not observed between mean normalized gene expression and prion permissibility (p > 0.05).
3.4 Transcript levels of Mock and Utah-inoculated clones combined
The data for the Mock and Utah-inoculated clones were combined to increase the sample sizes. For FSTL1, SPARC, and SDC4, the effect of inoculation was significant (p < 0.05), and the interaction of inoculation status and clone were also significant for SDC4 and SPARC (p < 0.05), so these genes were not included in this aspect of the study. However, data could be combined for AXL, FN1, BIRC5, and PRNP. For BIRC5 (p < 0.001) and PRNP (p < 0.01), the effect of clone remained statistically significant with two-way ANOVA, frequently involving the same original clones as described in sections 3.2 and 3.3, but with some modifications to be addressed. When data was combined for BIRC5, the significant variation in gene expression of clones 439Early and 439Late, compared to clone 440 were no longer significant (p > 0.05). In PRNP, although clone 439Early lacked a two-fold change in expression to other clones, collective data indicated that 439Early is weakly but significantly (p < 0.05) increased relative to 438. No significant correlations were identified between genes (p > 0.05) or with prion permissibility (p > 0.05).
Since inoculation status was influential to gene expression for FSTL1, SPARC, and SDC4 with two-way ANOVA, the ratio of normalized mean expression derived from both treatment groups (i.e., Mock-inoculated/Utah-inoculated) was constructed to evaluate if inoculation with scrapie prions influenced transcript levels of these three genes (Figure 3).
For SPARC, transcript levels were approximately three-fold higher for clones 438 and 434 (p < 0.05) as compared to clone 439Late subsequent to inoculation with scrapie prions.
3.5 Gene expression values relative to PRNP in Mock-inoculated clones
Although normalized PRNP expression was determined to lack two-fold changes in expression and displayed absent to minimal, significant variation amongst the clones in this, as well as another previous study [28], the influence of target genes relative to PRNP expression of other genes was evaluated to determine if target gene/PRNP ratios were significant relative to prion permissibility. Raw target gene/PRNP ratios were calculated using the relative quantity of expression (i.e., non-normalized) for target genes and PRNP, per clone. Fold changes in expression and correlations with prion permissibility were statistically evaluated as previously described in section 3.2. Statistically significant differences in expression were observed for BIRC5/PRNP (p < 0.001), SPARC/PRNP (p < 0.001), FSTL1/PRNP (p < 0.01) and SDC4/PRNP (p < 0.01) (Figure 4). In BIRC5/PRNP, differential gene expression of 439Late and 439Early were significantly decreased compared to clones 438 (p < 0.001 and p < 0.01), and 441 (p < 0.01 and p < 0.05), and this decreased expression was a five-to-eight-fold change for 439Late and 439Early compared to clone 438. Earlier in the study, significant fold changes with BIRC5 in 439Early and 439Late were also described for both treatment groups, independent of PRNP expression. Moreover, for BIRC5/PRNP, clone 439Late had significantly diminished expression compared to clones 434 (p < 0.01) and 440 (p < 0.05). SPARC/PRNP had significant decreased expression for clones 439Late and 439Early compared to clones 438 (p < 0.001 and p < 0.01) and 434 (p < 0.001 and p < 0.05), and a two-to-three-fold change in expression existed in clones 439Late and 439Early compared to clone 438. Also, for SPARC/PRNP, clones 440 and 441 had significantly diminished expression compared to 438 (p < 0.01 and p < 0.05), and clone 440 was significantly decreased from 434 (p < 0.05), but these changes were slightly less than two-fold relative to clone 438. In FSTL1, clone 439Late had significantly diminished expression from clones 438 (p < 0.05) and 434 (p < 0.05), but these differences were less than two-fold. For SDC4/PRNP expression, 439Late was significantly decreased compared to clones 438 (p < 0.01), 440 (p < 0.05), and 434 (p < 0.05), but these differences were also less than two-fold relative to clone 438. Relative to PRNP expression, significant differential gene expression remained absent for AXL and FN1.
Although an approximately three-fold change in BIRC5/PRNP expression was present in clone 439Late compared to 439Early, this difference was not significant (p > 0.05). No significant differences in expression of clones 439Early versus 439Late were present for any of the remaining gene ratios.
No significant correlations to prion permissibility were identified amongst all gene ratios.
3.6 Gene expression values relative to PRNP in Utah-inoculated clones
Fold change data and correlations to prion permissibility were constructed and statistically assessed as previously described in section 3.5. Statistically significant differences in expression were observed for AXL/PRNP (p < 0.05) (Figure 5). For AXL/PRNP, clone 441 had significantly diminished expression compared to clone 438 (p < 0.05), but this change was less than two-fold. Clones of BIRC5/PRNP lacked significant differential expression, but an
approximately three-to-six-fold change in expression was noted in clones 439Late and 439Early relative to clone 438. In FN1/PRNP, there was also an absence of significant differential expression amongst the clones, but two-fold changes in expression were present in clones 439Late, 439Early, 434, and 441, relative to clone 438, and mean expression was highest in the least permissive clone (i.e., 438). In FSTL/PRNP, no significant differential expression was observed amongst the clones, but a two-fold change was present in clone 439Late relative to clone 438. Significant differential gene expression was lacking for SDC4/PRNP and SPARC/PRNP.
For BIRC5/PRNP, clone 439Late had a two-fold higher transcript level relative to 439Early, but this change lacked statistical significance. No additional significant differences in expression were present in clone 439Early versus 439Late for any of the remaining gene ratios.
No significant correlations with prion permissibility were detected amongst the gene ratios.
3.7 Gene expression values relative to PRNP in Mock and Utah-inoculated clones
Target gene/PRNP expression ratios were pooled from Mock and Utah-inoculated clones for AXL/PRNP, BIRC5/PRNP and FN1/PRNP, and the impact of clone and inoculation status with expression, as well as correlations to prion permissibility were evaluated using similar procedures as described in section 3.4 (Figure 6). Neither inoculation status, nor the
interaction of clone with inoculation status were influential to expression of genes selected for this portion of the study, so treatment groups were pooled to increase sample sizes. For AXL/PRNP, the effect of clone was significant (p < 0.01) for clones 438 and 441. In BIRC5/PRNP, although the effect of clone remained significant to gene expression (p < 0.001) with many of the same clone comparisons depicted in section 3.5, further significant differential comparisons included clones 439Early and 434 (p < 0.01), clones 441 and 439Late (p < 0.01), and clones 440 and 438 (p < 0.05). Furthermore, BIRC5/PRNP differential expression became insignificant between clones 440 and 439Late (p > 0.05). Of interest, when FN1/PRNP was assessed jointly for both treatment groups, the effect of clone became significant (p < 0.01), and expression of clone 438 was significantly increased compared to 439Late (p < 0.05), 441 (p < 0.01) and 439Early (p < 0.05), with clones 441 and 439Early previously characterized as highly permissive to scrapie prion infection.
Additionally, BIRC5/PRNP (r = -0.7187, adjusted p-value < 0.05), and FN1/PRNP (r = -0.7841, adjusted p-value < 0.01) were strongly and negatively correlated with prion permissibility (Table 1).
Table 1. Significant correlations between target genes/PRNP and prion permissibility.

a Data was derived from at least 3 independent replicates. A p-value < 0.05 was considered significant. Ratio expression was scaled to 438.
b Adjusted p-value were derived from the formula (Pearson’s p-value x (total number of p-values/p-value rank)). The value obtained from this calculation was compared to the previous adjusted p-value, and the smaller of the two was recorded. This value is also referred to as the Benjamini-Hochburg p-value
(www.biostathandbook.com/multiplecomparisons.html).