Cancer Research Prevention Award  Metabolism
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online

Cancer Research 68, 34, January 1, 2008. doi: 10.1158/0008-5472.CAN-07-0988
© 2008 American Association for Cancer Research

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stearman, R. S.
Right arrow Articles by Geraci, M. W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Stearman, R. S.
Right arrow Articles by Geraci, M. W.

Molecular Biology, Pathobiology, and Genetics

A Macrophage Gene Expression Signature Defines a Field Effect in the Lung Tumor Microenvironment

Robert S. Stearman1, Lori Dwyer-Nield2, Michael C. Grady1, Alvin M. Malkinson2 and Mark W. Geraci1

Departments of 1 Medicine/Pulmonary Sciences and Critical Care Medicine, and 2 Pharmaceutical Sciences, University of Colorado Health Sciences Center, Denver, Colorado

Requests for reprints: Robert S. Stearman, Box C272, Room BB 3B10, 4200 East Ninth Avenue, Denver, CO 80262. Phone: 303-315-2317; Fax: 303-315-5632; E-mail: Robert.Stearman{at}uchsc.edu.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
One area of intensive investigation is to understand complex cellular and signaling interactions in the tumor microenvironment. Using a novel, although straightforward, microarray approach, we defined a gene expression signature from the lung tumor microenvironment in the murine A/J-urethane model of human lung adenocarcinoma. The tumor microenvironment is reflected by the composition of the cell types present and alterations in mRNA levels, resulting in a "Field Effect" around the tumor. The genes composing the Field Effect expression signature include proteases and their inhibitors, inflammation markers, and immune signaling molecules. By several criteria, the Field Effect expression signature can be attributed to the macrophage lineage, suggesting a qualitative change in the expression pattern of tumor-associated macrophages (TAM) observed in lung tumors. The protein expression levels for a number of Field Effect genes were verified by Western blot analysis of lung homogenates, and for their expression in macrophages and parenchymal cells outside of the tumors by immunohistochemistry. In addition, the Field Effect expression signature was used to classify bronchoalveolar lavage (BAL) cells from tumor-bearing or age-matched control mice. Using a variety of statistical measures, the Field Effect expression signature correctly classified the BAL cells >94% of the time. Finally, the protein levels for several Field Effect genes were higher in cell-free BAL fluid, indicating they may be secreted by the TAMs. This work suggests that TAMs generate a unique gene expression signature within the tumor microenvironment, and this signature could potentially be used for identifying lung cancer from BAL cells and/or fluid. [Cancer Res 2008;68(1):34–43]


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Lung cancer continues to be the number one cause of cancer mortality in the United States, with >210,000 new cases and 160,000 deaths estimated for 2007 (1). Tobacco smoke is the most widespread risk factor accounting for lung cancer in 85% to 90% of cases (2, 3). Environmental risks include occupational exposures to asbestos, uranium, and other mining activities, second-hand tobacco smoke, and household radon. In a study of Iowan women with (cases) or without (controls) lung cancer, the greatest risk factors were current or exsmoker status (adjusted odds ratio, 13.5), and an elevated odds ratio was present even 25 years after smoking cessation (4). Recently released data (1) covering 1975 to 2002, showed a relatively constant 5-year survival rate of ~15% after lung cancer diagnosis, independent of race or gender, with a higher 5-year survival rate (49%) if early stage disease is discovered before detectable lymph node or metastasis involvement. In contrast, colorectal cancer had an overall 5-year survival rate of 64% over the same study period. Although lung cancer is the most common cancer affecting both sexes, it had the third lowest 5-year survival rate, exceeding only liver and pancreatic cancers, which account for ~50,000 cases per year combined.

Survival studies indicate that early detection of lung cancer as stage I disease is the most significant factor in providing the highest 5-year survival rate (1, 5, 6). Unfortunately, there are no effective low-cost options available for general annual screening of the tobacco smoking population, estimated at 20% to 25% of the adults in the United States (1). Two recent studies examining the utility of spiral computed tomography as an annual screening method had conflicting results, although they used different clinical outcome measures (7, 8). Both studies did agree that spiral computed tomography does increase early detection and treatment of lung cancer in the easily identifiable at-risk population (a smoking history and/or occupational exposures; age >40 years). The significance of improved survival rate versus decreased mortality must await the results from two randomized controlled trials currently under way by a Dutch-Belgian collaboration and the National Cancer Institute.

One promise of gene expression microarray studies in lung cancer has been to develop novel approaches to diagnostics for disease classification and patient stratification into high probability treatment outcome groups (912). This overall strategy has found considerable success in the treatment of breast cancer, moving from research tool to clinical diagnostics, to identify best treatment outcome and recurrence using a 70-gene panel (MammaPrint; ref. 13). Recent reports suggest a similar approach is feasible using lung cancer tissue to identify high- and low-risk groups for recurrence and survival (1416), although these approaches do not address early detection methodologies.

As a first step toward identifying biomarkers of early lung cancer, we used the well-characterized A/J mouse-urethane model of human adenocarcinoma (12, 17). We reasoned that histologically appearing, normal lung tissue from urethane-treated mice would have an altered gene expression profile because of proximity to the tumor microenvironment, when compared with normal lung tissue from age-matched untreated animals. The altered gene expression could be derived from lung parenchymal cells and/or immune cells infiltrating and responding to the tumor. In lung cancer, the tumor microenvironment has been extensively studied, especially with regard to a wide variety of chemokine (C-C motif) signaling pathways after activation of nuclear factor-{kappa}B (1820). Chemotactic factors produced by tumors recruit blood-derived monocytes, and within the tumor microenvironment, they mature to tumor-associated macrophages (TAM; refs. 2123).

TAMs play a central role in many tumor-stromal interactions including angiogenesis, extracellular matrix remodeling, and invasion/metastasis (21, 24, 25). The degree of macrophage infiltration into human non–small cell lung cancers is correlated with the microvessel counts and inversely correlated with survival (21). In addition, alveolar-derived macrophages isolated by bronchoalveolar lavage (BAL) from lung cancer patients have altered in vitro characteristics (26). In transgenic mouse models, elevating the number of pulmonary macrophages overexpressing FGF10 in type 2 and Clara cells induced spontaneous development of lung tumors in the absence of exogenous carcinogens (27). Finally, chromosomal sites containing genes that regulate macrophage behavior map to lung cancer susceptibility sites (28). TAMs receive signals from diverse cell types within the tumor microenvironment and release signals that affect a variety of different cell types, including epithelial and endothelial cells, fibroblasts, and lymphocytes (24, 29). TAMs alter epithelial cell phenotype by stimulating proliferation, inhibiting apoptosis, and/or changing their morphology/differentiation status so that initiated epithelial cells adapt a more motile and less sessile configuration (30). Macrophages are phenotypically heterogeneous, and extremes of this continuum are called M1 (classic activation) and M2 (alternative activation; refs. 23, 29). M1 macrophages are stimulated by IFN{gamma} and microbial products (lipopolysaccharide), are bactericidal, and express inducible nitric oxide synthase (iNOS). The M2 state is induced by the cytokines interleukin (IL)-4 and IL-13, produces polyamines, lymphocyte suppression, and arginase expression. In A/J mice, alveolar macrophage phenotype changes from M0 (expressing small or negligible amounts of arginase or iNOS) in the absence of tumors, to M2 (arginase expressing) in adenoma-bearing mice, to M1 (iNOS expressing) in mice with carcinomas (22).

In this study, we compared the gene expression profiles of normal-appearing lung tissue adjacent to tumors with age-matched, untreated controls, and we identified a set of 46 genes specifically altered by the tumor microenvironment. Supporting the observed biology of TAMs, these genes are highly expressed in the monocyte-macrophage cell lineage. The increased protein expression for a number of these candidates was verified by immunohistochemistry of lung cancer tissue. Importantly, these genes had strong predictive value when used to classify expression data from BAL-derived macrophages from urethane-treated A/J mice in an independent study. These results suggest that gene expression information contained within samples from outside the tumor itself (surrogate tissues) are informative for predicting tumor status, and that TAMs play a central role in lung tumorigenesis.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Expression levels of specific genes verified by quantitative real-time PCR or Western blot analyses. The following genes were verified by the methods described below and are listed here using their standardized gene symbol notation3 in alphabetical order for reference: acid phosphatase 5 (Acp5; tartrate resistant), C-C motif ligand 6 (Ccl6), CD68 antigen (Cd68), chitinase 3-like 1 (Chi3l1), cathepsin D (Ctsd), cathepsin K (Ctsk), cathepsin Z (Ctsz), coagulation factor VII (F7), lipoprotein lipase (Lpl), leucine-rich {alpha}-2-glycoprotein 1 (Lrg1), lymphocyte antigen 75 (Ly75), nephroblastoma overexpressed gene (Nov), prostaglandin-endoperoxide synthase 1 [Ptgs1; cyclooxygenase (COX)1], signal-regulatory protein {alpha} (previously designated Ptpns1; Sirpa), secreted phosphoprotein 1 (osteopontin; Spp1).

Murine model of human lung cancer. As previously reported (12), male A/J mice (6–8 weeks old; The Jackson Laboratory) were given a single 1-mg/g i.p. injection of urethane or saline vehicle (age-matched untreated controls). At two time points (early, 24–26 weeks; late, 42 weeks) after urethane or saline treatment, mice were sacrificed by lethal pentobarbital injection, and their lungs excised for examination and molecular characterization. Tumors were dissected from neighboring tissue using a dissection microscope, and the adjacent tissue was placed in RNAlater (Ambion) for microarray analysis or in 20 mmol/L HEPES (pH 7.4), 2 mmol/L EDTA, 2 mmol/L EGTA, 1 mmol/L DTT, 10% glycerol, 5 µg/mL aprotinin, 10 µmol/L leupeptin, and 1 mmol/L 4-(2-aminoethly) benzenesulfonyfluoride for protein analysis. Lungs from age-matched untreated controls underwent similar processing.

RNA isolation and microarray analysis. RNA was isolated from lung tissues or BAL derived cells (below) using Trizol Reagent (Invitrogen) and RNeasy kits (Qiagen), followed by quality control measures using UV spectra characteristics (NanoDrop) and Bioanalyzer (Agilent) for size and integrity of the total RNA (12). Total RNA (2–5 µg from lung tissues or 0.25–1 µg from BAL cells) was converted to fragmented biotin-labeled cRNA by the Affymetrix kit reagents according to the manufacturer (Affymetrix). Labeled cRNAs were hybridized to MG-U74Av2 (lung tissue samples) or MOE430 2.0 (BAL samples) Affymetrix microarrays were washed, developed, and scanned using Fluidics and Scanner Workstations. The lung tissue microarray dataset consisted of four possible different samples: normal lung tissues or tissues adjacent to tumor from two different time points, 24 to 26 (early) or 42 (late) weeks after urethane injection. The number of replicate mice used for lung tissue microarrays was 5 early normal, 8 early adjacent, 7 late normal, and 7 late adjacent. Data quality characteristics of the MG-U74Av2 microarrays were as previously reported (12). The adjacent tissue datasets were deposited at Gene Expression Omnibus4 (31) under the series accession GSE2514 (12). Microarrays completed from BAL-derived cells were from 7 control and 11 urethane-treated A/J mice. Two control BAL RNAs were combined (C6_C13 BAL) to give sufficient total RNA for labeling, and 1 control sample (C36; indicated with an asterisk) was later identified as from a mouse with a spontaneous lung tumor. The MOE430 2.0 microarrays had an average scaling factor of 5.0 (SD = 1.0), average background intensity of 77 (SD = 64), and average percent present called 43% (SD = 5.5%). Affymetrix probe IDs from the various microarrays used in this analysis were annotated and cross-referenced using NetAffx (32). The new microarray datasets described in this work, including .cel and .exp files, have been deposited at the Gene Expression Omnibus under the series accession GSE7269. Processing and statistical analysis of the microarray datasets was completed using BRB-ArrayTools (v3.4.0a; 6/2006) freely available from Dr. Rich Simon (33).5 The BRB-ArrayTools suite makes extensive use of parametric and permutation analyses, as well as estimation of significance by False Discovery Rate. Classifiers are crossvalidated by a "leave one out" strategy by several different methods. More limited sets of probeIDs intensity data were analyzed using Excel XP and GraphPad Prism 4. Cluster and heatmap diagrams were generated in BRB-ArrayTools, which implements the algorithms according to Eisen et al. (34), using centered genes, centered correlations, and average linkage settings.

Measurement of RNA levels by quantitative real-time PCR. All reagent kits, primer/probe sets, and the GeneAmp 5700 Instrument were from Applied Biosystems, Inc. and used according to the manufacturer's instructions. Briefly, 100 ng of total RNA was converted to cDNA using the High Capacity cDNA Reverse Transcription kit in 100 µL reaction and diluted 2-fold with water. Five microliters of cDNA was used per primer/probe (1x concentration) set reaction in TaqMan Universal PCR MasterMix (20 µL final volume) and run for 40 cycles under standard conditions. Mouse β-actin endogenous control primer/probe set was used on all samples and a plasmid DNA dilution standard curve validated amplification efficiencies (r = 0.99; >5 replicates). Ct values are averages of RNA derived from 5 age-matched or 12 urethane-treated A/J mouse BAL cells, each sample run in duplicate. {Delta}Ct value of each gene was calculated relative to its β-actin (lower {Delta}Ct values indicate more abundant transcripts) and fold increase observed in urethane treated versus controls as 2 ({Delta}{Delta}Ct).

Protein biomarker analysis and immunohistochemistry. Protein extract preparation, Western methodology, and antibodies used are described in the Supplementary Method. For immunohistochemistry, A/J mouse lungs were perfused by instilling saline through the pulmonary artery to clear the blood from the lungs. Lungs were inflated with 10% buffered formalin before formalin fixation. Fixed lungs were paraffin embedded and 4 µm sections were cut and affixed to glass microscope slides. Tissue sections were rehydrated, endogenous peroxidases quenched in 1% H2O2 in methanol, and subjected to antigen retrieval in 100 mmol/L sodium citrate buffer in a steamer. After blocking in PBS containing 10% of the appropriate blocking serum, sections were incubated with primary antibody (Supplementary Method) for 90 min in a humidified chamber at 37°C. Proteins were visualized using Vectastain kits according to the manufacturer's instructions (Vector Laboratories) and diaminobenzidine.

BAL. The tracheas were cannulated with an 18-gauge angiocatheter, and 1 mL PBS+0.6 mmol/L EDTA was instilled and aspirated back into the syringe as described previously (35). For cell isolation, lungs were rinsed 5 times with 1 mL PBS, the rinses combined and centrifuged at 500 x g for 10 min. The rinses were then aspirated off, leaving ~500 µL per sample, and the cells transferred to a 1.5 mL microfuge tube, pelleted at 16,000 x g for 1 min, and directly lysed in 400 µL Trizol Reagent (Invitrogen) for RNA purification. For protein analysis, BAL fluid (~600 µL/mouse) was centrifuged at 16,000 x g to remove cells before being concentrated in a Microcon (Millipore) concentrator with a 3 kDa cutoff to ~50 µL/sample. Protein concentration was determined and samples were diluted into protein-loading buffer. Insufficient protein amounts of BAL samples from vehicle-treated mice did not allow comparable SDS-PAGE analyses, so only the samples from 24- and 42-week urethane-treated animals were analyzed for Ctsd, Ctsz, Chi3l1, Spp1, and Sirpa content. Four separate replicate mice were used to produce BAL fluid protein samples. The BAL protein yield from the 24-week urethane-treated mice was substantially lower than that from the 42-week urethane-treated mice. Immunoblot analysis was performed as described above, except that blots were stripped with 1% Tris-HCl (pH 6.7), 0.8% β-mercaptoethanol, and 2% SDS at 50°C for 30 min before reuse.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
We and others have used microarray technology to reveal expression differences between lung cancer and its surrounding tissue, with the goals of tumor classification, therapeutic responsiveness, recurrence, and murine modeling of this disease state (1012, 36, 37). As an extension of the analysis of the A/J-urethane murine model of human lung cancer (12), our approach here was to ascertain the influence of the tumor microenvironment on gene expression in adjacent tissue by comparison to age-, sex-, and environmental exposure–matched normal control A/J mice (Supplementary Fig. S1). Given the extensive literature on the interrelationships between the roles of inflammation processes, immune and stromal cells, and tumorigenesis (19, 21, 28), we hypothesized that a novel gene expression signature would be found in what is histologically appearing normal lung tissue as a result of the influence of the tumor microenvironment (herein called the "Field Effect"). The tumor microenvironment could produce a specific gene expression signature due to immune cell infiltration, endothelial cell alterations (e.g., angiogenesis), and/or the variety of cellular responsiveness to signaling pathways. In addition, a Field Effect signature may aid in biomarker discovery from a more distantly removed surrogate tissue source (38) such as BAL.

Developing an expression signature of the tumor microenvironment. Our microarray dataset consisted of normal or adjacent to tumor A/J mouse lung tissues from two different time points, 24 to 26 or 42 weeks after urethane injection (early versus late time points). Using a four-way classifier algorithm in BRB-ArrayTools (testing for significance in four categories: normal, adjacent, early, and late), a total of 204 probeIDs were identified with parametric P values of ≤0.0001 (see supplementary Excel Spreadsheet for complete signal intensity and annotation data). On the MG-U74Av2 microarray, with ~12,000 probeIDs, <2 probeIDs would be expected by random chance at this level of significance. As shown in a cluster dendrogram (Fig. 1A ), the 204 probeID classifiers generate a tree diagram with 2 main branches [early (yellow) versus late (black) samples; Age] subdivided by proximity to lung tumors [normal (green) versus adjacent (blue); Field]. Two criteria were jointly used to designate whether a probeID was classifying samples based on their proximity to lung tumor (Field) or because the animals in this study were used at two time points (Age). First, on a probeID basis, was the Student's t test for significant differences in the log2 signal intensities with a P value of <0.001, when comparing both potential pairs of samples (normal versus adjacent or early versus late). Second, which t test comparison was smaller (i.e., more statistically significant), the Field or Age comparison. Using both comparisons, the 204 probeIDs were divided into 48 Field probeIDs (24%; see Table 1 for listing of Field Effect genes), 139 Age probeIDs (68%), and 17 not clearly separated by these criteria into either category (8%). The Field probeIDs had significant absolute signal intensities (>64 arbitrary units) and each signal intensity distribution of the probeID was separable by sample type (Supplementary Fig. S2). If the Field Effect genes are reflective of the tumor microenvironment, the magnitude of the difference in signal intensities (adjacent-normal) should increase as the tumor becomes more aggressive at the later time point. A {Delta}{Delta}Signal Intensity graph (Supplementary Fig. S3) shows that 25 of the 43 up-regulated genes in the tumor microenvironment were significantly increased, whereas 3 of the 5 down-regulated genes were significantly decreased at the late time point.


Figure 1
View larger version (54K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 1. Microarray analysis of the A/J-urethane lung cancer model identifies Field Effect classifiers. A, RNAs from murine A/J-urethane lung tissues from either adjacent (A or A42) or age-matched normal controls (N or N42) were analyzed using MG-U74Av2 microarrays (Affymetrix) and the data were classified using a four-way comparison for age (24 or 42 wk after urethane treatment) or source of tissue (adjacent or normal control). The resulting 204 probeID classifier was used for supervised clustering of the datasets and displayed in a dendrogram format. Color coding bars are used to indicate tissue source (green, normal; blue, adjacent) and age (yellow, early, 24 wk postinjection; black, late, 42 wk postinjection). The datasets are distributed into two main branches by age (early versus late), and then each branch subdivided by location of tissue (adjacent versus normal). B, the microarray datasets were analyzed by supervised clustering using the Field Effect probeIDs (n48; adjacent versus normal) identified from the four-way comparison classifier (n204). The data are displayed in a dendrogram format with color coding as in Fig. 1A. Two main branches clearly divide the microarray datasets by adjacent and normal lung tissues. C, the microarray datasets were analyzed by supervised clustering using the Age Effect probeIDs (n139; early, 24 wk postinjection; late, 42 wk postinjection) identified from the four-way comparison classifier (n204). The data are displayed in a dendrogram format with color coding as in Fig. 1A. Two main branches are found, clearly dividing the microarray datasets by early andlate lung tissues. D, gene expression heatmap displays the degree of regulation for the Field Effect probeIDs with the probeIDs down-regulated in the adjacent tissue listed as the top five rows followed by the up-regulated probeIDs. Maximum color saturation was set at ~3-fold change in gene expression. The data are displayed above the heatmap as a dendrogram with the color coding as in Fig. 1A.

 

View this table:
[in this window]
[in a new window]

 
Table 1. Field effect genes

 
The Field and Age t test comparisons gave unambiguous categorization, as the ratio of t tests (e.g., Field t test P value/Age t test P value) were clearly different (median Field probeIDs t test ratio, 30,000; median Age probeIDs t test ratio, 78,000). To visualize the significance of these two sets of probeIDs, supervised clustering of the microarray datasets clearly identified the samples appropriately [Fig. 1B (Field) and C (Age)]. The expression level heatmap of the Field probeIDs shows the majority of genes (90%) are up-regulated (red) in the adjacent tissues, suggesting the presence of a novel gene expression signature (Fig. 1D).

Assessment of cell lineage characteristics within the tumor microenvironment. Two in silico analyses were done to ascertain if the Field genes can be attributed to a likely candidate for cellular origin. The murine Field probeIDs were converted to their appropriate human ortholog probeIDs for querying the SymAtlas database of normal human tissue expression (39). Each probeID was scored for its expression level by tissue type and marked as "overexpression" if the signal intensity of the tissue was significantly higher than its median level (Supplementary Fig. S4). A tissue was scored as overexpressing a gene when the signal intensity was greater than the median of all tissues plus thrice the SD of the measured signal intensities. The Field probeIDs were highly expressed in whole blood, specifically peripheral blood CD14+ monocytes and bone marrow CD33+ myeloid cells, the precursors to tissue matured macrophages. As a further test, the murine probeID orthologs of the human immune cell type specific genes, characterized by Du et al. (40), were used for supervised clustering of the A/J normal and adjacent microarray datasets (Supplementary Fig. S5). In no case were the major immune cell gene clusters able to distinguish the normal from adjacent samples. These results suggest a unique myeloid->monocyte->macrophage lineage for the cells found in the adjacent lung tissue as responsible for the gene expression signature.

Validation of increased protein levels of the Field genes in the tumor microenvironment. The protein level of five different Field genes (Ctsd, Ctsz, Chi3l1, Spp1, and Sirpa) was verified by Western blot analysis of normal and adjacent A/J lung homogenates (Fig. 2 ). In each case, the protein levels were higher in lung homogenates from adjacent tissue from urethane-treated A/J mice compared with untreated normal controls. Our previous work on Ptgs1 (COX1), another Field gene, had increased protein in lung homogenates from adjacent tissue (41).


Figure 2
View larger version (57K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 2. Western blot analysis for protein levels in whole homogenates. Whole homogenates from either normal or adjacent lung tissues were measured for the protein levels of five different Field genes (Ctsd, Ctsz, Chi3l1, Spp1, and Sirpa). In each case, the protein levels were higher in the adjacent tissues compared with the controls. The expression of Ctsz, Chi3l1, Spp1, and Sirpa were significantly increased (P < 0.05) as determined by densitometry measurements of the exposed films.

 
Immunohistochemical analyses for the same five proteins was carried out to determine if the up-regulated expression observed in lung homogenates can be attributed to specific cell types within the tumor microenvironment (Fig. 3 ). Generally, the overall intensity of immunohistochemical staining increased in macrophages at 42 weeks after urethane injection, along with their number and size in the adjacent tissue. Ctsd was present in all macrophages and RBC, as well as occasional Clara cells in all groups. Ctsz was not found in RBC or Clara cells, and was not expressed in all macrophages, but increased macrophage expression was seen in 42-week urethane-treated mouse lung. Chi3l1 protein was highly expressed in the large macrophages found in 42-week urethane-treated mouse lung but present at lower levels in Type 2 and Clara cells. Spp1 was found in macrophages and this expression did not depend on prior urethane-treatment, but its overall level increased in macrophages at 42 weeks after urethane treatment. Sirpa was expressed in Clara cells and macrophages in urethane-treated mice but was absent in these cell types in age-matched controls. Ptgs1 was previously shown to be expressed in Clara and Type 2 cells, and macrophages in adjacent tissue (41). At lower magnification (Supplementary Fig. S6), the tumor-stromal interface can be visualized. Generally, there is dark staining within the individual macrophages surrounding the tumors. The tumors show little staining for Cstz and Sirpa, some light staining for Cstd and Chi3l, but increased staining for Spp1. Interestingly, Spp1 is the only gene in the Field Effect classifier that was also found in our orthologous gene expression comparison of human lung adenocarcinoma and the A/J-urethane murine model (12). This suggests that Spp1 has an increasing expression gradient going from age-matched normal to adjacent to tumor tissue. Thus, 47 of the 48 Field Effect probeIDs are unique classifiers of the tumor microenvironment and are not the result of tumor cell contamination in the adjacent tissue or phagocytosis of tumor cells by the macrophages.


Figure 3
View larger version (83K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3. Immunohistochemical analysis of Field Effect proteins in fixed lung tissues. Formalin-fixed lung tissues from normal (42-wk untreated controls) and urethane-treated (24 or 42 wk after injection) A/J mice were stained for immunohistochemistry analyses of five different Field gene products (Ctsd, Ctsz, Chi3l1, Spp1, and Sirpa). Colored arrows, pointing toward representative cells or locations, show cell type identifications. Intensity of the brownish gray staining was generally highest in 42-wk urethane-treated animals. Twenty-four–week untreated animals were similar to 42-wk normal controls (not shown). Magnification bar, 40 µm.

 
Field Effect classifiers correctly predict BAL samples from tumor-bearing mice. If the Field genes are highly expressed in the vicinity of lung cancer due to infiltrating macrophages, we reasoned that the same genes would be useful in classifying microarray datasets derived from the cells recovered by BAL. As previously found (22), urethane-treated mice have a significant increase in total cells recovered by BAL compared with controls (3.7-fold; P < 0.0001). However, the cell type distribution was indistinguishable between urethane-treated and control mice, with macrophages the major cell type recovered (~95%) followed by lymphocytes (~5%) and neutrophils (~1%). RNA was prepared from cells recovered by BAL from 6 normal untreated controls and 11 urethane-treated A/J mice. These datasets were used to generate a BAL-specific classifier probeID lists (168 probeIDs; P < 0.0001), and a supervised cluster dendrogram using this list is shown in Fig. 4A . Each BAL sample is closely grouped in either tumor-bearing (T; red) or control (C; green) branches of the tree diagram with only sample C36 misidentified (magenta asterisk). On inspection of the formalin-fixed lungs from this mouse, it clearly had developed a spontaneous tumor, something the A/J strain is known to do as it ages (17, 42), suggesting why it was more closely grouped with the BAL samples from treated mice. The same BAL datasets were clustered with supervision using the independently derived Field probeIDs (49 probeIDs; Fig. 4B) producing a similar degree of accuracy as the BAL-derived classifiers. Using the prediction algorithm within BRB-ArrayTools, the Field probeIDs correctly predicted 94% to 100% of the BAL samples using a variety of statistical approaches (compound covariate, nearest neighbors, nearest centroid, and support vector machine; ref. 33). Interestingly, the Field gene and BAL-derived classifiers had only three genes in common: Acp5, Lrg1, and the expressed sequence tag (EST) 1100001G20Rik (Fig. 4C). Although the majority of the Field genes were not selected as BAL classifiers, most had significant absolute signal intensities (>64 arbitrary units) and nonoverlapping signal intensity distributions in the BAL microarray datasets (Supplementary Fig. S7).


Figure 4
View larger version (14K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 4. Field genes correctly classify microarray datasets from cells recovered by BAL. A, RNAs from cells in BAL fluid from urethane-treated or age-matched normal controls were analyzed using MOE430 2.0 microarrays (Affymetrix). The datasets were analyzed using a two-way comparison [tumor-bearing (T) versus control (C)], generating a set of BAL-specific classifier probeID lists (168 probeIDs; P < 0.0001). The microarray datasets are displayed after supervised clustering using the BAL-specific probeIDs with color-coding bars for control (green) or tumor-bearing (red) animals. Sample C36, from a control animal, nearest neighbors treated BAL cells; the formalin-fixed lung of C36 was examined and found to have a spontaneous tumor (magenta asterisks). B, microarray datasets from BAL fluid cells are displayed after supervised clustering using the Field Effect probeIDs (n49). The Field Effect probeIDs correctly identify the origin of the BAL cells, from control or urethane-treated animals, as well as the BAL-specific classifiers. Color coding and labeling are as in A. C, Venn diagram representing the overlap of the BAL-specific classifiers and Field Effect genes showing the Field genes are distinct (46 of 49; 94%) from the BAL classifiers. Only Acp5, Lrg1, and the EST 1100001G20Rik are shared between the two classifier lists.

 
The BAL samples were correctly classified using the Field Effect genes (derived from lung tissue samples), although the majority of these genes did not overlap with the BAL-specific classifiers. For this reason, we verified the BAL microarray signal intensities for a number of different Field genes by quantitative real-time PCR (qRT-PCR) analysis (Fig. 5A and B ). Theoretically, the Field gene classifiers could give accurate predictions of the BAL samples yet have opposite expression levels (e.g., up-regulated Field genes are down-regulated in BAL RNA). The mRNA levels for the 11 genes tested were significantly up-regulated in BAL cells from treated animals compared with control BAL cells, showing 1.5- to 21-fold increases (P ≤ 0.01). The qRT-PCR results further support the in silico assessment (Supplementary Figs. S4 and S5) that the Field Effect signature is attributable to macrophages because 95% of the recovered cells in BAL are macrophages. In addition, Western blot analysis showed increased protein levels in cell-free BAL fluid for the same proteins measured in lung homogenates (Fig. 5C). This suggests that macrophages infiltrating the tumor microenvironment have a unique Field gene expression signature reflected in both the mRNAs of the BAL cells and the secreted proteins in BAL of urethane-treated A/J mice.


Figure 5
View larger version (26K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 5. Validation of increased expression of Field Effect genes in BAL. A, mRNA levels for 11 different Field genes were determined using qRT-PCR with primer/probe pairs. {Delta}Ct values were measured relative to β-actin as the endogenous control mRNA in BAL from 5 normal and 11 treated animals. Smaller {Delta}Ct values indicate higher expression relative to β-actin in the BAL samples. Fold increase (tumor BAL/control BAL) was calculated by the {Delta}{Delta}Ct with the t test P value of <0.02 in all cases. Columns, mean; bars, the SE. B, the protein levels for five different Field genes (Ctsd, Ctsz, Chi3l1, Spp1, and Sirpa) were determined by Western blot analysis of cell-free BAL fluid. In each case, the protein levels were higher in the BAL fluid from treated animals compared with the controls. The expression of Ctsd, Chi3l1, and Spp1 was significantly increased (P < 0.05) as determined by densitometry measurements of the exposed films.

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
As with many cardiovascular diseases, lung cancer risk can easily be assigned to a specific patient population, with the overwhelming risk factors being cigarette smoke inhalation and environmental exposures (1). Unlike cardiovascular diseases, there are few clinically valid methods for identifying individuals with early stage disease, when it is most curable by surgical resection and chemotherapy. The controversial use of spiral computed tomography in asymptomatic smokers (7, 8) also suffers from high monetary cost and relatively poor rate of new patient identification during the extensive follow-up period. An alternative approach uses laser-induced fluorescence endoscope bronchoscopy to localize early stage bronchial lesions (43). Although screening studies continue to show the survival advantage of early diagnosis, no strategy other than direct biopsy exists for the evaluation of peripheral nodules. Our studies were designed to investigate the gene expression signature of the lung tumor microenvironment and evaluate the utility of BAL expression profiling.

In our previous work, we presented evidence using a microarray approach that the murine A/J mouse-urethane system is an excellent animal model of human adenocarcinoma due to the high degree of conservation in the altered gene expression patterns due to the cancer (12). We have now extended the study of the A/J-urethane model to ask if the tumor microenvironment develops a unique gene expression signature within the histologically normal tissue (adjacent to tumor) compared with untreated age-matched controls (normal). The 48 Field probeIDs represents 46 different genes with broadly functional groupings in matrix remodeling (cathepsins D, K, S, and Z; proteinase inhibitors cystatin B and {alpha}1-antitrypsin; and integrin {alpha} X), macrophage markers (CD68 and macrophage-expressed gene 1), and cellular immune functions (Ly75; Ccl6; CD200; colony-stimulating factor 2 receptor, β 1, and low affinity (granulocyte-macrophage); suppressor of cytokine signaling 3; and signal-regulatory protein-{alpha}). In addition, Ptgs1 was found up-regulated within the Field genes, a well-characterized marker for the role of inflammation in cancer (12, 35).

A number of recent reports have identified macrophage infiltration into the lung tumor microenvironment as an important consequence of multiple signaling pathways, including C-C motif and angiogenesis pathways (19, 21, 24). Our data argue strongly for the infiltration of a specific macrophage phenotype into the lung tumor microenvironment. These macrophages produce a gene expression signature (Field Effect genes) attributed to their presence in the tumor microenvironment, perhaps due to functional differences from the normally resident lung macrophages. The "strength" of the Field genes was shown by their predictive abilities in correctly classifying cells recovered by saline BAL of urethane or control A/J mouse lungs. In addition, a number of secreted proteins, predicted to be overexpressed with in the tumor microenvironment, were increased in cell-free BAL fluid. Because there is a high degree of similarity between the A/J-urethane model and human lung adenocarcinoma, we speculate that the Field genes could have predictive value in analyzing cells and/or fluid from BAL of cigarette-smoking patients. Ideally, a unique signature found by gene expression or proteomic analysis will be identified in high-risk asymptomatic individuals with early stage lung cancer.

The influence of immune cell infiltration on metastatic potential and patient survival has been examined in patients with colorectal cancer (44, 45). Instead of using a microarray approach, large-scale flow cytometry (48 CD markers and 10 surface markers/receptors) and qRT-PCR measurements (20 genes) were used to create a panel of cellular immune markers. Differing from our lung cancer results, their studies focused on the adaptive immune response, specifically the effector T cells, which were important in colorectal cancer, rather than inflammatory or immunosuppressive molecules. They found that T-cell markers of migration, activation, and differentiation were increased in colorectal tumors without signs of metastatic invasion. The density, type, and location of immune cells within the colorectal tumors, gave a better predictor of patient survival than the current histologic methods to stage their cancers and was reproduced in two additional patient populations.

The tumor microenvironment is clearly complex in terms of cell types within the milieu, regulation of signaling pathways, and tissue remodeling. In selected cases, potential novel therapeutic targets have been identified by analysis of the tumor microenvironment. Strikingly, the single largest functional group within the lung cancer Field genes is proteases (cathepsins D, K, S, and Z) and their inhibitors (cystatin B and {alpha}1-antitrypsin). The role of the cathepsin protease family in pancreatic cancer has been studied using both transgenic mouse models and pharmacologic agents (46, 47). Their initial findings, using microarray analysis of the RIP1-Tag2 murine pancreatic cancer model, showed higher expression of a number of cathepsins in tumor tissue. In addition, cathepsins are up-regulated in human pancreatic cancer and in a number of other tumor types (48, 49). Small molecule cathepsin inhibitors and specific cathepsin gene knockouts caused decreased pancreatic tumor formation, impaired angiogenesis, and lowered metastatic potential. Although cathepsins are typically lysosomal, in pancreatic cancer, they are found extracellularly at the interface between normal and tumor cells, implicating them in extracellular matrix remodeling and angiogenesis. We found elevated protein levels for cathepsins D and S in BAL fluid from urethane-treated A/J mice, suggesting these proteases are also present in the extracellular environment of lung adenocarcinoma.

In summary, many studies suggest that secreted chemotactic and differentiating factors affect TAMs in the tumor microenvironment (19, 23, 50). TAMs have a distinct phenotype that promotes inflammation, angiogenesis, and matrix remodeling, while suppressing the adaptive immune response. Our microarray analysis has identified a unique gene expression signature (Field genes) from TAMs in the tumor microenvironment that is distinct from normal tissue. The Field genes identified in TAMs were able to correctly classify BAL-derived cells from urethane-treated or control animals, which is one step further removed from direct measurement of lung biopsy tissue. Overexpression of several proteins from the Field gene list was validated in both lung homogenates and cell-free BAL fluid. Four different cathepsins were found as Field genes suggesting their potential role(s) in lung tumorigenesis. Future studies could examine the utility of the Field genes as human classifiers using BAL cells and fluid from control and lung cancer patients. In addition, this work suggests new approaches in murine preclinical studies, investigating the relationship between the cathepsin family members and lung cancer using pharmacologic agents and/or transgenic mouse models.


    Acknowledgments
 
Grant support: FAMRI Clinical Innovator Award (R.S. Stearman), R01 HL 72340 (M.W. Geraci), R01 CA 96133 (A.M. Malkison and M.W. Geraci), R01 CA33497 (A.M. Malkinson), and P50 CA 58187 (M.W. Geraci).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    Footnotes
 
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).

3 http://www.ncbi.nlm.nih.gov/gene Back

4 http://www.ncbi.nlm.nih.gov/geo/ Back

5 http://linus.nci.nih.gov/BRB-ArrayTools.html Back

Received 3/19/07. Revised 8/ 7/07. Accepted 10/ 8/07.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. American Cancer Society. Cancer Facts and Figures 2007. Available from: http://www.cancer.org/docroot/stt/stt_0.asp.
  2. Ruano-Ravina A, Figueiras A, Barros-Dios JM. Lung cancer and related risk factors: an update of the literature. Public Health 2003;117:149–56.[CrossRef][Medline]
  3. Ruano-Ravina A, Figueiras A, Barreiro-Carracedo MA, Barros-Dios J. Occupation and smoking as risk factors for lung cancer: a population-based case-control study. Am J Ind Med 2003;43:149–55.[CrossRef][Medline]
  4. Neuberger JS, Mahnken JD, Mayo MS, Field RW. Risk factors for lung cancer in Iowa women: implications for prevention. Cancer Detect Prev 2006;30:158–67.[CrossRef][Medline]
  5. Mountain CF. Revisions in the international system for staging lung cancer. Chest 1997;111:1710–7.[CrossRef][Medline]
  6. Mountain CF. Staging classification of lung cancer. A critical evaluation. Clin Chest Med 2002;23:103–21.[CrossRef][Medline]
  7. Bach PB, Jett JR, Pastorino U, Tockman MS, Swensen SJ, Begg CB. Computed tomography screening and lung cancer outcomes. JAMA 2007;297:953–61.[Abstract/Free Full Text]
  8. Henschke CI, Yankelevitz DF, Libby DM, Pasmantier MW, Smith JP, Miettinen OS. Survival of patients with stage I lung cancer detected on CT screening. N Engl J Med 2006;355:1763–71.[Abstract/Free Full Text]
  9. Beer DG, Kardia SL, Huang CC, et al. Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med 2002;8:816–24.[Medline]
  10. Hayes DN, Monti S, Parmigiani G, et al. Gene expression profiling reveals reproducible human lung adenocarcinoma subtypes in multiple independent patient cohorts. J Clin Oncol 2006;24:5079–90.[Abstract/Free Full Text]
  11. Meyerson M, Franklin WA, Kelley MJ. Molecular classification and molecular genetics of human lung cancers. Semin Oncol 2004;31:4–19.[Medline]
  12. Stearman RS, Dwyer-Nield L, Zerbe L, et al. Analysis of orthologous gene expression between human pulmonary adenocarcinoma and a carcinogen-induced murine model. Am J Pathol 2005;167:1763–75.[Abstract/Free Full Text]
  13. Buyse M, Loi S, van't Veer L, et al. Validation and clinical utility of a 70-gene prognostic signature for women with node-negative breast cancer. J Natl Cancer Inst 2006;98:1183–92.[Abstract/Free Full Text]
  14. Chen HY, Yu SL, Chen CH, et al. A five-gene signature and clinical outcome in non-small-cell lung cancer. N Engl J Med 2007;356:11–20.[Abstract/Free Full Text]
  15. Lu Y, Lemon W, Liu PY, et al. A gene expression signature predicts survival of patients with stage I non-small cell lung cancer. PLoS Med 2006;3:e467.[CrossRef][Medline]
  16. Potti A, Mukherjee S, Petersen R, et al. A genomic strategy to refine prognosis in early-stage non-small-cell lung cancer. N Engl J Med 2006;355:570–80.[Abstract/Free Full Text]
  17. Malkinson AM. Primary lung tumors in mice as an aid for understanding, preventing, and treating human adenocarcinoma of the lung. Lung Cancer 2001;32:265–79.[CrossRef][Medline]
  18. Aggarwal BB, Takada Y, Shishodia S, et al. Nuclear transcription factor NF-{kappa}B: role in biology and medicine. Indian J Exp Biol 2004;42:341–53.[Medline]
  19. Strieter RM, Belperio JA, Burdick MD, Sharma S, Dubinett SM, Keane MP. CXC chemokines: angiogenesis, immunoangiostasis, and metastases in lung cancer. Ann N Y Acad Sci 2004;1028:351–60.[CrossRef][Medline]
  20. Yamamoto Y, Gaynor RB. I{kappa}B kinases: key regulators of the NF-{kappa}B pathway. Trends Biochem Sci 2004;29:72–9.[CrossRef][Medline]
  21. Chen JJ, Lin YC, Yao PL, et al. Tumor-associated macrophages: the double-edged sword in cancer progression. J Clin Oncol 2005;23:953–64.[Abstract/Free Full Text]
  22. Redente EF, Orlicky DJ, Bouchard RJ, Malkinson AM. Tumor signaling to the bone marrow changes the phenotype of monocytes and pulmonary macrophages during urethane-induced primary lung tumorigenesis in A/J mice. Am J Pathol 2007;170:693–708.[Abstract/Free Full Text]
  23. Sica A, Schioppa T, Mantovani A, Allavena P. Tumour-associated macrophages are a distinct M2 polarised population promoting tumour progression: potential targets of anti-cancer therapy. Eur J Cancer 2006;42:717–27.[CrossRef][Medline]
  24. Mantovani A, Schioppa T, Porta C, Allavena P, Sica A. Role of tumor-associated macrophages in tumor progression and invasion. Cancer Metastasis Rev 2006;25:315–22.[CrossRef][Medline]
  25. Strieter RM, Belperio JA, Phillips RJ, Keane MP. CXC chemokines in angiogenesis of cancer. Semin Cancer Biol 2004;14:195–200.[CrossRef][Medline]
  26. Pouniotis DS, Plebanski M, Apostolopoulos V, McDonald CF. Alveolar macrophage function is altered in patients with lung cancer. Clin Exp Immunol 2006;143:363–72.[CrossRef][Medline]
  27. Clark JC, Tichelaar JW, Wert SE, et al. FGF-10 disrupts lung morphogenesis and causes pulmonary adenomas in vivo. Am J Physiol Lung Cell Mol Physiol 2001;280:L705–15.[Abstract/Free Full Text]
  28. Malkinson AM. Role of inflammation in mouse lung tumorigenesis: a review. Exp Lung Res 2005;31:57–82.[Medline]
  29. Mantovani A, Sica A, Sozzani S, Allavena P, Vecchi A, Locati M. The chemokine system in diverse forms of macrophage activation and polarization. Trends Immunol 2004;25:677–86.[CrossRef][Medline]
  30. Lewis CE, Pollard JW. Distinct role of macrophages in different tumor microenvironments. Cancer Res 2006;66:605–12.[Abstract/Free Full Text]
  31. Barrett T, Suzek TO, Troup DB, et al. NCBI GEO: mining millions of expression profiles-database and tools. Nucleic Acids Res 2005;33:D562–6.[Abstract/Free Full Text]
  32. Liu G, Loraine AE, Shigeta R, et al. NetAffx: Affymetrix probesets and annotations. Nucleic Acids Res 2003;31:82–6.[Abstract/Free Full Text]
  33. McShane LM, Radmacher MD, Freidlin B, Yu R, Li MC, Simon R. Methods for assessing reproducibility of clustering patterns observed in analyses of microarray data. Bioinformatics 2002;18:1462–9.[Abstract/Free Full Text]
  34. Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 1998;95:14863–8.[Abstract/Free Full Text]
  35. Bauer AK, Dwyer-Nield LD, Hankin JA, Murphy RC, Malkinson AM. The lung tumor promoter, butylated hydroxytoluene (BHT), causes chronic inflammation in promotion-sensitive BALB/cByJ mice but not in promotion-resistant CXB4 mice. Toxicology 2001;169:1–15.[CrossRef][Medline]
  36. Kakiuchi S, Daigo Y, Ishikawa N, et al. Prediction of sensitivity of advanced non-small cell lung cancers to gefitinib (Iressa, ZD1839). Hum Mol Genet 2004;13:3029–43.[Abstract/Free Full Text]
  37. von Eyben FE. Epidermal growth factor receptor inhibition and non-small cell lung cancer. Crit Rev Clin Lab Sci 2006;43:291–323.[CrossRef][Medline]
  38. Stearman RS, Keith RL, Dwyer-Nield L, et al. Expression analysis in lung cancer: the search for biomarkers. Am J Respir Cell Mol Biol 2004;31:S39–44.
  39. Su AI, Cooke MP, Ching KA, et al. Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci U S A 2002;99:4465–70.[Abstract/Free Full Text]
  40. Du X, Tang Y, Xu H, et al. Genomic profiles for human peripheral blood T cells, B cells, natural killer cells, monocytes, and polymorphonuclear cells: comparisons to ischemic stroke, migraine, and Tourette syndrome. Genomics 2006;87:693–703.[CrossRef][Medline]
  41. Bauer AK, Dwyer-Nield LD, Malkinson AM. High cyclooxygenase 1 (COX-1) and cyclooxygenase 2 (COX-2) contents in mouse lung tumors. Carcinogenesis 2000;21:543–50.[Abstract/Free Full Text]
  42. Shimkin MB, Stoner GD. Lung tumors in mice: application to carcinogenesis bioassay. Adv Cancer Res 1975;21:1–58.[Medline]
  43. Hirsch FR, Franklin WA, Gazdar AF, Bunn PA, Jr. Early detection of lung cancer: clinical perspectives of recent advances in biology and radiology. Clin Cancer Res 2001;7:5–22.[Abstract/Free Full Text]
  44. Galon J, Costes A, Sanchez-Cabo F, et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science 2006;313:1960–4.[Abstract/Free Full Text]
  45. Pages F, Berger A, Camus M, et al. Effector memory T cells, early metastasis, and survival in colorectal cancer. N Engl J Med 2005;353:2654–66.[Abstract/Free Full Text]
  46. Gocheva V, Zeng W, Ke D, et al. Distinct roles for cysteine cathepsin genes in multistage tumorigenesis. Genes Dev 2006;20:543–56.[Abstract/Free Full Text]
  47. Joyce JA, Baruch A, Chehade K, et al. Cathepsin cysteine proteases are effectors of invasive growth and angiogenesis during multistage tumorigenesis. Cancer Cell 2004;5:443–53.[CrossRef][Medline]
  48. Roshy S, Sloane BF, Moin K. Pericellular cathepsin B and malignant progression. Cancer Metastasis Rev 2003;22:271–86.[CrossRef][Medline]
  49. Turk V, Turk B, Guncar G, Turk D, Kos J. Lysosomal cathepsins: structure, role in antigen processing and presentation, and cancer. Adv Enzyme Regul 2002;42:285–303.[CrossRef][Medline]
  50. Stout RD, Jiang C, Matta B, Tietzel I, Watkins SK, Suttles J. Macrophages sequentially change their functional phenotype in response to changes in microenvironmental influences. J Immunol 2005;175:342–9.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
CarcinogenesisHome page
F. Colotta, P. Allavena, A. Sica, C. Garlanda, and A. Mantovani
Cancer-related inflammation, the seventh hallmark of cancer: links to genetic instability
Carcinogenesis, July 1, 2009; 30(7): 1073 - 1081.
[Abstract] [Full Text] [PDF]


Home page
BloodHome page
Z. Li, Y. Liu, S. Tuve, Y. Xun, X. Fan, L. Min, Q. Feng, N. Kiviat, H.-P. Kiem, M. L. Disis, et al.
Toward a stem cell gene therapy for breast cancer
Blood, May 28, 2009; 113(22): 5423 - 5433.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
S. Dubey and C. A. Powell
Update in Lung Cancer 2008
Am. J. Respir. Crit. Care Med., May 15, 2009; 179(10): 860 - 868.
[Full Text] [PDF]


Home page
BloodHome page
H. Chen, R. A. Campbell, Y. Chang, M. Li, C. S. Wang, J. Li, E. Sanchez, M. Share, J. Steinberg, A. Berenson, et al.
Pleiotrophin produced by multiple myeloma induces transdifferentiation of monocytes into vascular endothelial cells: a novel mechanism of tumor-induced vasculogenesis
Blood, February 26, 2009; 113(9): 1992 - 2002.
[Abstract] [Full Text] [PDF]


Home page
Proc Am Thorac SocHome page
T. Walser, X. Cui, J. Yanagawa, J. M. Lee, E. Heinrich, G. Lee, S. Sharma, and S. M. Dubinett
Smoking and Lung Cancer: The Role of Inflammation
Proceedings of the ATS, December 1, 2008; 5(8): 811 - 815.
[Abstract] [Full Text] [PDF]


Home page
Cancer Prevention ResearchHome page
K. Steiling, J. Ryan, J. S. Brody, and A. Spira
The Field of Tissue Injury in the Lung and Airway
Cancer Prevention Research, November 1, 2008; 1(6): 396 - 403.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stearman, R. S.
Right arrow Articles by Geraci, M. W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Stearman, R. S.
Right arrow Articles by Geraci, M. W.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online