Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

Immunological microenvironment in prostate cancer: high mast cell densities are associated with favorable tumor characteristics and good prognosis

BACKGROUND: Number of intratumoral mast cells predicts survival in various cancers. The prognostic significance of such mast cells in surgically treated prostate cancer is unknown. METHODS: Mast cell densities were determined in prostate cancer samples of more than 2,300 hormone-naïve patients using a tissue microarray format in correlation with clinical follow-up data. Mast cells were visualized immunohistochemically (c-kit). All patients were homogeneously treated by radical prostatectomy at a single institution. RESULTS: Mast cells were present in 95.9% of the tumor samples. Median mast cell number on the tissue spot was 9 (range: 0-90; median density: 31 mast cells/mm(2)). High mast cell densities were significantly associated with more favorable tumors having lower preoperative prostate-specific antigen (P = 0.0021), Gleason score (P < 0.0001) and tumor stage (P < 0.0001) than tumors with low mast cell densities. Prostate-specific antigen recurrence-free survival significantly (P = 0.0001) decreased with decline of mast cell density showing poorest outcome for patients without intratumoral mast cells. In multivariate analysis mast cell density narrowly missed to add independent prognostic information (P = 0.0815) for prostate-specific antigen recurrence. CONCLUSION: High intratumoral mast cell density is associated with favorable tumor characteristics and good prognosis in prostate cancer. This finding is consistent with a role of mast cells in the immunological host-defense reaction on prostate cancer. Triggering mast cell activity might expand immunotherapeutic strategies in prostate cancer.

Erythropoietin or Darbepoetin for patients with cancer--meta-analysis based on individual patient data

BACKGROUND: Erythropoiesis-stimulating agents (ESAs) reduce anemia in cancer patients and may improve quality of life, but there are concerns that ESAs might increase mortality. OBJECTIVES: Our objectives were to examine the effect of ESAs and identify factors that modify the effects of ESAs on overall survival, progression free survival, thromboembolic and cardiovascular events as well as need for transfusions and other important safety and efficacy outcomes in cancer patients. SEARCH STRATEGY: We searched the Cochrane Library, Medline, Embase and conference proceedings for eligible trials. Manufacturers of ESAs were contacted to identify additional trials. SELECTION CRITERIA: We included randomized controlled trials comparing epoetin or darbepoetin plus red blood cell transfusions (as necessary) versus red blood cell transfusions (as necessary) alone, to prevent or treat anemia in adult or pediatric cancer patients with or without concurrent antineoplastic therapy. DATA COLLECTION AND ANALYSIS: We performed a meta-analysis of randomized controlled trials comparing epoetin alpha, epoetin beta or darbepoetin alpha plus red blood cell transfusions versus transfusion alone, for prophylaxis or therapy of anemia while or after receiving anti-cancer treatment. Patient-level data were obtained and analyzed by independent statisticians at two academic departments, using fixed-effects and random-effects meta-analysis. Analyses were according to the intention-to-treat principle. Primary endpoints were on study mortality and overall survival during the longest available follow-up, regardless of anticancer treatment, and in patients receiving chemotherapy. Tests for interactions were used to identify differences in effects of ESAs on mortality across pre-specified subgroups. The present review reports only the results for the primary endpoint. MAIN RESULTS: A total of 13933 cancer patients from 53 trials were analyzed, 1530 patients died on-study and 4993 overall. ESAs increased on study mortality (combined hazard ratio [cHR] 1.17; 95% CI 1.06-1.30) and worsened overall survival (cHR 1.06; 95% CI 1.00-1.12), with little heterogeneity between trials (I(2) 0%, p=0.87 and I(2) 7.1%, p=0.33, respectively). Thirty-eight trials enrolled 10441 patients receiving chemotherapy. The cHR for on study mortality was 1.10 (95% CI 0.98-1.24) and 1.04; 95% CI 0.97-1.11) for overall survival. There was little evidence for a difference between trials of patients receiving different cancer treatments (P for interaction=0.42). AUTHORS' CONCLUSIONS: ESA treatment in cancer patients increased on study mortality and worsened overall survival. For patients undergoing chemotherapy the increase was less pronounced, but an adverse effect could not be excluded.

A very rare cancer in Down syndrome: medulloblastoma. Epidemiological data from 13 countries

Persons with Down syndrome (DS) uniquely have an increased frequency of leukemias but a decreased total frequency of solid tumors. The distribution and frequency of specific types of brain tumors have never been studied in DS. We evaluated the frequency of primary neural cell embryonal tumors and gliomas in a large international data set. The observed number of children with DS having a medulloblastoma, central nervous system primitive neuroectodermal tumor (CNS-PNET) or glial tumor was compared to the expected number. Data were collected from cancer registries or brain tumor registries in 13 countries of Europe, America, Asia and Oceania. The number of DS children with each category of tumor was treated as a Poisson variable with mean equal to 0.000884 times the total number of registrations in that category. Among 8,043 neural cell embryonal tumors (6,882 medulloblastomas and 1,161 CNS-PNETs), only one patient with medulloblastoma had DS, while 7.11 children in total and 6.08 with medulloblastoma were expected to have DS. (p 0.016 and 0.0066 respectively). Among 13,797 children with glioma, 10 had DS, whereas 12.2 were expected. Children with DS appear to be specifically protected against primary neural cell embryonal tumors of the CNS, whereas gliomas occur at the same frequency as in the general population. A similar protection against neuroblastoma, the principal extracranial neural cell embryonal tumor, has been observed in children with DS. Additional genetic material on the supernumerary chromosome 21 may protect against embryonal neural cell tumor development.

Distinct microRNA Expression Profile in Prostate Cancer Patients with Early Clinical Failure and the Impact of let-7 as Prognostic Marker in High-Risk Prostate Cancer

Background The identification of additional prognostic markers to improve risk stratification and to avoid overtreatment is one of the most urgent clinical needs in prostate cancer (PCa). MicroRNAs, being important regulators of gene expression, are promising biomarkers in various cancer entities, though the impact as prognostic predictors in PCa is poorly understood. The aim of this study was to identify specific miRNAs as potential prognostic markers in high-risk PCa and to validate their clinical impact. Methodology and Principal Findings We performed miRNA-microarray analysis in a high-risk PCa study group selected by their clinical outcome (clinical progression free survival (CPFS) vs. clinical failure (CF)). We identified seven candidate miRNAs (let-7a/b/c, miR-515-3p/5p, -181b, -146b, and -361) that showed differential expression between both groups. Further qRT-PCR analysis revealed down-regulation of members of the let-7 family in the majority of a large, well-characterized high-risk PCa cohort (n = 98). Expression of let-7a/b/and -c was correlated to clinical outcome parameters of this group. While let-7a showed no association or correlation with clinical relevant data, let-7b and let-7c were associated with CF in PCa patients and functioned partially as independent prognostic marker. Validation of the data using an independent high-risk study cohort revealed that let-7b, but not let-7c, has impact as an independent prognostic marker for BCR and CF. Furthermore, we identified HMGA1, a non-histone protein, as a new target of let-7b and found correlation of let-7b down-regulation with HMGA1 over-expression in primary PCa samples. Conclusion Our findings define a distinct miRNA expression profile in PCa cases with early CF and identified let-7b as prognostic biomarker in high-risk PCa. This study highlights the importance of let-7b as tumor suppressor miRNA in high-risk PCa and presents a basis to improve individual therapy for high-risk PCa patients.

Loss of TRAIL-receptors is a recurrent feature in pancreatic cancer and determines the prognosis of patients with no nodal metastasis after surgery

INTRODUCTION Agonistic antibodies targeting TRAIL-receptors 1 and 2 (TRAIL-R1 and TRAIL-R2) are being developed as a novel therapeutic approach in cancer therapy including pancreatic cancer. However, the cellular distribution of these receptors in primary pancreatic cancer samples has not been sufficiently investigated and no study has yet addressed the issue of their prognostic significance in this tumor entity. AIMS AND METHODS Applying tissue microarray (TMA) analysis, we performed an immunohistochemical assessment of TRAIL-receptors in surgical samples from 84 consecutive patients affected by pancreatic adenocarcinoma and in 26 additional selected specimens from patients with no lymph nodes metastasis at the time of surgery. The prognostic significance of membrane staining and staining intensity for TRAIL-receptors was evaluated. RESULTS The fraction of pancreatic cancer samples with positive membrane staining for TRAIL-R1 and TRAIL-R2 was lower than that of cells from surrounding non-tumor tissues (TRAIL-R1: p<0.001, TRAIL-R2: p = 0.006). In addition, subgroup analyses showed that loss of membrane staining for TRAIL-R2 was associated with poorer prognosis in patients without nodal metastases (multivariate Cox regression analysis, Hazard Ratio: 0.44 [95% confidence interval: 0.22-0.87]; p = 0.019). In contrast, analysis of decoy receptors TRAIL-R3 and -R4 in tumor samples showed an exclusively cytoplasmatic staining pattern and no prognostic relevance. CONCLUSION This is a first report on the prognostic significance of TRAIL-receptors expression in pancreatic cancer showing that TRAIL-R2 might represent a prognostic marker for patients with early stage disease. In addition, our data suggest that loss of membrane-bound TRAIL-receptors could represent a molecular mechanism for therapeutic failure upon administration of TRAIL-receptors-targeting antibodies in pancreatic cancer. This hypothesis should be evaluated in future clinical trials.

Investigation of IL23 (p19, p40) and IL23R highlights nuclear IL23p19 expression as a marker of indolent tumor features and favorable outcome in colorectal cancer patients

Objective: IL23 is involved in chronic inflammation but its role in cancer progression is not fully elucidated. Here we characterize IL23 subunits p40, p19 and IL23 receptor (IL23R) in the normal-adenoma-carcinomametastasis cascade of colorectal cancers and their relationship to clinicopathological and outcome data. Method: Immunohistochemistry for IL23R, IL12p40, IL23 and IL23p19 (monoclonal) was performed on a multi-punch tissue microarray (n=213 patients). Expression differences between normal-adenomas-cancerslymph nodes were evaluated. Correlation with clinicopathological and outcome data was undertaken. Results were validated on an independent cohort (n=341 patients). Results: An increased expression from normal-adenoma-cancer was observed (p<0.0001; all) followed by a marked reduction in lymph nodes (p<0.0001; all). Cytoplasmic and/or membranous staining of all markers was unrelated to outcome. Nuclear IL23p19 staining occurred in 23.1%and was associated with smaller tumor diameter (p=0.0333), early pT (p=0.0213), early TNM (p=0.0186), absence of vascular (p=0.0124) and lymphatic invasion (p=0.01493) and favorable survival (univariate (p=0.014) and multivariable (p=0.0321) analysis). All IL23p19 positive patients were free of distant metastasis (p=0.0146). Survival and metastasis results could be validated in Cohort 2. Conclusion: The presence of nuclear IL23p19 is related to indolent tumor features and favorable outcome supporting a more ‘protective’ role of this protein in colorectal cancer progression

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

DOWNREGULATION OF PAX2 SUPPRESSES OVARIAN CANCER CELL GROWTH

PAX2 is one of nine PAX genes regulating tissue development and cellular differentiation in embryos. PAX2 promotes cell proliferation, oncogenic transformation, cell-lineage specification, migration, and survival. Unattenuated PAX2 has been found in several cancer types. We therefore sought to elucidate the role of PAX2 in ovarian carcinomas. We found that PAX2 was expressed in low-grade serous, clear cell, endometrioid and mucinous cell ovarian carcinomas, which are relatively chemoresistant compared to high grade serous ovarian carcinomas. Four ovarian cancer cell lines, RMUGL (mucinous), TOV21G (clear cell), MDAH-2774 (endometrioid) and IGROV1 (endometrioid), which express high-levels of PAX2, were used to study the function of PAX2. Lentiviral shRNAs targeting PAX2 were used to knock down PAX2 expression in these cell lines. Cellular proliferation and motility assays subsequently showed that PAX2 stable knockdown had slower growth and migration rates. Microarray gene expression profile analysis further identified genes that were affected by PAX2 including the tumor suppressor gene G0S2. Reverse phase protein array (RPPA) data showed that PAX2 knockdown affected several genes that are involved in apoptosis, which supports the fact that downregulation of PAX2 in PAX2-expressing ovarian cancer cells inhibits cell growth. We hypothesize that this growth inhibition is due to upregulation of the tumor suppressor gene G0S2 via induction of apoptosis. PAX2 represents a potential therapeutic target for chemoresistant PAX2-expressing ovarian carcinomas.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

LMW-E MEDIATES MAMMARY TUMORIGENESIS BY DEREGULATING ACINAR MORPHOGENESIS & GENERATING CANCER STEM CELLS

Cyclin E is the regulatory subunit of the cyclin E/CDK2 complex that mediates the G1-S phase transition. N-terminal cleavage of cyclin E by elastase in breast cancer generates two low molecular weight (LMW) isoforms that exhibit both enhanced kinase activity and resistance to p21 and p27 inhibition compared to fulllength cyclin E. Clinically, approximately 27% of breast cancer patients overexpress LMW-E and associate with poor survival. Therefore, we hypothesize that LMW-E disrupts normal mammary acinar morphogenesis and serves as the initial route into breast tumor development. We first demonstrate that LMW-E overexpression in non-tumorigenic hMECs is sufficient to induce tumor formation in athymic mice significantly more than overexpression of full-length cyclin E and requires CDK2- associated kinase activity. Further in vivo passaging of these tumors augments LMW-E expression and tumorigenic potential. When subjected to acinar morphogenesis in vitro, LMW-E mediates significant morphological disruption by generating hyperproliferative and multi-acinar complexes. Proteomic analysis of patient tissues and tumor cells with high LMW-E expression reveals that the activation of the b-Raf-ERK1/2-mTOR pathway in concert with high LMW-E expression predicts poor patient survival. Combination treatment using roscovitine (CDK inhibitor) plus either rapamycin (mTOR inhibitor) or sorafenib (b-raf inhibitor) effectively prevented aberrant acinar formation in LMW-E-expressing cells by inducing the G1/S cell cycle arrest. In addition, the LMW-E-expressing tumor cells exhibit phenotypes characteristic of the EMT and enhanced cellular invasiveness. These tumor cells also enrich for cells with CSC phenotypes such as increased CD44hi/CD24lo population, enhanced mammosphere formation, and upregulation of ALDH expression and enzymatic activity. Furthermore, the CD44hi/CD24lo population also shows positive correlation with LMW-E expression in both the tumor cell line model and breast cancer patient samples (p<0.0001 & p=0.0435, respectively). Combination treatment using doxorubicin and salinomycin demonstrates synergistic cytotoxic effects in cells with LMW-E expression but not in those with full-length cyclin E expression. Finally, ProtoArray microarray identifies Hbo1 as a novel substrate of the cyclin E/CDK2 complex and its overexpression results in enrichment for CSCs. Collectively, these data emphasize the strong oncogenic potential of LMW-E in mammary tumorigenesis and suggest possible therapeutic strategies to treat breast cancer patients with high LMW-E expression.

HIV-related malignancies: A community-based study using a linkage of cancer registry and HIV registry data

Background. Increased incidence of cancer is documented in immunosuppressed transplant patients. Likewise, as survival increases for persons infected with the Human Immunodeficiency Virus (HIV), we expect their incidence of cancer to increase. The objective of this study was to examine the current gender specific spectrum of cancer in an HIV infected cohort (especially malignancies not currently associated with Acquired Immunodeficiency Syndrome (AIDS)) in relation to the general population.^ Methods. Cancer incidence data was collected for residents of Harris County, Texas who were diagnosed with a malignancy between 1975 and 1994. This data was linked to HIV/AIDS registry data to identify malignancies in an HIV infected cohort of 14,986 persons. A standardized incidence ratio (SIR) analysis was used to compare incidence of cancer in this cohort to that in the general population. Risk factors such as mode of HIV infection, age, race and gender, were evaluated for contribution to the development of cancer within the HIV cohort, using Cox regression techniques.^ Findings. Of those in the HIV infected cohort, 2289 persons (15%) were identified as having one or more malignancies. The linkage identified 29.5% of these malignancies (males 28.7% females 60.9%). HIV infected men and women had incidences of cancer that were 16.7 (16.1, 17.3) and 2.9 (2.3, 3.7) times that expected for the general population of Harris County, Texas, adjusting for age. Significant SIR's were observed for the AIDS-defining malignancies of Kaposi's sarcoma, non-Hodgkin's lymphoma, primary lymphoma of the brain and cancer of the cervix. Additionally, significant SIR's for non-melanotic skin cancer in males, 6.9 (4.8, 9.5) and colon cancer in females, 4.0 (1.1, 10.2) were detected. Among the HIV infected cohort, race/ethnicity of White (relative risk 2.4 with 95% confidence intervals 2.0, 2.8) or Spanish Surname, 2.2 (1.9, 2.7) and an infection route of male to male sex, with, 3.0 (1.9, 4.9) or without, 3.4 (2.1, 5.5) intravenous drug use, increased the risk of having a diagnosis of an incident cancer.^ Interpretation. There appears to be an increased risk of developing cancer if infected with the HIV. In addition to the malignancies routinely associated with HIV infection, there appears to be an increased risk of being diagnosed with non-melanotic skin cancer in males and colon cancer in females. ^

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Effects of erythropoiesis-stimulating agents on fatigue- and anaemia-related symptoms in cancer patients: systematic review and meta-analyses of published and unpublished data.

Background:Erythropoiesis-stimulating agents (ESAs) reduce the need for red blood cell transfusions; however, they increase the risk of thromboembolic events and mortality. The impact of ESAs on quality of life (QoL) is controversial and led to different recommendations of medical societies and authorities in the USA and Europe. We aimed to critically evaluate and quantify the effects of ESAs on QoL in cancer patients.Methods:We included data from randomised controlled trials (RCTs) on the effects of ESAs on QoL in cancer patients. Randomised controlled trials were identified by searching electronic data bases and other sources up to January 2011. To reduce publication and outcome reporting biases, we included unreported results from clinical study reports. We conducted meta-analyses on fatigue- and anaemia-related symptoms measured with the Functional Assessment of Cancer Therapy-Fatigue (FACT-F) and FACT-Anaemia (FACT-An) subscales (primary outcomes) or other validated instruments.Results:We identified 58 eligible RCTs. Clinical study reports were available for 27% (4 out of 15) of the investigator-initiated trials and 95% (41 out of 43) of the industry-initiated trials. We excluded 21 RTCs as we could not use their QoL data for meta-analyses, either because of incomplete reporting (17 RCTs) or because of premature closure of the trial (4 RCTs). We included 37 RCTs with 10 581 patients; 21 RCTs were placebo controlled. Chemotherapy was given in 27 of the 37 RCTs. The median baseline haemoglobin (Hb) level was 10.1 g dl(-1); in 8 studies ESAs were stopped at Hb levels below 13 g dl(-1) and in 27 above 13 g dl(-1). For FACT-F, the mean difference (MD) was 2.41 (95% confidence interval (95% CI) 1.39-3.43; P<0.0001; 23 studies, n=6108) in all cancer patients and 2.81 (95% CI 1.73-3.90; P<0.0001; 19 RCTs, n=4697) in patients receiving chemotherapy, which was below the threshold (⩾3) for a clinically important difference (CID). Erythropoiesis-stimulating agents had a positive effect on anaemia-related symptoms (MD 4.09; 95% CI 2.37-5.80; P=0.001; 14 studies, n=2765) in all cancer patients and 4.50 (95% CI 2.55-6.45; P<0.0001; 11 RCTs, n=2436) in patients receiving chemotherapy, which was above the threshold (⩾4) for a CID. Of note, this effect persisted when we restricted the analysis to placebo-controlled RCTs in patients receiving chemotherapy. There was some evidence that the MDs for FACT-F were above the threshold for a CID in RCTs including cancer patients receiving chemotherapy with Hb levels below 12 g dl(-1) at baseline and in RCTs stopping ESAs at Hb levels above 13 g dl(-1). However, these findings for FACT-F were not confirmed when we restricted the analysis to placebo-controlled RCTs in patients receiving chemotherapy.Conclusions:In cancer patients, particularly those receiving chemotherapy, we found that ESAs provide a small but clinically important improvement in anaemia-related symptoms (FACT-An). For fatigue-related symptoms (FACT-F), the overall effect did not reach the threshold for a CID.British Journal of Cancer advance online publication, 17 April 2014; doi:10.1038/bjc.2014.171 www.bjcancer.com.