* Department of Biochemistry and Molecular Biology, Medical School Duluth, University of Minnesota, Duluth, Minnesota 55812; Department of Veterinary and Biomedical Science, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota 55108;
Department of Mathematics and Statistics, College of Science and Engineering, University of Minnesota, Duluth, Minnesota 55812; and
Department of Molecular Biosciences, School of Veterinary Medicine, University of California, Davis, California 95616
1 To whom correspondence should be addressed at Department of Biochemistry and Molecular Biology, Medical School Duluth, University of Minnesota, Duluth, MN 55812. E-mail address: jregal{at}d.umn.edu.
Received June 28, 2005; accepted August 30, 2005
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: rodent; eosinophils; allergy; lung; nitric oxide.
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In general, it is assumed that mechanisms and mediators of the effector phase of asthma are common to all allergens, including low-molecular-weight occupational allergens such as TMA. Zimmermann et al. (2003) used two different antigens, OVA and Aspergillus, to elicit the asthmatic phenotype in the mouse and found that expression of 291 genes in the lung changed in common for both allergens. They described these as "signature genes" for allergic airway inflammation and have pursued their role in the effector phase of asthma. However, 236 and 205 genes were unique to sensitization and challenge with Aspergillus and OVA, respectively, suggesting that mechanisms may differ with the allergen. Zimmermann et al. (2003)
attributed these differences in gene expression between the two allergens to differences in exposure regimens (i.e., acute exposure to OVA and chronic exposure to Aspergillus). Clearly some differences in mechanisms of asthma may be related to differences in exposure conditions to the allergen. For example, in mouse models of asthma using OVA as the allergen, differences in the pathophysiology of the response have been attributed in part to variation in the dose and intensity of antigen administration, routes of sensitization and challenge, as well as the inbred strain of animal used in the mouse asthma models (Humbles et al., 2002
; Ma et al., 2002
; Regal, 2004
). Other evidence in guinea pig asthma models suggests that the role of the complement system and arachidonate metabolites may differ depending on whether the inciting allergen is OVA or TMA (Arakawa et al., 1993
; Fraser et al., 1995
; Regal and Fraser, 1996
). Also, in humans, variability of the asthma phenotype suggests that different mechanisms may result in asthma symptoms of different severity and character (Bel, 2004
).
We hypothesized that different allergens would evoke unique patterns of gene expression suggestive of unique effector mechanisms leading to the asthma phenotype. To test this hypothesis, gene expression was examined in the effector phase in lungs from Balb/c mice sensitized and challenged with two different allergens, TMA and OVA, evaluating eosinophilia as the asthma endpoint. To insure that expression pattern differences reflected mechanistic differences due to the allergen and not the exposure regimen (Humbles et al., 2002; Ma et al., 2002
; Regal, 2004
; Zimmermann et al., 2003
), identical acute sensitization and challenge regimens were used. Doses of TMA and OVA were chosen which evoked similar changes in eosinophil infiltration into the mouse lung. OVA-induced changes in gene expression in the present study confirmed increases in many genes reported previously, including arginase (Daheshia et al., 2002
; Scheerens et al., 2002
; Zhu et al., 2004
; Zimmermann et al., 2003
; Zuhdi Alimam et al., 2000
). In addition, TMA sensitization and challenge revealed a unique gene expression pattern compared to OVA. Further investigation focused on differences in arginase and changes in expression of genes involved in arginine metabolism and nitric oxide production. These data suggest unique effector mechanisms involving arginase exist for OVA- and TMA-induced asthma, and optimal therapy may need to target different pathways for asthma induced by different allergens.
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Animals were sensitized with OVA (0.3%) and TMA (3%) as previously described (Regal et al., 2001). Briefly, animals received OVA and TMA suspensions in corn oil intradermally on days 1 and 3, and intratracheally administered OVA or TMA-MSA (30 µg) on day 12. Animals were challenged on day 19 with 400 µg OVA, TMA-MSA, or MSA intratracheally to evoke the asthma phenotype. Seventy-two hours after intratracheal challenge, animals were anesthetized with pentobarbital and bled by cardiac puncture or retroorbital bleed. Lungs were removed for RNA isolation and analysis of eosinophil peroxidase (Regal et al., 2001
) as a measure of lung eosinophilia. Lung lobes for RNA isolation were flash frozen in liquid nitrogen or added to RNAlater (Ambion). Lobes for eosinophil peroxidase (EPO) were processed as previously described (Regal et al., 2001
). In a subset of these animals, a lung lobe was also removed for measurement of myeloperoxidase as an indicator of lung neutrophilia. In a separate group of animals (n = 19), after anesthesia with pentobarbital, lungs were lavaged with 2 x 0.9 ml PBS, and bronchoalveolar lavage (BAL) cells were identified and counted as previously described (Regal et al., 2001
).
Preparation of biotinylated target RNA for Affymetrix GeneChip analysis.
Procedures described in the Affymetrix GeneChip Expression Analysis Manual were used. Total RNA was isolated from lung lobes of six to eight animals in each of the four treatment groups and hybridized at one animal per chip. Lung lobes for RNA were thawed, homogenized with a Kontes Pellet Pestle motor and a 1.5-ml pellet pestle in RLT buffer (RNeasy Qiagen), and then further disrupted by passing through 19-, 23- and finally 25-gauge needles. The Qiagen protocol was followed to prepare the total RNA from the homogenized lung tissue, and spectrophotometric measurements were used to quantify yield. First- and second-strand cDNAs were synthesized from 15 µg of total RNA using a HPLC-purified T7-(dT)24 primer as recommended by Affymetrix and the SuperScript Choice kit from Invitrogen Life Technologies. The double-stranded cDNA was further purified using phase-lock gel phenol/chloroform extraction. Size distribution and yield of the cDNA were determined on a 1% agarose gel prior to the synthesis of biotin-labeled cRNA. The synthesis of biotin-labeled cRNA was as described in the Affymetrix manual, using the Bioarray High Yield RNA Transcript labeling ENZO Kit. Products of the in vitro transcription (from both before and after purification with the Qiagen RNeasy columns) were analyzed by gel electrophoresis. The cRNA was fragmented at 94°C for 35 min in Fragmentation buffer (Affymetrix) and checked by gel electrophoresis prior to delivery to the Biomedical Image Processing Laboratory of the University of Minnesota for hybridization with the murine Affymetrix GeneChip U74Av2.
Quality control, normalization, and statistical analysis of gene expression data.
Microarray Suite 5.0 was used for initial processing of the data. For each chip, we examined the chip image for gross flaws such as bubbles or streaks and for the alignment and intensity of both the corner signals and "landing lights" on the border of the chip. The Chip Report was inspected for each chip and the 5'/3' ratios, percent present calls, background level, and presence of the spike controls were compared to Affymetrix standards. In addition, dChip (Li and Wong, 2001) was used to identify outlier chips. GeneTraffic (Iobion) was used to normalize the data using Robust Multichip Analysis (RMA; Irizarry et al., 2003
). RMA-normalized, log-scale intensities from Gene Traffic were exported to SAS for identification of candidate genes by statistical analysis. Data were log transformed to insure a more normal distribution for statistical tests to be valid. Because of observed unequal variances in the gene expression data, significant treatment effects were tested using Satterthwaite's method for error degrees of freedom with PROC MIXED in SAS. In addition to a treatment effect, the model also included an effect to adjust for the date the chips were processed at the central array facility. The RMA method for normalization is designed to adjust for some differences in intensities, but there still remained a significant date effect after RMA normalization. Hence, we included the date effect in our ANOVA model (Iobion and Irizarry, personal communication).
Confirmation of selected microarray results using quantitative real-time PCR (qRT-PCR).
Total RNA was isolated according to the manufacturer's protocol (RNeasy, Qiagen) from lung tissue processed as described above. Lung samples were chosen from animals used for the gene expression studies, as well as from animals that had lung EPO measurements only. Two µg of total RNA were reverse transcribed using random hexamer primers (pf(N)6, Amersham Pharmacia Biotech Inc.) in a 20-µl reaction volume with 4 units of reverse transcriptase (Omniscript RT Kit, Qiagen), 0.4 units of RNase OUT Ribonuclease Inhibitor (Invitrogen Life Technologies), and a final concentration of 500 µM dNTPs in 1x reverse transcriptase buffer. The reverse transcription reaction was held at 37°C for 1 h. Subsequently, 0.5 µl of the 20-µl reverse transcription reaction for each animal was used in quantitative real-time PCR reactions (qRT-PCR) with each of seven primer sets (Integrated DNA Technologies, Inc., Table 1). A standard 15-µl qRT-PCR reaction mix contained 0.5 µl template (i.e., reverse transcription reaction product) and 7.5 µl SyberGreen Master Mix (Applied Biosystems). The qRT-PCR reactions were amplified in a 96-well plate format on a 7900HT Sequence Detection System (Applied Biosystems). A separate plate was used for each primer and contained experimental samples in duplicate for each animal. In addition to a no-template control, a control of reverse transcription reaction from standard lung RNA (Ambion) with G6pt1 (glucose-6-phosphate transporter protein 1) primers and a third control of reverse transcription reaction from standard lung RNA with test primers were included on each plate. The PCR profile was as follows: 50°C for 2 min, 10 min at 95°C, and then 40 cycles of 15 s at 95°C and 1 minute extensions at 60°C.
|
Arginase assay.
Arginase activity was measured by the method of Özer (1985). Briefly, lung samples were homogenized in 0.05 M Tris and 0.5% Triton-X 100 and then centrifuged 15 min at >900 x g (3000 rpm) in a refrigerated centrifuge (Jouan) to pellet debris. The supernatants were removed, and the MnCl2 concentration was adjusted to 1 mM. The samples were then heated to 55°C for 10 min to activate arginase. Assay conditions were as reported by Özer (1985)
with the following changes: 200 µl reactions per well were run in microtiter plates, maintaining reagent concentrations as reported in Özer (1985)
. NADH was substituted for NADPH, and 0.89 mM ADP was present in each reaction well. Serial dilutions of each sample were assayed with and without L-arginine substrate to control for background NADH consumption in the tissue extracts. A standard curve was determined from serial dilutions of bovine liver arginase (Calzyme).
Statistical analysis of cellular infiltration data, qRT-PCR and arginase activity.
All data were log transformed to insure a more normal distribution of the data and the valid application of statistical tests. Figures show the geometric mean +1 standard error, with significant comparisons indicated by an asterisk. Statistical significance was defined as p < 0.05. Throughout this study, four groups of animals were considered: (1) OVA sensitized, MSA challenged (OVAc), (2) OVA sensitized, OVA challenged (OVA), (3) TMA sensitized, MSA challenged (TMAc), (4) TMA sensitized, TMA-MSA challenged (TMA). Because of unequal variances, significances of treatment effects were tested using Satterthwaite's method for error degrees of freedom with PROC MIXED in SAS (SAS Institute Inc., Cary, NC). The following comparisons were made: OVAc versus OVA; TMAc versus TMA; difference in OVAc versus OVA to the difference in TMAc versus TMA; OVAc versus TMAc; and OVA versus TMA.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
A separate group of animals was anesthetized and lungs lavaged to assess cellular infiltration into the airspace with this sensitization and challenge regimen. Total white blood cells in the BAL were significantly increased over control in the TMA-sensitized and -challenged animals, but not in the OVA-sensitized and -challenged animals (data not shown). The total number of white blood cells in TMAc compared to OVAc did not differ. The change in white blood cells seen with TMA was primarily due to increases in eosinophils and macrophages. As seen in Figure 2, eosinophils were significantly increased in the bronchoalveolar lavage of the TMA group, but not in the OVA group, compared to their respective controls (lower brackets). The percent of eosinophils in the lavage increased from 10.5 to 23.5% for OVAc to OVA, and from 1.5% to 28.3% for TMAc to TMA. The eosinophil numbers in the OVAc group were significantly higher than the TMAc group, indicating a difference in the number of eosinophils between animals sensitized with OVA versus TMA at the time of challenge. Eosinophils in naïve animals were similar to the TMAc treatment group (data not shown), indicating that the elevated eosinophils in OVAc as compared to TMAc may have been due to residual effects of the intratracheal sensitization dose for OVAc. The number of macrophages in the BAL was also significantly increased in the TMA group compared to TMAc, but no significant increase was detected comparing OVA to OVAc. This represents a decrease in the percent of macrophages from 75% to 65% for OVAc to OVA and from 86% to 63% for TMAc to TMA. No treatment effect was seen for neutrophils (data not shown).
|
|
|
A subset of the uniquely up-or down-regulated genes is categorized in Table 3. The OVA animals display an up-regulation of numerous chemokines and receptors that show no change or decreased expression in TMA-treated animals. Ccl11 (eotaxin 1) is a mediator of eosinophil infiltration in OVA-induced asthma (Ganzalo et al., 1996). Eotaxin 2, a second eosinophil chemoattractant, is not represented on the U74Av2 chip. Two other chemoattractants up-regulated in both OVA- and TMA-induced asthma include the lectin YM1 (Chi3l3), which is identical to the lymphocyte-derived eosinophil chemotactic factor (ECF-L; Owhashi et al., 2000
), and AMCase, which has 67% protein identity to mouse YM1/ECF-L/Chi3l3 (Owhashi et al., 2000
). It should also be noted that many chemokines and receptors were coordinately up-regulated for OVA and TMA including Ccl8, Ccr5, Ccl9, Ccr1, and Ccl6 (Table 2 and Supplementary Table 4).
|
The greatest difference in expression of cell adhesion molecules was seen for the gene intelectin (three-fold difference between OVA and TMA). This gene is expressed in intestinal Paneth cells and is up-regulated in response to Trichinella infection (Pemberton et al., 2004). Its up-regulation in OVA-induced asthma has been previously described (Zimmermann et al., 2003
), but its role in asthma is unknown.
The complement component C5 has been described as a susceptibility locus for the development of airway hyperresponsiveness in OVA models of asthma (Hawlisch et al., 2004). In addition, C3a, a product of complement system activation, has been implicated as a critical mediator of airway hyperresponsiveness in OVA-induced asthma (Hawlisch et al., 2004
). Thus, we looked for differential expression of components of the complement system. As seen in Table 3, the complement genes Properdin and Daf2 were up-regulated in OVA-induced asthma but showed no change in TMA-treated animals. However, other complement genes C1qa, C1qb, C1qc, and C3 were coordinately up-regulated for both OVA- and TMA-induced asthma (Table 2 and Supplementary Table 4). Thus, expression of complement system genes does not fit a consistent pattern.
A very limited number of genes (3 + 1) demonstrated oppositional regulation in OVA- versus TMA-induced asthma (Fig. 3). Three genes (Sprr2a, pIgR and Fv-4) were up-regulated in OVA and down-regulated in TMA, whereas cytosolic acyl-CoA thioesterase 1 (Cte1) was significantly up-regulated in TMA and down-regulated in OVA. Sprr2a is a member of a family of small proline-rich proteins which form the precursor to the cornified envelope (CE) in stratified squamous epithelial cells and may be important in airway remodeling (Zimmermann et al., 2005). pIgR is the polymeric immunoglobulin receptor whose function includes the transport of IgA and IgM class polymeric immunoglobulins across the glandular and mucosal epithelia (reviewed in Mostov, 1994
). Sprr2a and pIgR are both STAT6-regulated proteins (Schjerven et al., 2000
; Zimmermann et al., 2005
). The differential regulation of Sprr2a and pIgR with OVA and TMA suggests coregulators interacting with STAT6 vary with allergen.
Confirmation of Microarray Analysis by qRT-PCR
To confirm the microarray results, a subset of genes was analyzed by qRT-PCR using the primers shown in Table 1. Lungs from six to eight animals were used for qRT-PCR determinations and included at least three animals in each treatment group that were also analyzed by microarrays. Genes with a variety of levels of expression were chosen, including genes involved in arginine metabolism. The gene G6pt1 (glucose-6-phosphate transporter protein 1) was selected as the housekeeping control by selecting the probe set from the microarray data with the smallest treatment effect (data not shown). qRT-PCR confirmed that G6pt1 expression does not significantly change with treatment (data not shown). Data from all primer sets were normalized to G6pt1 expression. The correlation (r) between microarray and qRT-PCR data for animals used in both assays ranged from 0.68 to 0.97 and was statistically significant (p < 0.05) for Gob-5, Ccl8, Chi3l3, AMCase, Arg1, Gatm, and Ccl11. In the case of Vcam1, Arg2, and Pglyrp1, the correlation was not significant. This lack of significance may be explained by the low expression levels observed for these genes in the microarray analysis, as Etienne et al. (2004) noted a reduced correlation between microarray and RT-PCR data with reduced expression levels, increasing number of absent calls, and increased distance between probes.
The relative RNA quantity for each animal determined by qRT-PCR was plotted by treatment for the selected genes (Fig. 4). In Panel A, statistical analysis indicated a significant and similar increase in Chi3l3, Ccl8, and Gob-5 for both allergens, consistent with the microarray data. In Panel B, statistical analysis of the qRT-PCR data indicated a significant increase in AMCase, Arg1, and Ccl11 for each allergen, with the increase for OVA/OVAc being significantly greater than that seen for TMA/TMAc. Gatm expression levels, shown in Panel C, significantly increased with OVA/OVAc, but not with TMA/TMAc. For each of the genes in Panels A, B, and C, the experimental/control ratio from microarray analysis showed the same relative changes as seen with qRT-PCR. In general, our data are consistent with Yuen et al. (2002) demonstrating qRT-PCR estimates of mRNA levels are consistently larger than those measured by microarray. No significant differences were detected by qRT-PCR analysis with Vcam1, Arg2, and Pglyrp1 (data not shown).
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The lung contains numerous cell types: endothelial, epithelial, dendritic, white blood cells, nerve cells, etc. Thus, gene expression changes noted in this study could be the result of changes in any one cell type, or changes in the percent of a single cell type, or a combination thereof. All of the cell types have been implicated as being important in the asthmatic response. Studying each individual cell type in vitro yields useful information, as does studying isolated airways and their reactivity. However, the interaction of cell types and mediators ultimately leads to the response in the whole lung with eosinophilia and airway obstruction. For example, gene expression changes in a purified macrophage population may be quite different from gene expression changes in a macrophage residing in the lung due to interactions with adjacent cell populations within the lung tissue. Thus, it is important to consider the whole lung and the complex interactions occurring in vivo, rather than just a single cell type.
In compiling a list of candidate genes from the array data, we utilized both a test of statistical significance and an expression ratio cutoff. Fold-change cutoffs used alone tend to result in a larger fraction of false positives due to outlier responses. After using a q-value cutoff of 0.1, we further narrowed our candidate gene list by declaring that an experimental-to-control ratio of 1.2 or greater (or 1/1.2 and lower) was necessary for us to consider the gene as potentially biologically significant. This further reduced our candidate gene list to 249 genes. A two-fold change was initially an arbitrary change reported as significant in many early microarray studies. However, this standard has changed considerably with advances in (a) the normalization of data (e.g., RMA) and (b) improved statistical analysis of array data. Our choice of 1.2 as an expression ratio cutoff with biological significance was based on a review of the literature and the fact that we were using a whole tissue where important changes in gene expression in a single cell type may be diluted. Because we were using a mixed tissue, we wanted to insure that any gene expression changes that might have been diluted due to the variety of cell types present were also still considered. Therefore, we used an experimental/control ratio of 1.2 that was both statistically significant and broad. Studies of Jin et al. (2001) and Pierce et al. (2003)
demonstrated that gene expression changes as low as 1.2-fold were biologically significant in a Drosophilia model and a colorectal cancer cell line, respectively. The actual biological significance of any change in gene expression can only be determined by critically testing focused hypotheses using more classical methods, including inhibitors and measurements of protein.
Our study demonstrated that many genes are up-regulated to a greater extent with TMA than OVA (e.g., Gob-5), just as many genes are up-regulated more with OVA than TMA (e.g., arginase). This suggests that differences in gene expression noted for OVA versus TMA are not likely due to differences in the dose of allergen or the extent of the allergic reaction. Additional unpublished studies using a different OVA sensitization and challenge regimen resulted in maximum lung EPO values comparable to the TMA treatment group in this study. In these unpublished studies with OVA as the allergen and a lower maximum lung EPO response, arginase up-regulation was also evident (experimental/control ratio of 10.3), suggesting that the extent of the eosinophilia does not determine the up-regulation of arginase. One would expect genes important in cellular infiltration to correlate in some manner with that biological endpoint. Clearly the differentially expressed genes are not all critical mediators of the eosinophil infiltration component of the asthmatic response, and the extent of up- or down-regulation for any given gene does not necessarily correlate with the magnitude of the allergic endpoint being measured. Our long-term goal is to link allergen-specific gene expression patterns with allergen-specific effector pathways or profiles of asthma symptoms. The differences in gene expression could reflect either differences in effector pathways leading to the same profile of asthma symptoms or differences in effector pathways leading to a different profile of asthma symptoms, i.e., asthma heterogeneity.
A single time point (72 h after challenge) was chosen to compare gene expression changes in the effector phase of asthma induced by the two different antigens, OVA and TMA. This 72 h time point was chosen because, in preliminary studies examining lung eosinophilia at 48 and 72 h after a single 400 µg TMA-MSA or OVA challenge, the change in lung EPO was as great or greater at 72 h than 48 h for both allergens. Similar time courses in cell infiltration are expected with these identical sensitization and challenge regimens, and the differences in gene expression are not attributable to different routes of administration of the antigens, the presence of different adjuvants, or different mouse strains. The phenotypes are similar in that the change in eosinophilia in the lung tissue is no different for the two allergens (Fig. 1). Unlike some more vigorous sensitization and challenge regimens, our method of sensitization and challenge with OVA and TMA did not result in any increase in neutrophils in the lung or BAL. However, we saw a significant increase in BAL macrophages and eosinophils with TMA sensitization and challenge, but not OVA sensitization and challenge (Fig. 2). We have not investigated whether gene expression profiles in TMA-induced asthma suggest greater macrophage chemotaxis and/or extravasation or instead result from greater expression of local myeloid progenitor cells. In a study by Tomkinson et al. (2001) using a single intranasal OVA challenge in sensitized animals, tissue eosinophilia peaked 48 h after OVA challenge, and BAL eosinophils peaked 4896 h after challenge and were sustained for at least 8 days. Thus, there is clear precedent for differences in tissue eosinophilia and BAL eosinophilia at any particular time point. In addition, studies of Pope et al. (2005)
have suggested that eotaxin 1 and 2 are mediators important in determining eosinophil accumulation in the airspace versus the tissue. Thus, the difference in BAL eosinophilia noted with the two allergens may reflect different effector pathways contributing to the heterogeneity of the allergic response with the two allergens. Clearly, further studies are indicated to rigorously examine the time course of changes in gene expression with the time course of each asthma symptom, e.g., airway hyperresponsiveness, mucus production, BAL and tissue eosinophilia, etc., for each allergen.
Of the 249 genes differentially regulated for OVA and TMA, we chose to focus on genes that were uniquely up- or down- regulated in arginine metabolism. Previous studies of others had implicated arginase as an important step in asthma pathogenesis (Zimmermann et al., 2003). A simplified diagram of our hypothetical working model of arginase action in the cell and potential fates of L-arginine in OVA- versus TMA-induced asthma is shown in Figure 6. Our gene expression results predict that OVA-induced asthma would increase levels of ornithine and perhaps the downstream products of spermine and polyamines due to increased expression of arginase 1, arginase 2, and Gatm and decreased production of NO. L-ornithine is a precursor element for production of collagen and polyamines, both of which are important for airway remodeling. This is consistent with the up-regulation of other remodeling genes noted in OVA-induced asthma but not TMA-induced asthma (Table 3). The decreased NO production predicted in the OVA model could be explained by a reduction in L-arginine availability due to increased activities and/or expression of arginase 1, arginase 2, and Gatm. Lower concentrations of L-arginine would ultimately favor the competitive inhibition of NOS by endogenously formed asymmetric dimethylarginine (ADMA) and N-monomethylarginine (NMMA). In contrast, the animals with TMA-induced asthma present a gene expression pattern favoring increased production of NO through the stable expression of arginase 1 and arginase 2 and the up-regulation of Ddah2, which will decrease the levels of competitive inhibitors of NOS (ADMA and NMMA). These data are consistent with results in a guinea pig model of TMA-induced asthma demonstrating an up-regulation of Ca2+-independent iNOS activity in the lung tissue and increased nitrite levels in BAL (Yan et al., 1995
). Whether the hypothetical NO production results in bronchodilation in the airway or increased reactive nitrogen species with protein nitration, inflammation, and cytotoxicity is unknown. The literature on NO in the lung is extensive and conflicting, and the end result of NO production depends on the concentration of NOS substrates and the site of NO generation (Ricciardolo et al., 2004
). Certainly, continued studies beyond gene expression are required to adequately test the hypotheses presented in this working model.
|
Remodeling and secretion may comprise a larger component of the response in OVA-induced asthma than in TMA-induced asthma, as suggested by the difference in gene expression profiles for Sprr2a, Sprr2b, RELM ß, scinderin, and Gob-4 between OVA and TMA and as noted in the working model in Figure 6. Sprr2a and RELM ß are produced in goblet cells (He et al., 2003; Komiya et al., 1999
) and contain STAT binding sites, consistent with many other genes involved in airway inflammation. Several other genes expressed in goblet cells (Gob-5, Muc5AC, and Tff2) have been previously described in OVA-induced asthma (Nakanishi et al., 2001
; Nikolaidis et al., 2003
; Zuhdi Alimam et al., 2000
) and were up-regulated in response to both OVA and TMA in our study.
The largest difference in expression detected between OVA and TMA was noted with Sprr2a, a member of a family of small proline-rich proteins which form the precursor to the cornified envelope (CE) in stratified squamous epithelial cells. This protein has been reported as STAT6-dependent for both OVA- and Aspergillus-induced asthma models (Zimmermann et al., 2004) and may be important in airway remodeling (Hohl et al., 1995
; Zimmermann et al., 2005
). However, in our TMA-induced model of asthma, we do not observe an increase in Sprr2a expression. This difference indicates that the regulation or timing of regulation for Sprr2a in TMA-treated animals differs from OVA-treated animals. If the regulation of Sprr2a is STAT6-dependent, then there may be coregulators of Sprr2a that differ between these asthma models. Further analysis is necessary to determine what if any coregulators exist for Sprr2a.
A group of 291 genes, proposed as "asthma signature genes" by Zimmermann et al. (2003), increased in expression to OVA- and Aspergillus-induced airway inflammation and included the genes Sprr2a, RELM ß, arginase 1, Scinderin, cathepsin B, and Ccl11. In the present model of OVA- and TMA-induced asthma, this particular subset of "asthma signature genes" were either increased in OVA- but not TMA-induced asthma (Sprr2a, RELM beta, Scinderin, cathepsin B) or elevated in both but to significantly different levels (Arg1, Ccl11). The disparities in gene expression between our study and that of Zimmermann et al. (2003)
may be due in part to different treatment regimens, samples sizes, or statistical analysis. Zimmermann et al. (2003)
used an acute and chronic treatment regimen for OVA and Aspergillus, respectively, whereas we used the same acute treatment regimen for both OVA and TMA. In addition, the sample sizes and statistical analysis differed in the two studies, with Zimmermann et al. (2003)
using sample sizes of two to four for treatment groups and selection criteria of p < 0.05 and fold change >2. In contrast, our study used an n of 6 or more for each treatment group as recommended by Cui and Churchill (2003)
and statistical analysis included computation of q-values to guard against false positives.
A cross-comparison of the results of our study and Zimmermann et al. (2003) results in a set of 62 genes that were up-regulated in response to all 3 allergens; OVA, TMA, and Aspergillus (Supplementary Table 2). Five of these genes, Ccl8, Chi3l3, Gob-5, Mmp12, and Tff2, have a 2-fold or greater increase in expression with all three allergens (this study and Zimmermann et al., 2003
). Together these data suggest that the effector mechanism of TMA-induced asthma may converge with the OVA and Aspergillus effector pathways at a point downstream of the asthma signature genes described by Zimmermann et al. (2003)
. Only continued studies with additional allergens will determine whether up-regulation of this set of 62 genes is indicative of asthma and suggestive of a common final effector pathway for all allergens. Nonetheless, TMA-induced allergic airway inflammation does not completely fit the model of asthma signature genes presented by Zimmermann et al. (2003)
.
Asthma is a complex lung disorder involving multiple cell types and signaling from within and beyond the lung. Although the phenotype of different allergens may appear similar, the mechanisms by which the asthmatic phenotype is presented appear to be quite different. This may have significant impact on potential therapeutic interventions in asthma. We have used microarray analysis and subsequent confirmation by qRT-PCR and an enzymatic assay to show that two occupational allergens, OVA and TMA, present different gene expression profiles, although their respective biological phenotypes are similar. Specifically, L-arginine utilization may be quite different in the two models due to contrasting levels of arginase activity. Such differences may result in alternate signaling pathways and cellular involvement producing the same constellation of asthma symptoms.
![]() |
SUPPLEMENTARY DATA |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
ACKNOWLEDGMENTS |
---|
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Bel, E. H. (2004). Clinical phenotypes of asthma. Curr. Opin. Pulm. Med. 10, 4450.[CrossRef][ISI][Medline]
Boeniger, M. F., Lummus, Z. L., Biagini, R. E., Bernstein, D. I., Swanson, M. C., Reed, C., and Massoudi, M. (2001). Exposure to protein aeroallergens in egg processing facilities. Appl. Occup. Environ. Health 16, 660670.
Cui, X., and Churchill, G. A. (2003). Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 4, 210219.[CrossRef][Medline]
Daheshia, M., Tian, N., Connolly, T., Drawid, A., Wu, Q., Bienvenu, J. G., Cavallo, J., Jupp, R., De Sanctis, G. T., and Minnich, A. (2002). Molecular characterization of antigen-induced lung inflammation in a murine model of asthma. Ann. N.Y. Acad. Sci. 975, 148159.
Donnelly, L. E., and Barnes, P. J. (2004). Acidic mammalian chitinasea potential target for asthma therapy. Trends Pharmacol. Sci. 25, 509511.[CrossRef][ISI][Medline]
Etienne, W., Meyer, M. H., Peppers, J., and Meyer, R. A., Jr. (2004). Comparison of mRNA gene expression by RT-PCR and DNA microarray. Biotechniques 36, 618626.[ISI][Medline]
Fraser, D. G., Regal, J. F., and Arndt, M. L. (1995). Trimellitic anhydride-induced allergic response in the lung: Role of the complement system in cellular changes. J. Pharmacol. Exp. Ther. 273, 793801.[Abstract]
Fulkerson, P. C., Zimmermann, N., Brandt, E. B., Muntel, E. E., Doepker, M. P., Kavanaugh, J. L., Mishra, A., Witte, D. P., Zhang, H., Farber, J. M., et al. (2004). Negative regulation of eosinophil recruitment to the lung by the chemokine monokine induced by IFN- (Mig, CXCL9). Proc. Natl. Acad. Sci. U.S.A. 101, 19871992.
Ganzalo, J. A., Jia, G. Q., Aguirre, V., Friend, D., Coyle, A. J., Jenkins, N. A., Lin, G. S., Katz, H., Lichtman, A., Copeland, N., et al. (1996). Mouse eotaxin expression parallels eosinophil accumulation during lung allergic inflammation but it is not restricted to a Th2-type response. Immunity 4, 114.[CrossRef][ISI][Medline]
Hawlisch, H., Wills-Karp, M., Karp, C. L., and Köhl, J. (2004). The anaphylatoxins bridge innate and adaptive immune responses in allergic asthma. Mol. Immunol. 41, 123131.[CrossRef][ISI][Medline]
He, W., Wang, M. L., Jiang, H. Q., Steppan, C. M., Shin, M. E., Thurnheer, M. C., Cebra, J. J., Lazar, M. A., and Wu, G. D. (2003). Bacterial colonization leads to the colonic secretion of RELMbeta/FIZZ2, a novel goblet cell-specific protein. Gastroenterology 125, 13881397.[CrossRef][ISI][Medline]
Hohl, D., de Viragh, P. A., Amiguet-Barras, F., Gibbs, S., Backendorf, C., and Huber, M. (1995). The small proline-rich proteins constitute a multigene family of differentially regulated cornified cell envelope precursor proteins. Dermatology 104, 902909.
Humbles, A. A., Lu, B., Friend, D. S., Okinaga, S., Lora, J., Al-Garawi, A., Martin, T. R., Gerard, N. P., and Gerard, C. (2002). The murine CCR3 receptor regulates both the role of eosinophils and mast cells in allergen-induced airway inflammation and hyperresponsiveness. Proc. Natl. Acad. Sci. U.S.A. 99, 14791484.
Irizarry, R. A., Bolstad, B. M., Collin, F., Cope, L. M., Hobbs, B., and Speed, T. P. (2003). Summaries of affymetrix genechip probe level data. Nucleic Acids Res. 31, e15.
Jin, W., Riley, R. M., Wolfinger, R. D., White, K. P., Passador-Gurgel, G., and Gibson, G. (2001). The contributions of sex, genotype and age to transcriptional variance in Drosophila melanogaster. Nat. Genet. 29, 389395.[CrossRef][ISI][Medline]
Komiya, T., Tanigawa, Y., and Hirohashi, S. (1999). Cloning of the gene gob-4, which is expressed in intestinal goblet cells in mice. Biochim. Biophys. Acta 1444, 434438.[ISI][Medline]
Li, C., and Wong, W. H. (2001). Model-based analysis of oligonucleotide arrays: Expression index computation and outlier detection. Proc. Natl. Acad. Sci. U.S.A. 98, 3136.
Ma, W., Bryce, P. J., Humbles, A. A., Laouini, D., Yalcindag, A., Alenius, H., Friend, S. S., Oettgen, H. C., Gerard, C., and Geha, R. S. (2002). CCR3 is essential for skin eosinophilia and airway hyperresponsiveness in a murine model of allergic skin inflammation. J. Clin. Invest. 109, 621628.
Martin-Verdeaux, S., Pombo, I., Iannascoli, B., Roa, M., Varin-Blank, N., Rivera, J., and Blank, U. (2002). Evidence of a role for Munc182 and microtubules in mast cell granule exocytosis. J. Cell Sci. 116, 325334.[CrossRef][ISI]
Mostov, K. E. (1994). Transepithelial transport of immunoglobulins. Ann. Rev. Immunol. 12, 6384.[CrossRef][ISI][Medline]
Nakanishi, A, Morita, S., Iwashita, H., Sagiya, Y., Ashida, Y., Shirafuji, H., Fujisawa, Y., Nishimura, O., and Fujino, M. (2001). Role of gob-5 in mucus overproduction and airway hyperresponsiveness in asthma. Proc. Natl. Acad. Sci. U.S.A. 98, 51755180.
Nikolaidis, N. M., Zimmermann, N., King, N. E., Mishra, A., Pope, S. M., Finkelman, F. D., and Rothenberg, M. E. (2003). Trefoil factor-2 is an allergen-induced gene regulated by Th2 cytokines and STAT6 in the lung. Am. J. Respir. Cell Mol. Biol. 29, 458464.
Owhashi, M., Arita, H., and Hayai, N. (2000). Identification of a novel eosinophil chemotactic cytokine (ECF-L) as a chitinase family protein. J. Biol. Chem. 275, 12791286.
Özer, N. (1985). A new enzyme-coupled spectrophotometric method for the determination of arginase activity. Biochem. Med. 33, 367371.[CrossRef][ISI][Medline]
Pemberton, A. D., Knight, P. A., Gamble, J., Colledge, W. H., Lee, J. K., Pierce, M., and Miller, H. R. (2004). Innate BALB/c enteric epithelial responses to Trichinella spiralis: Inducible expression of a novel goblet cell lectin, intelectin-2, and its natural deletion in C57BL/10 mice. J. Immunol. 173, 18941901.
Pierce, M., Wang, C., Stump, M., and Kamb, A. (2003). Overexpression of the ß-catenin binding domain of cadherin selectively kills colorectal cancer cells. Int. J. Cancer 107, 229237.[CrossRef][ISI][Medline]
Pope, S. M., Fulkerson, P. C., Blanchard, C., Akei, H. S., Nikolaidis, N. M., Zimmermann, N., Molkentin, J. D., and Rothenberg, M. E. (2005). Identification of a cooperative mechanism involving interleukin-13 and eotaxin-2 in experimental allergic lung inflammation. J. Biol. Chem. 280, 1395213961.
Regal, J. F. (2004). Immunologic effector mechanisms in animal models of occupational asthma. J. Immunotoxicol. 1, 2538.[CrossRef]
Regal, J. F., and Fraser, D. G. (1996). Systemic complement system depletion does not inhibit cellular accumulation in antihistamine pretreated allergic guinea pig lung. Int. Arch. Allergy Appl. Immunol. 109, 150160.
Regal, J. F., Mohrman, M. E., and Sailstad, D. M. (2001). Trimellitic anhydride-induced eosinophilia in a mouse model of occupational asthma. Toxicol. Appl. Pharmacol. 175, 234242.[CrossRef][ISI][Medline]
Ricciardolo, F. L. M., Sterk, P. J., Gaston, B., and Folkerts, G. (2004). Nitric oxide in health and disease of the respiratory system. Physiol. Rev. 84, 731765.
Scheerens, J., van Gessel, S. B., Nijkamp, F. P., and Folkerts, G. (2002). Eotaxin protein levels and airway pathology in a mouse model for allergic asthma. Eur. J. Pharmacol. 453, 111117.[CrossRef][ISI][Medline]
Schjerven, H., Brandtzaeg, P., and Johansen, F. E. (2000). Mechanism of IL-4-mediated up-regulation of the polymeric Ig receptor: Role of STAT6 in cell type-specific delayed transcriptional response. J. Immunol. 165, 38983906.
Storey, J. D., and Tibshirani, R. (2003). Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. U.S.A. 100, 94409445.
Stütz, A. M., Pickart, L. A., Trifilieff, A., Baumruker, T., Prieschl-Strassmayr, E., and Woisetschläger, M. (2003). The Th2 cell cytokines IL-4 and IL-13 regulate found in inflammatory zone 1/resistin-like molecule gene expression by a STAT6 and CCAAT/enhancer-binding protein-dependent mechanism. J. Immunol. 170, 17891796.
Tomkinson, A., Ciellewicz, G., Duez, C., Larson, K. A., Lee, J. J., and Gelfand, E. W. (2001). Temporal association between airway hyperresponsiveness and airway eosinophilia in ovalbumin-sensitized mice. Am. J. Respir. Crit. Care Med. 163, 721730.
Van Kampen, V., Merget, R., and Baur, X. (2000). Occupational airway sensitizers: An overview on the respective literature. Am. J. Ind. Med. 38, 164218.[CrossRef][ISI][Medline]
Yan, Z.-Q., Hannsson, G. K., Skoogh, B.-E., and Lötvall, J. O. (1995). Induction of nitric oxide synthase in a model of allergic occupational asthma. Allergy 50, 760764.[ISI][Medline]
Yuen, T., Wurmbach, E., Pfeffer, R. L., Ebersole, B. J., and Sealfon, S. C. (2002). Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res. 30, e48.
Zeiss, C. R., Mitchell, J. H., Peenen, V., Harris, J., and Levitz, D. (1990). A twelve year clinical and immunologic evaluation of workers involved in the manufacture of trimellitic anhydride (TMA). Allergy Proc. 11, 7177.[ISI][Medline]
Zeiss, C. R., Patterson, R., Pruzansky, J. J., Miller, M. M., Rosenberg, M., and Levitz, D. (1977). Trimellitic anhydride-induced airway syndromes: Clinical and immunologic studies. J. Allergy Clin. Immunol. 60, 96103.[CrossRef][ISI][Medline]
Zhu, Z., Zheng, T., Homer, R. J., Kim, Y. K., Chen, N. Y., Cohn, L., Hamid, Q., and Elias, J. A. (2004). Acidic mammalian chitinase in asthmatic Th2 inflammation and IL-13 pathway activation. Science 304, 16781682.
Zimmermann, N., Doepker, M. P., Witte, D. P., Stringer, K. F., Fulkerson, P. C., Pope, S. M., Brandt, E. B., Mishra, A., King, N. E., Nikolaidis, N. M., et al. (2005). Expression and regulation of small proline-rich protein 2 in allergic inflammation. Am. J. Respir. Cell Mol. Biol. 32, 428435.
Zimmermann, N., King, N. E., Laporte, J., Yang, M., Mishra, A., Pope, S. M., Muntel, E. E., Witte, D. P., Pegg, A. A., Foster, P. S., et al. (2003). Dissection of experimental asthma with DNA microarray analysis identifies arginase in asthma pathogenesis. J. Clin. Invest. 111, 18631874.
Zimmermann, N., Mishara, A., King, N. E., Fulkerson, P. C., Doepker, M. P., Nikolaidis, N. M., Kindinger, L. E., Moulton, E. A., Aronow, B. J., and Rothenberg, M. E. (2004). Transcript signatures in experimental asthma: Identification of STAT6-dependent and -independent pathways. J. Immunol. 172, 18151824.
Zuhdi Alimam, M., Piazza, F. M., Selby, D. M., Letwin, N., Huang, L., and Rose, M. C. (2000). Muc-5/5AC mucin messenger RNA and protein expression is a marker of goblet cell metaplasia in murine airways. Am. J. Respir. Cell Mol. Biol. 22, 253260.
|