Physiol. Genomics 27: 245-263, 2006.
First published August 1, 2006; doi:10.1152/physiolgenomics.00261.2005

1094-8341/06 $8.00
Received 19 October 2005;
accepted in final form 26 July 2006.
Physiological Genomics 27:245-263 (2006)
1094-8341/06 $8.00 © 2006 American Physiological Society
Temporal gene expression profiles of target-ablated olfactory epithelium in mice with disrupted expression of scavenger receptor A: impact on macrophages
M. L. Getchell1,2,
H. Li3,
R. A. Vaishnav4,
A. S. Borders4,
J. Witta5,
N. Subhedar4,
W. de Villiers5,
A. J. Stromberg3 and
T. V. Getchell2,4
1 Department of Anatomy and Neurobiology, University of Kentucky College of Medicine, Lexington, Kentucky
2 Sanders-Brown Center on Aging, University of Kentucky College of Medicine, Lexington, Kentucky
3 Department of Statistics, University of Kentucky College of Medicine, Lexington, Kentucky
4 Department of Physiology, University of Kentucky College of Medicine, Lexington, Kentucky
5 Department of Internal Medicine, University of Kentucky College of Medicine, Lexington, Kentucky
 |
ABSTRACT
|
|---|
Target ablation [removal of the olfactory bulb (OBX)] induces apoptotic death of olfactory sensory neurons (OSNs) and an immune response in which activation and recruitment of macrophages (m
s) into the olfactory epithelium (OE) occupy a central role. M
s phagocytose apoptotic neurons and secrete cytokines/growth factors that regulate subsequent progenitor cell proliferation and neurogenesis. Scavenger receptor A (SR-A) is a pattern recognition receptor that mediates binding of m
s to apoptotic cells and other relevant immune response functions. The aim of this study was to determine the impact of the absence of SR-A on the immune response to OBX. The immune response to OBX was evaluated in mice in which functional expression of the m
scavenger receptor (MSR) was eliminated by gene disruption (MSR/) and wild-type (wt) mice of the same genetic background. OBX induced significant apoptotic death of mature OSNs in the two strains. However, subsequent m
infiltration and activation and progenitor cell proliferation were significantly reduced in MSR/ vs. wt mice. Gene expression profiling at short intervals after OBX demonstrated significant differences in temporal patterns of expression of several gene categories, including immune response genes. Many immune response genes that showed different temporal patterns of expression are related to m
function, including cytokine and chemokine secretion, phagocytosis, and m
maturation and activation. These studies suggest that impairment of the immune response to OBX in the OE of MSR/ mice most likely resulted from decreased m
adhesion and subsequent reduced infiltration and activation, with a resultant decrease in neurogenesis.
neurogenesis; microarray; phagocytosis; chemokines; cytokines; dendritic cells
 |
INTRODUCTION
|
|---|
APOPTOSIS OF MATURE olfactory sensory neurons (OSNs) in the olfactory epithelium (OE) and their replacement by neurogenesis occur throughout the lifespan of vertebrates. To study the regulation of apoptosis and neurogenesis in the OE, investigators have utilized a model in which the olfactory bulb, which contains the synaptic targets of OSNs, is surgically removed [olfactory bulbectomy (OBX)], severing the OSN axons. This results in the synchronous apoptosis of mature OSNs (e.g., 12); activation, infiltration, and recruitment of macrophages (m
s) into the OE (e.g., 27); phagocytosis of apoptotic OSNs (e.g., 74); proliferation of OSN progenitors (e.g., 9) leading to neurogenesis; and epithelial remodeling (e.g., 25). Our laboratory and others have utilized OBX to investigate the multiple roles played by m
s in maintaining epithelial homeostasis, including phagocytosis of apoptotic OSNs (74), production of chemokines that regulate recruitment of monocytes/ m
s to the OE (26, 27, 50), and the secretion of cytokines and growth factors that stimulate neurogenesis and epithelial remodeling (25, 61).
Scavenger receptors are pattern-recognition receptors that contribute to tissue homeostasis and immune response (29, 75). Scavenger receptor A (SR-A, CD204), encoded by the murine gene Msr1 (macrophage scavenger receptor 1), is expressed primarily by mature m
s (but not monocytes) and dendritic cells (DCs), as well as some endothelial and vascular smooth muscle cells (13, 64). SR-A has been implicated in multiple processes, including the uptake of modified low-density lipoprotein, bacteria (e.g., 59, 64), and, of particular relevance to this study, apoptotic cells (64, 76); adhesion to cells such as activated B lymphocytes (81) and to extracellular matrix (e.g., 20, 30, 47, 65); and major histocompatibility complex (MHC) class I-restricted antigen presentation (58, 62). Its expression has been detected in models of neuronal degeneration (5, 60) and in microglia associated with ß-amyloid plaques in the brain (7, 11).
Phagocytosis, a key process in the immune response to OBX, is a highly complex, multistep, coordinated activity (e.g., 15, 51, 55, 67, 72). Apoptotic cells initiate phagocytosis by releasing "find-me" signals, which attract and activate m
s (e.g., 32, 52). Activated resident m
s secrete chemokines and cytokines that perpetuate and amplify the response. Chemokines are potent chemoattractants that are responsible for the recruitment of monocytes into tissues (e.g., 43), their subsequent activation and maturation (56), and the activation of microvascular endothelial cells to facilitate phagocyte recruitment (37). Cytokines also regulate monocyte maturation and m
activation into specific phenotypes (56) and stimulate local endothelial cell activation (45). Following infiltration into the tissue, activated phagocytes identify apoptotic cells through their display of "eat-me" or engulfment signals that bind to engulfment receptors on the phagocyte, sometimes in association with bridging molecules that link the two cells (e.g., Refs. 14, 17). This triggers signaling that results in the formation of the phagosome. The phagosome encloses the apoptotic cell and/or its fragments and ultimately processes them for disposal (22) or antigen presentation (33). Phagocytosis culminates in the release from m
s of anti-inflammatory mediators that suppress the development of tissue-destructive immune processes and of endothelial and epithelial growth and survival factors that promote neurogenesis and epithelial remodeling (16, 28).
The availability of a murine model in which exon 4 of the gene encoding SR-A was disrupted, blocking trimerization of the gene product into functional integral membrane receptors of both type I and type II SR-A (73), made it possible to investigate the potential involvement of SR-A in the immune response of the OE to OBX. MSR/ mice had normal blood levels of plasma cholesterol, triglycerides, and lipoproteins when maintained on a standard diet (73). Additionally, during their embryonic development, there was no impairment in the clearance of apoptotic cells due to the participation of and compensation by other cell types and receptors (46).
The aims of this study were 1) to determine if there were differences in the cellular responses to OBX, particularly m
infiltration into the OE and progenitor cell proliferation, between wild-type (wt) C57BL/6 and MSR/ mice, 2) to identify categories of genes in which significantly different temporal patterns of expression in the OE occurred in response to OBX in the two strains, and 3) to describe in detail the temporal pattern of expression of genes related to the innate immune response whose expression in the OE was significantly up- or downregulated by OBX in wt mice and to compare their expression patterns in the two strains. By comparing the cellular and molecular responses to OBX in the OE of wt and MSR/ mice, we identified key genes and pathways involved the innate immune response to OBX and identified a role for SR-A in OE homeostasis.
 |
MATERIALS AND METHODS
|
|---|
Animals and tissue preparation.
We obtained 36 6-wk-old male mice of the wt C57BL/6 strain from the Jackson Laboratory (Bar Harbor, ME). We obtained 36 6-wk-old male mice of the MSR/ strain, in which exon 4 of the Msr1 gene was disrupted (73), that were backcrossed onto a C57BL/6 background (2) from a breeding colony maintained at University of Kentucky. MSR/ mice were genotyped according to the protocol of Babaev et al. (2) to confirm the absence of SR-A. All mice were maintained in the facilities of the Department of Laboratory Animal Research on a 12:12-h light/dark cycle and were given food and water ad libitum. All protocols were implemented in accordance with National Institutes of Health guidelines and approved by the University of Kentucky Institutional Animal Care and Use Committee.
We used 21 mice of each strain for cellular studies to determine their olfactory phenotype and to investigate their response to bilateral OBX. Fifteen mice of each strain underwent bilateral OBX; three of the remaining six mice of each strain underwent sham bilateral OBX (0 h sham OBX), and three of each strain were euthanized without any surgical manipulation. The latter six mice were perfused, and tissues were removed and processed as described below; the tissues from these mice were used for histological and immunohistochemical procedures. No differences in the cellular parameters described in RESULTS were noted between the latter two groups; only the data from the 0-h sham OBX mice are reported. Both bilateral sham OBX and bilateral OBX were performed according to this laboratory's protocols as previously described (26, 61). All mice whose tissues were used for histological and immunohistochemical studies were injected with 5-bromo-2'-deoxyuridine (BrdU, 100 mg/kg ip; Sigma Chemical, St. Louis, MO) in sterile saline 2 h before perfusion. For tissue collection, three mice of each strain with no surgical manipulation and three mice of each strain per time point were euthanized by CO2 asphyxiation at 0 (sham), 2, 8, 16, or 48 h or 3 days post-OBX. These mice were perfused transcardially with 0.1 M phosphate-buffered saline (PBS, pH 7.4) followed by 3% paraformaldehyde. Tissues were postfixed for 2 h in 3% paraformaldehyde, cryoprotected, embedded, frozen, and sectioned as previously described (26, 61).
We used 15 mice of each strain at 0 h (sham) and 2, 8, 16, and 48 h post-OBX, 3/time point, for microarray analysis and real-time RT-PCR studies. They were euthanized by CO2 asphyxiation, and olfactory-nasal mucosa and liver were rapidly removed, weighed, frozen in liquid N2, and stored at 80°C.
All reagents used in animal preparation, tissue collection, and in situ hybridization protocols were molecular biology grade. Surgical instruments and disposables were RNase-free, and surgery and tissue collection and handling were performed with RNase-free techniques.
Histology and immunohistochemistry.
Cresyl violet staining was performed on defatted, rehydrated frozen sections with 1% aqueous cresyl violet (Sigma). DNA staining was performed on rehydrated frozen sections with a 1:105 dilution of bisbenzimide (Hoechst 33258, Sigma) in PBS.
A standard indirect immunofluorescence technique was used as previously described (61) to visualize immunoreactivity for the mature OSN marker olfactory marker protein (OMP; a generous gift of Dr. Frank Margolis) and macrophages (CD68 and F4/80). The height of the OE on the nasal septum was measured from the basement membrane to the top of the visible ciliary layer as previously described (61). The numbers of CD68+ and F4/80+ macrophages and BrdU+ proliferative basal cells were counted in the same region along 69 lengths of OE basement membrane, each measuring 200 µm in length as previously described (61). Statistical analysis was performed using two-tailed unpaired t-tests (GraphPad InStat v. 3.06; GraphPad Software, San Diego, CA).
RNA isolation.
Frozen tissue was pulverized in TRI reagent (Sigma-Aldrich) and processed through a QIAshredder column (Qiagen, Valencia, CA), and total RNA was purified using the Qiagen RNeasy Mini-Kit according to the manufacturer's protocol. Total RNA yield and purity were assessed spectrophotometrically and with the model 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA); all samples had A260/A280 ratios of 1.92.2 and showed two sharp peaks corresponding to 18S and 28S RNA on Bioanalyzer electropherograms.
High-density oligonucleotide arrays.
Affymetrix Murine Genome U74Av2 GeneChips (Affymetrix, Santa Clara, CA) were processed at the University of Kentucky Microarray Core Facility. Thirty GeneChips were used. Hybridization of unpooled RNA, normalization, and identification of present genes were performed as previously described (50). The microarray data set was deposited into the Gene Expression Omnibus database; the series accession number is GSE3455.
Statistical analyses of microarray data.
Mean hybridization signals were analyzed using a 2 x 5 ANOVA (overall P < 0.05), with strain as one factor and time post-OBX [0, 2, 8, 16, or 48 hours post-OBX (hpO)] as the other. The false discovery rate (FDR) correction (6, 66) at 1% was applied, reducing the effective P value to <3.2E-03. Identification of genes whose temporal pattern of expression differed significantly in the two strains was determined by using
new = 3.2E-03, determined by the previous FDR procedure as the cutoff point. The protected Fisher's least significant difference test procedure at the level of
new was adopted to control the family-wise error rate in identifying the subsequent interaction effect between strain and time, main effects, and/or simple effects. Data were analyzed using Excel (Microsoft, Redmond, WA), SAS (SAS Institute, Cary, NC), and SigmaPlot (SPSS, Chicago, IL). To determine if functional categories as defined by the Gene Ontology (http://www.geneontology.org) were overrepresented among genes with highly significant differences (P < 3.2E-03) in expression levels in the two strains atparticular post-OBX time points, we used the Expression Analysis Systematic Explorer (EASE) software (36). EASE analysis identified categories of genes with significantly different temporal patterns (EASE scores
0.05) in the two strains.
Real-time RT-PCR.
Real-time PCR was performed on aliquots of cDNA reverse-transcribed as previously described (50) from the same total RNA that was used for microarray hybridization. Primers for each gene, which spanned at least one intron, were: Il1rn, forward 5'-GCAAGCCTTCAGAATCTGGGATAC-3', reverse 5'-CTCAGAGCGGATGAAGGTAAAGCG-3' (48), annealing temperature, 58°C; Cxcl1, forward, 5'-TGCACCCAAACCGAAGTCATAG-3', reverse, 5'-GTGGTTGACACTTAGTGGTCTC-3' (80), annealing temperature, 52°C; Scarb1, forward, 5'-TTTCAGCAGGATCCATCTGGTGGA-3', reverse, 5'-AGTTCATGGGGATCCCAGTGAC-3' (70), annealing temperature, 58°C; H2-D1, forward, 5'-CTCTCCTGTGACATCCAGAGC-3', reverse, 5'-TTGGCTATGGAAGGGAACAC-3' (21), annealing temperature, 55°C. Reaction parameters were as follows: 2 min at 50°C; 15 min at 95°C; 40 cycles of 15 s each at 94°C, 30 s at the appropriate annealing temperature, and 1 min at 72°C. Controls and data analysis were performed as previously described (50).
In situ hybridization.
35S-UTP-labeled riboprobes were generated using a Maxiscript kit (Ambion, Austin, TX) and a fragment of the cDNA of mouse Scarb1 (GenBank accession number BC004656) corresponding to nucleotides 6611043 of the coding sequence as template. Riboprobes were synthesized, sections were hybridized with probes, and reactions were visualized as previously described (78).
 |
RESULTS
|
|---|
Cellular changes following target ablation.
The cellular organization of the OE of MSR/ mice (Fig. 1A) appeared indistinguishable from that of wt mice (e.g., 50, 61). OSNs, identified by immunohistochemistry with OMP (Fig. 1B), comprised the lower 3/4 of the OE, with a thin two- to three-cell high band of progenitor cells at the base of the epithelium. As in wt mice (61), small numbers of resident m
s expressing the pan-m
marker CD68 occurred primarily in the lamina propria beneath the OE; very rare intraepithelial resident m
s were also observed in both wt (Fig. 1C) and MSR/ mice. OE height was comparable in sham wt and MSR/ mice (Fig. 2A). At 3 days after OBX, OE height was reduced by 57% in wt mice and 68% in MSR/ mice and was significantly thinner (P < 1.0E-03) in both strains compared with their respective shams (Figs. 1D, 2A), due primarily to the loss of most of the OMP-immunoreactive OSNs (Fig. 1E); the height of the OE was 26% less in MSR/ than in wt mice following OBX (P < 1.0E-03). The numbers of infiltrating m
s increased significantly in both strains, particularly in the OE (Fig. 1F), but there were notable differences. While the numbers of CD68+ m
s were initially comparable and increased significantly (P < 1.0E-03) following OBX in both wt and MSR/ mice, significantly fewer (P = 0.012) infiltrated the OE in MSR/ mice (2.9-fold increase) vs. wt mice (4.8-fold increase, Fig. 2B); thus there were 48% fewer macrophages in the OE of MSR/ vs. wt mice after OBX. The numbers of activated m
s, identified by immunoreactivity for F4/80, increased significantly in wt mice (P = 0.023, 3.3-fold), but there was no increase in their numbers in MSR/ mice (Fig. 2C). The number of proliferative BrdU+ progenitor cells was significantly higher initially in MSR/ sham vs. wt sham mice (P = 2.0E-03, 5.1 fold); following OBX, the number increased significantly in both wt and MSR/ mice but was significantly lower in MSR/ vs. wt mice (P < 1.0E-03, 1.9-fold). Another difference between the two strains was the apparent persistence of apoptotic bodies in the OE of MSR/ vs. wt mice (Fig. 3, A and B) as indicated by nuclear staining with bisbenzimide. There appeared to be more individual large apoptotic bodies and more numerous clusters of small apoptotic bodies, particularly in the turbinates, in MSR/ vs. wt mice at 3 days post-OBX. These results demonstrated substantial differences in OBX-induced cellular responses, including decreased numbers of recruited and activated m
s, reduced progenitor cell proliferation, and an apparent defect in the clearance of apoptotic OSNs in the OE of MSR/ vs. wt mice.

View larger version (94K):
[in this window]
[in a new window]
|
Fig. 1. Effects of olfactory bulbectomy (OBX) on olfactory epithelium (OE) of macrophage scavenger receptor (MSR)/ mice. A and D: cresyl violet staining; B and E: olfactory marker protein (OMP) immunohistochemistry (IHC); C and F: CD68 IHC. A, D: in OE of 0-h sham MSR/ mice (A), pattern of organization was qualitatively the same as in wild-type (wt) mice; lightly stained sustentacular cell nuclei lay just below epithelial surface, and a thick layer of darkly stained nuclei of olfactory sensory neurons (OSNs) occupied most of the OE height; duct of Bowman's gland (d) penetrated basement membrane. In OE of MSR/ mice at 3 days post-OBX (D), most OSN nuclei were gone, resulting in large decrease in OE height; apoptotic nucleus (arrow) was present in OSN layer. Sustentacular cell nuclear layer remained intact. B, E: in 0-h sham MSR/ mice (B), closely packed OMP+ OSN cell bodies (arrows) occupied most of the lower 3/4 of OE, with dendrites projecting to OE surface and a thin band of unstained progenitor cells at base. In MSR/ mice at 3 days post-OBX (E), a large reduction in the number of OMP+ cells (arrows) resulted in decreased OE height. C, F: in OE of 0-h sham MSR/ mice (C), an occasional intraepithelial CD68+ macrophage (arrow) was observed. In OE of MSR/ mice at 3 days post-OBX (F), intraepithelial CD68+ macrophages (arrows) were more numerous. Scale bar (in D), 15 µm for AF.
|
|

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 2. Comparison of effects of OBX in wt and MSR/ mice. A: OE height, which was equivalent in 0-h shams of strains (open bars), was significantly reduced compared with respective shams at 3 days post-OBX (striped bars) in both strains (P < 0.001). OE height in MSR/ mice was significantly less than in wt mice at 3 days post-OBX (P < 0.001). B: numbers of CD68+ macrophages recruited to the OE at 3 days post-OBX increased significantly in wt and MSR/ mice (P < 0.001), but significantly fewer were recruited in MSR/ mice (P < 0.01). C: the numbers of F4/80+ activated macrophages recruited to the OE at 3 days post-OBX increased significantly in wt mice (P < 0.05), but there was no change in their numbers in MSR/ mice. D: 5 times more 5-bromo-2'-deoxyuridine (BrdU)+ progenitor cells occurred in OE of 0-h sham MSR/ vs. wt mice (P < 0.05). Numbers of BrdU+ cells increased significantly in both strains (P < 0.001) at 3 days post-OBX, but increase in proliferating progenitor cells was significantly less (P < 0.001) in MSR/ mice. *P < 0.05, **P < 0.01, ***P < 0.001.
|
|

View larger version (149K):
[in this window]
[in a new window]
|
Fig. 3. Apoptotic bodies in the OE at 3 days post-OBX. DNA staining demonstrated apoptotic bodies (arrows) in the OE of both wt mice (A) and MSR/ mice (B). In MSR/ mice, single apoptotic bodies (arrows) appeared larger and more numerous, and clusters of small apoptotic bodies (brackets) were more numerous. Scale bar (in A), 15 µm for A and B.
|
|
Global and interactive gene analyses.
The flow chart (Fig. 4) outlines the analyses applied to the hybridization signals on the GeneChips to identify genes in the OE that had significantly different patterns of expression at at least one time point following OBX in wt and MSR/ mice. Among the probe sets for known genes, an Absolute Call analysis identified 6,458 (70%) as present on at least one chip. A 2 (mouse strains) x 5 (0-h sham OBX and 2, 8, 16, and 48 hpO time points) ANOVA performed on the Present probe sets identified 4,006 significantly (P < 0.05) regulated genes. The FDR correction set to a stringent 1% was performed to evaluate the approximate number of type 1 errors represented as false positives, which reduced the number of significantly regulated probe sets by
52% to 2,092 genes with P values of <3.2E-03.

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 4. Flow chart demonstrating statistical analysis of microarray data. Probe sets/genes included in analysis are shown in bold. FDR, false discovery rate; NSD, no significant difference; hpO, hours postolfactory bulbectomy. See text for details.
|
|
A major aim of this study was to identify genes regulated by OBX as a result of OSN apoptosis and the subsequent innate immune response, neurogenesis, and epithelial remodeling that exhibited different temporal patterns of expression in MSR/ and wt mice. We separated the significantly regulated genes into two groups based on their expression levels in the 0-h sham OBX mice of each strain. For
52% of the genes, there was no significant difference (NSD, P > 0.05) between the magnitudes of the mean hybridization signals at 0 h in the two strains;
48% exhibited significant differences (P < 0.05) between the 0-h shams (Fig. 4). Among the 1,095 genes with NSD at 0 h, 293 genes (27%) had a highly significant difference (P values ranging from 7.8E-13 to 3.2E-3) between the magnitudes of their mean hybridization signals in the two strains at at least one time point post-OBX; that is, these genes differed significantly in their temporal pattern of expression in the two strains following OBX. We further classified these genes by identifying the time point at which the first statistically significant difference in expression between the two strains occurred. Seventy-six genes showed their first highly significant difference in the two strains at 2 h, 25 genes at 8 h, 56 genes at 16 h, and 136 genes at 48 h following target ablation. For example, there was NSD at 0 h in the magnitude of expression of Cp (ceruloplasmin, Fig. 5A), an IL (interleukin)-1ß-inducible gene expressed by m
s and glia that is upregulated following retinal and brain injury (10, 44, 49), in the wt and MSR/ strains. In wt mice, Cp was significantly upregulated at 16 hpO (P = 3.4E-04, 2.2-fold) and at 48 hpO (P = 7.8E-05, 3.8-fold). The first statistically significant difference between the expression of Cp in the two strains occurred at 16 hpO (P = 1.0E-03) when the magnitude and direction of expression in the two strains diverged sharply and the gene was expressed at a 1.7-fold higher level in wt than in MSR/ mice. At 48 hpO, Cp was again expressed at a significantly higher level in wt than in MSR/ mice (P = 4.2E-03, 2.5-fold). Of the remaining genes with NSD at 0 h, 337 (31%) had a significant difference (P values ranging from 3.3E-3 to 0.049) between the magnitudes of their mean hybridization signals in the two strains at at least one time point post-OBX. The remaining genes (465 genes, 42%) were main-effect genes; that is, their patterns of expression were parallel in the two strains, indicating that these genes responded to OBX with similar changes in their magnitude and direction of expression over time. For example, there was NSD in the magnitude of expression at 0 h of Bmp6 (bone morphogenetic protein 6, Fig. 5B), a member of the transforming growth factor (TGF)-ß superfamily of growth factors that is expressed by mature OSNs (63). In wt mice, Bmp6 was significantly downregulated at both 16 hpO (P = 0.013, 1.5-fold) and at 48 hpO (P = 0.017, 2.0-fold), and there was a parallel downregulation of expression at each time point in both wt and MSR/ mice, showing the main effect of treatment at all post-OBX time points.
Among the 997 genes for which there were significant differences (P < 0.05) at 0 h between the two strains, 283 genes (28%) exhibited highly significant differences between strain and time (P values ranging from 2.2E-09 to 3.1E-3). One hundred eighty-two genes showed their first statistically significant difference in expression at 2 h, 15 genes at 8 h, 40 genes at 16 h, and 46 genes at 48 h following target ablation. Of the remaining genes with a significant difference at 0 h, 369 (37%) had a significant difference (P values ranging from 3.3E-3 to 0.049) between the magnitudes of their mean hybridization signals in the two strains at at least one time point post-OBX. The remaining genes (345 genes, 35%) were main-effect genes.
In summary, 2,092 genes were differentially regulated (1% FDR, P = < 0.0032) in the two strains, of which about half exhibited either NSD or a significant difference in their initial level of expression in the two strains. About 28% of the differentially regulated genes displayed a temporal pattern of expression following OBX that showed highly significant differences (P < 3.2E-03) between wt and MSR/ mice.
Categorical gene analysis.
The genes used in Fig. 5 to illustrate the statistical categories of genes with a significant interaction between strain and time and main-effect genes suggest that these statistical categories cut across functional categories. Functional overrepresentation in the genes with NSD at 0 hpO and a significant difference at 0 hpO were analyzed separately using the EASE program and are listed in Supplemental Table 1. (The online version of this article contains supplemental data.) Genes with an SD at 0 hpO included intracellular genes related to ribosomes and mitochondrial electron transport. These categories are consistent with the biological and molecular categories of protein biosynthesis and RNA binding, and nucleoside triphosphate metabolism, hydrogen ion transporter activity, and oxidoredutase/antioxidant activity, respectively. The mitochondrial electron transport chain and hydrogen ion transporter activity/oxidoreductase activity categories contain many of the same genes, including those in the NADH dehydrogenase (ubiquinone), ubiquinol cytochrome c reductase, and cytochrome c oxidase families. In the hydrogen ion transporter activity category, 18/20 genes that are expressed exclusively in the mitochondrial electron transport chain were expressed at significantly higher levels in MSR/ mice, and in the antioxidant category, 7/11 genes, including peroxiredoxins 14, were expressed at higher levels, indicating the occurrence of oxidative stress in the OE of MSR/ mice. Another overrepresented category, regulation of apoptosis, contained primarily proapoptotic genes such as caspases 2, 3, and 9 that were expressed at significantly lower levels in MSR/ mice; however, apoptotic nuclei were not detected in the OEs of sham mice of either strain. At 2 hpO, genes with significant interactions in both groups (NSD or significant difference at 0 hpO) were primarily intracellular. Among genes with NSD at 0 hpO with a significant interaction at 2 hpO were genes associated with the endoplasmic reticulum in the secretory pathway (e.g., Rab6, Ssr1). The hydrolase activity and GTPase molecular function categories contained five of the same genes, including Rab6 as well as two genes of the Ras family and their related G proteins that are involved in cytoskeletal modulation. Among genes with a significant difference at 0 hpO and a significant interaction at 2 hpO were several transcription factors in the ribonucleoprotein complex/nucleotide/nucleic acid metabolism categories and several mitochondrial complex I NADH dehydrogenase flavoproteins in the NADH dehydrogenase/electron transporter categories. These results indicate that OBX-induced gene transcription and energy-dependent processes were regulated differently in the two strains. Of particular relevance to this study was the overrepresentation of genes with a significant difference at 0 hpO and a significant difference in expression levels at 2 hpO that are associated with the regulation of apoptosis. Among the genes in this category are Cebpb (CCAAT/enhancer binding protein-ß), Pdcd6 (programmed cell death 6), Rnf7 (ring finger protein 7), and Traf3 (TNF receptor-associated factor 3). Other apoptosis genes that were previously reported to be regulated by OBX in wt mice (25), such as Casp2 and Bax (see Apoptosis genes below), also show this pattern. No categories related to apoptosis were overrepresented at later time points, suggesting that the regulation of apoptosis genes by OBX was similar in the two strains by 1648 hpO.
Other genes of particular relevance to this study also appear in the overrepresented categories. Il1rn (IL-1 receptor antagonist) is in the growth factor activity category, probably relating to its role in the regulation of insulin and lipase levels, among genes with NSD at 0 hpO and a significant interaction at 8 hpO. Among the genes with NSD at 0 hpO that had a significant interaction at 16 hpO, Grn (granulin) encodes a growth-promoting protein that interacts with the product of Slpi (secretory leukocyte protease inhibitor) to promote epithelial wound healing. Among the genes with NSD at 0 hpO that had a significant interaction at 48 hpO, genes related to the extracellular space were significantly overrepresented, most likely related to the overrepresentation of genes in the immune response, chemotaxis, and chemokine activity biological/molecular categories that also had significant EASE scores at this time point. Genes with related functions in other over-represented categories at 48 hpO included Il1r2 (IL-1 receptor, type II) and Il1rl1 (IL-1 receptor-like 1) in the growth factor binding category; and Slpi in the endopeptidase/protease inhibitor activity category.
Although EASE analysis provides insight into functional categories of genes that were regulated differently in the two strains, the overrepresentation of a particular functional category does not provide information on the relative magnitude of the expression levels or direction of these changes. Therefore, we examined the detailed temporal expression patterns of genes whose expression level in wt mice was significantly regulated (P < 0.05) at one or more time points after OBX and that were classified in the following categories: 1) genes expressed in the OE primarily by OSNs; 2) genes previously associated with apoptotic pathways involved in OSN cell death; and 3) because SR-A is expressed primarily by m
s and endothelial cells, genes associated with cytokine and chemokine secretion and response, diapedesis, and phagocytosis that are relevant to OE degeneration and remodeling.
Olfactory-related genes.
As we previously reported (25), genes expressed only or primarily by OSNs showed a characteristic pattern of expression following OBX in wt mice, with significantly lower levels of gene expression beginning at 16 hpO for many genes (e.g., Bmp6, Fig. 5B; Omp, Fig. 6A) and for all by 48 hpO [e.g., Uchl1 (ubiquitin carboxy-terminal hydrolase L1) or PGP9.5, a neuronal marker; Fig. 6B]; this pattern indicates the loss of OSNs by apoptosis. In general, in MSR/ mice, these genes exhibited a sharp decrease in mean hybridization signals at 2 hpO followed by an increase to near 0-h sham levels at 8 hpO; at 16 hpO, the level of expression was greater in MSR/ than in wt mice but not significantly greater than in 0 hpO sham MSR/ mice. At 48 hpO, there was no significant difference between the two strains in the levels of expression for most of these genes. This pattern occurred in 15/19 genes expressed primarily by OSNs that we previously identified as regulated by OBX. Nine of the 19 OSN genes (Supplemental Table 2) were main-effect genes.

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 6. Expression patterns in response to OBX for genes expressed primarily by olfactory sensory neurons (OSNs; A, B) and genes involved in apoptosis (C, D). See Fig. 5 for symbols and text for details. A, B: mean hybridization signals for genes encoding olfactory marker protein (Omp, A) and ubiquitin carboxy-terminal hydrolase L1 (PGP9.5, Uchl1, B) were significantly decreased in wt mice at 16 and/or 48 hpO compared with 0-h shams, indicating OSN apoptosis. Both genes showed significantly different patterns of expression in MSR/ mice. C: mean hybridization signals for genes encoding BH3-interacting domain death agonist (Bid, C) and Bcl2-associated X protein (Bax, D) were significantly increased in wt mice at 16 and 48 hpO compared with 0-h shams and were expressed at similar levels at these time points in wt and MSR/ mice, indicating apoptotic progression in OSNs.
|
|
Apoptosis genes.
Because OBX causes apoptosis primarily of mature OSNs, we anticipated that relatively few apoptosis genes would exhibit significant differences in temporal patterns of expression in the two strains. Nine of the 17 apoptosis genes (Supplemental Table 3) were main-effect genes, including all of the caspases that were regulated by OBX (25). For example, the proapoptotic mediator Casp3 (caspase 3) was a main-effect gene with a significant difference at 0 h, and a significant increase in expression at 8 hpO (P = 0.012, 1.3-fold), 16 hpO (P = 2.0E-04, 1.5-fold), and 48 hpO (P = 6.0E-04, 1.4-fold) in wt mice, with a highly similar pattern in MSR/ mice. Another main-effect gene, proapoptotic Bid (BH3 interacting domain death agonist; Fig. 6C), had NSD at 0 h and significant increases in expression at 16 hpO (P = 8.5E-04, 1.8-fold) and 48 hpO (8.0E-06, 2.4-fold) in wt mice, with a highly similar pattern of expression in MSR/ mice. A downstream target of Bid, Ctsd (cathepsin D), which was significantly upregulated in wt mice at 48 hpO (P = 4.8E-03, 1.5-fold), was also a main-effect gene. Among apoptosis genes that showed different patterns of expression in the two strains was Bax (Bcl2-associated X protein, Fig. 6D), a proapoptotic regulator that significantly increased in expression in wt mice at 16 hpO (P = 0.012, 1.2-fold) and 48 hpO (P = 1.8E-03, 1.3-fold) and Casp2 (caspase 2), a negative regulator of apoptosis in neurons; notably, the levels of expression of these genes at 16 and 48 hpO was equivalent in the two strains. Apoptosis genes with highly significant differences in expression levels at 16 and/or 48 hpO were all expressed at lower levels in MSR/ mice than in wt mice.
Cytokines.
Our microarray results identified a potential candidate for an OSN-derived cytokine that may initiate the local immune response to OBX: IL-1ß. At 2 hpO, there was a significant upregulation of Il1b (IL-1ß, Fig. 7A) in both wt (P = 0.013, 2.1-fold,) and MSR/ mice (P = 0.033, 2.5-fold). This was consistent with the upregulation of the IL-1ß-inducible Ccl2 [chemokine (C-C motif) ligand 2, monocyte chemoattractant protein 1 (MCP-1)] mRNA by 8 hpO (see Chemokines below) and protein by 16 hpO (27). In contrast to the early parallel pattern of expression of Il1b in the two strains, there was a significant difference at 48 hpO (P = 9.8E-06, 8.5-fold), with an increase in expression in wt mice and a decrease in MSR/ mice. This was also consistent with the pattern of regulation of Ccl2 (see Chemokines below). Il1r2 (Fig. 7B), an IL-1 decoy receptor whose expression is upregulated concomitantly with that of IL-1ß, showed a similar early parallel pattern of expression in the two strains, followed by a significant difference in expression levels at 48 hpO (P = 1.4E-04, 3.8-fold), with an increase in expression in wt mice and a decrease in MSR/ mice. The expression pattern of Il1a (IL-1
, Fig. 7C) differed from that of Il1b in two ways: there was no upregulation in either strain at 2 hpO, and, at 48 hpO, the level of Il1a expression increased significantly in both strains (wt mice, P = 6.3E-04, 4.3-fold; MSR/ mice, P = 0.013, 2.1-fold), although it was expressed at a significantly lower level than in wt mice (P = 3.6E-05, 2.0-fold). In wt mice, Il1rn (IL-1 receptor antagonist, Fig. 7D) was significantly upregulated at 48 hpO (P = 2.1E-04; 5.2-fold); in MSR/ mice, the gene was significantly upregulated at 8 hpO (P = 0.011, 3.4-fold), differing significantly from its expression level in wt mice (4.2E-04, 3.2-fold). At 48 hpO, there was no upregulation of expression as seen in the wt mice, resulting in a significantly lower level of expression in the MSR/ mice (P = 4.5E-03, 14.3-fold). The two main-effect IL-1-related genes, Tollip (Toll interacting protein) and Tom1 (target of myb1 homolog, Supplemental Table 4), which interact with each other and are negative regulators of IL-1 signaling, were significantly downregulated at 48 hpO in wt mice (Tollip, P = 0.023, 1.4-fold; Tom1, P = 0.015, 1.4-fold).

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 7. Expression patterns in response to OBX for interleukin (IL)-1-related genes. See Fig. 5 for symbols and text for details. A: mean hybridization signals for Il1b (IL-1ß) increased significantly at 2 and 48 hpO in wt mice; in MSR/ mice, mean hybridization signals increase at 2 hpO but not at 48 hpO. B: mean hybridization signals for Il1r2 (IL-1 receptor, type II) increased in parallel in the 2 strains at 2 hpO, but the significant increase in wt mice at 48 hpO did not occur in MSR/ mice. C: the expression pattern of Il1a (IL-1 ) was different from that of Il1b in that there was no significant increase in expression until 48 hpO in wt mice. In MSR/ mice, mean hybridization signal increased but to a significantly lower level than in wt mice. D: mean hybridization signals for Il1rn (IL-1 receptor antagonist) increased significantly at 48 hpO in wt mice; in MSR/ mice, mean hybridization signals increase significantly over those of wt mice at 8 hpO but at 48 hpO decreased to significantly lower level.
|
|
Other m
-derived cytokines whose gene expression patterns exhibited significant differences between strain and time included members of the neuropoietic cytokine, interferon-
/ß, and complement families, and several related transcription factors and signaling molecules (Supplemental Table 4), and one TGF-ß-related anti-inflammatory cytokine that is associated with phagocytosis (Tgfbi, TGF-ß induced; Supplemental Table 6). Of the 25 cytokine and cytokine-related genes with significant differences in expression patterns between the two strains at one or more time points post-OBX, most were significantly upregulated in wt mice at 16 hpO (6 genes) and 48 hpO (11 genes). In MSR/ mice, most of these genes were not upregulated at 48 hpO, with near-sham levels of expression of six of eight IL-1-related genes, five of eight neuropoietic cytokine-related genes, one of three complement-related genes, two of three signaling molecules, and the TGF-ß-related gene.
Among genes in the neuropoietic or IL-6-type cytokine family with patterns of expression that showed significantly different patterns of expression in the two strains, Lifr (leukemia inhibitory factor receptor, Fig. 8A) was significantly upregulated at 48 hpO in wt mice (P = 0.019, 3.0-fold), whereas in MSR/ mice, its level of expression was not increased over 0-h sham levels, resulting in a significantly lower level of expression than in wt mice (P = 3.3E-04, 3.2-fold). For Il6 (IL-6, Fig. 8B), the pattern of expression in the two strains first differed significantly at 2 hpO (P = 0.026, 4.4-fold), when its expression was significantly upregulated in MSR/ mice but not wt mice. In wt mice, Il6 was upregulated significantly at 48 hpO (P = 1.2E-03, 3.0-fold), but in MSR/ mice, it was expressed at 0-h sham levels that were significantly lower than in wt mice (P = 0.011, 26.7-fold). The expression pattern of Osmr (oncostatin M receptor, Fig. 8C) differed from that of both Lifr and Il6 in that it was upregulated (significantly in MSR/ mice) at 8 hpO in both strains (wt, 2.8-fold; MSR/, P < 1.0E-04, 2.8-fold). The pattern of expression between the strains differed at 48 hpO, when in wt mice, Osmr expression was significantly upregulated (P = 0.012, 3.9-fold) to a level significantly higher than in MSR/ mice (P = 1.5E-03, 2.4-fold). The receptors through which these cytokines signal share a common element, gp130, encoded by Il6st (IL-6 signal transducer), which was upregulated significantly at 48 hpO in wt mice (P = 3.4E-03, 1.2-fold) and was a main-effect gene. Socs3 (suppressor of cytokine signaling 3, Fig. 8D), a regulator of signaling through gp130-containing receptors, was expressed in the three strains in a pattern resembling that of Osmr. At 8 hpO, there was a significant increase in expression over 0-h sham levels in both wt mice (P = 7.4E-03, 1.6-fold) and MSR/ mice (P = 0.011, 2.0-fold). At 48 hpO, there was a significant increase in expression in wt mice (P = 6.2E-03, 2.3-fold) that was significantly higher than in MSR/ mice (P = 2.9E-05, 3.1-fold).

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 8. Expression patterns in response to OBX for neuropoietic cytokine-related genes. See Fig. 5 for symbols and text for details. A: mean hybridization signal for Lifr (leukemia inhibitory factor receptor) increased significantly at 48 hpO in wt mice; in MSR/ mice, there was no significant increase in expression at 48 hpO. B: mean hybridization signals for Il6 (interleukin 6) increased significantly at 48 hpO in wt mice; they first differed significantly in the 2 strains at 2 hpO when its expression was significantly upregulated in MSR/ but not wt mice, and at 48 hpO when its expression was significantly upregulated in wt but not MSR/ mice. C: expression pattern of Osmr (oncostatin M receptor) differed from Lifr and Il6, being upregulated (significantly in MSR/ mice) at 8 hpO in both strains, but at 48 hpO, mean hybridization signals were significantly upregulated in wt but not MSR/ mice. D: similarly, Socs3 (suppressor of cytokine signaling 3) expression was significantly upregulated at 8 hpO in both strains, but only in wt mice at 48 hpO.
|
|
Chemokines.
The expression patterns of four CC chemokines and three CXCL chemokines differed significantly in wt and MSR/ mice (Supplemental Table 5). Ccl2 (Fig. 9A) showed a biphasic expression pattern in wt mice, with significant upregulation at 8 hpO (P = 1.6E-04, 14.7-fold), 16 hpO (P = 9.1E-03, 9.4-fold increase), and 48 hpO (P = 5.0E-04, 23.7-fold). In MSR/ mice, the pattern of expression was highly similar at 8 and 16 hpO, but at 48 hpO there was a significant difference in expression levels between the 2 strains (P = 2.2E-05, 6.1-fold). In wt mice, the other CC chemokines with highly significant differences in expression patterns, Ccl3 [chemokine (C-C motif) ligand 3, macrophage inflammatory protein (MIP)-1
] and Ccl6 [chemokine (C-C motif) ligand 6, c10] showed no change in expression levels until significant upregulation at 48 hpO (Ccl3, P = 3.4E-03, 5.2-fold; Ccl6, P = 0.011, 1.4-fold). In MSR/ mice, the expression of Ccl3 and Ccl6 remained at 0-h sham levels at all time points, resulting in significant differences in expression levels in the two strains (Ccl3, P = 8.1E-06, 3.8-fold; Ccl6, P = 0.018, 2.0-fold). The CXC chemokines showed similar patterns of expression to those of the CC chemokines. In wt mice, the expression levels of Cxcl1 [chemokine (C-X-C motif) ligand 1, Gro1 oncogene] and Cxcl2 [chemokine (C-X-C motif) ligand 2, Gro2, Fig. 9C] were significantly upregulated at only 48 hpO (Cxcl1, P = 2.7E-03, 7.5-fold; Cxcl2, P = 2.4E-03, 51.5-fold). In MSR / mice, the level of expression of Cxcl1 at 2 hpO was significantly greater than that in wt mice (P = 0.031, 5.7-fold), and for both, there were significant differences at 48 hpO with the wt strain (Cxcl1, P = 6.6E-03, 16.0-fold; Cxcl2, P = 1.9E-07, 42.7-fold).

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 9. Expression patterns in response to OBX for chemokine genes and vascular adhesion-related genes. See Fig. 5 for symbols and text for details. AC: chemokine genes. Mean hybridization signals for the chemokine Ccl2 [chemokine (C-C motif) ligand 2, monocyte chemoattractant protein (MCP)-1; A] increased significantly at 8 hpO in both strains and remained significantly elevated at 16 and 48 hpO in wt mice; in MSR/ mice, expression decreased at 16 and 48 hpO. Mean hybridization signals for chemokines Ccl3 [chemokine (C-C motif) ligand 3, macrophage inflammatory protein (MIP)-1 ; B] and Cxcl2 [chemokine (C-X-C motif) ligand 2, MIP-2; C] show significant increases in wt mice only at 48 hpO and no significant changes in MSR/ mice, resulting in significantly different expression patterns for all 3 chemokines only at 48 hpO. D and E: endothelial cell adhesion molecules Selp (selectin, platelet; D) and Vcam1 (vascular cell adhesion molecule 1, E) were significantly upregulated at 48 hpO in wt mice but not in MSR/ mice, resulting in significantly different expression patterns between the 2 strains. F: similarly, mean hybridization signals for S100a8 (S100 calcium binding protein A8), associated with diapedesis-related cytoskeletal remodeling, similarly increased significantly in wt mice at 48 hpO but not in MSR/ mice, resulting in significantly different expression patterns in the 2 strains.
|
|
Diapedesis-related genes.
Several genes encoding key adhesion molecules and their receptors as well as cytoskeletal components and associated signaling pathways associated with the diapedesis of recruited monocytes (10, 76) showed significantly different patterns of expression after OBX in wt and MSR/ mice (Supplemental Table 5). For example, Selp (P-selectin, Fig. 9D) was significantly upregulated only in wt mice at 48 hpO, when there was a significant increase in expression (P = 0.011, 8.6-fold); the lack of upregulation in MSR/ mice resulted in a significant difference in expression levels between the strains (P = 3.0E-05, 10.3-fold). Two genes encoding monocyte adhesion molecules were regulated by OBX, including Itgb2 (integrin-ß2, LFA-1) whose significant upregulation at 48 hpO in wt mice (P = 4.6E-03, 1.7-fold) was previously reported (25), were main-effect genes. In wt mice, Vcam1 (vascular cell adhesion molecule 1, Fig. 9E), an endothelial adhesion molecule, was significantly upregulated at 48 hpO (P = 4.0E-04, 2.5-fold) as previously reported (25); in MSR/ mice, its expression level remained at 0-h sham levels, resulting in a significant difference in expression between the two strains (P = 1.2E-04, 1.9-fold). The receptor for VCAM-1, Icam1 (intercellular adhesion molecule 1), was significantly upregulated at 48 hpO in wt mice (P = 1.1E-03, 3.1-fold); in MSR/ mice, Icam1 was expressed at significantly lower levels at both 16 hpO (P = 0.024, 2.8-fold) and 48 hpO (P = 9.0E-04, 3.4-fold).
Genes related to cytoskeletal remodeling were also differentially expressed in the two strains, including numerous genes encoding components of the actin cytoskeleton and signaling molecules associated with its regulation (not shown, see DISCUSSION); these patterns were consistent with the lower expression levels of Vcam1 and Icam1, whose interactions with their integrin ligands regulate cytoskeletal remodeling, in MSR/ vs. wt mice. The S100 calcium-binding proteins S100a8 (Fig. 9F) and S100a9, which regulate the interaction between microtubule and actin cytoskeletal remodeling pathways, were also differentially expressed in the two strains. In wt mice, both genes showed no significant changes in expression levels until 48 hpO, when both were significantly upregulated (S100a8, P = 2.1E-03, 4.0-fold; S100a9, P = 4.5E-03, 3.8-fold); both genes remained at 0-h sham levels in MSR/ mice, significantly lower than in wt mice (S100a8, P = 9.7E-03, 3.8-fold; S100a9, P = 0.016, 4.3-fold).
Phagocytosis: engulfment of apoptotic cells by macrophages.
Numerous phagocytosis-related genes are regulated by OBX (Supplemental Table 6). An enzyme that regulates the exposure of phosphatidylserine groups on the surface of apoptotic cells, Plscr1 (phospholipid scramblase 1, Fig. 10A), was significantly upregulated at 8 hpO (P = 1.3E-04, 2.0-fold) in wt mice, and its expression remained significantly elevated through 48 hpO (16 hpO, P = 8.7E-03, 1.6-fold; 48 hpO, P = 4.7E-04, 2.6-fold). In MSR/ mice, Plscr1 was expressed at significantly higher levels at 0 h than in wt mice; its expression at 8 hpO was significantly upregulated (P = 9.4E-03, 1.6-fold) as in wt mice, but by 48 hpO, it was expressed at a significantly lower level than in wt mice (P = 4.9E-03, 1.9-fold). Another eat-me signal on apoptotic cells, Anxa1 (annexin A1), was expressed at comparable levels in the shams of both strains (Fig. 10B
). In wt mice, Anxa1 was significantly upregulated at 8 hpO (P = 2.1E-03, 1.9-fold) and 16 hpO (P = 0.018, 1.6-fold), with a parallel pattern of expression in MSR/ mice. However, at 48 hpO, it was expressed at a significantly lower level in MSR/ mice (P = 2.1E-06, 2.0-fold).

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 10. Expression patterns in response to OBX for engulfment signal genes and scavenger receptors. A, B: "eat me" signals. Plscr1 (phosholipid scramblase 1, A) was significantly upregulated at 8, 16, and 48 hpO in wt mice; it was expressed at significantly higher level in 0 h sham MSR/ mice and was upregulated as in wt mice at 8 and 16 hpO but decreased significantly at 48 hpO. Similarly, Anxa1 (annexin AI, B) was significantly upregulated at 8, 16, and 48 hpO in wt mice, with similar upregulation at 8 and 16 hpO in MSR/ mice, but decreased significantly at 48 hpO. C, D: scavenger receptors. Mean hybridization signals for Scarb1 (scavenger receptor class B, member 1; C) increased significantly at 48 hpO in wt mice. Expression pattern was significantly different in MSR/ mice, with significantly lower expression at 2 hpO and significantly higher expression at 48 hpO. The scavenger receptor Cd68 (CD68, macrosialin; D) was significantly upregulated at 48 hpO in wt mice, with similar pattern of expression in MSR/ mice.
|
|

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 11. Expression patterns in response to OBX for phagosome genes, anti-inflammatory mediators, and macrophage (m )-specific genes. A, B: phagosome genes. OBX upregulated expression of Capg (gelsolin-like actin filament capping protein, A) in wt mice at 16 and 48 hpO. In MSR/ mice, Capg was not upregulated at 48 hpO. H2-D1 (B), which encodes a major compatibility class I antigen, was significantly upregulated at 16 and 48 hpO in wt mice, but was downregulated at these time points in MSR/ mice. Anti-inflammatory mediator gene Tgfbi (TGF-ß induced, C) was upregulated in wt mice at 16 and 48 hpO; its pattern of expression in MSR/ mice was significantly different at both 16 and 48 hpO. Slpi (secretory leukocyte protease inhibitor, D), another anti-inflammatory mediator gene, was significantly upregulated by OBX at 48 hpO in wt but not MSR/ mice. E, F: m genes. Arg1 (arginase 1, E), an indicator of alternative m activation, was significantly upregulated at 48 hpO in wt mice but not in MSR/ mice. M maturation indicator Chi3l1 (chitinase 3-like 1, F) was significantly upregulated in wt mice at 48 hpO; its pattern of expression in MSR/ mice was significantly different at both 16 and 48 hpO.
|
|
Among molecules that bridge apoptotic cells to phagocytes, Mfge8 (milk fat globule-EGF factor 8 protein) was a main-effect gene that was significantly upregulated in both strains at 48 hpO (wt mice, P = 0.017, 1.3-fold; MSR/ mice, P = 6.9E-03, 1.3-fold). The three bridging molecules regulated by OBX that had significantly different patterns of expression in the two strains differed significantly at 48 hpO, consistent with expression patterns in eat-me signals.
Engulfment receptors on m
s include several classes of scavenger receptors. In wt mice, Msr1 was upregulated at 48 hpO with marginal significance (P = 0.050, 1.6-fold); the upregulation was validated by real-time RT-PCR (P < 0.001, 25.0-fold; Fig. 12), indicating the participation of SR-A in the engulfment of apoptotic OSNs. Several other scavenger receptor genes were significantly upregulated in wt mice, including Cd36 (Cd36), a class B scavenger receptor, at 16 hpO, Lgals3bp (galectin-3 binding protein) and Stab1 (stabilin 1) at 16 and 48 hpO, and Scarb1 (scavenger receptor class B, member 1, Fig. 10C) and CD68 (Cd68, Fig. 10D), a class D scavenger receptor, at 48 hpO. Among the scavenger receptors that were regulated by OBX in wt mice, the expression level of only Scarb1 (Fig. 10C) was significantly upregulated in MSR/ compared with wt mice at 48 hpO (P = 0.022, 1.3-fold), suggesting that SR-B1 partially compensated for the absence of SR-A. We confirmed the upregulated expression of Scarb1 at 48 hpO with in situ hybridization (Fig. 13). Cd68 (Fig. 10D), which was significantly upregulated in both wt and MSR/ mice at 48 hpO (wt mice, P = 5.9E-03, 1.8-fold; MSR/ mice, P = 8.6E-03, 1.6-fold), was a main-effect gene.

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 12. Real-time RT-PCR analysis of Msr1 expression in olfactory epithelium. In wt mice, Msr1 was expressed at relatively low level in 0-hpO shams (open bars) that was strongly and significantly upregulated by OBX at 48 hpO (striped bars). In contrast, Msr1 expression was extremely low in MSR/ mice and was not significantly upregulated at 48 hpO, resulting in a significant difference in expression levels in the 2 strains. ***P = 0.001.
|
|

View larger version (125K):
[in this window]
[in a new window]
|
Fig. 13. Validation of Scarb1 upregulation by OBX in OE of MSR/ mice with in situ hybridization. A: section of nasal tissue from 0-h sham stained with Giemsa stain. OE is located between parallel rows of asterisks. B: autoradiograph of same section shows faint hybridization signals for Scarb1 over OE. C: section of nasal tissue at 48 hpO stained with Giemsa stain. OE thickness between parallel rows of asterisks is noticeably reduced. D: autoradiograph of same section shows strong hybridization signals for Scarb1 over OE. L, lumen of nasal cavity; S, nasal septum. Scale bar (in A), 50 µm for AD.
|
|
The genes encoding numerous phagosomal proteins were regulated by OBX and many showed significant differences in their expression patterns in the two strains. Among these were several actin- and tubulin-associated genes associated with cytoskeletal rearrangement. For example, Capg (gelsolin-like actin filament capping protein, Fig. 11A), was significantly upregulated in wt mice at 16 hpO (P = 0.013, 1.8-fold) and 48 hpO (P = 2.1E-03, 2.8-fold). In MSR/ mice, Capg was expressed at a significantly lower level at 48 hpO than in wt mice (P = 0.012, 2.3-fold). Similarly, Arhgdia (Rho GDP dissociation inhibitor-
) was also significantly upregulated in wt mice at 16 hpO (P = 7.4E-03, 1.3-fold) and 48 hpO (P < 1.0E-04, 1.6-fold); in MSR/ mice, Arhgdia was expressed at significantly lower levels at both 16 and 48 hpO (16 hpO, P = 0.046, 1.2-fold; 48 hpO, P = 0.033, 1.3-fold). H2-D1 (histocompatibility 2, D region locus 1, Fig. 11B), an antigen of the major histocompatibility class I associated with phagosomes, and Lgals3 (galectin-3), which localizes in forming and mature phagosomes, showed similar patterns of expression. H2-D1 was significantly upregulated at 16 hpO (P = 7.4E-03, 2.0-fold) and 48 hpO (P = 0.012, 2.0-fold) in wt mice. In MSR/ mice, it was expressed at significantly lower levels at 16 hpO (P = 1.1E-04, 2.4-fold) and 48 hpO (P = 8.8E-03, 1.7-fold).
Secretion of TGF-ß, the prototypical anti-inflammatory mediator released when m
s engulf apoptotic cells, was indicated by the upregulated expression of at least two genes induced by TGF-ß (e.g., Tgfbi, Fig. 11C; the growth factor-related Tgfbi4, not shown). In wt mice, Tgfbi, a mediator of epithelial and endothelial cell adhesion, was significantly upregulated at 16 hpO (P = 0.012, 1.9-fold) and 48 hpO (P = 6.5E-04, 3.3-fold). In contrast, in MSR/ mice, there was no upregulation of the gene, with resultant significant differences in expression compared with wt mice at both time points (16 hpO, P = 8.0E-03, 1.6-fold; 48 hpO, P = 1.7E-03, 3.0-fold). Slpi (Fig. 11D), a serine protease inhibitor with anti-inflammatory activity whose secretion is induced by the binding of m
s to apoptotic cells, was significantly upregulated over its levels in shams in both wt and MSR/ mice at 48 hpO (wt, P = 1.2E-03, 10.7-fold; MSR/, P = 2.2E-03, 11.2-fold), but there was a statistically significant difference between the levels of expression in the two strains (P = 4.8E-05, 2.2-fold), consistent with reduced m
activation in the MSR/ mice. Fpr1, formyl peptide receptor 1, a receptor for the anti-inflammatory lipoxin A4 that also interacts with annexin I to mediate anti-inflammatory activity, was significantly upregulated at 16 hpO (P = 1.7E-02, 2.4-fold) and at 48 hpO (P = 1.7E-03, 5.2 fold) in wt mice, but its expression level did not change in MSR/ mice, resulting in significant differences at both time points (16 hpO, P = 0.034, 3.0-fold; 48 hpO, P = 0.013, 2.9-fold).
Other macrophage genes.
In addition to chemokines, cytokines, and phagocytosis mediators, m
s synthesize other characteristic products whose gene expression levels may indicate their maturation or activation state; those regulated by OBX are listed in Supplemental Table 7. For example, Ch3l1 (chitinase 3-like 1, Fig. 11F), a gene induced as m
s near maturation from monocytes, was significantly upregulated at 48 hpO in wt mice (P = 4.1E-03, 3.5-fold); in MSR/ mice, the gene was expressed at significantly lower levels at both 16 hpO (P = 9.7E-03, 2.3-fold) and 48 hpO (P = 9.1E-03, 4.8-fold). Other genes induced as m
s mature, such as Cd53 and S100a9, were also expressed at significantly lower levels in MSR/ compared with wt mice. Activation state may also be indicated by gene expression levels. M
s metabolize arginine via two pathways: by using inducible nitric oxide synthase (NOS2) to yield the highly reactive NO as well as L-citrulline (classical activation), and by using arginase 1 to yield ornithine, a precursor of growth-promoting factors, and urea (alternative activation). In both wt and MSR/ mice, Nos2 was "always absent" at all time points. In contrast, Arg1 (Fig. 11E) was upregulated in wt mice at 48 hpO (P = 0.014, 7.6-fold), indicating that OE m
s exhibited an alternatively rather than classically activated phenotype in response to OBX (50). In MSR/ mice, Arg1 was expressed at significantly lower levels than in wt mice at 48 hpO (P = 4.9E-05, 44.6-fold), which, in the absence of Nos2 upregulation, suggests lack of m
activation rather than classical activation. Clecsf8 (C-type lectin, superfamily member 8), an endocytic receptor expressed primarily by activated m
s, showed the same pattern of regulation and of differences between the two strains.
DC genes.
M
s and DCs both express SR-A, are found at mucosal surfaces, phagocytose apoptotic cells, and express an overlapping repertoire of activation markers, cytokines, chemokines, and other function-related molecules; therefore, it is possible that some differences in the response to OBX in wt and MSR/ mice were due to differences in DC function. To evaluate this possibility, the expression profiles of two genes that are markers for DCs were analyzed. Ly75 (lymphocyte antigen 75, DC205) and Adam19 (a disintegrin and metallopeptidase domain 19) were expressed in the OE and were regulated by OBX (Supplemental Table 8). In wt mice, Ly75, a marker of immature DCs, was significantly upregulated at 48 hpO (P = 5.5E-03, 1.5-fold); in MSR/ mice, Ly75 was not regulated by OBX and at 48 hpO was expressed at significantly lower levels than in wt mice (P = 0.037, 1.4-fold), suggesting that the absence of SR-A affected DC function in response to OBX. Adam19 was significantly upregulated at 16 hpO in wt mice (P = 0.043, 1.4-fold) and was a main-effect gene.
Secretory immune system genes.
Activated B cells that produce polymeric IgA are components of the secretory immune system in the olfactory mucosa (23, 57). The expression of Pigr (polymeric immunoglobulin receptor), the transporter of polymeric IgA in Bowman's glands, was significantly upregulated in wt mice at 48 hpO (P = 4.9E-03, 2-fold); in MSR/ mice, it was expressed at a significantly higher level at 2 hpO (P = 0.034, 1.8-fold) and at a significantly lower level at 48 hpO (P = 0.030, 1.8-fold). Igj (immunoglobulin joining chain), a B lymphocyte product that links molecules of IgA into pentamers, was a main-effect gene that was upregulated in wt mice at 2 hpO (P = 0.025, 3.8-fold).
Molecular validation using real-time RT-PCR.
Using real-time RT-PCR to determine the expression levels of Msr1, we confirmed its significant upregulation by OBX in wt mice at 48 hpO (Fig. 12). The expression of Msr1 was negligible in the MSR/ mice, with no significant change in expression as a result of OBX.
In addition, we used real-time RT-PCR to determine the direction and magnitude of changes in expression levels for the interleukin gene Il1rn (Fig. 7D), the chemokine gene Cxcl1, the scavenger receptor gene Scarb1 (Fig. 10C), and the phagosome-associated gene H2-D1 (Fig. 11B) at time points at which the microarray data indicated that there were statistically significant differences in expression in the two strains (Table 1). In all cases, the direction of change in the microarray data and the real-time RT-PCR data was the same. Although fold changes determined by real-time RT-PCR are considerably larger for Il1rn and Cxcl1 compared with those obtained from microarray data, these genes also had relatively large fold changes compared with other genes in the microarray data. Because the expression levels of Scarb1 in the wt and MSR/ mice as determined by real-time RT-PCR were not significantly different although the direction of change was the same, we validated the upregulated expression of Scarb1 mRNA with in situ hybridization.
Molecular validation using in situ hybridization.
The expression of Scarb1 mRNA was evaluated in the OE of 0 hpO sham and 48 hpO MSR/ mice. In sham MSR/ mice (Fig. 13A), there was a faint band of hybridization signals localized over the OE (Fig. 13B). At 48 hpO, the OE was noticeably thinner (Fig. 13C), and hybridization signals increased markedly (Fig. 13D).
There was a very sparse, unlocalized distribution of radioactivity when the sense probe was used in place of the antisense probe in the OE from sham MSR/ mice, and there was no change in the OE from 48 hpO MSR/ mice (not shown). The stronger si