Received date: October 24, 2013; Accepted date: November 19, 2013; Published date: November 21, 2013
Citation: Kadam S, Cleverly AL, Farmen M, Grondin J, Cox YI, et al. (2013) A Canonical Transforming Growth Factor Beta-Dependent Signaling Pathway is Present in Peripheral Blood Cells of Cancer Patients with Skeletal Metastasis. J Mol Biomark Diagn 4:153. doi:10.4172/2155-9929.1000153
Copyright: © 2013 Kadam S, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Molecular Biomarkers & Diagnosis
Background: Patients with advanced metastatic cancer often have elevated levels of circulating transforming growth factor beta (TGF-β) that is thought to stimulate receptor mediated signaling through phosphorylated SMAD transcription factors in peripheral blood mononuclear cells.
Methods: To identify this TGF-β-dependent gene expression profile in cancer patients, we first evaluated a multi-gene expression profile from ex-vivo treated PBMCs with TGF-β1 stimulation. Change in expression, when challenged with a specific TGF-β receptor type kinase inhibitor was then derived to be ligand and inhibitor specific. Once this profile was established, we examined its role in identifying the activation of a canonical TGF-β1 signaling pathway in patient samples.
Results: We discovered a 37 gene sub-set where the expression profile from ex-vivo treated PBMCs was significantly associated with SMAD phosphorylation in cancer patients. We found significant correlation between the ex-vivo derived expression signature and circulating levels of TGF-β1 in patient samples. Additionally, we report association between the expression profile and the presence of several plasma proteins in disease samples that are known to be concomitantly present with TGF-β- dependent pathway activation.
Conclusion: An expression profile for TGF-β1 cytokine mediated-signaling in cancer patients is identified, which may serve as a biomarker to measure the pharmacodynamic effect of TGF-β inhibitors during clinical drug development and as marker of disease diagnosis.
Peripheral blood mononuclear cells, Skeletal metastasis, Transforming growth factor beta, LY2157299 monohydrate, SMAD
TGF-β1: Transforming Growth Factor Beta One; TGF-βRI: Transforming Growth Factor Beta Receptor Type One; SMAD: Mothers Against Decapentaplegic Homolog; PBMCs: Peripheral Blood Mononuclear Cells; CPT: Cell Preparation Tubes; PBS: Phosphate Buffered Saline; MAIP: Multianalyte Immunoassay Panel
The transforming growth factor beta (TGF-β) pathway plays a central role in a wide range of biological processes essential to embryonic development and tissue homeostasis .The canonical signaling event induced by TGF-β superfamily members is initiated by ligand-mediated hetero-oligomerization of the serine-threonine kinase receptors, which phosphorylate signaling of Regulatory Smads (RSmads),such as Smad 2/3 for TGF-βs and Smad1/5/8 for Bone Morphogenetic Proteins (BMPs).The activated R-Smads then form a heterocomplex with a common Smad (Co-Smad or Smad4) and the complex is translocated into the nucleus, resulting in transcriptional regulation of abroad repertoire of target genes [2,3]. Other Inhibitory SMADs (I-SMADs) such as SMAD6 and SMAD7 generally inhibit BMP signaling and TGF-β/activin signaling respectively and serve as negative regulators of TGF-β dependent pathway signaling [4,5]. In addition to this canonical Smad-dependent signaling, TGF-β and related ligands are also capable of activating small GTPases, mitogenactivated protein kinases, and phosphatidylinositol-3-kinase, initiating a so-called Smad-independent, non-canonical response .
TheTGF-β superfamily of ligands is normally thought to contribute to both normal development and tumor suppression . However during the course of cancer progression, the normal tumor suppression function of TGF-β is often lost and its over-expression is associated with creating a local immunosuppressive environment that enhances tumor growth and exacerbates invasive and metastatic tumor cell behavior [2-4].
It has been previously reported that elevation in pSMAD levels in Peripheral Blood Mononuclear Cells (PBMCs) can be detected after ex-vivo stimulation with TGF-β1 and this activation signal might represent a useful biomarker for evaluation of TGF-β inhibitors [6-8]. Furthermore, this potential for cytokine induced phosphorylation of SMAD was found to be correlated with serum TGF-β1 levels [6-8]. Separately, we have previously reported on the specificity of LY2157299 monohydrate kinase and its potent activity as an inhibitor of TGF-β receptor type one (TGF-βRI) mediated signaling using several pathway related protein markers including the phosphorylation of SMADs in tumor xenografts . However, quantitative measurement of pSMAD from tumor derived TGF-β requires separation of PBMCs from whole blood under conditions that minimize the effect of inherent platelet derived TGF-β; a significant pre-analytic source of variability. Secondly, this baseline variability in PBMC pSMAD measurement makes it difficult to estimate the residual post-translation SMAD phosphorylation capacity during an external cytokine challenge. We examined whether TGF-β1 addition to PBMCs ex-vivo, was accompanied by changes in global gene expression and whether addition of a specific TGF-βRI kinase inhibitor produced a specific changes in this expression pattern. We then reasoned that measurement of this change in gene expression might provide a handle to evaluate the influence of circulating TGF-β1 in the blood of advanced cancer patients. Independent for TGF-β-related signaling, this approach of using microarray-based transcriptional profiling from peripheral blood, has been successfully applied to evaluate the diagnostic potential of transcript derived “signatures” in multiple sclerosis , chronic fatigue syndrome , sickle cell anemia , post-traumatic stress disorder , Huntington´s disease  and acute respiratory distress syndrome . Hence, while not novel, the approach can help with the dose determination of new TGF-β inhibitors.
In parallel, over the last decade expression profiling procedures have been optimized with regard to the development of optimal standardized procedures for blood collection and subsequent RNA extraction [16,17]. Further improvement in reproducibility was possible with the assessment of precision of the various components of variability associated with sample handling and array data processing steps . These advances have led to the recognition of expression signature based “intrinsic subtypes” in tissues and are now used to classify clinical specimens into phenotypic and prognostic categories [19-21].
We have investigated the role of PBMCs and whole blood, as potential “surrogates” for the development of biomarkers of TGF-β dependent activation. The ex-vivo procedure developed for detection of the TGF-β1 dependent stimulation of normal pSMAD in PBMCs was used to follow changes in global gene expression with LY2157299 monohydrate; a specific inhibitor of TGF-βRI kinase. The measurement of ex-vivo transcript level changes was performed using RNA microarray and the resulting expression profiles derived from the ex-vivo assay were used to characterize expression patterns in the whole blood of advanced cancer patients. The underlying hypothesis assumes that patients with advanced cancer with elevated levels of plasma TGF-β will produce a measurable change in the pattern of expression in PBMCs and mirror the ongoing biological processes during cytokine dependent cancer progression. We found that a TGF-β stimulated and LY2157299 reversed 37-gene expression profile, derived from the ex-vivo PBMC analysis, was also present in whole blood from cancer patients. Expression index of the composite genes from patient blood was significantly correlated with PF4 adjusted pSMAD in corresponding PBMCs and circulating levels of TGF-β1 and other co-expressed factors known to be associated with canonical pathway activation.
Cell lines, compounds and reagents
Recombinant human TGF-β1 was purchased from R&D Systems, (catalog #240-B, Minneapolis, MN, USA). TRIzol cell lysis and RNA preservation reagent (Invitrogen, Carlsbad, CA) was used according to manufacturer’s recommendations. Tempus tubes (Applied Biosystems, Carlsbad, CA) and PAX gene blood RNA tubes (PreAnalytiX/Qiagen, Feldbachstrasse CH-8634 Hombrechtikon, Switzerland) were used for collection of whole blood and RNA preservation and extraction was performed according to manufacturer´s recommendations. LY2157299 monohydrate is a small molecule inhibitor of the TGF-β1 receptor kinase [9,22-24]. Microarray analysis was performed using Affymetrix, HG-U133A human array chips (Affymetrix, Santa Clara, CA). Blood used to recover PBMCs for ex-vivo experiments was from normal subjects and was collected in Cell Preparation Tubes (CPT) from Becton Dickinson (Franklin Lakes, NJ, USA).
PBMC separation and ex-vivo TGF-β1 stimulation assay
CPTs obtained from Becton Dickinson (Franklin Lakes, NJ, USA), containing a polyester gel and sodium citrate anticoagulant were used for blood collection and stored upright at room temperature for at least 2 h prior to sample collection to reduce temperature gradient shock during blood draw. Upon blood collection, the tube was immediately centrifuged in a clinical swinging bucket centrifuge at 800 g for 25 min at room temperature. The PBMC layer just above the gradient gel was gently resuspended into the separated upper plasma by slowly inverting the tube. Two 1.5 ml fractions of the supernatant located above the polymer barrier were transferred into two 12×75 mm polypropylene tubes. The polypropylene tubes were centrifuged at 1,600 g for 15 min and the supernatant was carefully removed. The PBMC pellet was resuspended in an equivalent volume of Phosphate Buffered Saline (PBS) with 1% Bovine Serum Albumin (BSA) recentrifuged at 1,600 g and washed once in PBS with 1% BSA. Finally, pellets from two tubes per donor were resuspended in PBS with 1% BSA and pooled for treatment. A 200 μl aliquot of the 3.0 ml PBMC sample was used for each duplicate treatment. A stock was prepared and appropriate volumes were added to ligand containing tubes to achieve a final concentration of 1, 10, 50 and 100 pM. An equivalent volume of PBS with 1% BSA was added to the unstimulated PBMC control “vehicle” sample tube. Tubes were incubated at 37°C in a Dubnoff metabolic water bath (Precision Scientific, Chicago, IL, USA) for 4 hours with shaking at 60 rpm. Where stated LY2157299 monohydrate was added to achieve a 5 μM concentration and the reaction incubated at 37°C for 5 min prior to addition of TGF-β1. At the end of the reaction period, samples were centrifuged for 15 min at 1,600 g and 800 μl of TriZol reagent was added to each tube followed by vigorous vortexing to allow cellular lysis and stabilize released RNA. Samples were then stored at -80°C and recovered for RNA extraction as necessary.
Sample collection and processing
Adult patients diagnosed with histological or cytological advanced or metastatic cancer involving the bone; able to safely donate approximately 100 ml of blood over a 2 to 4-week period; and not receiving any chemotherapy or radiotherapy while in the study were included . A total number of 28 patients were included and wholeblood and plasma was used for this study. Ten healthy volunteers were used to draw blood for expression analysis, 5 each in Paxgene and Tempus tubes. Among cancer patients, the majority were females (75%). The tumor type were breast cancer (n=19, all females), multiple myeloma (n=4, 2 males/2 females), prostate cancer (n=4) and renal cell carcinoma (n=1, male). A written informed consent was obtained in accordance with federal, state and institutional guidelines. Details on this patient population have been published elsewhere . Whole blood and serum samples were collected from patients with skeletal metastasis that were enrolled in a nondrug interventional study to donate blood. Visits were scheduled between a minimum of 48 h and a maximum of 7 days apart. Patients had one venipuncture at each visit for a maximum of four visits. At visit one, blood samples were collected for TGF-β1, Platelet Factor 4 (PF4), Multianalyte Immunoassay Panel (MAIP) and SMAD assay, which assessed phosphorylated SMAD (pSMAD) and total SMAD (tSMAD). PBMCs were obtained for baseline and ex-vivo assays to determine pSMAD and tSMAD expression as previously described [25,26]. Briefly, blood was drawn into three 8-mlvacutainer CPT containing 1 ml of 0.1 M sodium citrate for the SMAD assay. Tubes were inverted five times and centrifuged at room temperature for 25 min in a VanGuard V6500 centrifuge (SmithKline Beecham, PA, USA) at 1,300 g and processed to measure pSMAD and tSMAD levels in PBMCs. PBMCs were stimulated with TGF-β1 and the treated cells (both stimulated or unstimulated) were lysed and stored frozen at -80°C prior to analysis of pSMAD and total SMAD by ELISA as described elsewhere . Patient samples for RNA extraction and microarray analysis were collected in tempus tubes (ABI, Foster City, CA). A minimum volume of 3 ml blood was drawn in the 10 ml tempus tube containing 6 ml Red Blood Cell (RBC) lysis and RNA preservative solution. The tubes were shaken vigorously and stored at -20°C until RNA extraction. Whole blood for testing in the Multianalyte Immunoassay Panel (MAIP) was collected in sodium citrate tubes and centrifuged at 1,500 g for 15 minutes at 2-6°C. The topmost plasma portion was removed and recentrifuged under the same conditions. The resultant plasma sample was stored at -80°C until analysis.
RNA microarray, multianalyte immunoassay panel (MAIP), pSMAD and TGF-β1 analysis
Total RNA for microarray analysis from PBMCs was isolated from TRIzol aliquots following manufacturer’s instructions, and purified using an RNeasy kit (Qiagen Inc., Valencia, CA). RNA from whole blood collected in Tempus tubes was extracted using globin mRNA reduction kit (GLOBIN NcClear, ABI, Foster City, CA). RNA integrity and yield were assessed by absorbance at 260 and 280 nm and analysis with the Agilent RNA 6000 Nano LabChip (Agilent Technologies, Inc., Palo Alto, CA). 5 μg of RNA with OD 260:280 ratios of 1.8 and clear 18S and 28S ribosomal RNA bands on the Agilent chip (RIN>6) was processed for application to Affymetrix U133A human array chip. TGF-β levels were measured using a commercially available ELISA (R&D Systems, Minneapolis, MN, catalog # DB100) per the manufacturer’s protocol. Briefly, plasma samples underwent initial acidification/neutralization steps immediately prior to analysis to activate latent TGF-β to an immunoreactive form, under conditions which have been shown to minimize re-association with TGF-β binding proteins . The ELISA assay utilized soluble TGF-βRII coated on the ELISA plate for antigen capture, with detection via a HRP-linked polyclonal specific for TGF-β1. Recombinant human TGF-β1 (rhTGF-β1) included in the kit was used as a standard. In addition to TGF-β1, separate aliquots from these same plasma samples were also batch analyzed for additional markers of platelet activation, PF4 and Beta-Thromboglobulin (BTG) using commercial enzyme immunoassays (Asserachrom® PF4 and Asserachrom® BTG, both from Diagnostica Stago, Asnieres-Sur-Siene, France). Plasma samples were tested at Rules-Based Medicine (Austin, TX) in their human antigen Multiple Antigenic Peptides (MAP) platform.
Microarray data collection and statistical data analysis
Gene expression values were generated using Affymetrix Microarray Analysis Suite MAS5 algorithm to adjust array signal normalization, background adjustment and probe set/signal summarization. Array data that passed acceptable QC metrics for percent present call, scale factor and false positives was log2 transformed. The log2 transformed data was then standardized at the probeset level to generate Z-scores: where X is the expression value, is the mean expression value, S is the standard deviation, and n represents the number of samples . Standardization was performed separately for each probe set derived from ex-vivo PBMC dataset and the patient sample dataset. The array profile of processed gene signal values were clustered and converted to heat maps using Spot fire Decision Suite. Array profiles were also integrated using the Cytoscape software for biological network visualization and the accompanying plug-in, BiNGO was used to map statistically significant gene ontology using a p-value of 0.01 . Expression levels are presented as normalized z-scores which are derived from averaged probe set signal values across probe sets and annotated genes altered under various experimental conditions and in the whole blood of disease samples. Since some genes are down regulated by TGF-β1 stimulation, the sign on the z-score is reversed before the average is taken. This composite value is referred to as the expression index .
Pearson correlation and partial correlation was used to assess associations. To determine association with expression index of LY2157299 monohydrate reversed genes and plasma based proteins we adjusted the response variable to account for platelet effects using PF4 levels as previously described . Differential expression between ex-vivo PBMC treated groups was assessed by random effects ANOVA (SAS Proc mixed). Significance of gene expression differences between TGF-β1 stimulated and non-stimulated PBMCs was made by calculating the False Discovery Rate (FDR) at less than 0.1. Reversal of TGF-β1 effects by LY2157299 monohydrate was judged by p-values less than 0.05.
Development of a TGF-β1-dependent expression signature with PBMCs from healthy donors
Previous studies have confirmed that cell lines exposed to TGF-β1 produce specific changes in the expression of genes associated with TGF-βRI mediated signaling [8,30,31]. More recently Allendorph et al.  described engineered ligands and their ability to functionally dissect the diversity of the transforming growth factor receptor super family activated response. We reasoned that normal PBMCs exposed to TGF-β1 might produce a similar measurable “canonical” change in global expression provided the collection and handling of whole blood and the separation of PBMCs was adequately controlled to minimize the release of platelet derived transforming growth factors. The addition of a specific inhibitor of TGF-βRI dependent pathway signaling in such an ex-vivo model system would then allow characterization of the specific cytokine/inhibitor dependent expression response. A schematic of the experimental overview is represented in Figure 1. Using PBMCs from healthy donors we determined global expression patterns at TGF-β1 levels ranging from 0-100 pM, with and without pre incubation with 5 μM LY2157299 monohydrate. Development of this ex-vivo PBMC stimulation assay, its ability to detect TGF-β1 dependent response as measured by phosporylation of SMAD proteins and the corresponding inhibition with LY2157299 monohydrate has been previously reported . Here we have examined corresponding changes in global gene expression and evaluated the potential for deriving a canonical expression profile that may serve as a diagnostic tool to detect the presence of TGF-β1 dependent activation in blood from advanced cancer subjects.
Figure 1: Schematic representation of samples and analysis performed on blood collected from normal subjects for ex-vivo analysis and cancer patients for measurement of global gene expression, pSMAD in PBMCs and plasma analytes. Blood collected in CTAD tubes was used to separate PBMCs for ex-vivo expression analysis and pSMAD measurement. Plasma from EDTA tubes was used for analysis of TGF-ß1, PF4 and multiple analytes using the Multi Analyte Immunoassay Panel (MIAP). The Tempus blood collection tube was used for whole blood microarray analysis from patient samples.
The expression pattern of PBMCs stimulated with TGF-β1 with and without LY2157299 monohydrate addition is represented in Figure 2. From left to right, the lanes (1-5) represent subgroups of samples with increasing TGF-β1 levels followed by a subgroup (lanes 6-8) where 5μM LY2157299 monohydrate was added 30 min prior to addition of the various levels of TGF-β1. Lane 9 represents baseline expression in normal PBMCs with 5μM LY2157299 monohydrate alone. Although samples without addition of the kinase inhibitor were exposed to TGF-β1 for 4 h instead of 3.5 h the effect of the additional 30 min on expression was insignificant. There were 375 differentially expressed probe sets corresponding to 325 genes upon TGF-β1 exposure (at all levels examined) that were reversed with LY2157299 monohydrate addition (LY-reversed). In addition 132 probe sets corresponding to 114 genes were differentially expressed with TGF-β1 addition they were not reversed with LY2157299 monohydrate and were not included in further analysis (data not shown). The expression index calculated for the LYreversed 375probe sets showed a clear difference with TGF-β1 added compared with vehicle alone. There was a 2-fold change in expression index median even at the lowest concentration of TGF-β1 (1pM) that was almost completely reversed with the inhibitor. Interestingly, both with cytokine and with and without inhibitor addition, the pattern of the change became more diverse, with increasing TGF-β1concentration as more genes were activated perhaps as a result of non-specific TGF- β1activity. This effect which underscores the sensitivity of TGF-β signaling and the specificity of the expression index measure is evident in the width of the box that represents the upper and lower quartile in the box plot, Figure 2, lanes 2, 5 and 8. At lowest level of TGF-β tested activation and inhibition with LY2157299 is very sensitive and specific respectively, Figure 2 lanes 2 and 6. In each case, when pre-incubated with the LY2157299 for 30 min the corresponding expression index of LY-reversed genes was reduced with TGF-β1 addition. As expected, baseline LY-reversed expression index values were lowest and least variable with LY2157299 addition suggesting control of endogenous signaling from platelet-derived cytokines (see vehicle control).
Figure 2: Distribution of the expression index values for 375 probe-sets in normal PBMC replicates treated with various levels of TGF-ß1 with and without 5 µm LY2157299 monohydrate As described in materials and methods Lane (1) Vehicle represents no treatment and lane (9) represents LY 2159299 without TGF-ß1 addition. The box within each range of data points represents the upper and lower median of replicate values for each treatment group and line within the box represents the median value. The vertical line extending from the box where available represents the upper and lower range of data points for each experimental condition.
Development of a clinically relevant signature for systemic activation of the TGF-βRI-dependent pathway
We examined the ex-vivo 325 LY-reversed gene set for association with a comprehensive gene set consisting of 2,749 genes compiled by Xu and Kapoun  that was derived from several cell-lines treated with TGF-β1. We found 113 common genes with a fold change range from 6.9 to -4.9 (data not shown). The expression of these genes was also altered in tumor specimens compared with corresponding normal tissue . We rationalized that perhaps the effect of TGF-β1 and altered expression in tissues might be similar in patient derived PBMCs where ligand was peripherally present. Since we had previously determined that ex-vivo PBMC activation with TGF-β1 was also accompanied with accumulation of pSMAD , we hypothesized that residual phosphorylation capacity of SMAD proteins in PBMCs collected from cancer patients when stimulated with TGF-β1 (ss) compared with unstimulated (us) would be reflected by the level of expression of TGF-β1 specific genes in the PBMC sample. We therefore examined the 325 LY-reversed gene-set for genes with significant association with response to ex-vivo pSMAD stimulation from patients with advanced cancer to see if there was a systemic signature of TGF-β1 activation.
The analytical measurement of ex-vivo pSMAD previously described  was derived from TGF-β1 stimulated patient PBMCs and expressed as a ratio of stimulated (ss) versus unstimulated (su) ELISA signal values. We have previously demonstrated  that the pSMAD (ss/su) ratio was a reflection of the residual PBMC SMAD phosphorylation capacity and was inversely associated with plasma TGF-β1 levels in these patients. Since we had collected whole blood from subjects and measured PBMC pSMAD stimulation, ex-vivo, using the same conditions as those used to derive the 325LY-reversed gene set, we were able to determined which subset from 325 LY-reversed genes also significantly correlated with plasma pSMAD (ss)/(su). We found that the expression index of 37 of 325 ex-vivo derived genes was significantly associated with patients’ PBMC pSMAD (ss/su) response ratio, p-value <0.0001 as shown in Figure 3 where the dotted lines represent the confidence intervals around the simulated straight line correlation. We next examined the 37-gene expression index for association with plasma TGF-β1 levels in patient samples and found that it was directly correlated with a p-value equal to 0.0014 (Figure 4). The statistically significant association of the ex-vivo derived 37-gene expression index with two separate known TGF-β pathway activation measurements, plasma TGF-β1 and PBMC pSMAD from clinical specimens further validates the 37-gene signature as clinical “surrogate” for TGF-β1 pathway activation.
Figure 3: Association of pSMAD2 (ratio of TGF-ß1 stimulated and unstimulated) in patient PBMCs and the 37-gene whole blood Expression Index: Negative Correlation plot of the ratio of PBMC derived phosphorylated SMAD2/3, stimulated (SS)/unstimulated (SU)with the ex-vivo derived 37-gene expression index from whole blood of 28 cancer patients. The solid line represents the simulated correlation and the dotted lines represent the confidence interval. The horizontal dotted line intersects the area where the correlation was highest.
Figure 4: Positive Association of PF4 adjusted plasma TGF-ß1 level with the exvivo derived 37-gene expression index from whole blood of 28 cancer patients. The simulated solid line represents the correlation and the surrounding dotted lines represent confidence interval. The horizontal dotted line intersects the area of highest correlation.
TGF-β1 dependent gene profile is a surrogate for pathway activation in clinical specimens
TGF-β responsive changes in systemic expression have been previously reported in PBMCs [8,30]. However, affected cells were examined after separation from whole blood and cultured without adding TGF-β raising the possibility of cellular re-equilibration and loss of the systemic effect resulting from chronic exposure to circulating TGF-β. We sought to examine whether there was evidence of systemic effect in patient blood by immediately accessing mRNA from affected cells upon immediate cell lysis and mRNA preservation. Since platelet-derived TGF-β can confound the affect of circulating TGF-β on WBC expression, we examined the effect of normal whole blood collected in a PaxGene tube where cells are stabilized in an RNA preserving buffer system without lysis and compared this with blood collected in a Tempus tube where cells are immediately lysed in a an RNA preservative containing lysis buffer. We then compared the distribution of the expression index values for the ex-vivo derived 325 gene LY-revered genes and the 37-gene pSMAD correlative TGF-β1 responsive signature among Tempus tube collected patient blood and normal blood collected in PaxGene and Tempus tubes. As shown in Figure 5, the overall expression index for normal blood was higher in the PaxGene tubes compared with Tempus indicating that perhaps increased baseline expression in TGF-β1 responsive genes within the period of time it took to collection and preserve cells in this tube type. When the expression signal values for 37 genes TGF-β1 responsive signature from 28 patient blood samples were analyzed together with ex vivo PBMC data from normal and stimulated samples the resulting hierarchical heat map showed three main clusters as shown in Figure 6. A cluster that separated with normal samples collected in Tempus tubes and ex vivo PBMC vehicle controls that included 11 patient samples which we have termed “pathway not activated”. A cluster of 7 patient samples and 3 ex vivo PBMC vehicle controls termed “partially activated” and cluster with 10 patient samples that included all of the ex vivo PBMC samples with TGF-β1 addition which was termed “pathway activated”. Interestingly, we observed that the TGF-β1 gene expression was always elevated with TGF-β1 stimulated normal PBMCs but only sporadically in patient samples (Figure 6).
Figure 5: Representation of the expression index values for 375 probe-sets annotated to 325 genes and the reduced 37-gene set in the whole blood from 28 cancer patients and 5 normal adults collected in PaxGene and Tempus RNA preservation tubes. Several replicated are represented in each experimental condition. The box within each range of data points represents the upper and lower median of replicate values for each treatment group and line within the box represents the median. The vertical line extending from the box where available represents the upper and lower range of data points for each experimental condition. The 37-gene subset provided a broader dynamic range of expression index values compared with 325-gene set.
Figure 6: Hierarchically clustered heat map of the 37-gene expression signals from ex-vivo TGF-ß1 stimulated PBMCs, vehicle controls (PBMCs with no addition) and whole blood from 28 cancer patients and normal subjects. Two major clusters were seen, a TGF-ßPI pathway activated block to the left with all the ex-vivo TGF-ß1 stimulated samples and a block to the right with majority of normal (Tempus) and vehicle controls designated as pathway not-activated. A central block with intermediate expression was designated partially activated based on signal intensities.
Analysis of cytokine stimulated serum markers using a Multianalyte Immunoassay Panel (MAIP)
Since a major function of TGF-β1 appears to be the regulation of chemokine expression, we also measured several peripheral proteins known to be co-regulated with TGF-β1. The values from the MAIP protein measurements were tested for correlation with the corresponding gene probe set signals present on the Affymetrix U133 gene chip and the results are depicted in Table 1. Several cytokines appeared to be significantly correlated, both positively and negatively in this analysis. Based on this finding, we asked whether the elevated TGF-β1 in plasma was triggering the TGF-βRI signaling pathway in patient PBMCs. Also given that the correlation of MAIP and gene expression might be related to the tumor histology and influence its development, we evaluated whether there were specific correlations between gene expression and certain proteins represented in the MAIP. Such a correlation between gene expression in the blood and plasma protein levels may help to find new markers of response to TGF-βRI inhibitors. Hence, the correlation between ENA-78 and CXCL5/215101_s_at, VEGF and VEGFA/210512_s_at, as well as the correlation between MMP-9 and MMP9/203936_s_at, are intriguing observations for evaluation as potential predictive biomarkers of response in patient plasma.
|Protein||Gene Symbol/Probeset||Spearman R||p value|
|Prostate Acid Phosphatase||ACPP/204393_s_at||0.451||0.031|
Table 1: Association of secreted plasma protein levels with their corresponding gene expression signals in patient whole blood (cut-off p<0.1, unadjusted for multiplicity).
Recent progress in molecular oncology has provided increased optimism that empirical use of anticancer agents will soon be transformed into a practice of “personalized” cancer medicine, where molecular profile together with biomarker guided assessment of cellular content in patient’s tumor, rather than tissue type and population-based risk, will determine the choice of therapeutic intervention [32,33]. It has also become clear, that tumor response is guided by cellular molecular circuitry, in which these changes occur and are further facilitated by the presence of cytokines in the tumor microenvironment. This concept of context-dependent tumor addiction has important implications for the design of molecularly targeted therapies .
Although TGF-β signaling is well known to play a key role in tumor progression and as a soluble mediator of immune regulation, the effects of the presence of this family of cytokines on peripheral cells in patient blood has not been clearly dissected. In cancer patients, secreted proteins are known to be released from various locations, such as tumor tissue, circulating blood cells and also from normal tissues, such as endothelial cells making it difficult to assess their functional effect in patients chronically exposed to circulating cytokines. Generally, only when tumors secrete either highly specific proteins or at high levels, does the assessment of their presence in the circulation become a meaningful laboratory measure for disease management. Paradoxically, TGF-β isoforms are well known to act both as tumor suppressors and as tumor promoters and their subtle presence in the cellular environment can produce diverse outcomes including immune suppression and cellular senescence . It has also been suggested that chronic exposure of tumors to TGF-β might have an effect unlike acute exposure related expression changes derived from ex-vivo cell line based experiments and as such may require the analysis of both ex-vivo and patient derived specimens for biomarkers of TGF-β response . Since whole blood is readily accessible compared with tumor biopsies especially in advanced cancer patients, it is an ideal sample for assessment of TGF-β dependent response and measurement of pharmacodynamic effect during the course of drug development. Additionally, PBMCs typically have a higher cellular turnover than resident tumor cells and thus a subset of the newly generated cells would be expected to be relatively naïve to the effects of continuous TGF-β1 exposure unlike the primary tumor. These factors make whole blood an appropriate “surrogate” substrate for biomarkers of a systemic effect of TGF-β1 induced RNA expression in cancer patients.
We used an ex-vivo assay with normal PBMCs suspended in a buffered matrix where the in-situ transcriptional capacity was fixed and then evaluated TGF-β1-dependent mRNA expression as biomarker of its signaling activity. We then tested the PBMC based assay to determine the effect of addition of LY2157299 monohydrate, a small molecule inhibitor of the TGF-βRI. First, we found that a specific set of transcripts was differentially regulated with and without addition of LY2157299 monohydrate prior to stimulation with TGF-β1. Second, the pattern of change with inhibitor added appeared to mirror the effect seen when increasing levels of TGF-β1 suggesting that perhaps the broader response seen at higher levels of TGF-β1 (100 pM), resulted from significant saturation of cellular receptors, overwhelming specific TGF-β1 dependent transcriptional capacity and allowing nonspecific TGF-β1 derived expression causing a dilution of the specific response seen when the ligand was added at lower levels. This can become important when detecting TGF-β1gene signaling where there is too much circulating ligand, consequently making it harder to detect changes with LY2157299 treatment. Third, the ex-vivo assay was able to detect pre-analytical variability, due to suboptimal sample handling, known to set off release of platelet derived cytokines producing higher background expression prior to addition of exogenous TGF-β1 and the resulting loss of LY2157999 monohydrate dependent reversal. This is also evident in the difference in the expression between the control (vehicle) and the inhibitor-only lanes Figure 2. Gene ontology analysis of the two gene profiles indicated that LY2157299 monohydrate specifically and significantly affected ligand induced genes in both directions where as changes in the non reversible gene set were less significant and nonspecific (data not shown). The non-reversed genes were annotated to know genes with cellular “housekeeping” function such as ACTB and GAPDH and thus, the resulting expression index might reflect cellular transcription capacity.
There are conflicting reports of the association of TGF-β1 dependent expression in disease tissue, circulating levels and its effect on peripheral cells. For example, it was recently reported that in multiple sclerosis patients, while there was decreased activation of TGF-β1 pathway signaling in whole blood and a reduction in the level of regulated genes, peripheral cells from these subjects (ex-vivo) retained an intact exogenous TGF-β1 response . Our data similarly shows a relationship between pSMAD response and the amount of TGF-β1 stimulation. Since the tSMAD for paired unstimulated and stimulated PBMCs is strongly correlated, adjusting pSMAD for tSMAD in these samples was not necessary. However, since TGF-β1 is a potent signaling molecule, the magnitude of transcriptional regulation from chronic exposure to TGF-β1 is likely to be tightly controlled due to pleiotropic effects of this ligand. This may also explain the lack of a stronger correlation of the ex-vivo derived, LY2157299 monohydrate reversible expression index with peripheral TGF-β1 levels in patient samples. However, when the reversible expression index was adjusted for PF4, there was good correlation with plasma TGF-β1 levels, again suggesting that perhaps a naïve subset of peripheral cells was TGF-β1 responsive. Also, the correlation between pSMAD levels and gene expression suggests that the systemic 37-gene TGF-β1 responsive signature will be particularly important in conditions where TGF-β signaling is partially activated. Additional support for the association of the LY2157299 monohydrate reversible expression index with TGF-β1 was provided with the detection of co-regulated plasma derived proteins in patients. We found good overall correlation with multiple known co-regulated markers of TGF-β1 activity listed in Table 1. Specifically ENA-78, the protein product of CXCL5, MMP9 and VEGF are previously reported in animal models and tumor tissue to be tightly associated and prognostic of disease outcome in breast cancer . Similarly, we found other known compensatory proteins such as EPO, TNF-β and IL-3 to be negatively correlated with plasma TGF-β. This observation of correlation between blood-based proteins and gene expression in PBMCs further underscores that TGF-β signaling is closely regulated in circulation. Additionally we noticed that based on microarray signal intensities the expression level of TGFB1 was higher in TGF- β1stimulated PBMCs but only rarely in patient blood. It is worth noting that while the ex-vivo system provided a measure of change in expression that was clearly associated with known cytokine dependent changes in patient blood, the statistically significant observations reported here and the 37-gene TGF-β1 expression signature came from 28 patient sample of which 19 were diagnosed with breast cancer. It is likely that length and type of standard of care treatment in patients may influence global expression in patient PBMCs. The unequivocal association of the 37-gene expression index as a “surrogate” gene expression signature for TGF-β1 activation in tumors will require addition confirmation with robust clinical sampling. As drug development shifts towards target-based therapies  and the molecular assessment of disease subtypes becomes common place, the methodology described here in combination with specific inhibitors may become useful in disease diagnosis and the selection of appropriate treatment. Since TGF-βRI signaling is known to be elevated in patients with bone metastasis , this study illustrates an approach using an expression profile from peripheral cells to potentially identify patients with an active TGF-β signaling pathway supporting the disease condition [40,41].
The authors thank colleagues in the Laboratory of Experimental Medicine for valuable technical input. Special thanks to Crystal Dotson for technical support and Lalit Verma for providing access to biological network maps and gene ontology tools. We also thank Dr. Smita Agarwal, for writing assistance provided with the manuscript.