Original Research ARTICLE
Goats without Prion Protein Display Enhanced Proinflammatory Pulmonary Signaling and Extracellular Matrix Remodeling upon Systemic Lipopolysaccharide Challenge
- 1Faculty of Veterinary Medicine, Department of Production Animal Clinical Sciences, Norwegian University of Life Sciences, Sandnes, Norway
- 2Faculty of Veterinary Medicine, Department of Basic Sciences and Aquatic Medicine, Norwegian University of Life Sciences, Oslo, Norway
A naturally occurring mutation in the PRNP gene of Norwegian dairy goats terminates synthesis of the cellular prion protein (PrPC), rendering homozygous goats (PRNPTer/Ter) devoid of the protein. Although PrPC has been extensively studied, particularly in the central nervous system, the biological role of PrPC remains incompletely understood. Here, we examined whether loss of PrPC affects the initial stage of lipopolysaccharide (LPS)-induced acute lung injury (ALI). Acute pulmonary inflammation was induced by intravenous injection of LPS (Escherichia coli O26:B6) in 16 goats (8 PRNPTer/Ter and 8 PRNP+/+). A control group of 10 goats (5 PRNPTer/Ter and 5 PRNP+/+) received sterile saline. Systemic LPS challenge induced sepsis-like clinical signs including tachypnea and respiratory distress. Microscopic examination of lungs revealed multifocal areas with alveolar hemorrhages, edema, neutrophil infiltration, and higher numbers of alveolar macrophages, with no significant differences between PRNP genotypes. A total of 432 (PRNP+/+) and 596 (PRNPTer/Ter) genes were differentially expressed compared with the saline control of the matching genotype. When assigned to gene ontology categories, biological processes involved in remodeling of the extracellular matrix (ECM), were exclusively enriched in PrPC-deficient goats. These genes included a range of collagen-encoding genes, and proteases such as matrix metalloproteinases (MMP1, MMP2, MMP14, ADAM15) and cathepsins. Several proinflammatory upstream regulators (TNF-α, interleukin-1β, IFN-γ) showed increased activation scores in goats devoid of PrPC. In conclusion, LPS challenge induced marked alterations in the lung tissue transcriptome that corresponded with histopathological and clinical findings in both genotypes. The increased activation of upstream inflammatory regulators and enrichment of ECM components could reflect increased inflammation in the absence of PrPC. Further studies are required to elucidate whether these alterations may affect the later reparative phase of ALI.
The cellular prion protein (PrPC) is an evolutionarily conserved protein abundantly expressed in the central nervous system, and also at moderate levels peripherally, such as in the immune system and lungs (1–4). Misfolding of PrPC is the essential step for the development of prion disorders in humans and animals (5). Although intensively studied for several decades, the normal functions of PrPC have proven difficult to identify. Not surprisingly, the main focus has been on the nervous system, and PrPC has been linked to neuroprotection, signal transduction, and myelin maintenance, among others (6, 7). More recently, the association of PrPC with inflammatory responses and the immune system has been investigated. Lack of PrPC enhances the inflammatory reaction in several models and organ systems (8–11), and it has been proposed that PrPC plays important roles in neuroimmune cross-talk (12, 13).
In a line of Norwegian dairy goats, a mutation in the PrPC-encoding gene (PRNP) terminates PrPC synthesis. Animals homozygous for this mutation (PRNPTer/Ter) are the only known non–transgenic mammals devoid of PrPC (14). The goats are postulated to be scrapie-resistant, and constitute a unique animal model for studying PrPC physiology. No major loss-of-function phenotypes have been observed under normal herd conditions, but a hematological shift in red blood cells (15) and a primed state of type I interferon-stimulated genes are present at rest in these animals (16). When exposed to systemic lipopolysaccharide (LPS), PrPC-deficient goats display a prolonged sickness behavior and a slightly different immunological profile that is skewed toward a type I interferon response (17).
Acute lung injury (ALI) and its most severe manifestation, acute respiratory distress syndrome, contribute to increased morbidity and mortality in critically ill patients (18, 19). The most common cause of ALI is sepsis caused by a systemic bacterial infection (20). Early events in this syndrome can be mimicked by administration of LPS, a component of the outer membrane of Gram-negative bacteria. LPS is recognized by the toll-like receptor 4 (TLR4) complex expressed by innate immune cells and also by other cells such as endothelial cells within the lung (21). Systemic administration of LPS initially induces apoptosis of the alveolar capillary endothelium (22) and epithelium (23), and neutrophils are rapidly recruited. The migration of neutrophils into the alveoli is dependent upon crossing the basement membrane, possibly through gaps in the basal laminae (24), and is facilitated by proteolytic degradation of the highly organized extracellular matrix (ECM) (25). Once at the site of injury, activated macrophages and infiltrating neutrophils release a variety of proinflammatory cytokines, enzymes, and reactive oxygen and nitrogen species. These molecules participate in maintaining the inflammatory state and defeating invading pathogens, but also inflict damage of the lung parenchyma and microvascular structure. This results in proteinaceous edema and hemorrhage, as seen in the exudative phase of ALI (26). To the best of our knowledge, no studies have used RNA sequencing to investigate the pulmonary inflammatory response in mammals upon systemic LPS challenge.
We hypothesized that PrPC-deficient goats are more susceptible to inflammatory stimuli, and therefore examined the responsiveness to LPS of PRNPTer/Ter goats, compared with normal goats (PRNP+/+). This study shows a full-scale transcriptional profile coupled to morphological alterations in caprine lungs, an organ highly susceptible to systemic inflammation and investigates whether loss of PrPC affects the initial stage of ALI.
Materials and Methods
A total of 26 Norwegian dairy goat kids, 13 PRNPTer/Ter and 13 PRNP+/+ animals, were included in the study. Details about housing conditions can be found in reference (17) and an overview of the study groups including treatment, animal number, age, weight, and gender can be found in Table S1 in Supplementary Material. The animal experiment was performed in compliance with ethical guidelines and approved by the Norwegian Animal Research Authority (ID 5827, 6903, and 7881) with reference to the Norwegian regulation on animal experimentation (FOR-2015-06-18-761).
Sixteen goats (8 PRNPTer/Ter and 8 PRNP+/+) received a dual dose of LPS (Escherichia coli O26:B6, L2654 Sigma-Aldrich, USA) intravenously with a 24-h time interval between doses; 0.1 µg/kg (day 1) and 0.05 µg/kg (day 2). The control group consisted of 10 goats (5 PRNPTer/Ter and 5 PRNP+/+) and were given corresponding volumes of sterile saline. Clinical examinations, including respiratory rate, were performed at 12 time points during the first 7 h of day 1 and at 9 time points after the second LPS injection. Additionally, episodes of coughing and panting were recorded. The animals were euthanized by an overdose of pentobarbital 5 h after the second LPS challenge.
Gross Pathology and Histopathology Score
After euthanasia, goats were fully necropsied. The amount of fluid in the thoracic cavity and airways, and the distribution and amount of hemorrhages, atelectasis, edema, and emphysema in the lungs were recorded. Additionally, macroscopic changes compatible with lungworm (Muellerius capillaris) were noted.
Tissues were collected at standardized areas dorsally and ventrally from the left caudal lung lobe, as well as from areas with gross pathology, mainly hemorrhages. Tissues were immersion fixed in 4% formaldehyde for 1 week and then dehydrated in graded ethanol and paraffin embedded. Paraffin sections (4 µm) were mounted on slides and stained with hematoxylin and eosin (HE). Martius-scarlet-blue (MSB) staining was performed on selected sections to evaluate fibrin deposition. The HE sections were examined by light microscopy and a blinded, semiquantitative evaluation was performed to assess lung pathology. Both the airways and the lung interstitium were evaluated, focusing on the alveolar lumen and septa. Evaluation parameters included hemorrhage, hyperemia, fibrin deposition, necrosis, and neutrophil and macrophage infiltration. Pathological changes were scored as follows: 0 = minimal; 1 = little/a few, 2 = moderate; 3 = severe, including half-step grading. A total lung histopathology score was calculated for each section by summarizing the values from the above mentioned criteria.
6H4 and SAF32
Frozen lung tissue from a control PRNP+/+ goat was used to evaluate the pulmonary distribution of PrPC. After dissection, lung tissues were frozen in isopentane, transferred to liquid nitrogen, and stored at −70°C. Sections (8 µm) were cut on a cryostat and mounted on Superfrost® Plus slides (Menzel-Gläser, Thermo Scientific). The sections were fixed in formol-calcium, followed by 5 min inhibition of endogenous peroxidase activity (DAKO Peroxidase Block, Agilent). After blocking in normal goat serum (1:50) for 10 min, the sections were incubated with the primary anti-PrP antibodies 6H4 (1:500, Prionics AG) or SAF32 (1:1,000, Cayman) for 1 h at room temperature. Further steps were performed with the Envision + System-HRP DAB (Dako, K4007), according to the manufacturer’s instruction. Antibody specificity was checked by performing the same procedure on lung sections from a PRNPTer/Ter goat.
Formalin-fixed, paraffin-embedded sections (4 µm) were mounted on Superfrost® Plus slides (Menzel-Gläser, Thermo Scientific) and dried for 30 min at 58°C. The sections were deparaffinized in xylene and rehydrated through decreasing concentrations of graded ethanol. The PT link system (100°C for 15 min in Tris/EDTA-buffer, pH 9.1) was used for epitope retrieval. Endogenous peroxidase activity was blocked by incubation in 3% H2O2 in methanol for 20 min at room temperature. Sections were then blocked in normal goat serum (1:50), diluted in phosphate-buffered saline (PBS) and incubated with the primary antibody S100A8 (1:3,200, Abcam) for 1.5 h at room temperature. Further steps were performed with EnVison + System-HRP AEC (Dako, K4009).
Washing between steps was performed with PBS for both procedures. Sections were counterstained in hematoxylin, and mounted using Aquatex medium (Merck). All runs included a negative control section in which the primary antibody was replaced with 1% PBS. The number of cells that were positive for the S100 calcium-binding protein A8 (S100A8) was quantified within the field of view of three representative areas in each animal/section (400× magnification, ø = 450 µm).
RNA Extraction, Quality Control, and Pooling
Tissue samples were collected at a standardized area dorsally from the left caudal lung lobe within 15 min after euthanasia, and immediately immersed in RNA later (Invitrogen) and stored at −70°C. RNA was extracted from approximately 30 mg of lung tissue using the RNeasy Lipid Tissue Mini Kit (Qiagen, 74804) according to the manufacturer’s instruction. High quality RNA (RNA integrity number > 8) from individuals was pooled in equal amounts per treatment and genotype, reaching a final amount of 15 μg, as described in Salvesen et al. (17). RNA quality data are summarized in Table S2 in Supplementary Material.
RNA Sequencing and Differential Expression Analyses
RNA samples were shipped on dry ice to Novogene (Hong Kong) for RNA sequencing. Paired-end 150 bp sequencing was performed on an Illumina HiSeq2000. Quality control summary and mapping status are given in Tables S3 and S4 in Supplementary Material. Normalization was performed by the trimmed mean of M values method, p-value estimation assuming Poisson distribution, and false discovery rate (FDR) estimation by Benjamini-Hochberg correction of multiple testing (27, 28). Differential gene expression analysis was performed using DEGSeq (28) with the following criteria: Log2 ratio ± 0.59 (fold change ± 1.5) and a FDR-adjusted q-value (q < 0.05). The genomic response to LPS in each PRNP genotype was assessed by comparing LPS-groups with the saline-treated control of the matching genotype. Fragments per kilobase of exon per million fragments mapped values, which take into account the effects of both sequencing depth and gene length, were used to estimate gene expression levels.
Ingenuity Pathway Analysis (IPA) and Gene Ontology (GO) Enrichment
Gene lists from Capra hircus were first converted to gene symbols using the NCBI batch converter tool.1 Pathway analysis was performed with both IPA (v. 36601845) and GO analysis with the online DAVID 6.8 software (29, 30). Because the C. hircus genome is not available in both analyses, Homo sapiens orthologs were used instead. The lists of differentially expressed genes (DEGs), containing gene identifiers and corresponding expression values, were uploaded into the IPA software (Qiagen). The top enriched canonical pathways, upstream regulators and gene networks were identified by the “core analysis” function, using a cut-off in absolute fold change of 1.5 and a FDR < 0.05. The IPA software uses the Fishers exact test for pathway enrichment and Z-scores for pathway activation. A Benjamini-Hochberg adjusted p-value < 0.05 was considered significant for GO analysis. DEGs of the two genotypes were analyzed separately.
Validation of RNA-seq by Real-time qPCR
Complementary DNA was generated using RNA from each individual sample by the QuantiTect Reverse Transcription Kit (Qiagen) according to manufacturer’s instructions. A non-reverse transcriptase control and no template control were included. The expression of nine representative target genes [serum amyloid A3.2 (SAA3.2), CXCL10, ITGAM, CD14, IFI6, GADPH, S100A9, PRNP, and interleukin 1B (IL1B)] was investigated by the Light cycler 480 qPCR system. Real-time reactions were prepared using SYBR Green PCR Master Mix including 10 µl cDNA (1:10) in each reaction, and qPCR was performed with standard cycling conditions: initial denaturation for 5 min at 95°C, followed by 40 amplification cycles (10 s at 95°C, 10 s at 60°C, and 15 s at 72°C) and construction of melting curves. A standard curve was generated for each target gene to obtain primer amplification efficiencies, correlations, and dynamic range. Normalization was performed against the ACTB reference gene, which showed no significant change by LPS-treatment in the RNA-seq data, confirming stability as internal control. Relative expression was calculated using the 2−ΔΔCq method as described in Ref. (31). Primer sequences are given in Table S5 in Supplementary Material.
Descriptive and Statistical Analyses
Data are presented as mean ± SEM. Graphical and statistical analyses were performed in GraphPad Prism 6 (GraphPad Software Inc., USA) and Microsoft Excel 2013. A one-way ANOVA, followed by Dunnetts’s post hoc test, was performed when comparing three or more groups. Differences in lung histopathology score between two groups were assessed by the Mann–Whitney U-test. Significance levels p < 0.05; p < 0.01; p < 0.001 are indicated.
Systemic LPS Challenge Induces Tachypnea and Respiratory Distress
After LPS challenge, all eight goats from both genotypes displayed sepsis-like clinical signs, such as sickness behavior, fever, and tachycardia. PrPC-deficient goats displayed prolonged sickness behavior, and hematological analysis showed a rapid and profound leukopenia in both genotypes as we have previously described (17). Respiratory rate increased significantly in both PRNP genotypes after the initial high dose of LPS (0.1 µg/kg), peaking 1 h post-LPS challenge, and remained elevated at all subsequent time points on day 1 (Figure 1A). The mean respiratory rate was higher in the PRNPTer/Ter group at all time points during the first day, albeit not statistically significant. Periods of heavy abdominal breathing, coughing, and panting were observed in seven out of eight animals of both genotypes receiving LPS, whereas none of the control animals displayed alterations in clinical parameters. Clinical signs of respiratory distress were less profound on day 2, when the dosage of LPS was halved.
Figure 1. Respiratory rate and lung gross pathology after systemic lipopolysaccharide (LPS) challenge. Goats were subjected to a dual LPS challenge; 0.1 and 0.05 µg/kg with a 24-h time interval between doses. (A) The diagram shows the respiratory rate (breaths/min) after LPS challenge. Values are mean ± SEM. Symbols * and + represent significantly increased respiratory rate compared with baseline samples (0 h) of the matching genotype. (B) Dorsal view of the right caudal lung lobe. Subplural ecchymosis (arrow) and petechial hemorrhages (arrowhead) display a multifocal distribution.
Systemic LPS Challenge Induces Circulatory Disturbances and Inflammation
Multifocal ecchymosis and petechial hemorrhages were observed in four of eight PRNP+/+ goats and five of eight PRNPTer/Ter goats receiving LPS. The majority of hemorrhages were subpleural and dorsally located in the caudal lung lobes (Figure 1B). Two PRNPTer/Ter animals had hydrothorax with 90 and 400 mL of serous fluid within the thoracic cavity. Both these goats had mild to moderate ventral atelectasis. Four animals (PRNPTer/Ter) had sparse to moderate amounts of froth and fluid within the trachea and main bronchi, and one animal had a parenchymal lung edema, characterized by a wet cut surface. Sparse lungworm-changes (M. capillaris) were equally distributed between the genotypes in both the LPS and control groups.
The lung histopathology scores (Figure 2A) were significantly higher in the LPS-groups than in the saline controls (Figure 2B). Multifocal microscopic hemorrhages were recorded in seven out of eight animals receiving LPS in both genotypes, located both subpleurally and deeper within the lung parenchyma. These areas were characterized by alveolar hemorrhage and edema, increased amount of alveolar macrophages, and infiltration of neutrophils (Figure 2C). Some hyaline membranes were observed lining the alveolar wall (Figure 2D). Focally, an increased number of neutrophils were present in the alveolar septa, accompanied by single cell necrosis and hyperemia (Figure 2E). In these same foci, MSB staining showed fibrin deposition along alveolar septa of a few LPS-treated animals (Figure 2F). Immunohistochemical labeling of S100 calcium binding protein A8 (S100A8) increased upon LPS challenge, and was predominantly expressed by neutrophils, and at lower levels, in alveolar macrophages (Figures 2G–I). The histopathological changes were more profound dorsally than ventrally. There were no significant differences between PRNP genotypes.
Figure 2. Histopathology after lipopolysaccharide (LPS) challenge. (A) A blinded scoring of lung histopathology was performed as described in methods. Sections from the dorsal caudal lung lobe (CLL) and ventral CLL were scored in all animals, in addition to sections from areas with gross pathology, mainly hemorrhages (GPA). Symbols represent significantly increased lung histopathology score compared with saline controls of the matching genotype. GPA is compared with the dorsal CLL of the matching control. (B) Normal lung tissue from a saline control goat. (C) Alveolar edema and hemorrhages accompanied by neutrophil infiltration and increased number of alveolar macrophages. (D) Hyaline septa (arrows) and increased number of alveolar macrophages (arrowheads). (E) Area with moderate atelectasis and hyperemia. Several single cell necroses (arrows) can be seen within the alveolar septa. (F) MSB staining shows deposition of fibrin within the septa (arrows) and small amounts in the alveoli. Sections (C–F) are from areas with pathology in LPS-treated animals. (G) Mean number of S100A8-positive cells quantified within a representative 400× field of view (FOV, ø = 450 µm). Symbols represent significantly increased levels compared with saline controls of the matching genotype. (H) Normal lung tissue from a control animal displaying a few S100A8-positive neutrophils within the alveolar septa. (I) After LPS challenge, increased numbers of S100A8-positive neutrophils are sequestered within the septal capillaries and some infiltrate the alveoli (arrows). Alveolar macrophages displayed a weak S100A8-signal (arrowheads). Insertions show a magnified neutrophil (N) with segmented nucleus and an alveolar macrophage (Mφ). Magnification: (B,C) 200×; (D–I) 400×. Values are mean ± SEM.
PRNP Expression and Cellular Distribution of PrPC
The level of PRNP expression in the lungs was similar to that of the choroid plexus, and five times lower than in the hippocampus (17). After LPS challenge, PRNP expression increased by 1.6-fold in the lungs of PRNP+/+ goats, albeit not significantly. Expression of PRNP is reduced to about 25% in PRNPTer/Ter goats (Figure 3A). The localization of PrPC in the lungs of non-treated animals was investigated by use of 6H4 and SAF32 antibodies. A granular distribution of PrPC was observed within bronchial and bronchiolar epithelial cells of PRNP+/+ goats (Figure 3B). In the alveoli, PrPC was primarily localized at the luminal side of the alveolar septa, probably associated with alveolar pneumocytes (Figure 3C). No PrPC was detected in PRNPTer/Ter goats, confirming the absence of synthesized protein as well as the specificity of the antibodies (Figure 3D). Staining with SAF32 showed similar results (Figure S1 in Supplementary material).
Figure 3. PRNP expression and cellular distribution of PrPC. (A) Comparison of PRNP expression in three tissues from Norwegian dairy goats before and after lipopolysaccharide challenge. Symbols represent significantly (p < 0.001) reduced expression in PRNPTer/Ter goats, compared with PRNP+/+ goats of the same treatment. (B) PrPC shows a granular distribution within bronchial epithelial cells. (C) In the alveoli, PrPC labeling is associated with the luminal surface of the alveolar septa, probably associated with alveolar pneumocytes. (D) No PrPC was detected in sections from PRNPTer/Ter goats. All sections are from representative areas stained with the 6H4 antibody. Magnification (B,C) 400×; (D) 100×.
The Lung Transcriptome Is Profoundly Altered by Systemic LPS Challenge
After LPS challenge, a total of 432 [upregulated: 269, downregulated: 163] and 596 (upregulated: 339, downregulated: 257) genes were differentially expressed compared to the saline control group in normal and PrPC-deficient goats, respectively (Data Sheet S1 in Supplementary Material). The Venn diagrams show that 146 upregulated and 161 downregulated genes exclusively belonged to the PrPC-deficient goats when the filtration criteria (log2 ratio 0.59, q < 0.05) were applied (Figure 4A). However, closer investigation of the transcriptional data revealed that 95% of the DEGs were regulated in the same direction in both genotypes, albeit some excluded by fold change or q-value. Genes upregulated by a log2 ratio ≥ 3 (fold change ≥ 8) after LPS challenge are listed in Table 1. The most upregulated gene in both PRNP genotypes was the acute phase protein SAA3 (51-fold, PRNPTer/Ter; 25-fold, PRNP+/+). Other highly upregulated genes were cytokines, chemokines, and metal-binding proteins such as metallothioneins and several S100 proteins. The most downregulated gene was CYP1A1 (24-fold, PRNPTer/Ter; 22-fold, PRNP+/+), a member of the cytochrome P450 superfamily (Data Sheet S1 in Supplementary Material).
Figure 4. Differentially expressed genes after LPS challenge and biological function analysis. (A) Differential expression analyses were performed by comparing the LPS-groups to the saline control group of the matching genotype (log2 ratio 0.59 and q-value < 0.05). Venn diagrams show genes that overlap between the two genotypes. (B) Shown are the top 15 overrepresented upstream regulators identified by ingenuity pathway analysis (IPA). (C) Venn diagram showing the distribution of genes characterized as collagen catabolic process (GO:0030574) in the two PRNP genotypes. (D) The top 15 affected pathways identified by IPA are shown.
Affected Pathways Are Linked to LPS Challenge and Indicate Enhanced Inflammation and Turnover of ECM in PrPC-Deficient Goats
To assess the function of the DEGs, enriched pathways, networks and upstream regulators were identified by IPA. DEGs were also assigned to GO classifications, including biological process, cellular compartment and molecular function. In general, there was a high degree of overlap in the transcriptional response to LPS in normal and PrPC-deficient goats. As expected, a considerable amount of the DEGs belonged to the innate immune response and inflammatory response (Data Sheet S2 in Supplementary Material), and LPS was identified as the most significant upstream regulator of both genotypes (Figure 4B). The associated p-values indicate more enrichment of genes involved in the LPS-response in PRNPTer/Ter goats, particularly those stimulated by the proinflammatory cytokines IL-1β, IFN-γ, and TNF-α (Figure 4B). Other notably affected pathways in both genotypes involved acute-phase response signaling, antigen presentation, and granulocyte adhesion and diapedesis (Figure 4D). Although some upstream regulators of interferons were more significantly enriched in PrPC-deficient goats, predicted activation of interferon signaling were lower than in the normal goats. The full IPA output can be found in Data Sheet S3 in Supplementary Material.
Several of the GO-terms, such as collagen catabolic process (GO:0030574), ECM disassembly (GO0022617), and ECM structural constituents (GO:0005201) were exclusively enriched in PrPC-deficient goats (Data Sheet S2 in Supplementary Material). The genes were primarily in an upregulated state, and included collagen-encoding genes, secreted metalloproteinases, such as matrix metalloproteinases (MMPs) and ADAM genes, as well as cathepsins primarily located within lysosomes (Figure 4C). In line with this, the top canonical pathway in PRNPTer/Ter goats was related to turnover of ECM (Figure 4D). Also, IPA analysis revealed a significantly enriched network involved in connective tissue disorders (−logP 39, Figure 5).
Figure 5. Network of genes related to connective tissue disorders in PRNPTer/Ter goats. An overrepresented network is shown exclusively for PRNPTer/Ter goats of differentially expressed genes and interacting partners related to connective tissue disorders. Red means upregulated, green means downregulated.
Validation of RNA-seq Data by qPCR
The expression of nine target genes (SAA3.2, CXCL10, ITGAM, CD14, IFI6, GADPH, S100A9, PRNP, and IL1B) was investigated by qPCR on individual RNA samples from lung tissue (Figure 6). A strong correlation (r = 0.975, p < 0.0001, Pearson correlation) was observed between log2-transformed expression values of RNA-seq and qPCR, in line with our previous reports from other tissues (16, 17). Differential expression was confirmed by qPCR in 34 out of 36 comparisons.
Figure 6. Validation of RNA-seq data by qPCR. Nine target genes representing both differentially expressed and non-differentially expressed genes were measured by qPCR. For each gene, the expression levels are normalized relative to the saline-treated PRNP+/+ group. RNA-seq symbols above bars represent a significant differential expression (log2 ratio ± 0.59, q < 0.05). Bars (qPCR) display mean expression ± SEM with symbols representing significant differences, as assessed by multiple t-tests (p < 0.05).
We performed an LPS challenge in goats with and without PrPC and characterized the clinical, morphological, and transcriptional response to acute pulmonary inflammation. Both genotypes displayed tachypnea, clinical signs of respiratory distress, and histopathological evidence of ALI. A significant enrichment of ECM turnover, along with an increased activation score of genes related to the proinflammatory LPS response was observed in PrPC-deficient goats.
In the current study, more than 70% of the DEGs classified as being involved in collagen catabolic processes belonged exclusively to PRNPTer/Ter goats, and the connective tissue disorder network was significantly enriched. The collagen-encoding genes were primarily upregulated, probably as a result of resident fibroblast activation mediated by transforming growth factor (TGF)-β (32, 33). In line with this, the upstream regulator TGF-β showed an increased activation score in PrPC-deficient goats. Indeed, PrPC has been reported to modulate B1 integrin signaling directly (34, 35), which is a key regulator of TGF-β-dependent activation of fibroblasts (36). Homeostasis of the ECM is of major importance for the physiological function of most organs, and is maintained by secretion of proteases and degradation of ECM components, as well as synthesis of collagens (37). In addition, both ECM constituents and MMPs may act as modulators of inflammation by regulating barrier function along with cytokine and chemokine activity (25, 38, 39). Even though collagen deposition is an essential step in limiting alveolar damage and initiating repair (40), elevated levels of procollagen (COL3) and proteases (MMP1/3) have also been reported as biomarkers of poor prognosis in human ALI patients (41, 42). Altogether, this indicates that ECM degradation and subsequent collagen synthesis occurred to a greater extent in the absence of PrPC following LPS challenge. Although no significant histopathological differences between genotypes were observed in the acute phase of ALI, abnormal collagen homeostasis may lead to morphological alterations in the later reparative phase.
When analyzed by IPA, we found that genes associated with proinflammatory upstream regulators, such as TNF-α, IL-1β, and IFN-γ, were more activated in PrPC-deficient goats. We also observed a lower activation score of the interferon response pathway in the PRNPTer/Ter group, but this reflects that the PrPC-deficient control goats are in a primed state with respect to interferon-responsive genes (16), rather than having a reduced activation potential. Interestingly, type I interferons and some of the downstream effector genes are considered as important regulators of lung inflammation (43). Because most PrPC is localized to the cell surface, it has been hypothesized that PrPC binds extracellular proteins and facilitates transmembrane signaling processes (44). For instance, PrPC desensitized cells to TNF-α by restricting the level of available TNF-receptor present at the plasma membrane of neuroectodermal cells (34). We have recently proposed that PrPC protects tissues by modulating certain signaling pathways of innate immunity (17). This is in agreement with murine studies reporting that lack of PrPC exacerbates inflammation in different organ systems (8, 9, 11). Thus, increased expression of MMPs and ECM constituents might demonstrate an aggravated inflammatory damage of the lung parenchyma in goats devoid of PrPC.
We observed moderate levels of PRNP mRNA in the lungs, in accordance with reports from sheep and cattle (3, 4, 45, 46). After LPS challenge, PRNP expression increased slightly (1.6-fold) in PRNP+/+ goats, indicating a response to inflammatory stimuli. We have earlier published the same tendency in the hippocampus and choroid plexus (17). IHC showed that PrPC was distributed in several cell types throughout the lungs, similar to what has been reported in a bovine study (45). As expected, the absence of PrPC protein was confirmed in PRNPTer/Ter goats.
Systemic LPS challenge was associated with a dramatic change in the lung transcriptome, in which 432 and 596 genes were differentially expressed in normal and PrPC-deficient goats, respectively. Besides the genotype-dependent profile in ECM homeostasis and inflammatory regulators, there was a considerable overlap in the molecular response to LPS. Not surprisingly, the majority of the genes were annotated as being involved in the inflammatory and innate immune response, and showed a comparable profile to those reported in other species (47). The most upregulated gene in both genotypes was the acute phase protein SAA3, an isoform of the SAA apolipoproteins. Indeed, SAA3 seems to be a hallmark of acute inflammation in a range of tissues, at least in ruminants (17, 48). Interestingly, SAA3 may act as a ligand for the TLR4-MD2-complex, and play multiple roles in modulating the innate immune response (49). Four S100 proteins (S100G, S100A8, S100A9, and S100A12) were also among the most prominently upregulated genes in both genotypes. These low-molecular weight proteins have a Ca2+-binding capacity, and as well as being involved in regulation of proliferation, differentiation and apoptosis (50), also have antimicrobial activity against a variety of microorganisms (51). Because S100A8/A9/A12 genes are constitutively expressed in neutrophils, as demonstrated by S100A8 IHC in the present study, infiltrating neutrophils probably contribute to the upregulation of these genes in our model (52). Although the exact functions of S100 proteins in lung inflammation remain currently unresolved, at least S100A8/A9 seems to reduce the development of ALI (53). The most downregulated gene, CYP1A1, plays a key role in the metabolism of xenobiotics and some endogenous substrates (54). Similarly, systemic LPS suppressed CYP1A1 expression in the murine liver (55) and lungs (56). The downregulation could result from the formation of reactive oxygen species (ROS) generated by LPS challenge, as ROS have been shown to suppress CYP1A1 expression in vitro (57). Conversely, LPS treatment induced CYP1A1 expression in cultured human dendritic cells and in the thymus (58), suggesting that tissue-specific differences may occur.
Acute inflammation is a balancing act; although it is an essential component of an appropriate response and recovery from an insult, it may also damage the tissue in which it occurs. The multifocal areas with alveolar hemorrhage, edema and infiltration of inflammatory cells shown histologically, indicate disruption of the alveolar/capillary basement membrane. This is probably due to the massive release of proteases and inflammatory mediators (59, 60). The observed hydrothorax in two PrPC-deficient goats most likely represents transudative pleural effusion resulting from compromised integrity of the lung capillaries. Neutrophils were found sequestered within the pulmonary capillaries and some migrated into the alveoli upon LPS challenge. Correspondingly, several genes mapped to neutrophil chemotaxis, such as S100A8/A9/A12 proteins and chemokines (CCL2, CCL5 and CCL15), and IPA confirmed the involvement of granulocyte diapedesis. In some areas, the normal septal structure was replaced by fibrinous exudate compatible with hyaline septa or debris from necrotic cells (26). Although the mean lung histopathology score was slightly higher in PRNPTer/Ter goats than in normal goats, the difference was not statistically significant. The magnitude of LPS-effects on inflammation probably mask any putative histopathological differences associated with the demonstrated genotype-dependent differences in gene transcription, and thus would be difficult to detect by semiquantitative evaluation.
The substantial activation of the lung transcriptome and histopathological evidence of lung injury corresponded with the clinical signs observed. The presence of pulmonary intravascular macrophages, specialized in removing circulating bacteria and endotoxins, probably means that small ruminants are particularly sensitive to systemic LPS exposure (61). A significant increase in respiratory rate was observed, peaking at 1 h, and remained elevated the first 7 h postexposure. Moreover, periods of evident respiratory distress were observed in nearly 90% of the animals receiving LPS, indicating temporarily reduced lung function. This is in agreement with a significant enrichment of genes characterized as hypoxia-responsive. Given that only subtle effects have been attributed to PrPC, it was not surprising that loss of PrPC did not affect the overall clinical recovery and outcome.
In conclusion, we demonstrated a potent activation of the lung transcriptome that corresponded with morphological and clinical alterations. Although there was a substantial overlap in the response to LPS in both genotypes, increased turnover of ECM and activation of upstream inflammatory regulators were present in PrPC-deficient goats. Future studies should investigate whether these alterations affect the later reparative phase of ALI.
Availability of Data and Material
Original BAM files containing sequencing data are available from the SRA database (SRA accession number SRP116648). Data and material supporting the conclusions are contained within the manuscript and supplementing material.
The animal experiment was performed in compliance with ethical guidelines and approved by the Norwegian Animal Research Authority (ID 5827, 6903, and 7881) with reference to the Norwegian regulation on animal experimentation (FOR-2015-06-18-761).
CE, MT, AE, MB, and ØS designed the study and ØS, MR, CE, AE, and MT performed the experiments. ØS carried out laboratory procedures, performed the statistical analysis, and drafted the manuscript. MR performed the PrPC immunohistochemistry. JK performed the ingenuity pathway analysis. All authors have critically read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors gratefully thank Siri Bjerkreim Hamre and Mari Katharina Aas Ådland for skillful laboratory work. The authors acknowledge Lucy Robertson for proofreading the manuscript.
The research was supported by The Research Council of Norway (grant ID 227386/E40).
The Supplementary Material for this article can be found online at http://www.frontiersin.org/article/10.3389/fimmu.2017.01722/full#supplementary-material.
LPS, lipopolysaccharide; qPCR, quantitative real-time polymerase chain reaction; DEGs, differentially expressed genes; GO, gene ontology; IPA, ingenuity pathway analysis; IHC, immunohistochemistry; IL, interleukin; IFN, interferon; TNF, tumor necrosis factor; TGF, transforming growth factor.
1. McLennan NF, Rennison KA, Bell JE, Ironside JW. In situ hybridization analysis of PrP mRNA in human CNS tissues. Neuropathol Appl Neurobiol (2001) 27(5):373–83. doi:10.1046/j.0305-1846.2001.00343.x
2. Tanji K, Saeki K, Matsumoto Y, Takeda M, Hirasawa K, Doi K, et al. Analysis of PrPC mRNA by in situ hybridization in brain, placenta, uterus and testis of rats. Intervirology (1995) 38(6):309–15. doi:10.1159/000150457
3. Tichopad A, Pfaffl MW, Didier A. Tissue-specific expression pattern of bovine prion gene: quantification using real-time RT-PCR. Mol Cell Probes (2003) 17(1):5–10. doi:10.1016/S0890-8508(02)00114-7
4. Moudjou M, Frobert Y, Grassi J, La Bonnardiére C. Cellular prion protein status in sheep: tissue-specific biochemical signatures. J Gen Virol (2001) 82(Pt 8):2017–24. doi:10.1099/0022-1317-82-8-2017
6. Hirsch TZ, Hernandez-Rapp J, Martin-Lannerée S, Launay JM, Mouillet-Richard S. PrPC signalling in neurons: from basics to clinical challenges. Biochimie (2014) 104:2–11. doi:10.1016/j.biochi.2014.06.009
9. Gourdain P, Ballerini C, Nicot AB, Carnaud C. Exacerbation of experimental autoimmune encephalomyelitis in prion protein (PrPC)-null mice: evidence for a critical role of the central nervous system. J Neuroinflammation (2012) 9:25. doi:10.1186/1742-2094-9-25
10. Ding T, Zhou X, Kouadir M, Shi F, Yang Y, Liu J, et al. Cellular prion protein participates in the regulation of inflammatory response and apoptosis in BV2 microglia during infection with Mycobacterium bovis. J Mol Neurosci (2013) 51(1):118–26. doi:10.1007/s12031-013-9962-2
11. Liu J, Zhao D, Liu C, Ding T, Yang L, Yin X, et al. Prion protein participates in the protection of mice from lipopolysaccharide infection by regulating the inflammatory process. J Mol Neurosci (2014) 55(1):279–87. doi:10.1007/s12031-014-0319-2
12. Bakkebø MK, Mouillet-Richard S, Espenes A, Goldmann W, Tatzelt J, Tranulis MA. The cellular prion protein: a player in immunological quiescence. Front Immunol (2015) 6:450. doi:10.3389/fimmu.2015.00450
13. Mariante RM, Nóbrega A, Martins RA, Areal RB, Bellio M, Linden R. Neuroimmunoendocrine regulation of the prion protein in neutrophils. J Biol Chem (2012) 287(42):35506–15. doi:10.1074/jbc.M112.394924
15. Reiten MR, Bakkebø MK, Brun-Hansen H, Lewandowska-Sabat AM, Olsaker I, Tranulis MA, et al. Hematological shift in goat kids naturally devoid of prion protein. Front Cell Dev Biol (2015) 3:44. doi:10.3389/fcell.2015.00044
16. Malachin G, Reiten MR, Salvesen Ø, Aanes H, Kamstra JH, Skovgaard K, et al. Loss of prion protein induces a primed state of type I interferon-responsive genes. PLoS One (2017) 12(6):e0179881. doi:10.1371/journal.pone.0179881
17. Salvesen Ø, Reiten MR, Espenes A, Bakkebø MK, Tranulis MA, Ersdal C. LPS-induced systemic inflammation reveals an immunomodulatory role for the prion protein at the blood-brain interface. J Neuroinflammation (2017) 14(1):106. doi:10.1186/s12974-017-0879-5
21. Andonegui G, Bonder CS, Green F, Mullaly SC, Zbytnuik L, Raharjo E, et al. Endothelium-derived toll-like receptor-4 is the key molecule in LPS-induced neutrophil sequestration into lungs. J Clin Invest (2003) 111(7):1011–20. doi:10.1172/JCI16510
22. Fujita M, Kuwano K, Kunitake R, Hagimoto N, Miyazaki H, Kaneko Y, et al. Endothelial cell apoptosis in lipopolysaccharide-induced lung injury in mice. Int Arch Allergy Immunol (1998) 117(3):202–8. doi:10.1159/000024011
24. Walker DC, Behzad AR, Chu F. Neutrophil migration through preexisting holes in the basal laminae of alveolar capillaries and epithelium during streptococcal pneumonia. Microvasc Res (1995) 50(3):397–416. doi:10.1006/mvre.1995.1067
28. Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics (2010) 26(1):136–8. doi:10.1093/bioinformatics/btp612
30. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res (2009) 37(1):1–13. doi:10.1093/nar/gkn923
33. Armstrong L, Thickett DR, Mansell JP, Ionescu M, Hoyle E, Billinghurst RC, et al. Changes in collagen turnover in early acute respiratory distress syndrome. Am J Respir Crit Care Med (1999) 160(6):1910–5. doi:10.1164/ajrccm.160.6.9811084
34. Ezpeleta J, Boudet-Devaud F, Pietri M, Baudry A, Baudouin V, Alleaume-Butaux A, et al. Protective role of cellular prion protein against TNFα-mediated inflammation through TACE alpha-secretase. Sci Rep (2017) 7(1):7671. doi:10.1038/s41598-017-08110-x
35. Loubet D, Dakowski C, Pietri M, Pradines E, Bernard S, Callebert J, et al. Neuritogenesis: the prion protein controls beta1 integrin signaling activity. FASEB J (2012) 26(2):678–90. doi:10.1096/fj.11-185579
36. Hayashida T. Integrins modulate cellular fibrogenesis at multiple levels; regulation of TGF-beta signaling. Endocr Metab Immune Disord Drug Targets (2010) 10(4):302–19. doi:10.2174/1871530311006040302
41. Chesnutt AN, Matthay MA, Tibayan FA, Clark JG. Early detection of type III procollagen peptide in acute lung injury. Pathogenetic and prognostic significance. Am J Respir Crit Care Med (1997) 156(3):840–5. doi:10.1164/ajrccm.156.3.9701124
42. Fligiel SE, Standiford T, Fligiel HM, Tashkin D, Strieter RM, Warner RL, et al. Matrix metalloproteinases and matrix metalloproteinase inhibitors in acute lung injury. Hum Pathol (2006) 37(4):422–30. doi:10.1016/j.humpath.2005.11.023
47. Pennings JL, Kimman TG, Janssen R. Identification of a common gene expression response in different lung inflammatory diseases in rodents and macaques. PLoS One (2008) 3(7):e2596. doi:10.1371/journal.pone.0002596
48. Salvesen Ø, Reiten MR, Heegaard PM, Tranulis MA, Espenes A, Skovgaard K, et al. Activation of innate immune genes in caprine blood leukocytes after systemic endotoxin challenge. BMC Vet Res (2016) 12(1):241. doi:10.1186/s12917-016-0870-x
49. Deguchi A, Tomita T, Omori T, Komatsu A, Ohto U, Takahashi S, et al. Serum amyloid A3 binds MD-2 to activate p38 and NF-κB pathways in a MyD88-dependent manner. J Immunol (2013) 191(4):1856–64. doi:10.4049/jimmunol.1201996
51. Hsu K, Champaiboon C, Guenther BD, Sorenson BS, Khammanivong A, Ross KF, et al. Anti-infective protective properties of S100 calgranulins. Antiinflamm Antiallergy Agents Med Chem (2009) 8(4):290–305. doi:10.2174/187152309789838975
52. Foell D, Wittkowski H, Vogl T, Roth J. S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukoc Biol (2007) 81(1):28–37. doi:10.1189/jlb.0306170
55. Ke S, Rabson AB, Germino JF, Gallo MA, Tian Y. Mechanism of suppression of cytochrome P-450 1A1 expression by tumor necrosis factor-alpha and lipopolysaccharide. J Biol Chem (2001) 276(43):39638–44. doi:10.1074/jbc.M106286200
56. Wu D, Li W, Lok P, Matsumura F, Vogel CF. AhR deficiency impairs expression of LPS-induced inflammatory genes in mice. Biochem Biophys Res Commun (2011) 410(2):358–63. doi:10.1016/j.bbrc.2011.06.018
57. Morel Y, Barouki R. Down-regulation of cytochrome P450 1A1 gene promoter by oxidative stress – critical contribution of nuclear factor 1. J Biol Chem (1998) 273(41):26969–76. doi:10.1074/jbc.273.41.26969
58. Vogel CFA, Khan EM, Leung PSC, Gershwin ME, Chang WLW, Wu DL, et al. Cross-talk between aryl hydrocarbon receptor and the inflammatory response: a role for nuclear factor κB. J Biol Chem (2014) 289(3):1866–75. doi:10.1074/jbc.M113.505578
59. Pugin J, Verghese G, Widmer MC, Matthay MA. The alveolar space is the site of intense inflammatory and profibrotic reactions in the early phase of acute respiratory distress syndrome. Crit Care Med (1999) 27(2):304–12. doi:10.1097/00003246-199902000-00036
Keywords: acute lung injury, cellular prion protein, systemic inflammation, lipopolysaccharide, innate immunity, RNA sequencing, extracellular matrix remodeling, histopathology
Citation: Salvesen Ø, Reiten MR, Kamstra JH, Bakkebø MK, Espenes A, Tranulis MA and Ersdal C (2017) Goats without Prion Protein Display Enhanced Proinflammatory Pulmonary Signaling and Extracellular Matrix Remodeling upon Systemic Lipopolysaccharide Challenge. Front. Immunol. 8:1722. doi: 10.3389/fimmu.2017.01722
Received: 04 September 2017; Accepted: 21 November 2017;
Published: 06 December 2017
Edited by:Rudolf Lucas, Medical College of Georgia, Augusta University, United States
Reviewed by:Julio Villena, Centro de Referencia para Lactobacilos, Argentina
Justin J. Greenlee, National Animal Disease Center (USDA ARS), United States
Copyright: © 2017 Salvesen, Reiten, Kamstra, Bakkebø, Espenes, Tranulis and Ersdal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Øyvind Salvesen, email@example.com