The Journal of Allergy and Clinical Immunology
Volume 116, Issue 2 , Pages 305-311, August 2005

Dissecting asthma using focused transgenic modeling and functional genomics

  • Douglas A. Kuperman, PhD

      Affiliations

    • From the Department of Medicine, Allergy-Immunology Division, Northwestern University Feinberg School of Medicine, Chicago
    • Department of Medicine
    • Lung Biology Center
  • ,
  • Christina C. Lewis, PhD

      Affiliations

    • Department of Medicine
    • Lung Biology Center
  • ,
  • Prescott G. Woodruff, MD

      Affiliations

    • Department of Medicine
    • Division of Pulmonary and Critical Care Medicine
    • Cardiovascular Research Institute
  • ,
  • Madeleine W. Rodriguez, BS

      Affiliations

    • Department of Medicine
    • Lung Biology Center
  • ,
  • Yee Hwa Yang, PhD

      Affiliations

    • Department of Medicine
    • Lung Biology Center
  • ,
  • Gregory M. Dolganov, PhD

      Affiliations

    • Department of Medicine
    • Division of Pulmonary and Critical Care Medicine
    • Cardiovascular Research Institute
  • ,
  • John V. Fahy, MD

      Affiliations

    • Department of Medicine
    • Division of Pulmonary and Critical Care Medicine
    • Cardiovascular Research Institute
  • ,
  • David J. Erle, MD

      Affiliations

    • Department of Medicine
    • Lung Biology Center
    • Division of Pulmonary and Critical Care Medicine
    • Cardiovascular Research Institute
    • Program in Immunology, University of California San Francisco School of Medicine
    • Corresponding Author InformationReprint requests: David J. Erle, MD, UCSF Box 2922, San Francisco, CA 94143-2922.

Received 6 January 2005; received in revised form 28 February 2005; accepted 9 March 2005. published online 03 May 2005.

Chicago, Ill, and San Francisco, Calif

Article Outline

Background

Asthma functional genomics studies are challenging because it is difficult to relate gene expression changes to specific disease mechanisms or pathophysiologic features. Use of simplified model systems might help to address this problem. One such model is the IL-13/Epi (IL-13–overexpressing transgenic mice with STAT6 expression limited to epithelial cells) focused transgenic mouse, which isolates the effects of a single mediator, IL-13, on a single cell type, the airway epithelial cell. These mice develop airway hyperreactivity and mucus overproduction but not airway inflammation.

Objective

To identify how effects of IL-13 on airway epithelial cells contribute to gene expression changes in murine asthma models and determine whether similar changes are seen in people with asthma.

Methods

We analyzed gene expression in ovalbumin allergic mice, IL-13–overexpressing mice, and IL-13/Epi mice with microarrays. We analyzed the expression of human orthologues of genes identified in the mouse studies in airway epithelial cells from subjects with asthma and control subjects.

Results

In comparison with the other 2 models, IL-13/Epi mice had a remarkably small subset of gene expression changes. Human orthologues of some genes identified as increased in the mouse models were more highly expressed in airway epithelial cells from subjects with asthma than in controls. These included calcium-activated chloride channel 1, 15-lipoxygenase, trefoil factor 2, and intelectin.

Conclusion

The combination of focused transgenic models, DNA microarray analyses, and translational studies provides a powerful approach for analyzing the contributions of specific mediators and cell types and for focusing attention on a limited number of genes associated with specific pathophysiologic aspects of asthma.

Key words: Asthma, IL-13, calcium-activated chloride channel, trefoil factor, 15-lipoxygenase, intelectin

Abbreviations used: GAPD, Glyceraldehyde-3-phosphate dehydrogenase, IL-13/Epi, IL-13–overexpressing transgenic mice with STAT6 expression limited to epithelial cells, tg–IL-13, IL-13–overexpressing transgenic mice, STAT6, Signal transducer and activator of transcription factor 6, UCSF, University of California San Francisco

 

Asthma results from a complex interplay between genetic and environmental factors.1 Microarray studies of lung tissue from subjects with asthma and from animals with experimental allergic asthma typically reveal hundreds of genes that are differentially expressed in comparison with normal lung.2, 3, 4 However, it has been difficult to relate gene expression changes to specific mediators, cell types, or pathophysiologic features.

Models that focus on one particular mediator or cell type or on specific pathophysiologic features might help overcome this obstacle. We developed a transgenic mouse model to analyze effects of a single cytokine, IL-13, on a single cell type, the nonciliated airway epithelial cell.5 IL-13 expression is increased in the airways of subjects with asthma,6, 7 and IL-13 is necessary and sufficient for the development of experimental asthma.8, 9, 10 IL-13 activates the signaling molecule signal transducer and activator of transcription factor 6 (STAT6) in many cell types.11 Overexpression of IL-13 in the airways of mice with normal STAT6 expression (tg–IL-13 mice) causes mucus overproduction, airway hyperreactivity, inflammation, fibrosis, and emphysema.10 To isolate the effects of IL-13 on airway epithelial cells, we produced mice that overexpress IL-13 in the airway and express STAT6 only in nonciliated airway epithelial cells. These IL-13/Epi mice developed mucus overproduction and airway hyperreactivity but not inflammation, fibrosis, or emphysema.5

Here we apply a genomics-based approach to identify gene expression changes in the IL-13/Epi focused model, the tg–IL-13 model, and a conventional allergic asthma model. Inclusion of the focused model allowed us to pinpoint a remarkably small number of gene expression changes that were consistently associated with airway hyperreactivity and mucus production. Furthermore, we used these results to guide translational studies that show the relevance of some of these gene expression changes in people with asthma. These experiments demonstrate that focused transgenic models combined with microarrays can lead to an improved understanding of the pathogenesis of complex diseases such as asthma.

Back to Article Outline

Methods 

Mice 

The University of California San Francisco (UCSF) Committee on Animal Research approved the use of mice for these experiments. Care and use of animals complied with the United States Public Health Service's Policy on Humane Care and Use of Laboratory Animals by Awardee Institutions (#3400-01). Three groups of transgenic mice were used in these experiments: (1) CC10-IL-13+, Stat6+/− mice (tg–IL-13 mice); (2) CC10–IL-13+ Stat6−/− mice; and (3) CC10–IL-13+, Stat6−/−, CC10-hStat6+ mice (IL-13/Epi mice). CC10 refers to the Clara cell specific promoter used to express IL-13 and human STAT6 (hSTAT6).12 The development and characterization of these transgenic mice have been previously described.5 Mice used here were backcrossed 5 times onto the Balb/c genetic background. tg–IL-13 mice have an intact Stat6 gene and IL-13–driven activation of STAT6 in a wide range of cells in the lung, resulting in airway inflammation, mucus overproduction, airway hyperreactivity, subepithelial fibrosis, and emphysema. IL-13+ Stat6−/− mice (used as negative controls) lack STAT6 and did not develop any detectable IL-13–induced lung pathology. IL-13/Epi mice also lack mouse STAT6 but express human STAT6 selectively in airway epithelial cells. hSTAT6 is functional in mice, and epithelial-restricted activation of Stat6 by IL-13 induced airway hyperreactivity and mucus production without airway inflammation, subepithelial fibrosis, or emphysema. Wild-type Balb/c mice 6 to 8 weeks old were used for the ovalbumin challenge model. There were 5 mice in each experimental and control group.

Antigen sensitization and challenge 

Mice were sensitized by intraperitoneal administration of 50 μg grade V ovalbumin mixed with adjuvant (10 mg aluminum potassium sulfate) 3 times at weekly intervals. Control mice received adjuvant alone. Beginning 1 week after the last injection, mice were challenged 3 times by intranasal administration of ovalbumin (1 mg in 50 μL PBS) at daily intervals. Control mice were challenged with PBS alone. Tissues were harvested for isolation of RNA 24 hours after the last challenge.

Isolation and labeling of RNA from mice 

Whole-lung RNAs were purified by using Trizol (Invitrogen, Carlsbad, Calif). Integrity of all RNA samples used in this study was confirmed with a model 2100 bioanalyzer (Agilent Technologies, Inc, Palo Alto, Calif). Cy3-labeled and Cy5-labeled lung cDNAs were prepared as described.13 To obtain samples enriched for airway epithelial cell RNA, the superior portion of the trachea was cannulated, and the trachea and proximal major bronchi were excised from the thorax and slowly perfused with 0.35 mL lysis buffer (RNeasy kit; Qiagen Inc, Valencia, Calif). RNA was isolated from perfusate according to the manufacturer's instructions. Because the amount of RNA obtained from tracheal perfusates was only ∼1 μg, a T7 RNA polymerase-based method13 was used to prepare Cy3-labeled and Cy5-labeled amplified cRNAs for array hybridizations.

Microarray analysis 

Lung gene expression was analyzed by hybridizing Cy5-labeled cDNA from mouse lungs (5 mice per group, each hybridized separately) along with Cy3-labeled reference lung cDNA pooled from wild-type mice. Tracheal perfusate samples were analyzed similarly, except that amplified cRNA from each mouse was compared with an amplified cRNA reference pool made by using tracheal perfusate samples from wild-type mice. DNA microarrays used in these experiments were produced by using the Operon Biotechnologies (Huntsville, Ala) Mouse Genome Oligo 2.0 set of 70-mer oligonucleotides, supplemented by some additional 70-mers. A MIAME-compliant description of the array experiments and the raw array data are available from Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo, accession number GSE1438).

We used an approach that allowed us to estimate differential gene expression between the various groups we studied on the basis of linear models, as previously described.14, 15 To determine whether there were significant differences in gene expression between groups, we calculated the odds ratio (probability of being differentially expressed/probability of not being differentially expressed). When the log2 of the odds ratios (known as the B-statistic) was greater than 0, we classified the gene as differentially expressed.16, 17 Hierarchical clustering, a method for grouping together genes with similar expression patterns, was performed by using Acuity 4.0 software (Axon Instruments, Union City, Calif). In addition, genes were classified on the basis of expression patterns to determine whether they were increased relative to controls in 1, 2, or all 3 of the experimental groups (ovalbumin, tg–IL-13, and IL-13/Epi) as follows. Seven pseudogene vectors were created to represent genes that were increased only in 1 of the 3 models ([0,0,1]; [0,1,0]; and [1,0,0]), genes increased equally in 2 models ([0,1,1]; [1,0,1]; [1,1,0]), and genes increased equally in all 3 models ([1,1,1]). Each differentially expressed gene was assigned a vector in 3-dimensional space according to the median log fold-change gene expression values determined for that gene in the 3 experimental groups. All vectors were scaled to unit length, and each gene was matched to the closest pseudogene, as determined by Euclidean distance.

Human subjects 

These studies were approved by the UCSF Committee on Human Research and conducted in compliance with the Declaration of Helsinki principles. Written informed consent was obtained from all subjects. All subjects were adult nonsmokers (<10 pack-year total smoking history with last cigarette >1 year before the study). Medical histories were obtained, physical examinations were performed, symptom questionnaires were collected, and spirometry was performed as described.18 Airway reactivity was measured by determining the PC20.19 The 28 healthy control subjects had no history of lung disease and were not hyperreactive (PC20>16 mg/mL). All 30 subjects with asthma had a previous physician diagnosis of asthma, were hyperreactive (PC20 < 8 mg/mL), and used only short-acting inhaled β-adrenergic–agonist medications for therapy. Individuals with an asthma exacerbation or respiratory infection within the previous 6 weeks or significant medical problems other than asthma and those using inhaled or systemic corticosteroids or leukotriene antagonists were excluded.

Bronchial epithelial brushings 

Bronchoscopy was performed, and bronchial brushings were obtained randomly from right or left lower lobe bronchial segments by using 4 disposable cytology brushes. The brushes were gently vortexed in sterile saline. Cells from all brushes were pooled, yielding a single sample for each subject. An aliquot was removed for cytocentrifugation, stained with Diff-Quik (Baxter, McGraw Park, Ill), and examined by light microscopy. On average, the bronchial brushings contained 97% epithelial cells. Total RNA was extracted by using the RNeasy kit. Thirty-four of the 58 human epithelial cell RNA samples analyzed for this study were analyzed separately for another study (Woodruff et al, unpublished data).

Cultured human bronchial epithelial cells 

Air-liquid interface cultures were established by using published protocols.20 Primary normal human bronchial epithelial cells (lot 3F1191; Cambrex Bio Science, Baltimore, Md) were seeded onto 12-mm–diameter Corning Transwells (Corning, NY) containing 0.4-μm pores (4-5 replicates per condition). Cells were submerged in bronchial epithelial growth media (Clonetics) until confluent (3 days), and then the apical media was removed and the basolateral media was changed to differentiation media (1:1 of Dulbecco modified Eagle media to bronchial epithelial growth media) for 8 days. After the differentiation period, IL-13 (0-10 ng/mL) was maintained in the basolateral media over the period of the next 4 days. At the end of the IL-13 exposure, total RNA was isolated by using the RNeasy kit.

Analysis of gene expression by real-time PCR 

Primers and a probe for mouse genes (see Table E1 in the Journal's Online Repository at www.mosby.com/jaci) and human genes (see Table E2 in the Online Repository at www.mosby.com/jaci) were designed by using Primer Express software (Perkin Elmer, Boston, Mass). First-strand cDNA synthesis and PCR were performed by using ABI Prizm 7700 or 7900 Sequence Detection Systems (Applied Biosystems, Foster City, Calif). For mouse lung and cultured human airway epithelial cell studies, cycle thresholds for each gene were normalized to those of glyceraldehyde-3-phosphate dehydrogenase (GAPD). For human bronchial brushing samples, a 2-step PCR approach21 was used, and transcript copy numbers were normalized on the basis of the geometric mean expression values of 3 housekeeping genes (GAPD, elongation factor 1α1, and cyclophilin A) as described.22 To identify genes that were differentially expressed between subjects with asthma and normal subjects, we performed a Wilcoxon rank-sum test and calculated the corresponding adjusted P value by using the Westfall and Young23 maxT algorithm implemented in Bioconductor's multtest package (Bioconductor; open source software for bioinformatics, www.bioconductor.org).24

Back to Article Outline

Results 

Microarray analysis of gene expression in mouse asthma models 

We used DNA microarrays to analyze gene expression in 3 murine models of asthma. The first model was composed of mice that were sensitized and challenged with ovalbumin (ovalbumin model). The second model was composed of mice with transgenic overexpression of IL-13 in the lung (tg–IL-13 model); these mice had an intact gene for STAT6, a key signaling molecule required for IL-13 activity. The third model was composed of IL-13 transgenic mice with STAT6 expression limited to nonciliated airway epithelial cells (IL-13/Epi model). In whole-lung samples, 805 gene transcripts were differentially expressed (increased or decreased) in at least 1 model (Fig 1). Of these, 509 genes were increased or decreased by ≥2-fold in at least 1 model. There were 583 gene expression changes in the ovalbumin model and 351 changes in the tg–IL-13 model. In contrast, only 18 changes were seen in the IL-13/Epi focused transgenic model. Comparison of the focused transgenic model to the other models allowed for the identification of a small subset of gene expression changes that are attributable to a specific mechanism and also are relevant to pathogenesis in complex systems.

  • View full-size image.
  • Fig 1. 

    Lung gene expression changes in 3 mouse asthma models. Differentially expressed genes were arranged by hierarchical clustering. Each column represents data from 1 of 5 individual mice in each group. Colors represent fold-change compared with the appropriate controls. The arrow indicates a small group of genes that were increased in all 3 models. Ova, Ovalbumin.

To aid in interpretation, we grouped transcripts according to the magnitude of change in expression across the 3 models. Transcripts were grouped by determining whether the expression changes were seen in 1, 2, or all 3 of the models (see Table E3 in the Online Repository at www.mosby.com/jaci). There were 17 transcripts with at least 2-fold increases in the IL-13/Epi model but only 3 transcripts with at least 2-fold decreases; therefore, we focused our subsequent analyses on increased genes. A total of 239 transcripts were increased in the ovalbumin model but not induced (or induced to a much lesser extent) in the tg–IL-13 or IL-13/Epi models (Fig 2, A). These changes are likely largely attributable to allergen-induced lymphocyte activation and to production of mediators other than IL-13. There were 247 genes increased similarly in the ovalbumin and tg–IL-13 models but not in the IL-13/Epi model (Fig 2, B). These allergen-induced genes are apparently increased because of IL-13 effects on cells other than nonciliated airway epithelial cells, and the large number of gene transcripts is consistent with the idea that effects of IL-13 on these cells are important in allergic inflammatory responses.8, 9 A total of 73 transcripts were increased in the tg–IL-13 model but not increased (or increased to a much lesser extent) in the other 2 models (Fig 2, C). These represent genes that are induced by prolonged high-level overexpression of IL-13 but not by acute allergen challenge and may be involved in the pathogenesis of subepithelial fibrosis and emphysema, pathologic features present in the tg–IL-13 model but not the other 2 models.

  • View full-size image.
  • Fig 2. 

    Gene expression patterns. Grouping revealed genes with increased expression in (A) the ovalbumin (Ova) model only, (B) the Ova and tg–IL-13 models, (C) the tg–IL-13 model only, and (D) all 3 models. The number of genes (left) and representative genes (right) is shown for each group. Phenotypic attributes of each model are shown at the bottom. AHR, Airway hyperreactivity.

We were especially interested in identifying genes that were increased in all 3 models. Analysis of whole-lung samples revealed 35 genes that had similar fold increases in all models (Fig 2, D). Because some epithelial gene expression changes might have been undetectable in whole lung samples, we also analyzed gene expression in tracheal perfusate samples enriched for airway epithelial cell RNA. In these samples, we detected 276 genes with altered expression (increased or decreased) in at least 1 model (see Table E4 in the Online Repository at www.mosby.com/jaci). Of these, 107 genes were changed by ≥2-fold. There were 43 genes changed in the ovalbumin model and 239 in the tg–IL-13 model. There were 76 genes changed in the IL-13/Epi model, but only 22 of these were changed by 2-fold or more. Some genes had similar expression changes in both analyses, but many other changes were identified only in lung samples or only in tracheal samples. In part, this reflects the fact that the tracheal sample was enriched for large airway epithelial cells, whereas the whole-lung sample contained smaller airway epithelial cells as well as many other cell types. In addition, we expect that there were differences in the sensitivity of the lung and tracheal sample analyses because we did not need to use RNA amplification for lung samples but did need to amplify tracheal samples. By combining microarray results from whole-lung and tracheal perfusate samples, we identified 18 genes that were increased at least 2-fold by allergen and by direct effects of IL-13 on airway epithelial cells (Table I). We also included trefoil factor 2 because previous PCR analysis5 suggested that the microarrays underestimated changes in expression of this gene.

Table I. Genes induced by allergen and by direct effects of IL-13 on airway epithelial cells
Fold-change by arraysFold-change by PCR
SymbolDescriptionOvatg–IL-13IL-13/EpiOvatg–IL-13IL-13/Epi
Clca3Chloride channel calcium activated 378.359.828.2820.31530.71448.2
Slc26a4Solute carrier family 26, member 49.929.820.018.740.110.3
RetnlaResistin-like α4.323.823.7NDNDND
ItlnIntelectin7.610.819.29.573.8146.2
NadsynNAD synthetase 117.212.56.40.70.50.8
Atp2a1Ca++ transporting ATPase15.38.03.90.40.10.9
AMCaseAcidic mammalian chitinase2.512.38.63.77.44.4
Muc5acMucin 5, subtypes A and C12.16.83.833.523.47.8
Agr2Anterior gradient 23.87.74.311.121.610.2
Tff1Trefoil factor 13.36.45.752.6592.2224.7
Reg3gRegenerating islet-derived 3γ6.44.13.2NDNDND
PigrPolymeric immunoglobulin receptor2.58.42.31.512.32.2
Muc5bMucin 5, subtype B4.05.12.83.67.413.2
D630002J15RikRIKEN cDNA D630002J15 gene2.84.32.0NDNDND
ScinScinderin2.82.93.021.921.834.5
Alox1515-lipoxygenase3.91.92.05.21.62.8
H2-Q7Histocompatibility 2, Q region locus 72.42.12.1NDNDND
1110001D15RikRIKEN cDNA 1110001D15 gene2.12.12.4NDNDND
Tff2Trefoil factor 24.54.51.83.58.64.9

Genes with a 2.0-fold or greater increase in expression in both the ovalbumin (Ova) allergic model and the IL-13/Epi model are included. Array fold-change values represent the larger median fold change from lung or tracheal perfusate compared with the appropriate control group.

ND indicates that we did not use PCR to determine fold change for these mouse genes, which do not have obvious human orthologues.

PCR validation of gene expression changes in mouse models 

In preparation for translational studies, we generated a validated list of genes with human orthologues that were increased in both ovalbumin and IL-13/Epi mice. Five transcripts listed in Table I did not have obvious human orthologues. Eleven of the remaining 14 transcripts were increased by at least 2-fold in both ovalbumin allergic and IL-13/Epi mice by PCR (Table I). Our mouse modeling approach therefore resulted in the selection of 11 genes that could be studied in people with asthma.

Airway epithelial gene expression in human subjects with and without asthma 

We used real-time PCR to analyze gene expression in airway epithelial cells from 30 subjects with mild to moderate asthma and 28 controls (Table II). There were substantial and highly significant increases in the expression of the calcium activated chloride channel Clca1 (orthologue of murine Clca3) and intelectin in subjects with asthma (Table III and Fig 3, A). There were smaller increases in expression of 15-lipoxygenase and trefoil factor 2 in subjects with asthma. One transcript, mucin 5b, was less abundant in subjects with asthma, and other transcripts examined were not significantly different. We were unable to consistently detect acidic mammalian chitinase transcripts in either group by using 2 different sets of PCR primers and probes.

Table II. Characteristics of human subjects
Control subjectsSubjects with asthmaP
Number2830
Age36 ± 838 ± 13NS
Sex16 F/12 M18 F/12 MNS
PC20 (mg/mL)60.8 ± 1.40.9 ± 1.3<.0001
FEV1/forced vital capacity (%)80.6 ± 1.471.1 ± 1.4<.0001
FEV1 (% predicted)106.5 ± 13.285.2 ± 13.9<.0001

Values represent means ± SDs.

Table III. Gene expression in airway epithelial cells from control subjects and subjects with asthma
Median gene copy number
SymbolDescriptionFold DifferenceControlAsthmaticAdjusted P
Clca1Chloride channel calcium activated 1216.50.06 × 1048.21 × 104.0001
ItlnIntelectin3.41.73 × 1056.12 × 105.0002
Alox1515-lipoxygenase1.32.28 × 1063.27 × 106.06
Tff2Trefoil factor 22.70.55 × 1041.00 × 104.06
Muc5acMucin 5, subtypes A and C1.51.67 × 1072.07 × 107NS
Tff1Trefoil factor 11.41.47 × 1052.17 × 105NS
Agr2Anterior gradient 21.10.90 × 1071.04 × 107NS
Slc26a4Solute carrier family 26, member 40.83.72 × 1043.35 × 104NS
ScinScinderin1.04.70 × 1054.99 × 105NS
Muc5bMucin 5, subtype B−5.01.84 × 1060.42 × 106.0001
AMCaseAcidic mammalian chitinaseNDND

Fold difference between subjects with asthma and control subjects, adjusted for age and sex.

Human Clca1 is orthologous to mouse Clca3.

ND indicates that this transcript was not detected in most samples from both groups.

  • View full-size image.
  • Fig 3. 

    Intelectin (Itln) expression in human airway epithelial cells. A, Intelectin gene expression by airway epithelial cells from control subjects and subjects with asthma. Medians and interquartile ranges are shown. B, IL-13 treatment of cultured human airway epithelial cells induced increased expression of intelectin transcripts.

IL-13 induces intelectin expression in cultured human airway epithelial cells 

One of the transcripts substantially increased in airway epithelial cells from subjects with asthma was intelectin, a pattern recognition molecule that has not been previously implicated in asthma. IL-13 treatment of cultured primary human bronchial epithelial cells led to a large increase in intelectin transcripts (Fig 3, B), demonstrating that IL-13 directly increases intelectin expression and that other cell types are not required.

Back to Article Outline

Discussion 

We combined focused transgenic modeling, functional genomics, and translational studies in human subjects to help understand important aspects of the complex pathogenesis of asthma. As expected, microarray analysis of an allergic model of asthma revealed hundreds of gene expression changes. Analysis of the tg–IL-13 model, where IL-13 acts on many different cell types to produce extensive lung pathology, also showed hundreds of changes, many of which were similar to those seen in the allergic model. By using the data from these 2 models, it would have been extremely challenging to determine how particular gene expression changes relate to specific disease mechanisms or to select a reasonable number of candidate genes for further study. By analyzing a third model (IL-13/Epi) that focused on the effects of a single mediator, IL-13, on a single cell type, the airway epithelial cell, we were able to generate a dramatically shorter list of differentially expressed genes. These expression changes were associated with airway hyperreactivity and mucus overproduction, which are features of all 3 models, but not with inflammation, fibrosis, or emphysema, which are absent in the IL-13/Epi model.

Our approach was designed to reduce the complexity inherent in microarray studies of disease pathogenesis. Other microarray studies have reduced complexity by analyzing homogeneous populations of cultured cells or cells isolated from tissues by laser capture microscopy or flow-cytometric sorting. Although those approaches can be useful, an important advantage of the focused transgenic model approach that we describe here is that it isolates the effects that a single mediator exerts on a single cell type in vivo.

We used translational studies to determine whether the gene expression changes first identified by using mouse models were also seen in human disease. Four of the orthologues identified in the models (the chloride channel Clca1, intelectin, 15-lipoxygenase, and trefoil factor 2) were increased in airway epithelial cells from people with asthma. In another study involving a subset of the human subjects used for this study, we found that expression of Clca1 transcripts and protein was increased in airway epithelial cells from subjects with asthma (Woodruff et al, unpublished data). Previous reports suggest that induction of the mClca3/hClca1 chloride channel in airway epithelial cells is important for mucus production.25, 26, 27 15-Lipoxygenase produces 15S-hydroxyeicosatetraenoic acid, which has been reported to trigger mucus secretion in dogs,28 promote contraction of human airway smooth muscle,29 and potentiate increases in allergen-induced early asthmatic responses in human beings.30 Trefoil factor 2 was shown to be increased in an allergic mouse model of asthma,31 and we now show that trefoil factor 2 is increased in airway epithelial cells from people with asthma. Trefoil factor 2 promotes migration of cultured bronchial epithelial cells32 and contributes to repair of injured intestinal epithelium.33 Intelectin is a recently described pattern recognition molecule not previously implicated in asthma. There are 2 related intelectins, Itln1 and Itln2. The probes and primers that we used were designed to recognize Itln1, but they might not distinguish between these 2 very closely related sequences. The molecular patterns recognized by intelectin include furanosides such as galactofuranose.34 Galactofuranosyl residues are present in bacterial and fungal cell walls and in protozoan parasites but not in mammalian cells. The possible contributions of intelectin to asthma pathogenesis require further exploration, but intelectin in the airway might alter the response of subjects with asthma to infection or colonization with bacterial or fungal pathogens.

We did not find significant increases in expression of the other human homologues that we examined. To some extent, this likely reflects difficulties inherent in detecting gene expression changes in limited samples (epithelial brushings) taken from human subjects with stable mild to moderate disease, as opposed to whole-lung analysis of mice with homogeneous genetic backgrounds housed in a controlled environment and subjected to a potent disease-inducing stimulus. One transcript, Muc5b, was significantly increased in each of the mouse models, but expression of both the transcript and the protein is decreased in people with asthma (Woodruff et al, unpublished data). Differences in Muc5b expression might relate to the fact that human Muc5b is expressed predominantly in submucosal glands.35 These glands are abundant in human airways but are absent in the mouse, with the exception of a single gland in the trachea. This is an example of the limits of mouse modeling. Our PCR assay was unable consistently to detect the transcript for another orthologue, AMCase, but a recent report demonstrated that airway epithelial production of AMCase is also increased in subjects with asthma.36 Including translational studies helped to clarify the relevance of the animal models to human disease and allowed us to draw novel inferences about the activity of a specific mechanism in human disease.

We conclude that the combination of focused transgenic models, DNA microarray analyses, and translational studies provides a powerful approach for analyzing the contributions of specific mediators and cell types and for focusing attention on a limited number of genes associated with specific pathophysiologic aspects of complex diseases like asthma.

Back to Article Outline

 

We thank Dean Sheppard and Andrea Barczak for their advice and Xiaozhu Huang, Louis Nguyenvu, Michael Salazar, and the staffs of the Sandler Center Animal Physiology and Microscopy Core and the UCSF National Heart, Lung, and Blood Institute Shared Microarray Facility for technical assistance.

Back to Article Outline

Appendix. Supplementary data 

Back to Article Outline

References 

  1. Bousquet J, Jeffery PK, Busse WW, Johnson M, Vignola AM. Asthma: from bronchoconstriction to airways inflammation and remodeling. Am J Respir Crit Care Med. 2000;161:1720–1745
  2. Laprise C, Sladek R, Ponton A, Bernier MC, Hudson TJ, Laviolette M. Functional classes of bronchial mucosa genes that are differentially expressed in asthma. BMC Genomics. 2004;5:21
  3. Zimmermann N, King NE, Laporte J, Yang M, Mishra A, Pope SM, et al. Dissection of experimental asthma with DNA microarray analysis identifies arginase in asthma pathogenesis. J Clin Invest. 2003;111:1863–1874
  4. Zou J, Young S, Zhu F, Gheyas F, Skeans S, Wan Y, et al. Microarray profile of differentially expressed genes in a monkey model of allergic asthma. Genome Biol. 2002;3:research0020.1-0020.13
  5. Kuperman DA, Huang X, Koth LL, Chang GH, Dolganov GM, Zhu Z, et al. Direct effects of interleukin-13 on epithelial cells cause airway hyperreactivity and mucus overproduction in asthma. Nat Med. 2002;8:885–889
  6. Huang SK, Xiao HQ, Kleine-Tebbe J, Paciotti G, Marsh DG, Lichtenstein LM, et al. IL-13 expression at the sites of allergen challenge in patients with asthma. J Immunol. 1995;155:2688–2694
  7. Humbert M, Durham SR, Kimmitt P, Powell N, Assoufi B, Pfister R, et al. Elevated expression of messenger ribonucleic acid encoding IL-13 in the bronchial mucosa of atopic and nonatopic subjects with asthma. J Allergy Clin Immunol. 1997;99:657–665
  8. Grunig G, Warnock M, Wakil AE, Venkayya R, Brombacher F, Rennick DM, et al. Requirement for IL-13 independently of IL-4 in experimental asthma. Science. 1998;282:2261–2263
  9. Wills-Karp M, Luyimbazi J, Xu X, Schofield B, Neben TY, Karp CL, et al. Interleukin-13: central mediator of allergic asthma. Science. 1998;282:2258–2261
  10. Zhu Z, Homer RJ, Wang Z, Chen Q, Geba GP, Wang J, et al. Pulmonary expression of interleukin-13 causes inflammation, mucus hypersecretion, subepithelial fibrosis, physiologic abnormalities, and eotaxin production. J Clin Invest. 1999;103:779–788
  11. Nelms K, Keegan AD, Zamorano J, Ryan JJ, Paul WE. The IL-4 receptor: signaling mechanisms and biologic functions. Annu Rev Immunol. 1999;17:701–738
  12. Stripp BR, Sawaya PL, Luse DS, Wikenheiser KA, Wert SE, Huffman JA, et al. cis-acting elements that confer lung epithelial cell expression of the CC10 gene. J Biol Chem. 1992;267:14703–14712
  13. Barczak A, Rodriguez MW, Hanspers K, Koth LL, Tai YC, Bolstad BM, et al. Spotted long oligonucleotide arrays for human gene expression analysis. Genome Res. 2003;13:1775–1785
  14. Diaz E, Yang YH, Ferreira T, Loh KC, Okazaki Y, Hayashizaki Y, et al. Analysis of gene expression in the developing mouse retina. Proc Natl Acad Sci U S A. 2003;100:5491–5496
  15. Jin W, Riley RM, Wolfinger RD, White KP, Passador-Gurgel G, Gibson G. The contributions of sex, genotype and age to transcriptional variance in Drosophila melanogaster. Nat Genet. 2001;29:389–395
  16. Lonnstedt I, Speed TP. Replicated microarray data. Statistica Sinica. 2002;12:31–46
  17. Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3;Available at: http://www.bepress.com/sagmb/vol3/iss1/art3/
  18. Claman DM, Boushey HA, Liu J, Wong H, Fahy JV. Analysis of induced sputum to examine the effects of prednisone on airway inflammation in asthmatic subjects. J Allergy Clin Immunol. 1994;94:861–869
  19. Chai H, Farr RS, Froehlich LA, Mathison DA, McLean JA, Rosenthal RR, et al. Standardization of bronchial inhalation challenge procedures. J Allergy Clin Immunol. 1975;56:323–327
  20. Atherton HC, Jones G, Danahay H. IL-13-induced changes in the goblet cell density of human bronchial epithelial cell cultures: MAP kinase and phosphatidylinositol 3-kinase regulation. Am J Physiol Lung Cell Mol Physiol. 2003;285:L730–L739
  21. Dolganov GM, Woodruff PG, Novikov AA, Zhang Y, Ferrando RE, Szubin R, et al. A novel method of gene transcript profiling in airway biopsy homogenates reveals increased expression of a Na+-K+-Cl- cotransporter (NKCC1) in asthmatic subjects. Genome Res. 2001;11:1473–1483
  22. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:research0034.1-0034.11
  23. Westfall PH, Young SS. Resampling-based multiple testing: examples and methods for p-value adjustment. New York: John Wiley & Sons; 1993;
  24. Ge Y, Dudoit S, Speed TP. Resampling-based multiple testing for microarray data hypothesis. Test. 2003;12:1–44
  25. Hoshino M, Morita S, Iwashita H, Sagiya Y, Nagi T, Nakanishi A, et al. Increased expression of the human Ca2+-activated Cl- channel 1 (CaCC1) gene in the asthmatic airway. Am J Respir Crit Care Med. 2002;165:1132–1136
  26. Nakanishi A, Morita S, Iwashita H, Sagiya Y, Ashida Y, Shirafuji H, et al. Role of gob-5 in mucus overproduction and airway hyperresponsiveness in asthma. Proc Natl Acad Sci U S A. 2001;98:5175–5180
  27. Zhou Y, Shapiro M, Dong Q, Louahed J, Weiss C, Wan S, et al. A calcium-activated chloride channel blocker inhibits goblet cell metaplasia and mucus overproduction. Novartis Found Symp. 2002;248:150–165discussion 65-70, 277-82
  28. Johnson HG, McNee ML, Sun FF. 15-Hydroxyeicosatetraenoic acid is a potent inflammatory mediator and agonist of canine tracheal mucus secretion. Am Rev Respir Dis. 1985;131:917–922
  29. Copas JL, Borgeat P, Gardiner PJ. The actions of 5-, 12-, and 15-HETE on tracheobronchial smooth muscle. Prostaglandins Leukot Med. 1982;8:105–114
  30. Lai CK, Polosa R, Holgate ST. Effect of 15-(s)-hydroxyeicosatetraenoic acid on allergen-induced asthmatic responses. Am Rev Respir Dis. 1990;141:1423–1427
  31. Nikolaidis NM, Zimmermann N, King NE, Mishra A, Pope SM, Finkelman FD, et al. Trefoil factor-2 is an allergen-induced gene regulated by Th2 cytokines and STAT6 in the lung. Am J Respir Cell Mol Biol. 2003;29:458–464
  32. Oertel M, Graness A, Thim L, Buhling F, Kalbacher H, Hoffmann W. Trefoil factor family-peptides promote migration of human bronchial epithelial cells: synergistic effect with epidermal growth factor. Am J Respir Cell Mol Biol. 2001;25:418–424
  33. Farrell JJ, Taupin D, Koh TJ, Chen D, Zhao CM, Podolsky DK, et al. TFF2/SP-deficient mice show decreased gastric proliferation, increased acid secretion, and increased susceptibility to NSAID injury. J Clin Invest. 2002;109:193–204
  34. Tsuji S, Uehori J, Matsumoto M, Suzuki Y, Matsuhisa A, Toyoshima K, et al. Human intelectin is a novel soluble lectin that recognizes galactofuranose in carbohydrate chains of bacterial cell wall. J Biol Chem. 2001;276:23456–23463
  35. Sharma P, Dudus L, Nielsen PA, Clausen H, Yankaskas JR, Hollingsworth MA, et al. MUC5B and MUC7 are differentially expressed in mucous and serous cells of submucosal glands in human bronchial airways. Am J Respir Cell Mol Biol. 1998;19:30–37
  36. Zhu Z, Zheng T, Homer RJ, Kim YK, Chen NY, Cohn L, et al. Acidic mammalian chitinase in asthmatic Th2 inflammation and IL-13 pathway activation. Science. 2004;304:1678–1682

 Supported by National Institute of Health grants HL56835 and HL72301 and by the UCSF Sandler Center for Basic Research in Asthma.Disclosure of potential conflict of interest: D. A. Kuperman, none disclosed. C. A. Lewis, none disclosed. P. G. Woodruff, none disclosed. M. W. Rodriguez, none disclosed. Y. H. Yang, none disclosed. G. M. Dolganov, none disclosed. J. V. Fahy, none disclosed. D. J. Erle, none disclosed.

PII: S0091-6749(05)00598-1

doi:10.1016/j.jaci.2005.03.024

The Journal of Allergy and Clinical Immunology
Volume 116, Issue 2 , Pages 305-311, August 2005