Abstract
Cancer-associated fibroblasts (CAF) have been reported to support tumor progression by a variety of mechanisms. However, their role in the progression of non–small cell lung cancer (NSCLC) remains poorly defined. In addition, the extent to which specific proteins secreted by CAFs contribute directly to tumor growth is unclear. To study the role of CAFs in NSCLCs, a cross-species functional characterization of mouse and human lung CAFs was conducted. CAFs supported the growth of lung cancer cells in vivo by secretion of soluble factors that directly stimulate the growth of tumor cells. Gene expression analysis comparing normal mouse lung fibroblasts and mouse lung CAFs identified multiple genes that correlate with the CAF phenotype. A gene signature of secreted genes upregulated in CAFs was an independent marker of poor survival in patients with NSCLC. This secreted gene signature was upregulated in normal lung fibroblasts after long-term exposure to tumor cells, showing that lung fibroblasts are “educated” by tumor cells to acquire a CAF-like phenotype. Functional studies identified important roles for CLCF1–CNTFR and interleukin (IL)-6–IL-6R signaling in promoting growth of NSCLCs. This study identifies novel soluble factors contributing to the CAF protumorigenic phenotype in NSCLCs and suggests new avenues for the development of therapeutic strategies. Cancer Res; 72(22); 5744–56. ©2012 AACR.
Introduction
Cancer is initiated and progresses within a microenvironment that is itself altered as a consequence of the tumorigenic process. Stromal cells in contact with cancer cells secrete growth factors and cytokines that may act directly by signaling to tumor cells or indirectly by recruiting other stromal components to promote tumor progression (1–3). An important aspect of this process is the expansion of cancer-associated fibroblasts (CAF; refs. 4–6). CAFs are a diverse population of stromal cells with distinct characteristics in different tumors and tissues. To date, they have been best characterized in human breast and prostate cancer (7, 8) as well as in several mouse cancer models (9–11). However, the exact origins of these cells and how they acquire their tumor-promoting capacities remains poorly understood. Several studies have identified key CAF-secreted molecules that play a role in tumor progression. For example, in human breast cancer, CAFs recruit endothelial cells and signal directly to tumor cells via the cytokine SDF-1 (12). In mouse dermal, cervical and islet cell cancers, CAFs secrete a proinflammatory array of cytokines in response to NF-κB activation (11). CAFs have also been reported to modulate the self-renewal and proliferation of cancer stem cells in vitro and in vivo (13, 14). These observations link CAFs to the regulation of multiple processes that control tumor development. Thus, CAF–tumor cell interactions play an important role in cancer progression and represent a relatively unexplored target for cancer therapy. Consequently, understanding the precise molecular mechanisms that lead to tumor promotion by CAFs is an important goal in cancer biology.
Lung cancer is one of the leading causes of cancer death in the world. Non–small cell lung cancer (NSCLC) accounts for most lung cancer cases and the overall 5-year survival of patients with this disease remains approximately 15% (15). Most NSCLC tumors are characterized by a significant desmoplastic/fibrotic component at the time of diagnosis (16, 17). Several studies have suggested a role for lung fibroblasts in NSCLCs using in vitro models (18–20). Recently, CAF-secreted VEGF has been suggested to promote tumor invasion (21). Whether other secreted proteins also play a role in lung CAF-induced tumor proliferation has not been reported. Here, we report a systematic cross-species study of the molecular phenotype of normal lung fibroblasts and CAFs and identify novel mediators of tumor–stroma cross-talk in NSCLCs.
Oncogenic Kras-driven mouse models of lung cancer recapitulate many features of human NSCLCs (22), including oncogene-specific gene signatures (23). In this study, normal mouse lung fibroblasts and CAFs isolated from a Kras-driven lung adenocarcinoma model were used to identify individual genes and gene signatures relevant to the pathogenesis of lung cancer. While both normal mouse lung fibroblasts and CAFs contribute to tumor progression, lung CAFs fostered tumor growth in vivo more efficiently than normal mouse lung fibroblasts. A gene signature specific to lung CAFs including multiple secreted proteins and genes related to inflammatory processes was identified and shown to be a marker of poor survival in patients with NSCLCs. Mouse and human lung fibroblasts upregulated genes in this gene signature after exposure to tumor cells, strongly suggesting that CAFs derive from normal resident lung fibroblasts. Finally, a critical and nonredundant functional role for one member of the CAF-secreted gene signature, CLCF1, was confirmed using both knockdown and overexpression studies. While interleukin (IL)-6 also contributes to tumor growth, signaling through the IL-6 receptor (IL-6R) is dispensable for tumor growth. In contrast, blockade of CLCF1–CNTFR paracrine signaling significantly decreases the growth of tumor cells in vivo. Thus, these studies underscore a previously uncharacterized role for CLCF1 in paracrine signaling in NSCLCs.
Materials and Methods
Reagents
Open reading frames for lacZ and clcf1 were obtained from Open Biosystems (Thermo Fisher). Il6 cDNA was synthesized from mouse CAF RNA with the DyNAmo cDNA Synthesis Kit (New England Biolabs). All open reading frames were PCR-amplified using reverse primers that lacked the stop codon for each gene, cloned into an entry D-TOPO vector, and transferred to a pL6/V5-DEST lentiviral vector using Gateway technology (Invitrogen). Mouse-raised antibodies to α-smooth muscle actin, vimentin, cytokeratin, β-actin (Sigma), and V5 (Invitrogen), and rabbit-raised antibodies to β-actin (LiCOR) were used.
Cell lines
LKR10 and LKR13 cells are mouse lung cancer cell lines derived by serial passage of minced lung adenocarcinoma tissues from 2 tumors isolated from separate lobes of the same KrasLA1 mouse (24). These cells were a gift of Dr. Julien Sage (Stanford University School of Medicine, Stanford, CA). A549 and H1299 cells were from American Type Culture Collection. LSZ2 cells were derived through xenograft passages from KrasLSLG12D mice and established as cell lines in vitro. All cell lines were grown in Dulbecco's' Modified Eagle's Media (DMEM) supplemented with 10% FBS and 1% penicillin/streptomycin. MLE12 cells were purchased from American Type Culture Collection and grown according to vendor's specifications. These cells were originally isolated by Wikenheiser and colleagues (25).
Isolation of normal lung fibroblasts and CAFs
Mouse fibroblasts were isolated as described (26). Cells could be grown in vitro for 4 to 6 passages before replicative senescence became apparent under these culture conditions. Human normal lung fibroblasts and CAFs were isolated from normal lung or lung adenocarcinomas obtained with approval from the Stanford Institutional Review Board (Stanford, CA). Tissues were digested to obtain single cells and grown in 96-well culture plates before expansion for different analyses. For immortalization of normal mouse lung fibroblasts, cells were transfected using a 1:10 ratio of pEmpty Puro and pSV40Tag vectors using Lipofectamine (Invitrogen) and following manufacturer's specifications. The vectors were a gift of Dr. Steve Artandi (Stanford University School of Medicine, Stanford, CA).
Retro and lentiviral infections
Virus was produced by transfection into 293FT cells of a retroviral vector expressing eGFP (27) using FuGENE (Roche) according to manufacturer's recommendations. A pLKO.1 vector expressing mCherry was constructed by replacing the puromycin cassette with mCherry cDNA (Clontech). pLKO and pL6 vectors were used to produce lentivirus by transfection into 293FT cells as previously described (28). Viruses were filtered and applied directly to cells for infection at a multiplicity of infection (MOI) < 1.
Cell proliferation assay
Cells were plated into 96-well plates and treated with MTT. Experiments were read on the indicated days according to manufacturer's instructions (Cell proliferation Kit I, Roche). Data were normalized to day 0 of experiment. The 3T3 assays were as described (29) except for normal mouse lung fibroblasts and CAFs, 200,000 cells were plated and counted every 3 days.
Co-injection experiments
Xenograft experiments were carried out in accordance with Stanford Institutional Animal Care and Use Committee guidelines. For experiments with mouse lung cancer cells, 5 × 103 LKR10, LKR13, or LSZ2 were resuspended, alone or in combination with 5 × 104 normal mouse lung fibroblasts or CAFs (ratio 1:10), in 100 mL of serum-free DMEM and injected subcutaneously into the 2 lower flanks of athymic Balb/cnu/nu mice (Charles River). Beginning 1 week after injection, tumor dimensions were measured every 3 to 4 days, and tumor volume was calculated using volume = π/6 × (length) × (width)2. For human experiments, 5 × 105 A549 or H1299 cells were injected alone or in combination with 5 × 106 normal mouse lung fibroblasts or iCAFs.
Conditioned media and coculture experiments
Conditioned media from normal fibroblasts/CAFs: 1 × 106 lung fibroblasts and CAFs were plated into 10 cm and media were changed on the following day. Fresh media (DMEM, 10% FBS, and 1% penicillin/streptomycin) was added, and cells were grown for another 2 days. Media was then removed and concentrated using Amicon Ultra columns (Millipore). Total protein was determined using the BCA Kit (Pierce) and indicated concentrations were used to treat LKR10 or LKR13 cells. Conditioned media from LKR10: 2 × 106 cells were plated and grown in fresh media for 5 days, and conditioned media was harvested; 1:1 ratio of filtered conditioned media plus fresh media was used to treat normal mouse lung fibroblasts for 36 hours. For coculture, mCherry–normal fibroblasts and LKR10 cells were mixed in a 1:2 ratio for 5 or 15 days. Then, mCherry-positive cells were sorted out for gene expression analysis. Normal mouse lung fibroblasts were cocultured with CAFs for 5 or 15 days, and mCherry-negative cells sorted for mRNA analysis. Human normal lung fibroblasts were cocultured with GFP-expressing A549 cells at a 1:2 ratio for 12 and 24 days, followed by sorting of GFP-negative cells for gene expression analysis. G10 and G13 cells were cocultured with normal mouse lung fibroblasts overexpressing clcf1, il6, and control cells at a 1:5 ratio for 7, 14, and 21 days.
Western blotting
Cells were scraped, washed, and lysed in NP-40 buffer supplemented with protease inhibitor cocktail (Roche), 25 mmol/L sodium fluoride and 1 mmol/L sodium orthovanadate (Sigma). Protein samples were resolved by SDS-PAGE and transferred to Immobilon-FL membranes (Millipore) and incubated in blocking buffer (Li-COR Biosciences) for 1 hour before addition of primary antibodies. Antibody incubations and washing steps were done according to Li-COR protocols. Membranes were scanned using an Odyssey Imager (Li-COR Biosciences).
Immunofluorescence and trichrome staining
Cells were trypsinized and grown on glass slides overnight. Cells were fixed using 4% paraformaldehyde, processed, and stained with indicated primary antibodies. AlexaFluor secondary antibodies were from Invitrogen. Bioquant software was used to montage the entire xenograft sections and calculate the percentage of trichrome staining (Bioquant Inc.). Images were taken using an Eclipse E800 microscope (Nikon).
Quantitative reverse transcriptase-PCR
RNA was isolated using TRIzol reagent (Invitrogen) following manufacturer's specifications and cleaned-up using Qiagen miniRNA columns (Qiagen). cDNA was synthesized with a DyNAmo cDNA Synthesis Kit (New England Biolabs) and quantitative reverse transcriptase-PCR (qRT-PCR) was carried out using SYBRGreen (Applied Biosystems; see Supplementary Information for primer sequences).
Gene expression analysis
The Robust Multichip Average (RMA) Express 1.0.4 program (30) was used for background adjustment and quantile RMA normalization of the 45,101 probe sets encoding mouse genome transcripts. Linear Models for Microarray data (LIMMA; ref. 31) was used for comparison between the groups on RMA-normalized signal intensities. After preprocessing steps, a moderated Student t test was used to select differentially expressed genes between groups (CAFs and normal mouse lung fibroblasts). P values were adjusted for multiple testing using the Benjamini–Hochberg method. Fold change was calculated by comparing gene expression levels between CAFs and normal mouse lung fibroblasts samples and expressed as the ratio between the averages of normalized intensities of the 2 groups. Only genes with fold change threshold greater than 2 and adjusted P < 0.075 were considered as differentially expressed. The overexpressed genes identified from LIMMA were analyzed using Ingenuity Pathway Analysis (IPA) software (32). For network analysis, we chose to include only direct interactions in IPA. We used Onto-Express Ensemble (33, 34) for pathway analysis and Gene Ontology (GO) studies of the overexpressed genes from LIMMA. P values were computed using the hypergeometric distribution and adjusted for multiple hypotheses using Benjamini–Hochberg correction.
Meta-analysis
We conducted meta-analysis of 5 human NSCLC data sets obtained using DNA microarrays. Samples included tumor and stroma. Details of the meta-analysis procedure will be described elsewhere (Khatri and colleagues, manuscript in preparation). Briefly, we estimated the effect size for each gene in each data set as Hedges' adjusted g (35). The study-specific effect sizes for each gene were then combined into a single meta-effect size using a linear combination of study-specific effect sizes, fi, where each study-specific effect size was weighted by inverse of the variance in the corresponding study (Eq. A). After computing meta-effect size, significant genes were identified using Z-statistic, and p-values were corrected for multiple hypotheses testing using Benjamini–Hochberg false discovery rate (FDR) correction.
Survival analysis
Survival analysis was conducted on both gene sets and individual genes. Four different datasets were used (36). Data were generated from RNA samples that included both tumor and stromal populations. Only samples at N0 and N1 stage and with time period less than 60 months for the study were analyzed. Each data set was normalized independently and standardization (mean centering), and scaling (dividing the centered value by the standard deviation) was conducted before uni- and multivariate analysis. High and low expressed groups were determined by taking the average of the standardized values across all samples for a specific gene. For gene sets, a summation of all the genes for a particular sample was calculated, and the average of the sum values was used to determine the high and low groups. Samples that had a value less than the calculated average were assigned as the low-risk group and the rest to high-risk group. Cox HRs were calculated for both uni- and multivariate analysis to estimate the relative risk along with confidence interval (CI) and log-rank tests.
Statistical analyses
A 2-tailed t test was used for comparisons of groups. Error bars correspond to either SD or SEM. Significant P values in the text correspond to <0.05 (*), <0.01 (**), or <0.001 (***).
Results
CAFs promote lung tumor development in vivo
To compare the effect of normal mouse lung fibroblasts and lung CAFs on lung cancer progression, we used cells derived from the KrasLA1 mouse model of lung adenocarcinoma (24). KrasLA1 lung tumors contained a significant fibrotic component, suggesting a role for CAFs in the progression of these tumors (Fig. 1A). CAFs isolated from tumors in this model and normal lung fibroblasts from littermate wild-type mice were expanded in vitro for approximately 7 to 10 days before gene expression analysis or in vivo experiments (see Materials and Methods). Expression of fibroblast-specific protein-1 (fsp-1; ref. 37) was detected in both normal mouse lung fibroblasts and CAFs, whereas markers of other cells types were very low or undetectable (Supplementary Fig. S1A). Both normal mouse lung fibroblasts and CAFs showed expression of mesenchymal markers [α-smooth muscle actin (α-sma) and vimentin] by immunofluorescence (Fig. 1B) whereas they did not express cytokeratins (Supplementary Fig. S1B). As expected (12, 21), expression levels of α-sma were higher in CAFs than in normal mouse lung fibroblasts. It has been suggested that CAFs may be derived from epithelial cells that have undergone an epithelium-to-mesenchyme transition (38). However, KrasLA1 CAFs did not show recombination of the Kras locus, which occurs in epithelial tumors in this model (Supplementary Fig. S1C), showing that they are not of tumor origin. Finally, neither normal mouse lung fibroblasts nor CAFs formed subcutaneous tumors when injected alone (data not shown).
Lung CAFs stimulate lung cancer growth in vivo. A, representative image of a lung adenocarcinoma section from a KrasLA1 mouse stained with Masson's trichrome. Scale bar, 300 μm. Bar graph represents the percentage of the stromal component in lung adenocarcinomas (n = 4). B, immunofluorescence on normal fibroblasts (NFs) and CAFs using antibodies against SMA and vimentin. Scale bars, 200 μm. C and D, LKR10 and LKR13 xenografts injected alone or with either NFs or CAFs. Endpoint is 22 days for LKR10 and 20 days for LKR13. Graphs show average tumor volume at indicated time points (n = 10 per group). Error bars indicate ±SEM. E, representative images of Masson's trichrome staining of xenografts derived from tumor cells alone (LKR13) or tumor cells plus NFs or CAFs. Bar graph shows percentage of stromal component for each type of tumor (n = 4 tumors per group). Scale bars, 300 μm. F, xenograft experiment of LKR13 injected alone or with either NFs or CAFs that have been passaged in vitro for >2 weeks. Endpoint is 24 days. Graphs show averaged tumor volume at indicated time points (n = 10 injections per group). Error bars indicate ±SEM. G and H, tumor cells LKR10 (G) and LKR13 (H) were treated with 2.5 μg/mL of conditioned media from NFs and CAFs, and cell viability was analyzed by MTT 72 hours after treatment. Results are representative of 3 independent experiments using conditioned media from 3 different NF and CAF lines. All P values are for a 2-tailed t test. *, P < 0.05; **, P < 0.01; and ***, P < 0.001.
Next, mouse lung cancer cell lines derived from the same KrasLA1 model (LKR10 and LKR13) were grafted alone or co-injected with either normal mouse lung fibroblasts or CAFs (ratio 1:10) into immunodeficient mice. Tumors harboring CAFs were significantly larger than tumors harboring only tumor cells. Interestingly, tumors harboring normal mouse lung fibroblasts were also larger than tumors lacking fibroblasts, although significantly smaller than those harboring CAFs (Fig. 1C and D). Normal mouse lung fibroblasts and CAFs had similar proliferative capacity in vitro (Supplementary Fig. S1D). In addition, the fibrotic component was similar in tumors that incorporated normal mouse lung fibroblasts or CAFs (Fig. 1E) and closely recapitulated that of the KrasLA1 tumors (Fig. 1A). Thus, the differential contribution to tumor growth of CAFs and normal mouse lung fibroblasts is not likely to be due to a difference in their proliferative potential in vivo. The ability of normal fibroblasts and CAFs to influence tumor growth in vivo was maintained after 2 weeks in culture (Fig. 1F), showing that the tumor-promoting capacity of normal mouse lung fibroblasts and CAFs is stable over time and independent of ongoing contact with tumor cells, at least after 2 weeks in culture. Therefore, resident lung fibroblasts contribute to tumorigenesis in vivo and CAFs provide a further advantage to tumor growth.
CAFs have been proposed to promote tumor growth through direct interaction with tumor cells or by recruitment of other cells to the tumor microenvironment (11, 12). To determine whether CAFs from KrasLA1 mice directly stimulate the growth of tumor cells through secreted factors, LKR10 and LKR13 cells were treated in vitro with conditioned media from normal mouse lung fibroblasts or CAFs. Proliferation of both lung cancer cell lines was significantly increased by conditioned media from CAFs when compared with that of normal mouse lung fibroblasts or control medium (Fig. 1G and H). This suggests that CAFs act, at least in part, by secreting soluble ligands that directly stimulate the growth of tumor cells.
CAFs are characterized by a secretory and inflammatory expression profile
Gene expression microarrays were used to identify changes between normal mouse lung fibroblasts and CAFs that may contribute to their different tumor-enhancing capacity. Unsupervised hierarchical clustering separated normal mouse lung fibroblasts and CAF gene expression into 2 distinct clusters (Supplementary Fig. S2A). One hundred sixty-four genes were upregulated in CAFs compared with normal mouse lung fibroblasts (CAF-HIGH), whereas only 11 genes were downregulated in CAFs (CAF-LOW; fold change > 2 or <0.5 and FDR < 0.075; Supplementary Fig. S2B and Supplementary Table S1). Interrogation of the molecular pathways represented in CAF-HIGH genes using Onto-Express (34) identified cytokine–cytokine receptor interaction (P = 5.04 × 10−4), apoptosis (P = 0.0086), and Jak-STAT signaling (P = 0.0087) as pathways enriched in CAF-HIGH genes (Supplementary Table S2). In addition, genes coding for proteins secreted into the extracellular space were significantly enriched in the CAF-HIGH gene list (corrected P = 0.0026 and 0.0058, respectively) as were membrane proteins (corrected P = 5.8 × 10−4; Supplementary Table S3). Thus, extracellular proteins secreted by CAFs may represent a key component of their phenotype. A heat map for these 15 secretome genes is shown in Fig. 2A. Differential expression of most of these secreted genes in CAFs was validated by qRT-PCR (Fig. 2B).
A secretory and inflammatory gene signature specific to CAFs. A, heat map of soluble ligands identified in the CAF gene signature. NFs, normal fibroblasts. B, mRNA expression analysis of the 15 genes identified in A. Results are shown as CAF expression (n = 5 samples) relative to average NF expression (n = 5 samples) for each gene indicated in the x-axis. *, P < 0.05; **, P < 0.01; and ***, P < 0.001. C, 2-dimensional representation of network 1 obtained using IPA software. Genes in red correspond to those with a fold change >2 in CAFs. Secreted proteins are highlighted by red circles. Lines indicate interaction between 2 genes and arrows show direction of a gene interaction when information is publicly available.
To identify whether the gene signature of CAFs was centered around specific gene networks, we used a database of curated information about gene interactions (see Materials and Methods). This identified 4 major gene networks based on the human orthologues of the CAF-HIGH genes (Supplementary Table S4). These networks were enriched for molecular pathways of inflammatory response and disease, cell death, cellular growth and development, and cell-to-cell signaling and interaction (P < 0.0001; Supplementary Table S5). Interestingly, the inflammatory response pathway was represented in 2 networks, highlighting the biologic relevance of this cellular function in CAFs. The 15 secretome genes were particularly enriched in “network 1,” which contained 10 of these secreted proteins (Fig. 2C). Thus, an inflammatory expression profile involving multiple secreted proteins is characteristic of lung CAFs. Collectively, these results identify secreted proteins that are part of an inflammatory phenotype in lung CAFs that might contribute to tumor development.
A CAF secreted-gene signature is a marker of poor survival in human NSCLCs
To probe the relevance of these findings to human NSCLC, we conducted cross-species gene expression analysis using publicly available human NSCLC gene expression data sets. As NSCLC samples often contain a significant percentage of tumor stroma, we reasoned that expression profile analyses from NSCLC biopsy samples could identify gene expression changes due to overexpression of stromal genes rather than genes expressed in tumor cells themselves. First, gene expression of the human orthologues of the mouse-defined secreted gene signature was assessed in 5 independent microarray data sets from normal human lung and NSCLC samples (36, 39–42). Several of these genes (CCL7, CCL13, CLCF1, EREG, IGFBP5, IL11, LIF, and NGF) were upregulated in NSCLCs compared with normal lung (Supplementary Fig. 3A and Supplementary Table S6). Next, the prognostic potential of the CAF gene signature was queried in a large expression profile data set of human NSCLC (36). In 2 independent data sets, univariate analysis revealed that high levels of the secreted CAF gene signature was associated with poor survival in patients with NSCLC (CAN/DF: log-rank test, P = 0.0002, HR = 3.16 ± 0.26; and UM: log-rank test, P = 0.002, HR = 2.93 ± 0.14; Fig. 3A and B). Survival prediction using this gene set was independent of stage and gender in a multivariate analysis (CAN/DF: log-rank test, P = 0.001, HR = 3.12 ± 0.26; and UM: log-rank test, P = 0.007, HR = 2.91 ± 0.14). These observations show that the secreted gene signature identified in mouse CAFs is highly relevant to human NSCLC and is a strong predictor of patient survival.
Survival prognosis of a mouse CAF-secretory gene signature in human NSCLCs. A and B, Kaplan–Meier plots of patients with NSCLCs from 2 different data sets, Dana Farber Cancer Institute (CAN/DF; Boston, MA), and University of Michigan (UM; Ann Arbor, MI), using high and low levels of a secretory gene signature corresponding to 15 secretome genes identified in mouse CAFs.
Tumor–stroma cross-talk directly alters the phenotype of resident lung fibroblasts to become CAFs
The mechanism by which normal lung fibroblasts are activated to become CAFs is poorly understood. For example, it is unclear whether tumor cells can directly activate normal mouse lung fibroblasts or whether other stromal components are first recruited into the tumor microenvironment leading to an indirect mechanism for normal mouse lung fibroblasts activation. Furthermore, it is unclear whether normal mouse lung fibroblasts are altered to become CAFs or whether CAFs perhaps arise from a rare subpopulation of either resident cells or cells mobilized from bone marrow–derived mesenchymal stem cells (37, 43–46). Therefore, we tested whether tumor cells could induce changes via secreted proteins that would directly lead to upregulation of CAF-secreted genes in normal mouse lung fibroblasts. Conditioned media from mouse tumor cells were used to treat immortalized normal fibroblasts (see Materials and Methods and Supplementary Fig. S4A), and the expression of the secretory gene signature was assessed by qRT-PCR 36 hours later. Normal mouse lung fibroblasts subjected to tumor cell conditioned media displayed increased expression (>1.5-fold) of mRNAs for 11 of 15 CAF genes identified by microarray analysis (Fig. 4A). Similar findings were observed in an independent normal mouse lung fibroblast line (Supplementary Fig. S4B). The upregulation of the secretory gene signature in normal mouse lung fibroblasts after a short exposure to tumor cell conditioned media suggests that CAFs may arise directly from reprogramming of resident normal mouse lung fibroblasts rather than requiring expansion of a rare subpopulation of cells within the stromal compartment. However, expression of these secreted proteins in normal mouse lung fibroblasts returned to basal levels after 72 hours of treatment (data not shown), suggesting that tumor cell conditioned media provide a transient effect that fails to permanently reprogram normal mouse lung fibroblasts toward a CAF phenotype.
Lung tumor cells and CAFs induce a CAF-like phenotype in normal fibroblasts (NFs). A, NFs were grown in the presence of conditioned media from NFs or LKR10 cells for 24 hours, and the expression levels of CAF soluble ligands were assessed by qRT-PCR. Conditioned media (CM) was harvested from cells after 5 days in culture. Results are relative to conditioned media from NFs and representative of 3 experiments. B, mCherry-NFs were cocultured in vitro with NFs or LKR10 cells for 5 (dark gray) or 15 (light gray) days, FACS sorted, and studied by qRT-PCR for the expression levels of CAF-secreted proteins. Experiments were carried out in triplicate and repeated at least twice. Results are relative to coculture with NFs. C, human NFs were cocultured with A549 cells expressing GFP for 24 days. NFs were then sorted out for expression analysis of indicated genes. Results are relative to human NFs cocultured with GFP-positive human NFs and representative of 2 independent experiments. D, NFs were cocultured in vitro with mCherry-NFs or mCherry-CAFs (ratio 1:2) for 5 (dark gray) and 15 (light gray) days and the 15 secretome genes were studied by qRT-PCR. E, immunostaining of xenografts using anti-mCherry and anti-vimentin antibodies. Arrowheads point at endogenous, host fibroblasts. Scale bar, 200 μm. F, averaged tumors from LKR13 cells injected alone, with NFs, or NFs passaged in vivo with G13 cells (induced CAFs or iCAFs; n = 8 mice per group). Endpoint is 20 days. G, cell proliferation assay of NFs and iCAFs measured by MTT analysis. H, same experiment as in F, using human NSCLC cell lines (A549) and one iCAF line (iCAF2). Endpoint is 34 days (A549). Error bars indicate ±SEM. All P values are for a 2-tailed t test. *, P < 0.05; **, P < 0.01; and ***, P < 0.001.
To determine whether direct cell-to-cell contact or continuous stimulation by tumor cells would provide stable upregulation of the CAF-secretory gene signature, normal mouse lung fibroblasts and tumor cells were cocultured in vitro, and the expression profile of the secreted gene signature was studied in these “educated” normal mouse lung fibroblasts. Lung fibroblasts were modified to express a fluorescent reporter (mCherry) to allow isolation by fluorescence-activated cell sorter (FACS) after coculture (Supplementary Fig. S4C). An upregulation of 9 of 15 secreted proteins, similar to what was observed after conditioned media treatment, was observed after 5 days of coculture. A further increase in the expression level of the majority of secreted proteins (13 of 15 genes) was observed after 15 days (Fig. 4B). These results were confirmed using a coculture system incorporating human normal lung fibroblasts and GFP-positive NSCLC cell lines (A549), indicating the relevance of tumor–normal fibroblast interactions to human lung cancer. Exposure of human normal fibroblasts to NSCLCs for 12 days enhanced the expression of 9 genes of the secretory gene signature >1.5-fold (data not shown). Upregulation of 12 of 15 members of the secretory gene signature was further increased after 24 days of coculture in vitro (Fig. 4C). These studies show the existence of functional cross-talk between tumor cells and fibroblasts and suggest that (i) the CAF phenotype is the consequence of direct and chronic stimulation of normal lung fibroblasts by tumor cells and (ii) CAFs can be derived directly from resident lung fibroblasts and do not necessarily arise from a distinct population of mesenchymal cells recruited from a distant site.
While the above results suggest that tumor cells directly contribute to normal fibroblast–to-CAF transition, CAFs may also establish homotypic interactions with resident normal lung fibroblasts and contribute to further recruitment of fibroblasts toward a CAF phenotype. Treatment of normal mouse lung fibroblasts with CAF conditioned media did not upregulate expression of secretory genes (data not shown) indicating that, unlike tumor cell conditioned media, CAF conditioned media is not sufficient to induce gene expression changes in normal mouse lung fibroblasts. However, coculture of normal mouse lung fibroblasts and CAFs led to upregulation of expression of 8 of 15 CAF genes in normal mouse lung fibroblasts (Fig. 4D and Supplementary Fig. S4D). Taken together, these results suggest that tumor cells initiate a transition of normal mouse lung fibroblasts to a proinflammatory CAF state that is further propagated to other normal mouse lung fibroblasts by already activated CAFs.
To analyze the relevance of these findings in vivo and to determine whether normal lung fibroblasts could be stably reprogrammed to a CAF-like phenotype, fluorescently labeled normal mouse lung fibroblasts (mCherry) and tumor cells (GFP) were co-injected to form xenografts (Supplementary Fig. S4E). Tumors from these co-injections yielded similar volumes to those observed with primary lung fibroblasts (Supplementary Fig. S4F) showing that manipulation of lung fibroblasts and tumor cells with reporter constructs did not result in significantly phenotypic changes. Double staining with anti-mCherry and anti-vimentin antibodies showed the presence of exogenous, donor-derived normal mouse lung fibroblasts within expanding tumors (Fig. 4E). Thus, lung fibroblasts expand within xenografts and are an active and sustained component of developing tumors in vivo. The in vitro data above suggested that tumor cells can progressively drive the normal fibroblast–to-CAF transition and that this transition is further promoted by CAFs acting directly on normal mouse lung fibroblasts. To study whether lung fibroblasts could acquire a stable CAF-like phenotype after in vivo passage with tumor cells, normal mouse lung fibroblasts were sorted from tumors and expanded in vitro for approximately 7 to 10 days. Next, tumor cells were injected alone, with parental normal lung fibroblasts, or with fibroblasts isolated from tumor xenografts (here called induced CAFs or iCAFs). Secondary xenografts incorporating iCAFs were larger than those including “uneducated” primary normal mouse lung fibroblasts and the difference closely paralleled that between primary normal mouse lung fibroblasts and CAFs (Fig. 4F). These differences were not due to variations in the proliferation rate between iCAFs and normal mouse lung fibroblasts (Fig. 4G). iCAFs isolated from secondary xenografts (siCAF) also recapitulated the tumor-promoting abilities shown by iCAFs from primary xenografts (Supplementary Fig. S4G). These results indicate that normal mouse lung fibroblasts can be “educated” to display characteristics of CAFs in vivo and that the long-term interaction of these fibroblasts with tumor cells provides a point of no return from a CAF-like phenotype even after in vitro expansion for more than 1 week.
Next, we tested the relevance of the CAF phenotype established in mouse fibroblasts to human lung cancer oncogenesis. The normal mouse lung fibroblasts and iCAFs described above were mixed with the human NSCLC cell lines A549 or H1299. A549 and H1299 cells grew significantly better when mixed with iCAFs (Fig. 4H and Supplementary Fig. S4H). Thus, the acquisition of a CAF-like phenotype in resident lung fibroblasts leads to the promotion of oncogenesis in vivo via the direct stimulation of tumor cell proliferation by secreted proteins. These findings support the relevance of studies of mouse-derived lung CAFs to human lung cancer.
The cytokines clcf1 and il6 contribute to tumor progression in vivo
As normal mouse lung fibroblasts can transition to iCAFs that are phenotypically similar to CAFs isolated directly from tumors, the most functionally relevant secreted genes that drive the CAF phenotype should be common to both iCAFs and CAFs. To identify these genes, the expression profile of the CAF-secretory gene signature was queried in normal mouse lung fibroblasts compared with iCAFs. Six genes (clcf1, cxcl14, ereg, gdnf, il6, and ngf) in the original CAF signature were consistently upregulated by at least 1.5-fold in iCAFs (Fig. 5A and Supplementary Fig. S5A). Thus, this subset of the secreted gene signature may be the most critical for the CAF phenotype. As a further functional test, normal fibroblasts and CAFs were isolated from primary human lung adenocarcinomas (Supplementary Fig. S5B). Two pairs of matched primary human normal mouse lung fibroblasts and CAFs were analyzed for the expression of these 6 genes. Notably, only CLCF1 and IL-6 were upregulated in human CAFs compared with their matched normal mouse lung fibroblasts (Fig. 5B and Supplementary Fig. S5C), suggesting that these 2 genes may play a central role in driving the protumorigenic phenotype of CAFs.
Clcf1- and il6-expressing normal fibroblasts (NF) enhance tumor growth in vivo. A, mRNA expression levels of clcf1, cxcl14, ereg, gdnf, il6, and ngf analyzed in 2 iCAF lines analyzed by qRT-PCR. Results are normalized against NF gene expression levels (solid line). B, expression levels of CLCF1 and IL-6 mRNA in human NFs and CAFs were assessed by qRT-PCR. C, mRNA expression levels analysis in NFs overexpressing clcf1 and il6. D, V5-tag protein detection by Western blotting of cells in C. Cells were treated with 10 μg/mL brefeldin A and 5 μg/mL monensin for 5 hours to inhibit the secretory pathway and facilitate protein detection by Western blotting. E and F, xenograft experiment of LKR10 and LKR13 injected with control lacZ-, clcf1-, or il6-overexpressing NFs. Graphs show averaged tumor volume at indicated time points (n = 10 injections per group). Endpoints are 23 days for LKR10 and 19 days for LKR13. Error bars indicate ±SEM. All P values are for a 2-tailed t test. *, P < 0.05. G, cell proliferation assay of NFs overexpressing lacZ, clcf1, and il6 at day 3 of experiment, as measured by MTT. H and I, mRNA expression levels of tumor cells (LKR10 and LKR13), primary CAFs, control NFs, and NFs overexpressing either clcf1 or il6 assessed by qRT-PCR. Results are relative to NFs_lacZ (black line).
To test the functional role of clcf1 and il6 in lung cancer in vivo, normal mouse lung fibroblasts were transduced with lentiviral vectors containing cDNAs coding for either of these 2 mRNAs (Fig. 5C and D). Fibroblasts overexpressing clcf1, il6, or control (lacZ) were co-injected with tumor cells. Co-injection with lung fibroblasts overexpressing clcf1 and il6 yielded significantly larger tumors (Fig. 5E and F). Similar results were obtained using an independent lung tumor cell line derived from the KrasLSLG12D mouse model (ref. 47; Supplementary Fig. S5D). Importantly, overexpression of clcf1 or il6 did not alter the proliferation rate of the corresponding fibroblasts (Fig. 5G). Thus, clcf1 and il6 secreted by CAFs play a direct role in stimulating tumor growth in vivo. This effect is likely due to direct stimulation of tumor cells by fibroblasts, as both clcf1 and il6 overexpressors increased the relative cell number of cancer cells when compared with control fibroblasts when cocultured in vitro (Supplementary Fig. S5E and S5F). No expression of either clcf1 or il6 was detected in tumor cells (Fig. 5H and I), suggesting that the primary source of both cytokines is derived from CAFs which establish a paracrine signaling axis to promote tumor cell growth.
Suppression of clcf1–cntfr paracrine signaling axis decreases lung cancer in vivo
Clcf1 and il6 act through a receptor complex containing specific α-subunit receptors, cntfr and il6r, which recognize and bind clcf1 and il6, respectively. These distinct α-subunits act together with the common β-subunit receptor gp130 to activate JAK-STAT signaling (48). Abrogating the ligand–receptor interaction via decreased expression of cntfr and il6r in tumor cells would thus be expected to halt tumor growth in vivo. To test this hypothesis, 2 independent short hairpin RNAs (shRNA) to cntfr and il6r were validated and used to infect LKR10 cells (Fig. 6A and Supplementary Fig. S6A). We co-injected CAFs along with cntfr-and il6r-deficient LKR10 cells or controls and assessed tumor volume. Tumor cells carrying single shRNAs directed against cntfr were significantly smaller than controls (Fig. 6B). A similar trend was observed with 2 independent il6r shRNAs; however the differences with regard to control cells were not statistically significant for both hairpins (Supplementary Fig. S6B). Importantly, loss of cntfr or il6r did not alter tumor cell growth in vitro, suggesting that the in vivo effect is likely due to an important role of this signaling pathway in mediating paracrine effects in vivo but not intrinsic tumor proliferation (Fig. 6C and Supplementary Fig. S6C). Taken together, the results of the overexpression and knockdown studies indicate that while both clcf1 and il6 stimulate tumor growth, il6–il6r signaling may be dispensable whereas loss of the clcf1–cntfr axis causes a consistent decrease in tumor formation, suggesting that the 2 signaling pathways are distinct and yet together constitute an important paracrine signaling network between CAFs and tumor cells in NSCLCs.
Decreased response to CAF-derived clcf1 in cntfr-depleted tumor cells. A, qRT-PCR analysis of LKR10 cells transduced with 2 independent shRNAs to cntfr or a control shRNA. B, tumor volume analysis of xenografts including CAFs and cntfr-depleted LKR10 cells. Error bars indicate ±SEM. All P values are for a 2-tailed t test. **, P < 0.01. Endpoint is 26 days. C, cell viability analysis of LKR10 cells expressing cntfr shRNAs assessed by MTT.
Discussion
Accumulating evidence points to a critical role for CAFs in tumor progression (4–6). As increased stromal proportion is a marker of poor survival in patients with NSCLCs (16, 17, 49), CAFs are likely to play an important role in lung cancer progression. We used a cross-species approach to identify a secretory gene signature in CAFs that promotes lung cancer growth and to query individual components of this signature for relevance to the phenotype of CAFs and to human NSCLCs. These studies link secretion of CLCF1 and IL-6 from CAFs to tumor progression in vivo. CLCF1 and IL-6 are members of the IL-6 family (50). CLCF1 (cardiotrophin-like cytokine factor 1, also known as NNT-1 or BSF-3) is a cytokine that has been implicated in the pathogenesis of multiple myeloma (51) but has not been reported to play a role in NSCLCs. While autocrine secretion of IL-6 has been proposed to occur in lung cancer cell lines carrying EGF receptor (EGFR) mutations (52), our results suggest that IL-6 is produced primarily by CAFs and not tumor cells. Indeed, the functional analysis presented here suggests that CLCF1 and IL-6 secretion by CAFs play a significant role in the contribution of this component of the tumor microenvironment to tumor growth.
CLCF1 and IL6 are cytokines that signal through multimeric complexes that include the gp130 receptor/β-subunit and ultimately trigger a signaling cascade that is mediated by the Jak-STAT and mitogen-activated protein kinase (MAPK) pathways (50). The specificity of the cellular responses elicited by these cytokines is governed by the expression patterns of different α-subunit receptors, CNTFR for CLCF1, and IL-6R for IL-6, which can be activated by each cytokine independently (50). Interestingly, overexpression of another IL-6 family member upregulated in mouse CAFs, IL11, which signals through an alternative α-subunit receptor (IL11RA), does not induce tumor growth in vivo in our model system (data not shown), suggesting that distinct receptor complexes may be activated in specific tumor types. Our loss-of-function experiments show a more prominent role of the CLCF1–CNTFR axis in our model and suggest that CLCF1 may be particularly relevant to NSCLC progression. Thus, these results provide rationale for the development of new therapies targeting CAFs and suggest that inhibition of the CLCF1–CNTFR axis may be a novel strategy.
In vitro and in vivo experiments in models of breast and skin carcinoma have reported that tumor cells can induce a CAF-like phenotype in normal fibroblasts, suggesting that direct tumor cell–fibroblast cross-talk plays an important role in tumor development (11, 45). In this study, we show that transient upregulation of some CAF-specific genes in normal mouse lung fibroblasts can be induced by tumor-derived secreted cues. However, the ability of normal mouse lung fibroblasts to fully upregulate a secretory CAF gene signature is only achieved through long-term direct exposure to tumor cells. Our results suggest that while normal mouse lung fibroblasts are responsive to tumor-derived signals, the generation of normal mouse lung fibroblasts populations stably expressing a CAF-like gene signature (and thus, recapitulating the CAF tumor-promoting phenotype) likely involves genetic and/or epigenetic changes resulting from chronic tumor cell-to–normal fibroblast interactions (53, 54). The fact that normal mouse lung fibroblasts passaged through a xenograft (iCAFs) are phenotypically indistinguishable from primary CAFs after 3 to 4 weeks of in vivo growth (despite intermediate expansion in vitro for >1 week) supports the notion that long-term interaction between tumor cells and normal mouse lung fibroblasts leads to the acquisition of a stable CAF phenotype. In addition, in vitro experiments described here also show a role for human NSCLCs in the reprogramming of human normal fibroblasts. Thus, tumor cells are critical in regulating the reprogramming of lung fibroblasts to CAFs in both mouse and human lung cancer. Further work will be required to elucidate the molecular mechanisms that induce this transition.
Cross-species gene expression analysis has been used successfully to identify mouse gene signatures relevant to human NSCLCs (23). By analyzing mouse stromal cells, we identified a gene signature predictive of survival in human NSCLCs, supporting the relevance of studies of stromal compartments in mice to the pathogenesis of human cancer. Previous work identified an 11-gene signature of human CAFs in NSCLC samples with prognostic relevance (55). Surprisingly, we found no significant overlap between our mouse CAF-HIGH gene set and the upregulated CAF genes uncovered in the former study (data not shown). This discrepancy is likely due to the heterogeneous genetic background of the different human samples and highlights the ongoing value of cross-species comparison as a tool for functional genomics studies.
In summary, the data presented in this work describe a novel CAF gene signature with significant relevance in human NSCLCs. We show that CAF-derived CLCF1 and IL-6 exert specific paracrine effects to promote tumor growth in vivo and may be particularly important components of the CAF secretome. In addition, we establish that direct interaction between tumor cells and lung fibroblasts is sufficient to alter the phenotype of the later to become CAFs. These observations highlight a prominent role for CAFs in NSCLCs and provide a framework for further mechanistic studies to identify relevant therapeutic targets. Finally, these studies underscore the need to further elucidate the mechanisms that lead to the conversion of lung fibroblasts to CAFs as this may provide a way to inhibit the “education” of normal mouse lung fibroblasts to a protumorigenic CAF state.
Disclosure of Potential Conflicts of Interest
A.J. Butte is a consultant/advisory board member for NuMedii, Personalis, Ansh Labs, Roche, Johnson and Johnson, and Samsung Advanced Institute of Technology. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: S. Vicent, E.A. Sweet-Cordero
Development of methodology: S. Vicent, E.A. Sweet-Cordero
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S. Vicent, L.C. Sayles, Y. Zheng, A.K. Gillespie, N. Clarke, Y. Xu, J. Shrager, C.D. Hoang, E.A. Sweet-Cordero
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Vicent, L.C. Sayles, D. Vaka, P. Khatri, O. Gevaert, R. Chen, S. Plevritis, A.J. Butte, E.A. Sweet-Cordero
Writing, review, and/or revision of the manuscript: S. Vicent, J. Shrager, C.D. Hoang, A.J. Butte, E.A. Sweet-Cordero
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S. Vicent, L.C. Sayles, J. Shrager, E.A. Sweet-Cordero
Study supervision: S. Vicent, J. Shrager, A.J. Butte, E.A. Sweet-Cordero
Grant Support
This work was supported by the Alfonso Martín Escudero Foundation and the Tobacco-Related Disease Research Program (TRDRP) of California (UCAEF SPO#44443; to S. Vicent), the National Cancer Institute (1R01CA138256-01; to P. Khatri, A.J. Butte, and E.A. Sweet-Cordero), the Center for Cancer Systems Biology (CCSB) at Stanford [funded by a grant from the National Cancer Institute (U54 CA149145) to O. Gevaert and S. Plevritis], the Fund for Scientific Research Flanders (FWO-Vlaanderen to O. Gevaert), and an American Cancer Society Research Scholar Award (RSG-08-132-01-CDD; to E.A. Sweet-Cordero).
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.
Acknowledgments
The authors thank Julien Sage, Laura Attardi, Steve Artandi, Monte Winslow, and all the members of the Sweet-Cordero laboratory for helpful discussions. They also thank Oihana Murillo for technical assistance with flow cytometric analyses.
Footnotes
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).
- Received March 26, 2012.
- Revision received August 15, 2012.
- Accepted August 17, 2012.
- ©2012 American Association for Cancer Research.