Gene expression signature of Gleason score is associated with prostate cancer outcomes in a radical prostatectomy cohort.

Prostate cancer (PCa) is a leading cause of cancer-related mortality worldwide. Gleason score (GS) is one of the best predictors of PCa aggressiveness, but additional tumor biomarkers may improve its prognostic accuracy. We developed a gene expression signature of GS to enhance the prediction of PCa outcomes. Elastic net was used to construct a gene expression signature by contrasting GS 8-10 vs. ≤6 tumors in The Cancer Genome Atlas (TCGA) dataset. The constructed signature was then evaluated for its ability to predict recurrence and metastatic-lethal (ML) progression in a Fred Hutchinson (FH) patient cohort (N=408; NRecurrence=109; NMLprogression=27). The expression signature included transcripts representing 49 genes. In the FH cohort, a 25% increase in the signature was associated with a hazard ratio (HR) of 1.51 (P=2.7×10-5) for recurrence. The signature's area under the curve (AUC) for predicting recurrence and ML progression was 0.68 and 0.76, respectively. Compared to a model with age at diagnosis, pathological stage and GS, the gene expression signature improved the AUC for recurrence (3%) and ML progression (6%). Higher levels of the signature were associated with increased expression of genes in cell cycle-related pathways and decreased expression of genes in androgen response, estrogen response, oxidative phosphorylation, and apoptosis. This gene expression signature based on GS may improve the prediction of recurrence as well as ML progression in PCa patients after radical prostatectomy.

Gene expression panel predicts metastatic-lethal prostate cancer outcomes in men diagnosed with clinically localized prostate cancer.

Prognostic biomarkers are needed to distinguish patients with clinically localized prostate cancer (PCa) who are at high risk of metastatic progression. The tumor transcriptome may reveal its aggressiveness potential and have utility for predicting adverse patient outcomes. Genomewide gene expression levels were measured in primary tumor samples of 383 patients in a population-based discovery cohort, and from an independent clinical validation dataset of 78 patients. Patients were followed for ≥ 5 years after radical prostatectomy to ascertain outcomes. Area under the receiver-operating characteristic curve (AUC), partial AUC (pAUC, 95% specificity), and P-value criteria were used to detect and validate the differentially expressed transcripts. Twenty-three differentially expressed transcripts in patients with metastatic-lethal compared with nonrecurrent PCa were validated (P < 0.05; false discovery rate < 0.20) in the independent dataset. The addition of each validated transcript to a model with Gleason score showed that 17 transcripts significantly improved the AUC (range: 0.83-0.88; all P-values < 0.05). These differentially expressed mRNAs represent genes with diverse cellular functions related to tumor aggressiveness. This study validated 23 gene transcripts for predicting metastatic-lethal PCa in patients surgically treated for clinically localized disease. Several of these mRNA biomarkers have clinical potential for identifying the subset of PCa patients with more aggressive tumors who would benefit from closer monitoring and adjuvant therapy.

Bevacizumab May Differentially Improve Ovarian Cancer Outcome in Patients with Proliferative and Mesenchymal Molecular Subtypes.

Purpose: Recent progress in understanding the molecular biology of epithelial ovarian cancer has not yet translated into individualized treatment for these women or improvements in their disease outcome. Gene expression has been utilized to identify distinct molecular subtypes, but there have been no reports investigating whether or not molecular subtyping is predictive of response to bevacizumab in ovarian cancer.Experimental Design: DASL gene expression arrays were performed on FFPE tissue from patients enrolled on the ICON7 trial. Patients were stratified into four TCGA molecular subtypes. Associations between molecular subtype and the efficacy of randomly assigned therapy with bevacizumab were assessed.Results: Molecular subtypes were assigned as follows: 122 immunoreactive (34%), 96 proliferative (27%), 73 differentiated (20%), and 68 mesenchymal (19%). In univariate analysis patients with tumors of proliferative subtype obtained the greatest benefit from bevacizumab with a median PFS improvement of 10.1 months [HR, 0.55 (95% CI, 0.34-0.90), P = 0.016]. For the mesenchymal subtype, bevacizumab conferred a nonsignificant improvement in PFS of 8.2 months [HR 0.78 (95% CI, 0.44-1.40), P = 0.41]. Bevacizumab conferred modest improvements in PFS for patients with immunoreactive subtype (3.8 months; P = 0.08) or differentiated subtype (3.7 months; P = 0.61). Multivariate analysis demonstrated significant PFS improvement in proliferative subtype patients only [HR, 0.45 (95% CI, 0.27-0.74), P = 0.0015].Conclusions: Ovarian carcinoma molecular subtypes with the poorest survival (proliferative and mesenchymal) derive a comparably greater benefit from treatment that includes bevacizumab. Validation of our findings in an independent cohort could enable the use of bevacizumab for those patients most likely to benefit, thereby reducing side effects and healthcare cost. Clin Cancer Res; 23(14); 3794-801. ©2017 AACR.

Novel Gene Expression Signature Predictive of Clinical Recurrence After Radical Prostatectomy in Early Stage Prostate Cancer Patients.

BACKGROUND: Current clinical tools have limited accuracy in differentiating patients with localized prostate cancer who are at risk of recurrence from patients with indolent disease. We aimed to identify a gene expression signature that jointly with clinical variables could improve upon the prediction of clinical recurrence after RP for patients with stage T2 PCa. METHODS: The study population includes consented patients who underwent a radical retropubic prostatectomy (RP) and bilateral pelvic lymph node dissection at the University of Southern California in the PSA-era (1988-2008). We used a nested case-control study of 187 organ-confined patients (pT2N0M0): 154 with no recurrence ("controls") and 33 with clinical recurrence ("cases"). RNA was obtained from laser capture microdissected malignant glands representative of the overall Gleason score of each patient. Whole genome gene expression profiles (29,000 transcripts) were obtained using the Whole Genome DASL HT platform (Illumina, Inc). A gene expression signature of PCa clinical recurrence was identified using stability selection with elastic net regularized logistic regression. Three existing datasets generated with the Affymetrix Human Exon 1.0ST array were used for validation: Mayo Clinic (MC, n = 545), Memorial Sloan Kettering Cancer Center (SKCC, n = 150), and Erasmus Medical Center (EMC, n = 48). The areas under the ROC curve (AUCs) were obtained using repeated fivefold cross-validation. RESULTS: A 28-gene expression signature was identified that jointly with key clinical variables (age, Gleason score, pre-operative PSA level, and operation year) was predictive of clinical recurrence (AUC of clinical variables only was 0.67, AUC of clinical variables, and 28-gene signature was 0.99). The AUC of this gene signature fitted in each of the external datasets jointly with clinical variables was 0.75 (0.72-0.77) (MC), 0.90 (0.86-0.94) (MSKCC), and 0.82 (0.74-0.91) (EMC), whereas the AUC for clinical variables only in each dataset was 0.72 (0.70-0.74), 0.86 (0.82-0.91), and 0.76 (0.67-0.85), respectively. CONCLUSIONS: We report a novel gene-expression based classifier identified using agnostic approaches from whole genome expression profiles that can improve upon the accuracy of clinical indicators to stratify early stage localized patients at risk of clinical recurrence after RP. Prostate 76:1239-1256, 2016. © 2016 Wiley Periodicals, Inc.

Genetic suppression reveals DNA repair-independent antagonism between BRCA1 and COBRA1 in mammary gland development.

The breast cancer susceptibility gene BRCA1 is well known for its function in double-strand break (DSB) DNA repair. While BRCA1 is also implicated in transcriptional regulation, the physiological significance remains unclear. COBRA1 (also known as NELF-B) is a BRCA1-binding protein that regulates RNA polymerase II (RNAPII) pausing and transcription elongation. Here we interrogate functional interaction between BRCA1 and COBRA1 during mouse mammary gland development. Tissue-specific deletion of Cobra1 reduces mammary epithelial compartments and blocks ductal morphogenesis, alveologenesis and lactogenesis, demonstrating a pivotal role of COBRA1 in adult tissue development. Remarkably, these developmental deficiencies due to Cobra1 knockout are largely rescued by additional loss of full-length Brca1. Furthermore, Brca1/Cobra1 double knockout restores developmental transcription at puberty, alters luminal epithelial homoeostasis, yet remains deficient in homologous recombination-based DSB repair. Thus our genetic suppression analysis uncovers a previously unappreciated, DNA repair-independent function of BRCA1 in antagonizing COBRA1-dependent transcription programme during mammary gland development.

Expression of cell cycle-regulated genes and prostate cancer prognosis in a population-based cohort.

BACKGROUND: Prostate cancer (PCa) is clinically and biologically heterogeneous, making it difficult to predict at detection whether it will take an indolent or aggressive disease course. Cell cycle-regulated genes may be more highly expressed in actively dividing cells, with transcript levels reflecting tumor growth rate. Here, we evaluated expression of cell cycle genes in relation to PCa outcomes in a population-based cohort. METHODS: Gene expression data were generated from tumor tissues obtained at radical prostatectomy for 383 population-based patients (12.3-years average follow-up). The overall mean and individual transcript levels of 30 selected cell cycle genes was compared between patients with no evidence of recurrence (73%) and those who recurred (27%) or died (7%) from PCa. RESULTS: The multivariate adjusted hazard ratio (HR) for a change from the 25th to 75th percentile of mean gene expression level (range 8.02-10.05) was 1.25 (95%CI 0.96-1.63; P = 0.10) for PCa recurrence risk, and did not vary substantially by Gleason score, TMPRSS2-ERG fusion status, or family history of PCa. For lethal PCa, the HR for a change (25th to 75th percentile) in mean gene expression level was 2.04 (95%CI 1.26-3.31; P = 0.004), adjusted for clinicopathological variables. The ROC curve for mean gene expression level alone (AUC = 0.740) did not perform as well as clinicopathological variables alone (AUC = 0.803) for predicting lethal PCa, and the addition of mean gene expression to clinicopathological variables did not substantially improve prediction (AUC = 0.827; P = 0.18). Higher TK1 expression was strongly associated with both recurrent (P = 6.7 × 10(-5)) and lethal (P = 6.4 × 10(-6)) PCa. CONCLUSIONS: Mean expression level for 30 selected cell cycle-regulated genes was unrelated to recurrence risk, but was associated with a twofold increase in risk of lethal PCa. However, gene expression had less discriminatory accuracy than clinical variables alone for predicting lethal events. Transcript levels for several genes in the panel were significantly overexpressed in lethal versus non-recurrent PCa.

CDK6 levels regulate quiescence exit in human hematopoietic stem cells.

Regulated blood production is achieved through the hierarchical organization of dormant hematopoietic stem cell (HSC) subsets that differ in self-renewal potential and division frequency, with long-term (LT)-HSCs dividing the least. The molecular mechanisms underlying this variability in HSC division kinetics are unknown. We report here that quiescence exit kinetics are differentially regulated within human HSC subsets through the expression level of CDK6. LT-HSCs lack CDK6 protein. Short-term (ST)-HSCs are also quiescent but contain high CDK6 protein levels that permit rapid cell cycle entry upon mitogenic stimulation. Enforced CDK6 expression in LT-HSCs shortens quiescence exit and confers competitive advantage without impacting function. Computational modeling suggests that this independent control of quiescence exit kinetics inherently limits LT-HSC divisions and preserves the HSC pool to ensure lifelong hematopoiesis. Thus, differential expression of CDK6 underlies heterogeneity in stem cell quiescence states that functionally regulates this highly regenerative system.

Transcriptome profiling of archived sectioned formalin-fixed paraffin-embedded (AS-FFPE) tissue for disease classification.

BACKGROUND: Archived tissues from previously completed prospective trials represent invaluable resource for biomarker development. However, such specimens are often stored as sections on glass slides, in which RNA is severely degraded due to prolonged air exposure. We evaluated whether a proportion of archived sectioned formalin-fixed paraffin-embedded (AS-FFPE) tissues yield transcriptome profiles comparable to freshly cut (FC) FFPE tissues, which can be used for retrospective class prediction analysis. METHODS: Genome-wide transcriptome profiles of 6 to 7-year-old AS-FFPE tissue sections (generated from 5 to 16-year-old blocks) of 83 hepatocellular carcinoma (HCC) and 47 liver cirrhosis samples were generated by using whole-genome DASL assay (Illumina) and digital transcript counting (nCounter) assay (NanoString), and gene signature-based prediction of HCC subclasses and prognosis was compared with previously generated FC-FFPE profiles from the same tissue blocks. RESULTS: RNA quality and assay reproducibility of AS-FFPE RNA were comparable to intermediate to poor quality FC-FFPE samples (RNA Integrity Number: up to 2.50, R-square for technical replicates: up to 0.93). Analyzable transcriptome profiles were obtained in 64 (77%) HCC and 36 (77%) cirrhosis samples. Statistically more confident predictions based on random resampling-based method (nearest template prediction) were obtained in 37 (58%) HCC and 13 (36%) cirrhosis samples. Predictions made in FC-FFPE profiles were reproduced in 36 (97%) HCC and 11 (85%) cirrhosis AS-FFPE profiles. nCounter assay was tested in 24 cirrhosis samples, which yielded confident prediction in 15 samples (63%), of which 10 samples (67%) showed concordant predictions with FC-FFPE profiles. CONCLUSIONS: AS-FFPE tissues yielded poorer quality RNA and transcriptome profiles compared to FC-FFPE tissues. Statistically more confident class prediction was feasible in 37 of 83 HCC samples and 13 of 47 cirrhosis samples. These results suggest that AS-FFPE tissues can be regarded as a resource for retrospective transcriptome-based class prediction analysis when they are the only available materials.

Reduced lymphoid lineage priming promotes human hematopoietic stem cell expansion.

The hematopoietic system sustains regeneration throughout life by balancing self-renewal and differentiation. To stay poised for mature blood production, hematopoietic stem cells (HSCs) maintain low-level expression of lineage-associated genes, a process termed lineage priming. Here, we modulated expression levels of Inhibitor of DNA binding (ID) proteins to ask whether lineage priming affects self-renewal of human HSCs. We found that lentiviral overexpression of ID proteins in cord blood HSCs biases myeloerythroid commitment at the expense of lymphoid differentiation. Conversely, reducing ID2 expression levels increases lymphoid potential. Mechanistically, ID2 inhibits the transcription factor E47 to attenuate B-lymphoid priming in HSCs and progenitors. Strikingly, ID2 overexpression also results in a 10-fold expansion of HSCs in serial limiting dilution assays, indicating that early lymphoid transcription factors antagonize human HSC self-renewal. The relationship between lineage priming and self-renewal can be exploited to increase expansion of transplantable human HSCs and points to broader implications for other stem cell populations.

APOBEC3B upregulation and genomic mutation patterns in serous ovarian carcinoma.

Ovarian cancer is a clinically and molecularly heterogeneous disease. The driving forces behind this variability are unknown. Here, we report wide variation in the expression of the DNA cytosine deaminase APOBEC3B, with elevated expression in the majority of ovarian cancer cell lines (three SDs above the mean of normal ovarian surface epithelial cells) and high-grade primary ovarian cancers. APOBEC3B is active in the nucleus of several ovarian cancer cell lines and elicits a biochemical preference for deamination of cytosines in 5'-TC dinucleotides. Importantly, examination of whole-genome sequence from 16 ovarian cancers reveals that APOBEC3B expression correlates with total mutation load as well as elevated levels of transversion mutations. In particular, high APOBEC3B expression correlates with C-to-A and C-to-G transversion mutations within 5'-TC dinucleotide motifs in early-stage high-grade serous ovarian cancer genomes, suggesting that APOBEC3B-catalyzed genomic uracil lesions are further processed by downstream DNA "repair" enzymes including error-prone translesion polymerases. These data identify a potential role for APOBEC3B in serous ovarian cancer genomic instability.

Regulation of adipose oestrogen output by mechanical stress.

Adipose stromal cells are the primary source of local oestrogens in adipose tissue, aberrant production of which promotes oestrogen receptor-positive breast cancer. Here we show that extracellular matrix compliance and cell contractility are two opposing determinants for oestrogen output of adipose stromal cells. Using synthetic extracellular matrix and elastomeric micropost arrays with tunable rigidity, we find that increasing matrix compliance induces transcription of aromatase, a rate-limiting enzyme in oestrogen biosynthesis. This mechanical cue is transduced sequentially by discoidin domain receptor 1, c-Jun N-terminal kinase 1, and phosphorylated JunB, which binds to and activates two breast cancer-associated aromatase promoters. In contrast, elevated cell contractility due to actin stress fibre formation dampens aromatase transcription. Mechanically stimulated stromal oestrogen production enhances oestrogen-dependent transcription in oestrogen receptor-positive tumour cells and promotes their growth. This novel mechanotransduction pathway underlies communications between extracellular matrix, stromal hormone output, and cancer cell growth within the same microenvironment.

The transcriptional architecture of early human hematopoiesis identifies multilevel control of lymphoid commitment.

Understanding how differentiation programs originate from the gene-expression 'landscape' of hematopoietic stem cells (HSCs) is crucial for the development of new clinical therapies. We mapped the transcriptional dynamics underlying the first steps of commitment by tracking transcriptome changes in human HSCs and eight early progenitor populations. We found that transcriptional programs were extensively shared, extended across lineage-potential boundaries and were not strictly lineage affiliated. Elements of stem, lymphoid and myeloid programs were retained in multilymphoid progenitors (MLPs), which reflected a hybrid transcriptional state. By functional single cell analysis, we found that the transcription factors Bcl-11A, Sox4 and TEAD1 (TEF1) governed transcriptional networks in MLPs, which led to B cell specification. Overall, we found that integrated transcriptome approaches can be used to identify previously unknown regulators of multipotency and show additional complexity in lymphoid commitment.

MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32.

BACKGROUND: Although microRNAs (miRNAs) are implicated in osteosarcoma biology and chemoresponse, miRNA prognostic models are still needed, particularly because prognosis is imperfectly correlated with chemoresponse. Formalin-fixed, paraffin-embedded tissue is a necessary resource for biomarker studies in this malignancy with limited frozen tissue availability. METHODS: We performed miRNA and mRNA microarray formalin-fixed, paraffin-embedded assays in 65 osteosarcoma biopsy and 26 paired post-chemotherapy resection specimens and used the only publicly available miRNA dataset, generated independently by another group, to externally validate our strongest findings (n = 29). We used supervised principal components analysis and logistic regression for survival and chemoresponse, and miRNA activity and target gene set analysis to study miRNA regulatory activity. RESULTS: Several miRNA-based models with as few as five miRNAs were prognostic independently of pathologically assessed chemoresponse (median recurrence-free survival: 59 months versus not-yet-reached; adjusted hazards ratio = 2.90; P = 0.036). The independent dataset supported the reproducibility of recurrence and survival findings. The prognostic value of the profile was independent of confounding by known prognostic variables, including chemoresponse, tumor location and metastasis at diagnosis. Model performance improved when chemoresponse was added as a covariate (median recurrence-free survival: 59 months versus not-yet-reached; hazard ratio = 3.91; P = 0.002). Most prognostic miRNAs were located at 14q32 - a locus already linked to osteosarcoma - and their gene targets display deregulation patterns associated with outcome. We also identified miRNA profiles predictive of chemoresponse (75% to 80% accuracy), which did not overlap with prognostic profiles. CONCLUSIONS: Formalin-fixed, paraffin-embedded tissue-derived miRNA patterns are a powerful prognostic tool for risk-stratified osteosarcoma management strategies. Combined miRNA and mRNA analysis supports a possible role of the 14q32 locus in osteosarcoma progression and outcome. Our study creates a paradigm for formalin-fixed, paraffin-embedded-based miRNA biomarker studies in cancer.

Highly parallel genome-wide expression analysis of single mammalian cells.

BACKGROUND: We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs. METHODOLOGY/PRINCIPAL FINDINGS: The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R(2)~0.76-0.80 between individual cells and R(2)~0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R(2)~0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input. CONCLUSIONS/SIGNIFICANCE: In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.

Angiogenic mRNA and microRNA gene expression signature predicts a novel subtype of serous ovarian cancer.

Ovarian cancer is the fifth leading cause of cancer death for women in the U.S. and the seventh most fatal worldwide. Although ovarian cancer is notable for its initial sensitivity to platinum-based therapies, the vast majority of patients eventually develop recurrent cancer and succumb to increasingly platinum-resistant disease. Modern, targeted cancer drugs intervene in cell signaling, and identifying key disease mechanisms and pathways would greatly advance our treatment abilities. In order to shed light on the molecular diversity of ovarian cancer, we performed comprehensive transcriptional profiling on 129 advanced stage, high grade serous ovarian cancers. We implemented a, re-sampling based version of the ISIS class discovery algorithm (rISIS: robust ISIS) and applied it to the entire set of ovarian cancer transcriptional profiles. rISIS identified a previously undescribed patient stratification, further supported by micro-RNA expression profiles, and gene set enrichment analysis found strong biological support for the stratification by extracellular matrix, cell adhesion, and angiogenesis genes. The corresponding "angiogenesis signature" was validated in ten published independent ovarian cancer gene expression datasets and is significantly associated with overall survival. The subtypes we have defined are of potential translational interest as they may be relevant for identifying patients who may benefit from the addition of anti-angiogenic therapies that are now being tested in clinical trials.

Non-housekeeping genes expressed in human trabecular meshwork cell cultures.

PURPOSE: To identify non-housekeeping genes definitively expressed in the human trabecular meshwork (TM). METHODS: Microarray gene expression data on TM cultured cells from four studies were compared. Genes that were queried in at least three studies and assessed to be expressed in at least three studies were considered definitively expressed genes of the human TM. Housekeeping genes were removed from this set of genes. The non-housekeeping TM gene profile was analyzed for pathway enrichment and microRNA targeting, using bioinformatics tools. The results were compared with results of previous non-array based studies. RESULTS: Nine hundred and sixty-two genes were identified as non-housekeeping TM expressed genes. Analysis of these by Kyoto Encyclopedia of Genes and Genomes led to identification of two enriched biologic pathways that achieved a highly significant Bonferroni p-value (p≤0.01): focal adhesion and extracellular matrix (ECM)-receptor interaction. Many of the genes were previously implicated in TM-related functions and the TM-associated disease glaucoma; however, some are novel. MicroRNAs known to be expressed in the trabecular meshwork were predicted to target some of the genes. Ten genes identified here, ALDH1A1 (aldehyde dehydrogenase 1 family, member A1), CDH11 (cadherin 11, type 2, OB-cadherin), CXCR7 (chemokine (C-X-C motif) receptor 7), CHI3L1 (chitinase 3-like 1), FGF2 (fibroblast growth factor 2), GNG11 (guanine nucleotide binding protein [G protein], gamma 11), IGFBP5 (insulin-like growth factor binding protein 5), PTPRM (protein tyrosine phosphatase, receptor type, M), RGS5 (regulator of G-protein signaling 5), and TUSC3 (tumor suppressor candidate 3), were also reported as TM expressed genes in three earlier non-microarray based studies. CONCLUSIONS: A transcriptome consisting of 962 non-housekeeping genes definitively expressed in the human TM was identified. Multiple genes and microRNAs are proposed for further study for a better understanding of TM physiology.

Gene expression profiling in formalin-fixed, paraffin-embedded tissues using the whole-genome DASL assay.

Here, we provide a detailed technical description of a gene expression assay (Whole-Genome DASL (WG-DASL)), which not only enables whole-genome transcriptional profiling of degraded material, such as formalin-fixed, paraffin-embedded tissues, but is also capable of generating robust profiles with low input intact RNA. The WG-DASL assay combines target-specific annealing, extension, and ligation events followed by universal PCR and labeling steps to generate highly multiplexed Cy3-labeled products. These short products, which are single-stranded, are directly hybridized to a whole-genome expression BeadChip (HumanRef-8) containing probe content corresponding to ∼24 K RefSeq transcripts. After washing and imaging, fluorescence emissions are quantitatively recorded for each probe using high-resolution confocal scanners and imaging software. GenomeStudio software allows direct analysis of mRNA expression data and provides results in standard file formats that can be readily exported and analyzed with most standard gene expression analysis software programs. This technology is particularly useful for genome-wide expression profiling in degraded, archived material, including limited quantities of clinical samples, such as microdissected and biopsied materials.

Gene-expression signature of vascular invasion in hepatocellular carcinoma.

BACKGROUND & AIMS: Vascular invasion is a major predictor of tumor recurrence after surgical treatments for hepatocellular carcinoma (HCC). While macroscopic vascular invasion can be detected by radiological techniques, pre-operative detection of microscopic vascular invasion, which complicates 30-40% of patients with early tumors, remains elusive. METHODS: A total of 214 patients with hepatocellular carcinoma who underwent resection were included in the study. By using genome-wide gene-expression profiling of 79 hepatitis C-related hepatocellular carcinoma samples (training set), a gene-expression signature associated with vascular invasion was defined. The signature was validated in formalin-fixed paraffin-embedded tissues obtained from an independent set of 135 patients with various etiologies. RESULTS: A 35-gene signature of vascular invasion was defined in the training set, predicting vascular invasion with an accuracy of 69%. The signature was independently associated with the presence of vascular invasion (OR 3.38, 95% CI 1.48-7.71, p=0.003) along with tumor size (diameter greater than 3 cm, OR 2.66, 95% CI 1.17-6.05, p=0.02). In the validation set, the signature discarded the presence of vascular invasion with a negative predictive value of 0.77, and significantly improved the diagnostic power of tumor size alone (p=0.045). CONCLUSIONS: The assessment of a gene-expression signature obtained from resected biopsied tumor specimens improved the diagnosis of vascular invasion beyond clinical variable-based prediction. The signature may aid in candidate selection for liver transplantation, and guide the design of clinical trials with experimental adjuvant therapies.

Combining clinical, pathology, and gene expression data to predict recurrence of hepatocellular carcinoma.

BACKGROUND & AIMS: In approximately 70% of patients with hepatocellular carcinoma (HCC) treated by resection or ablation, disease recurs within 5 years. Although gene expression signatures have been associated with outcome, there is no method to predict recurrence based on combined clinical, pathology, and genomic data (from tumor and cirrhotic tissue). We evaluated gene expression signatures associated with outcome in a large cohort of patients with early stage (Barcelona-Clinic Liver Cancer 0/A), single-nodule HCC and heterogeneity of signatures within tumor tissues. METHODS: We assessed 287 HCC patients undergoing resection and tested genome-wide expression platforms using tumor (n = 287) and adjacent nontumor, cirrhotic tissue (n = 226). We evaluated gene expression signatures with reported prognostic ability generated from tumor or cirrhotic tissue in 18 and 4 reports, respectively. In 15 additional patients, we profiled samples from the center and periphery of the tumor, to determine stability of signatures. Data analysis included Cox modeling and random survival forests to identify independent predictors of tumor recurrence. RESULTS: Gene expression signatures that were associated with aggressive HCC were clustered, as well as those associated with tumors of progenitor cell origin and those from nontumor, adjacent, cirrhotic tissues. On multivariate analysis, the tumor-associated signature G3-proliferation (hazard ratio [HR], 1.75; P = .003) and an adjacent poor-survival signature (HR, 1.74; P = .004) were independent predictors of HCC recurrence, along with satellites (HR, 1.66; P = .04). Samples from different sites in the same tumor nodule were reproducibly classified. CONCLUSIONS: We developed a composite prognostic model for HCC recurrence, based on gene expression patterns in tumor and adjacent tissues. These signatures predict early and overall recurrence in patients with HCC, and complement findings from clinical and pathology analyses.

Expression profiling of formalin-fixed paraffin-embedded primary breast tumors using cancer-specific and whole genome gene panels on the DASL® platform.

BACKGROUND: The cDNA-mediated Annealing, extension, Selection and Ligation (DASL) assay has become a suitable gene expression profiling system for degraded RNA from paraffin-embedded tissue. We examined assay characteristics and the performance of the DASL 502-gene Cancer Panel v1 (1.5K) and 24,526-gene panel (24K) platforms at differentiating nine human epidermal growth factor receptor 2- positive (HER2+) and 11 HER2-negative (HER2-) paraffin-embedded breast tumors. METHODS: Bland-Altman plots and Spearman correlations evaluated intra/inter-panel agreement of normalized expression values. Unequal-variance t-statistics tested for differences in expression levels between HER2 + and HER2 - tumors. Regulatory network analysis was performed using Metacore (GeneGo Inc., St. Joseph, MI). RESULTS: Technical replicate correlations ranged between 0.815-0.956 and 0.986-0.997 for the 1.5K and 24K panels, respectively. Inter-panel correlations of expression values for the common 498 genes across the two panels ranged between 0.485-0.573. Inter-panel correlations of expression values of 17 probes with base-pair sequence matches between the 1.5K and 24K panels ranged between 0.652-0.899. In both panels, erythroblastic leukemia viral oncogene homolog 2 (ERBB2) was the most differentially expressed gene between the HER2 + and HER2 - tumors and seven additional genes had p-values < 0.05 and log2 -fold changes > |0.5| in expression between HER2 + and HER2 - tumors: topoisomerase II alpha (TOP2A), cyclin a2 (CCNA2), v-fos fbj murine osteosarcoma viral oncogene homolog (FOS), wingless-type mmtv integration site family, member 5a (WNT5A), growth factor receptor-bound protein 7 (GRB7), cell division cycle 2 (CDC2), and baculoviral iap repeat-containing protein 5 (BIRC5). The top 52 discriminating probes from the 24K panel are enriched with genes belonging to the regulatory networks centered around v-myc avian myelocytomatosis viral oncogene homolog (MYC), tumor protein p53 (TP53), and estrogen receptor α (ESR1). Network analysis with a two-step extension also showed that the eight discriminating genes common to the 1.5K and 24K panels are functionally linked together through MYC, TP53, and ESR1. CONCLUSIONS: The relative RNA abundance obtained from two highly differing density gene panels are correlated with eight common genes differentiating HER2 + and HER2 - breast tumors. Network analyses demonstrated biological consistency between the 1.5K and 24K gene panels.