Gene expression analysis of bone metastasis and circulating tumor cells from metastatic castrate-resistant prostate cancer patients
Journal of Translational Medicinevolume 14, Article number: 72 (2016)
Characterization of genes linked to bone metastasis is critical for identification of novel prognostic or predictive biomarkers and potential therapeutic targets in metastatic castrate-resistant prostate cancer (mCRPC). Although bone marrow core biopsies (BMBx) can be obtained for gene profiling, the procedure itself is invasive and uncommon practice in mCRPC patients. Conversely, circulating tumor cells (CTCs), which are likely to stem from bone metastases, can be isolated from blood. The goals of this exploratory study were to establish a sensitive methodology to analyze gene expression in BMBx and CTCs, and to determine whether the presence or absence of detectable gene expression is concordant in matching samples from mCRPC patients.
The CellSearch® platform was used to enrich and enumerate CTCs. Low numbers of PC3 prostate cancer (PCa) cells were spiked into normal blood to assess cell recovery rate. RNA extracted from recovered PC3 cells was amplified using an Eberwine-based procedure to obtain antisense mRNA (aRNA), and assess the linearity of the RNA amplification method. In this pilot study, RNAs extracted from CTCs and PCa cells microdissected from formalin-fixed paraffin-embedded BMBx, were amplified to obtain aRNA and assess the expression of eight genes functionally relevant to PCa bone metastasis using RT-PCR.
RNAs were successfully extracted from as few as 1–5 PCa cells in blood samples. The relative expression levels of reference genes were maintained after RNA amplification. The integrity of the amplified RNA was also demonstrated by RT-PCR analysis using primer sets that target the 5′-end, middle, and 3′-end of reference mRNA. We found that in 21 out of 28 comparisons, the presence or absence of detectable gene expression in CTCs and PCa cells microdissected from single bone lesions of the same patients was concordant.
This exploratory analysis suggests that aRNA amplification through in vitro transcription may be useful as a method to detect gene expression in small numbers of CTCs and tumor cells microdissected from bone metastatic lesions. In some cases, gene expression in CTCs and BMBxs was not concordant, raising questions about using CTC gene expression to make clinical decisions.
More than 80 % of patients with metastatic castrate-resistant prostate cancer (mCRPC) present skeletal metastases, which invariably lead to an incurable disease for which we only have treatments that become palliative in nature [1, 2]. Thus, the identification of gene signatures pertaining to prostate cancer (PCa) bone metastasis is fundamental to the development of novel therapeutic targets and/or the identification of prognostic or predictive biomarkers.
The vast majority of gene profile analyses using biopsies from PCa patients have been performed in primary tumors or metastatic lesions other than bone [3–5], which might not provide clues about gene transcriptional changes occurring in skeletal metastasis. Although gene analysis of bone marrow core biopsies (BMBxs) would be of importance to identify the genetic make-up of PCa cells that colonize and interact with the bone microenvironment, BMBxs are not routinely performed in mCRPC patients due to the time-consuming and invasive nature of the procedure. Furthermore, in most cases, the combination of imaging studies and clinical information has proven to provide diagnostic accuracy for the assessment of suspected bone metastases . In contrast, circulating tumor cells (CTCs) can be obtained repeatedly and non-invasively through routine blood draws, and isolated through different techniques utilizing cell-surface antigens, or other physical, and/or biological properties of cancer cells [7–9].
CTCs are rare cancer cells transported through the peripheral circulation. In localized, non-metastatic cancers, CTCs are thought to emanate from primary tumors. CTCs may also be released from metastatic lesions and found in peripheral blood of patients with advanced cancer, thus providing an opportunity for “liquid biopsy” that may offer information on the evolution of the disease during treatment [9, 10]. The assessment of CTC numbers using the CellSearch® platform has been validated as a tool of clinical utility to monitor treatment response and predict survival in mCRPC patients .
Previous studies by others have demonstrated clearly the capacity of gene expression analyses of RNA obtained in large amounts from bone metastases derived from rapid autopsies of PCa patients [12, 13]. Here, we wished to determine the feasibility of gene expression analysis of both CTCs and BMBxs in living mCRPC patients with bone scan evidence of skeletal metastasis. Since CTCs in mCRPC patients are likely to derive from bone metastatic deposits, we also aimed to explore similarities/differences in gene expression in both CTCs and BMBxs obtained from the same patients.
PC3 and LNCaP PCa cells were obtained from ATCC (Manassas, VA) and provided by Dr. Leland Chung (Cedars-Sinai Medical Center, Los Angeles, CA), respectively. Both cell lines were grown in RPMI-1640 with 10 % heat-inactivated fetal bovine serum (FBS). Authentication of the human cell lines used here was verified through short tandem repeat profiling by the Research Technology Support Facility of Michigan State University.
Cell spike-in experiments
For validation studies, 80 % confluent PC3 and LNCaP monolayers were trypsinized, and diluted in complete culture medium to adjust their cell density to 50 cells/mL. Cell suspensions (100 μL) were plated into the wells of 96-well cell culture plates. Twenty-four hours later, cells in all wells were counted under the microscope, and those with exactly five cells were trypsinized and added to 7.5 mL of blood from healthy volunteer donors drawn into CellSave® preservative tubes (Janssen Diagnostics, LLC). CellSave® tubes containing blood spiked with cultured PCa cells, as well as matched non-spiked blood samples serving as negative controls, were processed in triplicate using the semiautomated CellSearch® Circulating Epithelial Cell Kit (Janssen) in the CellTracks® Autoprep® system (Janssen) for cell enumeration at the Biobanking and Correlative Sciences Core at Karmanos Cancer Institute. Briefly, epithelial cells present in peripheral blood were magnetically captured with a ferrofluid-coupled antibody targeting the epithelial cell adhesion molecule (EpCAM), then immunostained with allophycocyanin (APC)-labeled antibodies to the leukocyte marker CD45, phycoerythrin (PE)-labeled antibodies to epithelial markers cytokeratins (CKs) 8, 18, and 19, and stained with the nuclear stain 4,2-diamidino-2-phenylindole-dihydrochloride (DAPI). The sample was transferred automatically to a cartridge in a MagNest, and finally scanned with the semi-automated fluorescence optical system Cell-Tracks Analyzer II® (Janssen). Objects preselected and displayed by the system in a gallery were defined as CTCs and counted by a trained operator if they were round to oval in shape, 4 µm in size or larger, positive for the epithelial marker (CK-PE) and nuclear stain (DAPI) with at least 50 % overlap between the CK-PE-positive cytoplasm and the nucleus, and negative for the leukocyte marker (CD45-APC). In addition, 5 PC3 or LNCaP cells were spiked into 7.5 mL of blood from healthy volunteer donors drawn into fixative-free K3EDTA tubes (BD Vacutainer®, Becton–Dickinson) (three experiments performed by three different investigators), to validate the sensitivity of the RNA amplification method used by us (see description below). To this end, the spiked blood was immediately processed using the CellSearch® Profile Kit (Janssen) and the CellTracks® Autoprep® system, and then the tube employed containing a final volume of 900 µL CTC enriched sample was placed in a DynaMag™-15 magnet (Invitrogen) for 10 min. After carefully aspirating the liquid without disturbing the ferrofluid bead pellet concentrated on the tube wall, 1 mL of TRIzol® (Invitrogen) was added, and the tube vortexed to lyse cells and inactivate nucleases. The lysate was transferred to a 1.5 mL RNase-free Eppendorf tube, and stored at −80 °C until RNA isolation and antisense RNA (aRNA) amplification was performed.
Patient enrollment and eligibility
Written informed consent was obtained from 24 mCRPC patients enrolled to participate in the human protocol # 2011-060 (Principal Investigator: Dr. Michael L. Cher), and approved by Karmanos Cancer Institute and Wayne State University Institutional Review Board. Study inclusion criteria were: (a) prior diagnosis of mCRPC characterized by rising PSA level or clinical disease progression despite a castrate level of serum testosterone; (b) clinical decision to start a new anti-cancer systemic therapy, and no treatment with any investigational drug within 2 weeks prior to blood draw/tissue biopsy proposed herein; (c) metastatic deposit visible on an imaging study obtained within 8 weeks prior to initiation of new systemic therapy; (d) radiological evidence of accessibility to a metastatic deposit in bone by computed tomography (CT); (e) performance status of 0–2 by ECOG/Zubrod criteria; (f) Absolute neutrophil count ≥1500 mm3, hemoglobin ≥9.0 g/dL, platelets ≥100,000/mm3; (g) Prothrombin Time and Partial Thromboplastin Time should be < institutional upper limit of normal; (h) presence of one or more CTCs per 7.5 mL, assessed using the Veridex CellSearch® Profile Kit assay, as explained above. The characteristics of the patients enrolled in this study are described in Table 1.
Blood sample processing
Blood from each patient was initially collected in CellSave® preservative tubes and processed for CTC enumeration as described above. Patients who were positive for CTCs were subjected to a new blood draw into K3EDTA anticoagulant tubes used in the CellSearch® Profile Kit, for CTC enrichment and RNA extraction as described above. The TRIzol® lysate obtained was transferred to a 1.5 mL RNase-free Eppendorf tube, and stored at -80 °C until needed.
Procurement and processing of tumor tissue from bone metastases
On the same day blood was drawn for CTC enrichment and RNA extraction, patients then underwent CT-guided BMBx using a battery-powered drill and biopsy needle set. Four to six BMBxs were obtained from the iliac bone of each patient. One of the cores was placed in an RNase-free Eppendorf tube and flash frozen and transported in liquid nitrogen for storage at −80 °C until processing. The remaining BMBxs were fixed for 24 h in 4 % paraformaldehyde in diethylpyrocarbonate (DEPC)-treated PBS, and decalcified in 10 % EDTA in autoclaved DEPC-treated water, pH 7.0, with agitation at room temperature for 3 days. After decalcification, each BMBx was progressively dehydrated with increasing concentrations of ethanol, and immediately paraffin-embedded using Precision Cut Paraffin (Thermo Scientific). All aqueous solutions were prepared using DEPC-treated water.
Laser capture microdissection
Two adjacent 5-µm sections were obtained from the formalin-fixed paraffin-embedded (FFPE) BMBxs, using at all times RNase-free technique. Immunohistochemistry (IHC) for cytokeratin was performed on one of the slides for identification of PCa cells, as previously described [14, 15], as PCa cells are the only cytokeratin-positive cells expected to be found within the BMBx. The adjacent tissue section was mounted onto an RNase-free polyethylene naphthalate (PEN) membrane glass slide (Arcturus), which was previously sprayed with RNase AWAY™ (Thermo Scientific), washed twice with DEPC-treated, nuclease-free water (Fisher Scientific), and then exposed to UV for 30 min under a laminar flow hood to render the PEN membrane more hydrophilic and improve adherence of the specimen. After air drying for about 2 h under a laminar flow hood, the section mounted on the slide was deparaffinized, hydrated with decreasing graded alcohols made with DEPC-treated water, and rapidly stained with Harris hematoxylin and alcoholic Eosin Y solution (H&E), washed twice with 100 % ethanol and xylene, then held in xylene until initiation of the laser capture microdissection (LCM) session. The H&E-stained slide was air dried and loaded onto the ArcturusXT™ LCM System. Metastatic lesions identified in the FFPE sample were excised using both the infrared (IR) and ultraviolet (UV) microdissection lasers, and collected on individual CapSure® Macro LCM Caps (Arcturus). When needed, the adjacent section immunostained for cytokeratin was used for guidance during the LCM session to identify epithelial (PCa) cells within the bone marrow. Tumor material was deemed to be inadequate for LCM and RNA collection if the area identified as tumor was smaller than 9000 µm2, with 150 PCa cells on average. To collect total RNA from the microdissected tissue, the CapSure® Macro LCM Cap containing the captured material was immediately placed into a 0.5 mL RNase-free microcentrifuge tube filled with 50 µL of TRIzol® reagent, vortexed for 5 min to lyse the cells, and the lysate obtained stored at −80 °C until needed.
RNA isolation and amplification of small amounts of RNA
Cultured PCa cells and CTCs were lysed with TRIzol® reagent, as described above. Flash frozen BMBxs were placed in pre-chilled microtubes containing 2.8 mm ceramic beads and TRIzol® reagent using a Precellys®24 homogenizer (Peqlab LLC), 1 cycle × 30 s at 6000 rpm at 4 °C, similarly as described by others . High-quality total RNA (DNA-free) was purified from cultured cells, CTCs, LCM, and frozen BMBxs lysed with TRIzol® using Direct-zol™ RNA MiniPrep (Zymo Research), as per manufacturer’s instructions. For gene analysis of CTCs enriched from patients’ blood, spiked-in cultured cells, or tumor cells metastatic to bone recovered by LCM, RNAs were amplified using an antisense mRNA (aRNA) amplification system based on the Eberwine’s procedure , using the MessageAmp™ II aRNA Amplification Kit (Invitrogen). The protocol was performed as recommended by the manufacturer, except for the replacement of ArrayScript™ RT by Superscript® VILO™ Master Mix (Invitrogen) and the replacement of 10 × First Strand buffer from MessageAmp™ II aRNA Amplification Kit (Invitrogen) by 5 × First Strand buffer from Superscript® II Reverse Transcriptase (Invitrogen). This was done to increase the reverse transcription (RT) of total RNA into a first strand of complimentary DNA (cDNA). SuperScript® VILO™ Master Mix contains a recombinant ribonuclease inhibitor and SuperScript® III RT, which is among the best performing reverse transcriptases in terms of reproducibility and sensitivity for low copy RNA levels . After RNase treatment, a second strand cDNA is generated by DNA polymerase. The resulting double-stranded cDNA was then used as a template for T7-RNA polymerase for in vitro transcription (IVT) into aRNA (also known as complimentary RNA, cRNA) and amplification. The procedure was repeated in a second round of amplification when additional aRNA yield was needed. After IVT, double-stranded cDNAs were removed by treatment with DNase I, and the amplified RNA was purified.
Concentrations of total RNA and amplified aRNA obtained from cells spiked in blood, as well as from amplified aRNA obtained from CTCs and flash frozen or FFPE BMBxs, were quantified by absorbance measurements using an Epoch™ Microplate Spectrophotometer for micro-volume analysis (BioTek). Purity of total RNA and amplified aRNA was assessed through the ratio of the absorbance of the samples at 260 and 280 nm (A260/A280).
Assessment of amplified aRNA integrity
Integrity of mRNA is usually assessed by the ratio of the 28S:18S rRNA species shown in denaturing agarose gel electrophoresis, based on the assumption that the quality of rRNA (more than 80 % of total RNA) reflects that of underlying mRNA. However, this approach cannot be used to assess the integrity of amplified aRNA, which does not contain rRNA and derives from minute amounts of total RNA. Therefore, we assessed the integrity of aRNA with RT-PCR (see below) using two sets of primers to probe different positions (5′ end, middle, and 3′ end regions) of EpCAM and GAPDH transcripts, by verifying the presence of each respective amplicon on agarose gels (Fig. 2a), following a strategy similar to that described by Nolan .
Gene expression analysis
Total RNA or amplified aRNA was reverse transcribed into cDNA using iScript™ cDNA synthesis kit (Bio-Rad) according to manufacturer’s protocol. For RT-PCR, resultant cDNAs were used as a template in a PCR reaction using DreamTaq DNA polymerase (Life Technologies). Forward and reverse primers used are listed in Table 2. The following amplification conditions were used: an initial denaturation at 95 °C for 3 min, followed by 35–40 cycles (except for GAPDH, 25 cycles) of 95 °C for 30 s, 52–55 °C for 30 s and 72 °C for 1 min, followed by a final extension 72 °C for 5 min. PCR products were resolved on a 2 % agarose gel and visualized by ethidium bromide staining. DNA bands were visualized using a ChemiDoc XRS gel documentation system (Bio-Rad).
For reverse transcriptase quantitative real-time quantitative PCR (RT-qPCR), the Mastercycler RealPlex2 (Eppendorf) real-time PCR system and GoTaq qPCR Master Mix (Promega) were used. Thermal cycle parameters were as follows: initial activation at 95 °C for 2 min, 40 cycles of denaturation at 95 °C for 15 s, annealing 55 °C for 15 s, and extension at 72 °C for 30 s. The mean cycle threshold (Ct) for each gene was normalized to levels of the housekeeping gene GAPDH in the same sample. Relative fold changes in expression for each gene were calculated by the delta–delta-CT method .
As a proof of concept, we selected eight genes for analysis with RT-PCR in CTCs and LCM BMBxs based on their relevance to PCa bone metastasis: (a) EpCAM, which codes for a transmembrane epithelial glycoprotein  overexpressed in adenocarcinomas  and used in the CellSearch® system to enrich CTCs via immunomagnetic separation; (b) PSA (prostate-specific antigen), a well know biomarker for PCa screening found to be positive in mCRPC patients ; (c) BMP-7 (bone morphogenetic protein-7), a member of the transforming growth factor-beta (TGF-β) family that is usually expressed in osteoblastic bone metastases of PCa , and increased in CRPC patients ; (d) MMP-14 (a.k.a. MT1-MMP), a membrane-tethered matrix metalloproteinase that we found to be expressed by PCa cells in skeletal metastases, and contribute to bone remodeling and intraosseous tumor growth ; (e) TMPRSS2-ERG gene rearrangement due to chromosomal translocations that fuse the androgen-regulated TMPRSS2 promoter with the ETS family transcription factor ERG, which is expressed in about half of advanced PCas [27–29] and has been associated with an increased risk death in PCa patients [30–32]; (f) IL-6 (interleukin-6), which is highly expressed by PCa cells with aggressive phenotype , and has been associated with resistance to chemotherapy in CRPC  and bone remodeling in PCa bone metastases [35, 36]; (g) Vimentin, an intermediate filament protein expressed in mesenchymal cells frequently used as a marker of epithelial to mesenchymal transition (EMT) ; and (h) GAPDH, a housekeeping gene used as a control.
Data obtained using RT-qPCR are presented as mean values ± SD and analyzed using the ANOVA test, with p < 0.05 considered as statistically significant.
Validation of linear amplification of RNA for gene analysis in small number of PCa cells
To mimic the CTC number threshold currently used to predict overall survival (OS) in men with mCRPC [11, 21, 38–40] and establish the method for RNA amplification from limited quantities of total RNA for gene expression, we spiked blood from healthy male volunteers with five PC3 or 5 LNCaP cells. Using the CellSearch® Circulating Epithelial Cell Kit in the CellTracks® Autoprep® system, we found that the recovery rate of PC3 and LNCaP cells spiked into normal blood was 76.7 ± 6.2 % (n = 6) and 80.2 ± 3.0 % (n = 6) at the five and 50 cell level, respectively, confirming high efficiency of PCa cell recovery from blood. No epithelial cells were found in non-spiked matched blood used as negative controls. After we confirmed the cell recovery rate using the Circulating Epithelial Cell enrichment system, we then spiked five PC3 or LNCaP cells in new blood of healthy men drawn in K3EDTA tubes and processed using the CellSearch® Profile Kit, used for RNA analysis. Minute amounts of total RNA extracted from the PC3 and LNCaP cells recovered through that CellSearch® Profile Kit and the CellTracks® Autoprep® system were subjected to one and two round amplification using the MessageAmp™ II aRNA Amplification Kit. As shown in Table 3, this methodology was efficient enough to generate sufficient quantities of aRNA of good quality (A260/A280 ratio of 1.8–2.0) from a few PCa cells for gene expression analyses, such as RT-PCR, RT-qPCR, and cDNA microarrays if desired. To establish the nature of linear amplification, we compared relative expression levels in a set of genes with high, medium, or low expression levels, in unamplified RNA extracted from 2 × 105 PC3 cells (~6–9 µg of total RNA, equivalent to ~180–270 pg of mRNA) and one- or two-round amplified aRNA derived from five PC3 cells (Table 3). The relative abundance of transcripts for the different genes assessed by RT-qPCR was maintained, confirming the linear amplification of RNA (Fig. 1).
Gene expression analysis in CTCs and BMBx from mCRPC patients with bone metastasis
To check the purity and quality of the aRNA amplified from RNA extracted from CTCs and PCa cells microdissected from BMBxs, we designed primers to target different regions of EpCAM and GAPDH transcripts, epithelium-specific and housekeeping genes, respectively (Fig. 2a). The primers targeting the 5′ end and the middle regions of the transcripts produced the expected size of the resultant PCR products as effectively as those targeting the 3′ end (Fig. 2b), indicating that the integrity of the aRNA amplified from both CTCs and microdissected FFPE PCa bone metastases was good. Having confirmed that the quality of the amplified aRNA was adequate, we analyzed all the CTC and BMBx samples included in the study using RT-PCR.
Of 24 patients with bone scan/CT scan evidence of one or more skeletal metastases enrolled in the study, nine had no CTCs and therefore did not undergo bone biopsy. Table 1 summarizes the demographic, clinical, and tumor characteristics of the 15 patients that had one or more EpCAM-positive CTCs in their blood. One patient was eliminated from the study due to undetectable aRNA in CTCs even after two-rounds of amplification. In addition, another patient was excluded because his BMBx showed sarcomatoid tumor cells. Of the 13 patients remaining, a second blood sample was obtained and subjected to CTC enrichment using the CellSearch® Profile kit for RNA extraction, and subsequent aRNA amplification. RT-PCR analysis of aRNA amplified from CTCs revealed that the gene expression varied from one patient to another. As expected, EpCAM, which codes for the protein used for CTC enrichment, was expressed in all cases. PSA and IL-6 were also expressed by CTCs from all patients (Fig. 3a). Unexpectedly, the transcript for Vimentin was expressed in many of patients’ CTCs (Fig. 3a). Importantly, aRNA amplification from as little as one CTC (e.g., patient 04 M–S) was sufficient to perform gene analysis, consistent with our results with amplified aRNA obtained from five cultured PCa cells.
Of 13 patients with detectable CTCs, 7 (54 %) had cytokeratin-positive cells on immunohistochemical analysis of their BMBx. This yield is consistent with that reported in the literature [41, 42]. However, only four had a large enough area occupied by cancer cells for LCM and RNA collection. In those cases, we compared aRNA amplified from RNA extracted from ground flash frozen BMBx with that obtained from LCM FFPE sections of matching BMBx of the same patient. In the latter, areas exclusively containing PCa cells, as identified by H&E and confirmed by IHC for pan-CK, were UV laser microdissected from the bone specimens mounted on membrane glass slides, and captured with the help of a gentle IR laser, using the ArcturusXT™ LCM System (Fig. 3b). We observed that, although RNA yield was better in ground frozen bone core biopsies, the expression level of most of the genes specific to PCa was considerably diminished compared to amplified aRNA derived from PCa cells microdissected from FFPE BMBxs after equal number of PCR cycles and normalization to GAPDH (not shown). Conversely, genes commonly expressed by bone cells, such as RANK (receptor activator of nuclear factor κ-B), were highly expressed in aRNA derived from ground frozen BMBx, but minimally or not expressed in pure PCa cell populations acquired through LCM in FFPE BMBx (not shown). These results demonstrate that analysis of frozen bone samples is limited by contamination with normal bone cells.
Using at all times RNase-free instruments and technique, and aqueous reagents prepared in nuclease-free water, we were able to amplify good quality aRNA from the total RNA collected from the FFPE BMBxs (Fig. 2b). Since the aRNA yield obtained was good, we analyzed in BMBxs the same genes studied in CTCs (Fig. 3c). In order to explore whether or not there was concordance in presence or absence of gene expression as detected by conventional RT-PCR in CTCs and matched BMBxs, we lined up digitally the bands obtained for all the genes amplified in CTCs and BMBxs (same number of cycles) side by side (Additional file 1: Figure S1). We found that in 21 out of 28 comparisons, the presence or absence of detectable gene expression in CTCs and PCa cells microdissected from single bone metastasis of the same patients was concordant.
Precise analysis of human tissue is necessary for evaluation of gene expression. In this study, we wished to develop methods of gene expression analysis of CTCs and BMBxs in patients with mCRPC. We found that LCM, in conjunction with reliable methods of gene amplification, is an option for the specific isolation and molecular analysis of small number of PCa cells present in BMBxs. In recent years, LCM systems with IR and UV lasers combined with efficient software have been developed to isolate homogeneous cells precisely identified from heterogeneous tissues based on morphological criteria, and complemented by IHC phenotyping of the cell type of interest. In these LCM systems, the low-energy IR laser is fired to adhere spots within the selected tissue area to the PEN membrane of the glass slide, while the UV rapidly and precisely cuts out the zone defined to collect it in a collection cap . In addition to RNA analysis, different types of molecular analyses can be performed on cells procured by LCM, including DNA and proteomic analyses [44, 45].
In this exploratory study with a small cohort of mCRPC patients, we demonstrated the feasibility of amplifying aRNA for gene analysis in limited numbers of PCa cells microdissected from FFPE BMBxs. We show that using a modified Eberwine’s procedure, adequate yields of RNA can be obtained for gene expression profiling. Our results obtained under tightly controlled conditions show good quality aRNA not only in freshly isolated CTCs, but also in aRNA amplified from FFPE BMBxs despite several studies having reported RNA degradation in formalin-fixed samples [46–48]. To avoid RNA degradation, we found that critical measures such as rapid processing of BMBxs after collection, preparation of all aqueous solutions using DEPC-treated water, and decontamination of all materials and surfaces with ready-to-use surfactants that remove RNA and RNases from lab equipment, are mandatory throughout the procedure. Previous studies have reported LCM and DNA sequencing from a frozen bone metastasis of a mCRPC patient , or gene expression analysis of total RNA directly obtained from snap-frozen bone marrow biopsies largely replaced by tumor  or frozen bone metastatic cores isolated at autopsy [12, 13, 51, 52]. However, in addition to our earliest analysis in a few PCa bone metastasis , the present study is the only report that, to the best of our knowledge, uses LCM and aRNA amplification for gene analysis of limited number of tumor cells microdissected from FFPE BMBx obtained from living mCRPC patients. In spite of better RNA recovery and quality ascribed to frozen tissue [43, 44, 54], we found that its use as homogenate might be inadequate for gene expression analysis of PCa bone micrometastasis, due to dilution of PCa-specific genes by other genes expressed predominantly by bone cells. In that sense, the use of LCM to identify and capture precise cells from FFPE BMBxs seems to be more accurate to study gene expression after aRNA amplification. Limitations of this approach include the yield of bone biopsies (~50 %) and the need for sufficient tissue (~9000 µm2) for LCM.
Besides the clinical utility of the CTC count to monitor treatment response and predict survival in mCRPC patients , enriched CTCs can also be used for molecular studies that may provide important clues to understand the biology behind metastatic dissemination of PCa. Analysis of CTCs may also lead to the discovery of novel predictive and prognostic biomarkers or demonstrate targets for therapy. Many groups have succeeded in genomic and transcriptomic profiling of CTCs in patients with PCa [29, 55–62]; however, to our knowledge, none of them have compared molecular features of CTCs with those of bone micrometastasis in mCRPC patients. Because significantly higher CTC numbers are detected in mCRPC patients with bone metastasis relative to those without bone metastatic lesions , CTCs in patients are likely to derive from bone metastatic deposits. Thus, we hypothesized that gene signatures of skeletal metastasis may be mirrored in CTCs in each mCRPC patient. Our study suggests some concordance in presence or absence of gene expression in CTCs and single bone lesions of the same mCRPC patients. An obvious limitation of our study is the small sample size and the biopsy-targeting of single bone lesions; future studies with more patients and perhaps more biopsy sites in individual patients may provide additional clarification.
Data collected from autopsy studies in men who died of mCRPC have revealed substantial heterogeneity among tumor cells in bone metastases within the same patient [63–65]. Therefore, we surmise that presence or absence of detectable gene expression in CTCs, rather than gene expression levels, may provide a more clinically relevant overview of the genomic landscape in each mCRPC patient. Detectable gene expression in CTCs may suggest a druggable target; however, this result may ultimately be unreliable due to heterogeneity of target gene expression among separate metastatic deposits within individual patients. Similarly, absence of target gene expression in CTCs may not rule out target gene expression in some metastatic deposits. Therefore, this exploratory analysis raises concerns about relying on CTCs to predict response to therapy.
In this study we used CellSearch®, an EpCAM-based platform utilized for capturing of epithelial tumor cells from peripheral blood of mCRPC patients. Like other CTC-enriching technologies that rely on the expression of epithelial markers, the CellSearch® system does not have the capacity to isolate CTCs that might have gone through EMT. However, we found that the genes that code for the epithelial marker EpCAM and the mesenchymal marker vimentin were co-expressed in many of our patients’ CTC samples, suggesting a partial EMT. These results are in agreement with the findings of other groups that reported expression of EMT markers in EpCAM-positive CTCs [66–68]. We also found Vimentin to be expressed by most of the pure PCa tissues microdissected from BMBx. Our results are supported by IHC studies by Sethi et al. that show vimentin overexpression at the invasive front of bone metastases of PCa patients and E-cadherin within the center of the lesion , and those of Bryden et al. who reported lower E-cadherin expression in poorly differentiated bone metastases than in more differentiated ones . The expression of mesenchymal markers in PCa metastatic to bone could be ascribed to osseous factors that offer a selective growth advantage for tumor cells with a mesenchymal phenotype, as opposed to other tissues (e.g., lung) that promote mesenchymal-to-epithelial transition in colonizing cancer cells [71, 72].
To our knowledge, this exploratory study is the first to demonstrate the feasibility of simultaneous gene analysis in small numbers of CTCs and laser-captured, microdissected bone metastases in living mCRPC patients. Our data suggest that partial EMT might take place in bone metastatic PCa cells and in CTCs captured using EpCAM-based enrichment platforms. In this study, we found some concordance in gene expression profile of CTCs and PCa cells microdissected from single bone lesions of the same patients. This exploratory analysis raises concerns about relying on CTCs to predict response to therapy. Further investigations using a larger cohort of patients and new technologies capable to enrich the full spectrum of CTCs are needed to define more precisely the genomic landscaping of mCRPC through the use of these “liquid biopsies”.
Inoue T, Segawa T, Kamba T, Yoshimura K, Nakamura E, Nishiyama H, et al. Prevalence of skeletal complications and their impact on survival of hormone refractory prostate cancer patients in Japan. Urology. 2009;73(5):1104–9. doi:10.1016/j.urology.2008.07.062.
Berruti A, Tucci M, Mosca A, Tarabuzzi R, Gorzegno G, Terrone C, et al. Predictive factors for skeletal complications in hormone-refractory prostate cancer patients with metastatic bone disease. Br J Cancer. 2005;93(6):633–8. doi:10.1038/sj.bjc.6602767.
Dunn TA, Fedor HL, De Marzo AM, Luo J. Molecular profiling of indolent human prostate cancer: tackling technical challenges to achieve high-fidelity genome-wide data. Asian J Androl. 2012;14(3):385–92. doi:10.1038/aja.2011.147.
Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18(1):11–22. doi:10.1016/j.ccr.2010.05.026.
Tomlins SA, Mehra R, Rhodes DR, Cao X, Wang L, Dhanasekaran SM, et al. Integrative molecular concept modeling of prostate cancer progression. Nat Genet. 2007;39(1):41–51. doi:10.1038/ng1935.
Ciray I, Astrom G, Sundstrom C, Hagberg H, Ahlstrom H. Assessment of suspected bone metastases. CT with and without clinical information compared to CT-guided bone biopsy. Acta Radiol. 1997;38(5):890–5.
Alix-Panabieres C, Pantel K. Challenges in circulating tumour cell research. Nat Rev Cancer. 2014;14(9):623–31. doi:10.1038/nrc3820.
Zhe X, Cher ML, Bonfil RD. Circulating tumor cells: finding the needle in the haystack. Am J Cancer Res. 2011;1(6):740–51.
Miyamoto DT, Sequist LV, Lee RJ. Circulating tumour cells-monitoring treatment response in prostate cancer. Nat Rev Clin Oncol. 2014;11(7):401–12. doi:10.1038/nrclinonc.2014.82.
Alix-Panabieres C, Pantel K. Circulating tumor cells: liquid biopsy of cancer. Clin Chem. 2013;59(1):110–8. doi:10.1373/clinchem.2012.194258.
de Bono JS, Scher HI, Montgomery RB, Parker C, Miller MC, Tissing H, et al. Circulating tumor cells predict survival benefit from treatment in metastatic castration-resistant prostate cancer. Clin Cancer Res. 2008;14(19):6302–9. doi:10.1158/1078-0432.CCR-08-0872.
Morrissey C, True LD, Roudier MP, Coleman IM, Hawley S, Nelson PS, et al. Differential expression of angiogenesis associated genes in prostate cancer bone, liver and lymph node metastases. Clin Exp Metastasis. 2008;25(4):377–88. doi:10.1007/s10585-007-9116-4.
Mehra R, Kumar-Sinha C, Shankar S, Lonigro RJ, Jing X, Philips NE, et al. Characterization of bone metastases from rapid autopsies of prostate cancer patients. Clin Cancer Res. 2011;17(12):3924–32. doi:10.1158/1078-0432.CCR-10-3120.
Wiesner C, Nabha SM, Dos Santos EB, Yamamoto H, Meng H, Melchior SW, et al. C-kit and its ligand stem cell factor: potential contribution to prostate cancer bone metastasis. Neoplasia. 2008;10(9):996–1003.
Dong Z, Bonfil RD, Chinni S, Deng X, Trindade Filho JC, Bernardo M, et al. Matrix metalloproteinase activity and osteoclasts in experimental prostate cancer bone metastasis tissue. Am J Pathol. 2005;166(4):1173–86. doi:10.1016/S0002-9440(10)62337-1.
Carter LE, Kilroy G, Gimble JM, Floyd ZE. An improved method for isolation of RNA from bone. BMC Biotechnol. 2012;12:5. doi:10.1186/1472-6750-12-5.
Van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH. Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Natl Acad Sci USA. 1990;87(5):1663–7.
Okello JB, Rodriguez L, Poinar D, Bos K, Okwi AL, Bimenya GS, et al. Quantitative assessment of the sensitivity of various commercial reverse transcriptases based on armored HIV RNA. PLoS One. 2010;5(11):e13931. doi:10.1371/journal.pone.0013931.
Nolan T, Hands RE, Bustin SA. Quantification of mRNA using real-time RT-PCR. Nat Protoc. 2006;1(3):1559–82. doi:10.1038/nprot.2006.236.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. doi:10.1006/meth.2001.1262.
Litvinov SV, Velders MP, Bakker HA, Fleuren GJ, Warnaar SO. Ep-CAM: a human epithelial antigen is a homophilic cell-cell adhesion molecule. J Cell Biol. 1994;125(2):437–46.
Went PT, Lugli A, Meier S, Bundi M, Mirlacher M, Sauter G, et al. Frequent EpCam protein expression in human carcinomas. Hum Pathol. 2004;35(1):122–8.
Greene KL, Albertsen PC, Babaian RJ, Carter HB, Gann PH, Han M, et al. Prostate specific antigen best practice statement: 2009 update. J Urol. 2009;182(5):2232–41. doi:10.1016/j.juro.2009.07.093.
Masuda H, Fukabori Y, Nakano K, Takezawa Y, CSuzuki T, Yamanaka H. Increased expression of bone morphogenetic protein-7 in bone metastatic prostate cancer. Prostate. 2003;54(4):268–74. doi:10.1002/pros.10193.
Morrissey C, Brown LG, Pitts TE, Vessella RL, Corey E. Bone morphogenetic protein 7 is expressed in prostate cancer metastases and its effects on prostate tumor cells depend on cell phenotype and the tumor microenvironment. Neoplasia. 2010;12(2):192–205.
Bonfil RD, Dong Z, Trindade Filho JC, Sabbota A, Osenkowski P, Nabha S, et al. Prostate cancer-associated membrane type 1-matrix metalloproteinase: a pivotal role in bone response and intraosseous tumor growth. Am J Pathol. 2007;170(6):2100–11. doi:10.2353/ajpath.2007.060720.
Tomlins SA, Rhodes DR, Perner S, Dhanasekaran SM, Mehra R, Sun XW, et al. Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer. Science. 2005;310(5748):644–8. doi:10.1126/science.1117679.
Helgeson BE, Tomlins SA, Shah N, Laxman B, Cao Q, Prensner JR, et al. Characterization of TMPRSS2:eTV5 and SLC45A3:ETV5 gene fusions in prostate cancer. Cancer Res. 2008;68(1):73–80. doi:10.1158/0008-5472.CAN-07-5352.
Danila DC, Anand A, Sung CC, Heller G, Leversha MA, Cao L, et al. TMPRSS2-ERG status in circulating tumor cells as a predictive biomarker of sensitivity in castration-resistant prostate cancer patients treated with abiraterone acetate. Eur Urol. 2011;60(5):897–904. doi:10.1016/j.eururo.2011.07.011.
Demichelis F, Fall K, Perner S, Andren O, Schmidt F, Setlur SR, et al. TMPRSS2:ERG gene fusion associated with lethal prostate cancer in a watchful waiting cohort. Oncogene. 2007;26(31):4596–9. doi:10.1038/sj.onc.1210237.
Attard G, Clark J, Ambroisine L, Fisher G, Kovacs G, Flohr P, et al. Duplication of the fusion of TMPRSS2 to ERG sequences identifies fatal human prostate cancer. Oncogene. 2008;27(3):253–63. doi:10.1038/sj.onc.1210640.
John J, Powell K, Conley-Lacomb MK, Chinni SR. TMPRSS2-ERG fusion gene expression in prostate tumor cells and its clinical and biological significance in prostate cancer progression. J Cancer Sci Ther. 2008;4(4):94–101. doi:10.4172/1948-5956.1000119.
Wegiel B, Bjartell A, Culig Z, Persson JL. Interleukin-6 activates PI3 K/Akt pathway and regulates cyclin A1 to promote prostate cancer cell survival. Int J Cancer. 2008;122(7):1521–9. doi:10.1002/ijc.23261.
Ara T, Declerck YA. Interleukin-6 in bone metastasis and cancer progression. Eur J Cancer. 2010;46(7):1223–31. doi:10.1016/j.ejca.2010.02.026.
Domingo-Domenech J, Oliva C, Rovira A, Codony-Servat J, Bosch M, Filella X, et al. Interleukin 6, a nuclear factor-kappaB target, predicts resistance to docetaxel in hormone-independent prostate cancer and nuclear factor-kappaB inhibition by PS-1145 enhances docetaxel antitumor activity. Clin Cancer Res. 2006;12(18):5578–86. doi:10.1158/1078-0432.CCR-05-2767.
Morrissey C, Lai JS, Brown LG, Wang YC, Roudier MP, Coleman IM, et al. The expression of osteoclastogenesis-associated factors and osteoblast response to osteolytic prostate cancer cells. Prostate. 2010;70(4):412–24. doi:10.1002/pros.21075.
Mendez MG, Kojima S, Goldman RD. Vimentin induces changes in cell shape, motility, and adhesion during the epithelial to mesenchymal transition. FASEB J. 2010;24(6):1838–51. doi:10.1096/fj.09-151639.
Thalgott M, Rack B, Maurer T, Souvatzoglou M, Eiber M, Kress V, et al. Detection of circulating tumor cells in different stages of prostate cancer. J Cancer Res Clin Oncol. 2013;139(5):755–63. doi:10.1007/s00432-013-1377-5.
Scher HI, Jia X, de Bono JS, Fleisher M, Pienta KJ, Raghavan D, et al. Circulating tumour cells as prognostic markers in progressive, castration-resistant prostate cancer: a reanalysis of IMMC38 trial data. Lancet Oncol. 2009;10(3):233–9. doi:10.1016/S1470-2045(08)70340-1.
Olmos D, Arkenau HT, Ang JE, Ledaki I, Attard G, Carden CP, et al. Circulating tumour cell (CTC) counts as intermediate end points in castration-resistant prostate cancer (CRPC): a single-centre experience. Ann Oncol. 2009;20(1):27–33. doi:10.1093/annonc/mdn544.
Ross RW, Halabi S, Ou SS, Rajeshkumar BR, Woda BA, Vogelzang NJ, et al. Predictors of prostate cancer tissue acquisition by an undirected core bone marrow biopsy in metastatic castration-resistant prostate cancer–a Cancer and Leukemia Group B study. Clin Cancer Res. 2005;11(22):8109–13. doi:10.1158/1078-0432.CCR-05-1250.
Spritzer CE, Afonso PD, Vinson EN, Turnbull JD, Morris KK, Foye A, et al. Bone marrow biopsy: RNA isolation with expression profiling in men with metastatic castration-resistant prostate cancer–factors affecting diagnostic success. Radiology. 2013;269(3):816–23. doi:10.1148/radiol.13121782.
Bidarimath M, Edwards AK, Tayade C. Laser capture microdissection for gene expression analysis. Methods Mol Biol. 2015;1219:115–37. doi:10.1007/978-1-4939-1661-0_10.
Mukherjee S, Rodriguez-Canales J, Hanson J, Emmert-Buck MR, Tangrea MA, Prieto DA, et al. Proteomic analysis of frozen tissue samples using laser capture microdissection. Methods Mol Biol. 2013;1002:71–83. doi:10.1007/978-1-62703-360-2_6.
Domazet B, Maclennan GT, Lopez-Beltran A, Montironi R, Cheng L. Laser capture microdissection in the genomic and proteomic era: targeting the genetic basis of cancer. Int J Clin Exp Pathol. 2008;1(6):475–88.
Chung JY, Braunschweig T, Hewitt SM. Optimization of recovery of RNA from formalin-fixed, paraffin-embedded tissue. Diagn Mol Pathol. 2006;15(4):229–36. doi:10.1097/01.pdm.0000213468.91139.2d.
von Ahlfen S, Missel A, Bendrat K, Schlumpberger M. Determinants of RNA quality from FFPE samples. PLoS One. 2007;2(12):e1261. doi:10.1371/journal.pone.0001261.
Arzt L, Kothmaier H, Quehenberger F, Halbwedl I, Wagner K, Maierhofer T, et al. Evaluation of formalin-free tissue fixation for RNA and microRNA studies. Exp Mol Pathol. 2011;91(2):490–5. doi:10.1016/j.yexmp.2011.05.007.
Van Allen EM, Foye A, Wagle N, Kim W, Carter SL, McKenna A, et al. Successful whole-exome sequencing from a prostate cancer bone metastasis biopsy. Prostate Cancer Prostatic Dis. 2014;17(1):23–7. doi:10.1038/pcan.2013.37.
Stanbrough M, Bubley GJ, Ross K, Golub TR, Rubin MA, Penning TM, et al. Increased expression of genes converting adrenal androgens to testosterone in androgen-independent prostate cancer. Cancer Res. 2006;66(5):2815–25. doi:10.1158/0008-5472.CAN-05-4000.
Larson SR, Zhang X, Dumpit R, Coleman I, Lakely B, Roudier M, et al. Characterization of osteoblastic and osteolytic proteins in prostate cancer bone metastases. Prostate. 2013;73(9):932–40. doi:10.1002/pros.22639.
Zhang X, Coleman I, Brown LG, True LD, Kollath L, Lucas JM, et al. SRRM4 expression and the loss of REST activity may promote the emergence of the neuroendocrine phenotype in castration-resistant prostate cancer. Clin Cancer Res. 2015. doi:10.1158/1078-0432.CCR-15-0157.
Mainetti LE, Zhe X, Diedrich J, Saliganan AD, Cho WJ, Cher ML, et al. Bone-induced c-kit expression in prostate cancer: a driver of intraosseous tumor growth. Int J Cancer. 2015;136(1):11–20. doi:10.1002/ijc.28948.
Li J, Gregory SG, Garcia-Blanco MA, Armstrong AJ. Using circulating tumor cells to inform on prostate cancer biology and clinical utility. Crit Rev Clin Lab Sci. 2015;52(4):191–210. doi:10.3109/10408363.2015.1023430.
Giesing M, Driesel G, Molitor D, Suchy B. Molecular phenotyping of circulating tumour cells in patients with prostate cancer: prediction of distant metastases. BJU Int. 2012;110(11 Pt C):1202–11. doi:10.1111/j.1464-410X.2012.11534.x.
Attard G, Swennenhuis JF, Olmos D, Reid AH, Vickers E, A’Hern R, et al. Characterization of ERG, AR and PTEN gene status in circulating tumor cells from patients with castration-resistant prostate cancer. Cancer Res. 2009;69(7):2912–8. doi:10.1158/0008-5472.CAN-08-3667.
Miyamoto DT, Lee RJ, Stott SL, Ting DT, Wittner BS, Ulman M, et al. Androgen receptor signaling in circulating tumor cells as a marker of hormonally responsive prostate cancer. Cancer Discov. 2012;2(11):995–1003. doi:10.1158/2159-8290.CD-12-0222.
Cann GM, Gulzar ZG, Cooper S, Li R, Luo S, Tat M, et al. mRNA-Seq of single prostate cancer circulating tumor cells reveals recapitulation of gene expression and pathways found in prostate cancer. PLoS One. 2012;7(11):e49144. doi:10.1371/journal.pone.0049144.
Jiang Y, Palma JF, Agus DB, Wang Y, Gross ME. Detection of androgen receptor mutations in circulating tumor cells in castration-resistant prostate cancer. Clin Chem. 2010;56(9):1492–5. doi:10.1373/clinchem.2010.143297.
Leversha MA, Han J, Asgari Z, Danila DC, Lin O, Gonzalez-Espinoza R, et al. Fluorescence in situ hybridization analysis of circulating tumor cells in metastatic prostate cancer. Clin Cancer Res. 2009;15(6):2091–7. doi:10.1158/1078-0432.CCR-08-2036.
Magbanua MJ, Sosa EV, Scott JH, Simko J, Collins C, Pinkel D, et al. Isolation and genomic analysis of circulating tumor cells from castration resistant metastatic prostate cancer. BMC Cancer. 2012;12:78. doi:10.1186/1471-2407-12-78.
Lohr JG, Adalsteinsson VA, Cibulskis K, Choudhury AD, Rosenberg M, Cruz-Gordillo P, et al. Whole-exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer. Nat Biotechnol. 2014;32(5):479–84. doi:10.1038/nbt.2892.
Shah RB, Mehra R, Chinnaiyan AM, Shen R, Ghosh D, Zhou M, et al. Androgen-independent prostate cancer is a heterogeneous group of diseases: lessons from a rapid autopsy program. Cancer Res. 2004;64(24):9209–16. doi:10.1158/0008-5472.CAN-04-2442.
Suzuki H, Freije D, Nusskern DR, Okami K, Cairns P, Sidransky D, et al. Interfocal heterogeneity of PTEN/MMAC1 gene alterations in multiple metastatic prostate cancer tissues. Cancer Res. 1998;58(2):204–9.
Roudier MP, True LD, Higano CS, Vesselle H, Ellis W, Lange P, et al. Phenotypic heterogeneity of end-stage prostate carcinoma metastatic to bone. Hum Pathol. 2003;34(7):646–53.
Armstrong AJ, Marengo MS, Oltean S, Kemeny G, Bitting RL, Turnbull JD, et al. Circulating tumor cells from patients with advanced prostate and breast cancer display both epithelial and mesenchymal markers. Mol Cancer Res. 2011;9(8):997–1007. doi:10.1158/1541-7786.MCR-10-0490.
Chen CL, Mahalingam D, Osmulski P, Jadhav RR, Wang CM, Leach RJ, et al. Single-cell analysis of circulating tumor cells identifies cumulative expression patterns of EMT-related genes in metastatic prostate cancer. Prostate. 2013;73(8):813–26. doi:10.1002/pros.22625.
Pal SK, He M, Wilson T, Liu X, Zhang K, Carmichael C, et al. Detection and phenotyping of circulating tumor cells in high-risk localized prostate cancer. Clinical Genitourin Cancer. 2015;13(2):130–6. doi:10.1016/j.clgc.2014.08.014.
Sethi S, Macoska J, Chen W, Sarkar FH. Molecular signature of epithelial-mesenchymal transition (EMT) in human prostate cancer bone metastasis. Am J Transl Res. 2010;3(1):90–9.
Bryden AA, Freemont AJ, Clarke NW, George NJ. Paradoxical expression of E-cadherin in prostatic bone metastases. BJU Int. 1999;84(9):1032–4.
Turley EA, Veiseh M, Radisky DC, Bissell MJ. Mechanisms of disease: epithelial-mesenchymal transition–does cellular plasticity fuel neoplastic progression? Nat Clin Pract Oncol. 2008;5(5):280–90. doi:10.1038/ncponc1089.
Kyprianou N. ASK-ing EMT not to spread cancer. Proc Natl Acad Sci USA. 2010;107(7):2731–2. doi:10.1073/pnas.0914721107.
RDB and H-RCK conceived and designed the experiments. WJC, DSMO, AJN, and LEM performed the experiments. HAD performed image-guided bone marrow core biopsies. RDB, H-RCK, MLC, and EH analyzed the data. EH and MLC contributed patients’ biospecimens. RDB wrote the paper. All authors read and approved the final manuscript.
This work was supported by a Grant from The Sage Foundation (RDB and MLC, support for KCI Bone Metastasis Biobank), Department of Defense grant W81XWH-11-1-0050 (subcontract to HRCK, RDB and MLC), and Wayne State School of Medicine (start-up funds to RDB). The authors would like to thank the patients who consented to provide blood and bone biopsies that were used in this study.
The authors declare that they have no competing interests.