This article has Open Peer Review reports available.
BMP2 response pattern in human lung fibroblasts predicts outcome in lung adenocarcinomas
© Rajski et al.; licensee BioMed Central. 2015
Received: 5 October 2014
Accepted: 17 April 2015
Published: 29 April 2015
Bone morphogenetic proteins play important roles in development, morphogenesis and cancer. With this study we aimed to characterize the response of lung stromal fibroblasts to BMPs and their antagonists on a genome wide level and investigate its potential role in human lung adenocarcinomas.
We used an ex vivo culture model and measured gene expression changes in human lung fibroblasts after stimulation with BMPs and their antagonists using HEEBO microarrays. The in vitro data were correlated with in vivo observations in published expression datasets of human lung adenocarcinomas.
We have systematically analyzed the response to BMP2, BMP4, BMP7 and their antagonists, Gremlin and Noggin, to define common and specific gene expression patterns. A BMP2 induced gene expression signature was defined, which is specific for stromal fibroblasts. Gene expression profiles from lung adenocarcinoma biopsies were analyzed to determine the prognostic significance of the “Fibroblast specific BMP2 induced gene list”. This gene list successfully segregated patients with different prognostic outcome in 3 datasets. In a small dataset (Garber et al.) there was a strong trend for a worse prognosis of patients with adenocarcinomas of all stages over-expressing the “Fibroblast specific BMP2 induced gene list”. In two larger datasets with stage I adenocarcinomas we observed a significantly worse disease-free (p = 0.002, Lee et al. and p = 0.002, Bhattacharjee et al.) and overall survival (p = 0.0002).
The effects of BMPs and their antagonists are heterogeneous in different cell types. The gene expression pattern induced by BMP2 in primary lung fibroblasts may predict outcomes of patients with lung adenocarcinomas.
The bone morphogenetic proteins (BMPs) are secreted signaling molecules that play important roles in numerous aspects of metazoan development, including stem cell regulation, fate determination, body axis specification, cell polarity, proliferation and differentiation . Thus, alterations might play a role in cancer.
It has long been established that epithelial-mesenchymal interactions are essential for development and organogenesis  such as the skin, teeth, gut, and lungs. More recently, it has been recognized that cancer is a tissue-based disease of stromal cells co-evolving with carcinoma cells. Increasing evidence suggests that tumorigenesis is dependent upon contextual signals received from the stroma  and epithelial-mesenchymal interactions are playing an important role in the initiation and progression of tumorigenesis. As the tumor develops, the surrounding stroma provides tumor cells with growth factors and favorable matrix components that foster proliferation, migration and colonization of distant organs. In fact, the reciprocal interaction of stromal cells and carcinoma cells drives carcinogenesis. The most abundant mesenchymal cells in the stroma of most tissues and organs are fibroblasts. These cells secrete extracellular matrix components and signaling molecules that contribute to the establishment of customized microenvironments for epithelial cells and provide a specialized niche for tissue stem cells . Dissecting the molecular mechanisms that underlie cell-to-cell signaling and crosstalk between epithelial cells, stem cells, fibroblasts and other stromal cells in tissue microenvironments is essential to understand normal development and carcinogenesis . Here, we focus on one molecular element that may play a role in epithelial-mesenchymal interactions in development and cancer – the transcriptional response of fibroblasts to BMP-induced signals. Given the multifaceted roles of BMP signaling in these processes , its transcriptional effects on fibroblastic stromal cells are of great interest. However, in lung cancer, the genome-wide effects of the involved signaling molecules on stromal fibroblasts have not yet been systematically described. Therefore, to build a systematic foundation of knowledge, we focused on the effects of BMPs and their antagonists, Gremlin and Noggin, on transcription in normal human lung fibroblasts. Importantly, we also demonstrate that the BMP2 gene signature identified in the lung fibroblasts in vitro has significant prognostic relevance in human lung adenocarcinomas and may thus serve as novel prognostic cancer markers and therapeutic targets in future.
Primary human fibroblasts (CCL-171) and the human breast cancer cell lines MDA-MB-231 and T47D were obtained from the American Type Culture Collection (ATTC, Atlanta, Georgia, USA) on 2005 November 15. After resuscitation the cells were propagated in Dulbecco’s modified Eagle’s medium (D-MEM, Invitrogen, Carlsbad, USA) supplemented with 10% heat-inactivated FBS (Invitrogen, Carlsbad, CA, USA), 4.5 g/l glucose, 4 mM L-glutamine, 100 U/ml penicillin and 100 μg/ml streptomycin (Gibco, Carlsbad, CA, USA). The cells were maintained by regular passages when confluence was reached and used for the experiments within 3–4 months. The study was approved by the Ethikkommission beider Basel, Switzerland (approval No. 271/05).
BMP stimulation experiments
For the experiment, 30’000 cells/cm2 were seeded in 3 ml of 5% FBS D-MEM and incubated for 6 h to allow attachment. The cells were washed extensively with phosphate-buffered saline and starved for 48 h in fresh low-serum D-MEM supplemented with 0.2% FBS. The cells were starved to reduce the effects of any stimulation from the regular cell culture medium. The medium was subsequently replaced with fresh low-serum D-MEM with or without 200 ng/ml BMP2 (human recombinant in Escherichia coli; Sigma Aldrich, St. Louis, MO, USA), 24 ng/ml BMP4, 200 ng/ml BMP7, 240 ng/ml Noggin or 1 μg/ml Gremlin. The cells were stimulated for 24 h, and the RNA was harvested to test the effects of BMP and their antagonists on mRNA expression patterns.
RNA extraction and amplification
After aspirating the culture medium, the cell monolayer was washed once with phosphate-buffered saline. The cells were lysed in a buffer containing guanidine isothiocyanate (RLT buffer, QIAGEN, CA, USA). The total RNA was extracted with an RNeasy kit (QIAGEN, Valencia, CA, USA) according to the manufacturer’s instructions. The RNA concentration was measured with a NanoDrop system spectrophotometer (ND-1000 Spectrophotometer Technologies, Wilmington, NC, USA). The integrity of extracted RNA was assessed by electrophoresis in a 1% agarose gel in MOPS buffer. For mRNA amplification, an Amino Allyl MessageAmp™ II aRNA Amplification Kit was used (Ambion, Austin, TX, USA). The amplification of mRNA from 500 ng total RNA, purification of the cDNA, in vitro transcription and purification of aRNA were performed according to the manufacturer’s instructions. The integrity and quantity of the amplified RNA were verified as described above.
Gene expression analysis with HEEBO arrays
Global gene expression was assessed using the Human Exonic Evidence Based Oligonucleotide (HEEBO) microarrays produced at the Stanford Functional Genomic Facility (Stanford, CA, USA). The HEEBO microarrays consisted of 44,544 70-mer probes, which included the following: (a) constitutive exonic probes (30,718); (b) alternatively spliced/skipped exonic probes (8,441); (c) non-coding RNA probes (196); (d) BCR/TCR Genic/regional probes (372); (e) other probes (843); and (f) controls.
The complete details regarding the clones on the arrays can be found at the Stanford functional genomics facility website (https://microarray.org/sfgf/). For the microarray experiments, 8 μg of amplified RNA (aRNA) was mixed with the doping controls. The samples were vacuum dried, resolved in coupling buffer and labeled with Cy5 dye. The labeled samples were pooled with equal amounts of reverse color Cy3-labeled amplified reference RNA from Stratagene (Stratagene, Santa Rosa, CA, USA). The labeled aRNA was purified using an AminoAllyl MessageAmp™ II aRNA Amplification Kit (Ambion) according to the manufacturer’s instructions and fragmented using fragmentation reagents (Ambion). The fragmented probe was added to a hybridization buffer containing Cot/PolyA/tRNA (0.05 μg/uL each), 0.3% SDS and 3.3 × SSC supplemented with HEPES buffer. Following a denaturing step at 100° C, the probe was placed on the microarray for competitive hybridization. After 18 h, the slides with hybridized probe were sequentially washed, immediately dried in an ozone-free environment and scanned using an Axon Scanner 4100A (Axon Instruments, CA, USA). The data analysis and clustering Microarray fluorescent image analysis were performed using the software Genepix Pro version 5.0 22.214.171.124 (Axon Instruments, Union City, CA, USA). Any spots with obvious array artifacts or poor technical quality were manually removed from any further analysis. The raw data files were stored in the Stanford Microarray Database-Princeton University and are freely available (http://smd.princeton.edu/cgi-bin/publication/viewPublication.pl?pub_no=1095)  as well as in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE 66142 (http://www.ncbi.nlm.nih.gov/geo/).
The data were expressed as the log2 ratio of the fluorescence intensity of the sample and fluorescence intensity of the reference for each element on the array. A sequential data filtering procedure was applied to include only measurements fulfilling our quality requirements (i.e. data with a Cy3 channel or Cy5 channel mean intensity over median background intensity greater than 1.5 and a Pearson regression correlation of the pixels within the spot of greater than 0.6). The genes that did not meet these criteria for at least 80% of the measurements across the experimental samples were excluded from further analysis. The data were evaluated by unsupervised hierarchical clustering with the Cluster® software using Pearson correlation (non-centered metric) and average linkage and displayed using TreeView® software .
SAM is a statistical approach to identify genes with expression patterns that are significantly associated with specific characteristics of the sample sets . The Excel-SAM-package version 2.1. was obtained from the webpage of the Tibshirani Lab: http://statweb.stanford.edu/~tibs/SAM/. SAM analysis was applied to the BMP dataset to examine differences between stimulated and unstimulated cells and between cells exposed to different stimuli. A two-way, unpaired test or a multiclass test was performed to compare the two groups of interest. A 10-nearest neighbor imputation engine was applied to estimate missing data , and 100 permutations were performed to compute the expected values and calibrate false positives. The average amount of imputed values was 5.17 +/− 4.56 percent.
Correlation of gene expression changes upon exposure to different stimuli
To test for similarities between the expression changes of single genes upon stimulation with BMP2, BMP4, BMP7, Noggin and Gremlin we calculated the Pearson correlation coefficient using the R software (R Development Core Team. A language and environment for statistical computing. Available: http://www.R-project.org) and displayed the results as scatter plots.
Human cancer datasets
The dataset published by Garber and colleagues  contains global gene expression profiles for 67 human lung cancers derived from 56 patients. 41 of them were adenocarcinomas with survival data for 24 patients.
The dataset published by Lee and colleagues  contains global gene expression profiles of 138 human lung cancers, of which 63 were adenocarcinomas associated with survival data. The dataset published by Bhattacharjee  contains mRNA expression levels of 12,600 transcript sequences from 186 lung tumor samples, including 139 adenocarcinomas. Of these, 125 samples were associated with clinical data (some patients were in multiple runs). The Bhattacharjee dataset was obtained from the Broad Institute website, and the Garber dataset was acquired from the SMD webpage . A list of 156 unique clones comprising 67 genes building the “Fibroblast specific BMP2 induced gene list” was extracted from the three datasets. To avoid possible overweighting of clones from Unigene clusters that matched more than one probe on the array, the expression values derived from probes matched to the same Unigene cluster were averaged. Only genes that had greater than 80% data values and tumor samples from patients having complete clinical data were used. The resulting dataset was subjected to average linkage unsupervised clustering, and the results were displayed using TreeView software . All statistical tests were performed using the R statistical software (version 2.10.1). Survival curves were obtained using the Kaplan-Meier estimator and univariate Cox proportional hazards regression models were fitted (R package “survival”) R Development Core Team. A language and environment for statistical computing. Available: http://www.R-project.org. Accessed 2013 December 2. Hazard ratios and cox p-values are given to all the curves.
The patients in the Lee  and the Bhattacharjee datasets  were stratified based on a continuous score derived from the “Fibroblast specific BMP2 induced gene list”, as described previously. Briefly, the average expression level of the genes in the “Fibroblast specific BMP2 induced gene list” was calculated for each patient attributing a score. The patients were then divided into two groups, separating them by the median value of the continuous scores. Only patients with stage I cancers were included in the survival analysis. Overall survival was based on the number of deaths from any cause, and the patients were censored at the final follow up. Disease-specific survival analysis was based on the number of deaths caused by the disease, with patients censored at the last follow-up. In this analysis, patients who died from other causes were considered alive and censored. The survival statistics were calculated as described above.
GO::TermFinder uses a list of genes as input and determines whether those genes have gene ontology (GO) terms overrepresented in their combined set of annotations compared with what would be expected by chance for a randomly selected group of genes from the background population of all of the genes . In our analysis, to calculate the frequency of particular annotations we used the gene lists from specific clusters and the full parental gene lists from the same heatmap as a background file. For a SAM-derived signature, we used an input file for SAM analysis as the gene list. For p-value correction we used the Bonferroni correction method, as described by Boyle et al. .
BMP2 stimulation of lung stromal fibroblasts induces regulatory genes of the BMP signaling pathway
A notable feature of the transcriptional response was the induction of genes with direct roles in the regulation of BMP signaling either as inhibitors (e.g., BAMBI) or positive regulators (e.g., BMPR1A). Co-expression of both negative and positive regulators of BMP signaling provides a network of opposing control mechanisms that are likely to contribute to the robust and precise regulation of BMP target gene expression. Furthermore, we detected an up-regulation of BMP target genes such as SMAD7, ID1, ID2 and ID3. Other genes such as BMP4 and CXCL12 were systematically down-regulated. The down-regulation of CXCL12 indicates a tie to stem cell regulation . The “CCL-171-derived BMP2” signature as marked in Figure 1 (Additional file 2: Table S2), also contains genes that are known to be involved in developmental processes, including GATA6, DKK1, ID2, BAMBI, SERPINE1, and PTGIS.
The identification of transcription factors and other regulatory genes involved in the transcriptional response (e.g. HEY1, PHTF2, GATA6, and LMCD1) suggests that a more extensive reprogramming of transcription occurs downstream of the primary target genes; the mechanism of target gene activation is likely to be far more complex than direct activation through SMADs, involving a cascade of downstream transcription factors and regulatory proteins that play a role in further regulation of secondary BMP2 target genes. Among the transcription factors regulated by BMP signal transduction, it was interesting to find HEY1 and GATA6, which play important roles in pattern formation, morphogenesis and body-axis specification. The increased expression of TPM1, TPM2 and TNC  indicates a fundamental reprogramming of normal lung fibroblasts such that they exhibited a gene expression pattern that is typically found in carcinoma-associated fibroblasts.
For an unbiased assessment of which features are shared by members of the “CCL-171-derived BMP2” signature and to verify the significance of enrichment of a specific gene ontology, we used the GO::TermFinder tool . The analysis revealed that the “CCL-171-derived BMP2” signature as shown in Figure 1 (Additional file 2: Table S2) is significantly enriched for genes involved in biological processes such as the cellular response to chemical stimuli and the BMP signaling pathway with Bonferoni corrected p-values of 0.00017 and 0.006, respectively (Additional file 3: Table S3).
Furthermore to identify genes with a significant change in expression levels we performed a two class unpaired SAM analysis using a false discovery rate <1%. 76 significantly induced and 151 significantly repressed genes after BMP2 stimulation are shown in Additional file 4: Table S4. According to our expectations, all the significantly up-regulated genes are contained within the list of Additional file 1: Table S1. In addition 10 genes not comprised within Additional file 1: Table S1 were found to be significantly down-regulated with SAM analysis, namely: GNG11, ABHD5, PHF17, GAS1, ARSI, JUN, KLHL29, CHEK2, ST6GAL1, PDGFRA, H1F0, BTG1 (fold changes of these genes: 2.15-2.58 fold decrease).
A global picture of genes that are differentially expressed in response to BMP stimulation or inhibition in lung stromal fibroblasts
A comparison of the effect of BMP2 stimulation in lung stromal fibroblasts and breast cancer cells
To determine whether the response to BMP2 is a general effect or specific for mesenchymal cells such as the lung fibroblasts (CCL-171), we compared the gene expression profile of BMP2-stimulated CCL-171 cells with the expression profiles of the breast cancer cells MDA-MB-231 and T47D treated with BMP2. First, to systematically define a gene signature reflecting the common response to BMP2 shared by CCL-171, MDA-MB-231 and T47D cells, we used a two class SAM with block permutation to identify genes with a significant change using a false discovery rate of less than 0.8%. A common set of 80 genes that were significantly induced, and 200 genes that were significantly repressed in all three cell types were identified (Additional file 6: Table S6).
BMP4 stimulation of lung stromal fibroblasts induces a time-dependent gene expression signature
Prognostic significance of the “Fibroblast specific BMP2 induced gene list” in human lung adenocarcinomas
To further support this finding we analyzed the stage I lung adenocarcinomas published by Lee et al. . From 67 genes building “Fibroblast specific BMP2 induced gene list” 48 were present in Lee dataset. Patients stratification based on unsupervised hierarchical clustering divided patients into two groups. Patients carrying tumors with high expression levels of the “Fibroblast specific BMP2 induced gene list” showed a significantly higher risk of death than the patients with a low expression level (76% versus 38% after 24 months, hazard ratio (HR): 0.49, cox-p =0.002) (Figure 5B).
Gene expression signatures have been described in a number of studies as surrogate phenotypes to explore the relevance of in vitro findings in vivo in human cancer [20-22]. In this study, we have taken this approach to explore the effects of BMP stimulation of fibroblastic stroma in vitro on lung cancer. It has long been established that epithelial-mesenchymal interactions are essential for lung development and organogenesis . Given the multifaceted roles of BMP signaling in these processes [24,25], investigating its transcriptional effects on fibroblastic stromal cells might improve our understanding of lung carcinogenesis. In this study we provide a comprehensive genomic analysis of gene expression changes in lung fibroblasts in response to stimulation with BMP2, BMP4, BMP7 and their antagonists Gremlin and Noggin. Our results reveal alterations in the expression of a diverse spectrum of genes, reflecting changes in various components of key signaling pathways involved in tumorigenesis.
Stimulation of primary human lung fibroblasts with BMP2 significantly changes their gene expression pattern through signaling via the BMP pathway, as shown by changes in the expression of genes in the BMP pathway and BMP target genes, such as BMPR2A, BAMBI, SMAD7ID1, ID2, and ID3. Of specific interest is the induction of TPM1, TPM2, TNC and POSTN, suggesting a fundamental reprogramming of lung fibroblasts such that the cells exhibit a gene expression profile that is typical for carcinoma-associated fibroblasts. Of specific interest is also the expression of genes, which might support cancer stem cells by forming a stem cell niche. POSTN e.g. has recently been shown to be an essential component of the metastatic stem cell niche . Preparation of the metastasis location is an essential step in the cascade of effects leading to the development of fatal metastatic disease. It would be of great interest to further investigate and clarify the role of BMPs in the lung cancer stem cell niche. Furthermore as a result of the GO-term analysis, the mRNA expression levels of GATA6, ID2, AHR, IGFBP7, BAMBI, SERPINE1, SMAD6, TRNP1, and CYR61 increased in response to BMP2 stimulation, suggest that BMP2 stimulation of fibroblasts plays a role in the regulation of cell proliferation. GATA6, DKK1, ID2, BAMBI, SERPINE1, CYR61 and PTGIS indicate a positive regulation of a developmental process. The local regulation of diverse cell signaling pathways by BMP in fibroblasts could therefore have multifaceted consequences for tissue microenvironments in vivo, including the balance between proliferation and development. The interconnectivity of these regulatory systems may contribute to the robustness of stem cell niches and the precise spatial and temporal control of differentiation.
To clarify the in vivo relevance of our in vitro finding of an activation of the fibroblastic stroma, we correlated the expression of the “fibroblast specific BMP2 induced gene list” with the outcome of lung adenocarcinoma biopsies. Up-regulation of the “fibroblast specific BMP2 induced gene list” is associated with a significantly worse disease-free and overall survival as shown in two datasets of operated stage I lung adenocarcinoma patients, and one set of human lung adenocarcinoma biopsies.
Since segregation of patient samples using unsupervised hierarchical clustering might be unstable and dependent on missing values, we confirmed the finding using a more robust approach based on a score derived from the mean expression values (centroid) of all the genes in the signature. This suggests that the “fibroblast specific BMP2 induced gene list” is prognostic in lung adenocarcinomas. This hints to the in vivo biological relevance of our in vitro finding of an up-regulation of this gene list by BMP2 stimulation in lung fibroblasts. Of note this does not necessarily mean that any gene subsets or single genes of the “fibroblast specific BMP2 induced gene list” have the same prognostic power, neither it excludes existence of other, parallel prognostic signatures.
In line with our findings, BMP2 has been shown to be over-expressed in lung cancer . BMP2 serum levels have been shown to be higher in lung cancer patients than in normal healthy controls . There was a positive correlation of tumor burden and tumor stage. Not surprisingly BMP2 serum levels were correlated with patient outcome . In our analysis the effect of BMP2 stimulation expressed by the “fibroblast specific BMP2 induced gene list” is correlated with patient outcome even in the same tumor stage, namely stage I tumors. This hints to a prognostic relevance of the BMP2 effect in lung cancer, which is independent from stage and tumor burden. Further in vitro evidence suggests an essential role of BMP2 in lung cancer. Knockdown experiments of BMP2 in A549 lung adenocarcinoma cells had shown decreased growth and invasion potential by changes in the tumor cells gene expression pattern revealing therapeutic targets and strategies .
Our analysis of the changes in the expression pattern of the stromal fibroblasts upon BMP2 stimulation adds additional insights in tumorigenesis and potential therapeutic targets. Stromal-targeted therapy in lung cancer has not yet become a clinical reality due to several challenges. Identification of tumorigenic mechanisms and factors, as e.g. BMP signaling in stromal cells, is a prerequisite for the design of potential novel therapeutic strategies.
Bone morphogenetic proteins play important roles in development, morphogenesis and cancer. We showed that the effects of BMPs and their antagonists, gremlin and noggin, are heterogeneous in different cell types. The gene expression pattern induced by BMP2 in lung fibroblasts significantly correlated with the prognosis of patients with stage I lung adenocarcinomas.
We would like to thank Michael Fero, Elena Seraia and the staff of the Stanford Functional Genomics Facility for supplying us with the HEEBO microarrays that were used for this study. We thank Kathy Ball, Janos Demeter, John Matese and the staff from SMD for support with the microarray data storage and analysis tools. This work was supported in part by the Swiss National Science Foundation grant 320000-112794/1(http://www.snf.ch) and a grant from Oncosuisse (http://www.oncosuisse.ch).
- Bragdon B, Moseychuk O, Saldanha S, King D, Julian J, Nohe A. Bone morphogenetic proteins: a critical review. Cell Signal. 2011;23(4):609–20.View ArticlePubMedGoogle Scholar
- Langenfeld EM, Calvano SE, Abou-Nukta F, Lowry SF, Amenta P, Langenfeld J. The mature bone morphogenetic protein-2 is aberrantly expressed in non-small cell lung carcinomas and stimulates tumor growth of A549 cells. Carcinogenesis. 2003;24(9):1445–54.View ArticlePubMedGoogle Scholar
- Langenfeld E, Kong Y, Langenfeld J. Bone morphogenetic protein 2 stimulation of tumor growth involves the activation of Smad-1/5. Oncogene. 2005;25(5):685–92.View ArticleGoogle Scholar
- Sharpe P, Ferguson M. Mesenchymal influences on epithelial differentiation in developing systems. J Cell Sci.1988;1988(Supplement 10):195–230.Google Scholar
- Sonnenschein C, Soto AM. An updated review of environmental estrogen and androgen mimics and antagonists. J Steroid Biochem Mol Biol. 1998;65(1):143–50.View ArticlePubMedGoogle Scholar
- Zhang J, Niu C, Ye L, Huang H, He X, Tong W-G, et al. Identification of the haematopoietic stem cell niche and control of the niche size. Nature. 2003;425(6960):836–41.View ArticlePubMedGoogle Scholar
- Pietras K, Ã–stman A. Hallmarks of cancer: interactions with the tumor stroma. Exp Cell Res. 2010;316(8):1324–31.View ArticlePubMedGoogle Scholar
- Singh A, Morris RJ. The Yin and Yang of bone morphogenetic proteins in cancer. Cytokine Growth Factor Rev. 2010;21(4):299–313.View ArticlePubMedPubMed CentralGoogle Scholar
- Hubble J, Demeter J, Jin H, Mao M, Nitzberg M, Reddy T, et al. Implementation of genepattern within the stanford microarray database. Nucleic Acids Res. 2009;37 suppl 1:D898–901.View ArticlePubMedGoogle Scholar
- 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(25):14863–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, et al. Missing value estimation methods for DNA microarrays. Bioinformatics. 2001;17(6):520–5.View ArticlePubMedGoogle Scholar
- Garber ME, Troyanskaya OG, Schluens K, Petersen S, Thaesler Z, Pacyna-Gengelbach M, et al. Diversity of gene expression in adenocarcinoma of the lung. Proc Natl Acad Sci U S A. 2001;98(24):13784–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee E-S, Son D-S, Kim S-H, Lee J, Jo J, Han J, et al. Prediction of recurrence-free survival in postoperative nonâ€“small cell lung cancer patients by using an integrated model of clinical information and gene expression. Clin Cancer Res. 2008;14(22):7397–404.View ArticlePubMedGoogle Scholar
- Bhattacharjee A, Richards WG, Staunton J, Li C, Monti S, Vasa P, et al. Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc Natl Acad Sci. 2001;98(24):13790–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Calvi L, Adams G, Weibrecht K, Weber J, Olson D, Knight M, et al. Osteoblastic cells regulate the haematopoietic stem cell niche. Nature. 2003;425(6960):841–6.View ArticlePubMedGoogle Scholar
- Chiquet-Ehrismann R, Mackie EJ, Pearson CA, Sakakura T. Tenascin: an extracellular matrix protein involved in tissue interactions during fetal development and oncogenesis. Cell. 1986;47(1):131–9.View ArticlePubMedGoogle Scholar
- Malanchi I. Santamaria-MartÃnez A, Susanto E, Peng H, Lehr H-A, Delaloye J-F, Huelsken J: Interactions between cancer stem cells and their niche govern metastatic colonization. Nature. 2011;481(7379):85–9.View ArticlePubMedGoogle Scholar
- Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, et al. Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004;2(2), E7.View ArticlePubMedPubMed CentralGoogle Scholar
- Buess M, Nuyten DS, Hastie T, Nielsen T, Pesich R, Brown PO. Characterization of heterotypic interaction effects in vitro to deconvolute global gene expression profiles in cancer. Genome Biol. 2007;8(9):R191.View ArticlePubMedPubMed CentralGoogle Scholar
- Nuyten DS, Kreike B, Hart AA, Chi JT, Sneddon JB, Wessels LF, et al. Predicting a local recurrence after breast-conserving therapy by gene expression profiling. Breast Cancer Res. 2006;8(5):R62.View ArticlePubMedPubMed CentralGoogle Scholar
- Masters JR. Epithelial-mesenchymal interaction during lung development: the effect of mesenchymal mass. Dev Biol. 1976;51(1):98–108.View ArticlePubMedGoogle Scholar
- Sr B. Henderson R, Winnier G, Oikawa T, Hogan B: Evidence from normal expression and targeted misexpression that bone morphogenetic protein (Bmp-4) plays a role in mouse embryonic lung morphogenesis. Development. 1996;122(6):1693–702.Google Scholar
- Weaver M, Yingling JM, Dunn NR, Bellusci S, Hogan B. Bmp signaling regulates proximal-distal differentiation of endoderm in mouse lung development. Development. 1999;126(18):4005–15.PubMedGoogle Scholar
- Malanchi I, Santamaria-MartÃnez A, Susanto E, Peng H, Lehr H-A, Delaloye J-F, et al. Interactions between cancer stem cells and their niche govern metastatic colonization. Nature. 2012;481(7379):85–9.View ArticleGoogle Scholar
- Choi YJ, Kim ST, Park KH, Oh SC, Seo JH, Shin SW, et al. The serum bone morphogenetic protein-2 level in non-small-cell lung cancer patients. Med Oncol. 2012;29(2):582–8.View ArticlePubMedGoogle Scholar
- Fei Z-H, Yao C-Y, Yang X-L, Huang X-E, Ma S-L. Serum BMP-2 up-regulation as an indicator of poor survival in advanced non-small cell lung cancer patients. Asian Pac J Cancer Prev. 2013;14:5293–9.View ArticlePubMedGoogle Scholar
- Ye X, Niu X, Tang N, Xu Y, Li Z, Yu Y, et al. Adenovirus mediated knockdown of bone morphogenetic protein 2 inhibits human lung cancer growth and invasion in vitro and in vivo. Int J Immunopathol Pharmacol. 2011;25(4):967–76.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.