RNA expression analysis from formalin fixed paraffin embedded tissues

Formalin fixation and paraffin embedding (FFPE) is the most commonly used method worldwide for tissue storage. This method preserves the tissue integrity but causes extensive damage to nucleic acids stored within the tissue. As methods for measuring gene expression such as RT-PCR and microarray are adopted into clinical practice there is an increasing necessity to access the wealth of information locked in the Formalin fixation and paraffin embedding archives. This paper reviews the progress in this field and discusses the unique opportunities that exist for the application of these techniques in the development of personalized medicine.

Association of gene expression with sequential proliferation, differentiation and tumor formation in murine skin

Differential gene expression in two established initiation and promotion skin carcinogenesis models during promotion and tumor formation was determined by microarray technology with the purpose of distinguishing the genes more associated with neoplastic transformation from those linked with proliferation and differentiation. The first model utilized dimethylbenz[a]anthracene initiation and 12-O-tetradecanoylphorbol 13-acetate (TPA) promotion in the FVB/N mouse, and the second TPA promotion of the Tg.Ac mouse, which is endogenously initiated by virtue of an activated Ha-ras transgene. Comparison of gene expression profiles across the two models identified genes whose altered expression was associated with papilloma formation rather than TPA-induced proliferation and differentiation. DMBA suppressed TPA-induced differentiation which allowed identification of those genes associated more specifically with differentiation rather than proliferation. EASE (Expression Analysis Systemic Explorer) indicated a correlation between muscle-associated genes and skin differentiation, whereas genes involved with protein biosynthesis were strongly correlated with proliferation. For verification the altered expression of selected genes were confirmed by RT-PCR; Carbonic anhydrase 2, Thioredoxin 1 and Glutathione S-transferase omega 1 associated with papilloma formation and Enolase 3, Cystatin 6 and Filaggrin associated with TPA-induced proliferation and differentiation. In situ analysis located the papillomas Glutathione S-transferase omega 1 expression to the proliferating areas of the papillomas. Thus we have identified profiles of differential gene expression associated with the tumorigenesis and promotion stages for skin carcinogenesis in the mouse.

Immune activation in Retinal Aging: A Gene Expression Study.

PURPOSE
To investigate changes in gene expression during aging of the retina in the mouse.

METHODS
Total RNA was extracted from the neuroretina of young (3-month-old) and old (20-month-old) mice and processed for microarray analysis. Age-related, differentially expressed genes were assessed by the empiric Bayes shrinkagemoderated t-statistics method. Statistical significance was based on dual criteria of a ratio of change in gene expression >2 and a P < 0.01. Differential expression in 11 selected genes was further verified by real-time PCR. Functional pathways involved in retinal ageing were analyzed by an online software package (DAVID-2008) in differentially expressed gene lists. Age-related changes in differential expression in the identified retinal molecular pathways were further confirmed by immunohistochemical staining of retinal flat mounts and retinal cryosections.

RESULTS
With ageing of the retina, 298 genes were upregulated and 137 genes were downregulated. Functional annotation showed that genes linked to immune responses (Ir genes) and to tissue stress/injury responses (TS/I genes) were most likely to be modified by ageing. The Ir genes affected included those regulating leukocyte activation, chemotaxis, endocytosis, complement activation, phagocytosis, and myeloid cell differentiation, most of which were upregulated, with only a few downregulated. Increased microglial and complement activation in the aging retina was further confirmed by confocal microscopy of retinal tissues. The most strongly upregulated gene was the calcitonin receptor (Calcr; >40-fold in old versus young mice).

CONCLUSIONS
The results suggest that retinal ageing is accompanied by activation of gene sets, which are involved in local inflammatory responses. A modified form of low-grade chronic inflammation (para-inflammation) characterizes these aging changes and involves mainly the innate immune system. The marked upregulation of Calcr in ageing mice most likely reflects this chronic inflammatory/stress response, since calcitonin is a known systemic biomarker of inflammation/sepsis. © Association for Research in Vision and Ophthalmology.

Hyperacetylation in prostate cancer induces cell cycle aberrations, chromatin reorganization and altered gene expression profiles.

Histone acetylation is a fundamental mechanism in the regulation of local chromatin conformation and gene expression. Research has focused on the impact of altered epigenetic environments on the expression of specific genes and their pathways. However, changes in histone acetylation also have a global impact on the cell. In this study we used digital texture analysis to assess global chromatin patterns following treatment with trichostatin A (TSA) and have observed significant alterations in the condensation and distribution of higher-order chromatin, which were associated with altered gene expression profiles in both immortalised normal PNT1A prostate cell line and androgen-dependent prostate cancer cell line LNCaP. Furthermore, the extent of TSA-induced disruption was both cell cycle and cell line dependent. This was illustrated by the identification of sub-populations of prostate cancer cells expressing high levels of H3K9 acetylation in the G2/M phase of the cell cycle that were absent in normal cell populations. In addition, the analysis of enriched populations of G1 cells showed a global decondensation of chromatin exclusively in normal cells.

Thrombopoietin, flt3-ligand and c-kit-ligand modulate HOX gene expression in expanding cord blood CD133(+) cells

Haemopoietic stem/progenitor cell (HSPC) development is regulated by extrinsic and intrinsic stimuli. Extrinsic modulators include growth factors and cell adhesion molecules, whereas intrinsic regulation is achieved with many transcription factor families, of which the HOX gene products are known to be important in haemopoiesis. Umbilical cord blood CD133(+) HSPC proliferation potential was tested in liquid culture with 'TPOFLK' (thrombopoietin, flt-3 ligand and c-kit ligand, promoting HSPC survival and self-renewal), in comparison to 'K36EG' (c-kit-ligand, interleukins-3 and -6, erythropoietin and granulocyte colony-stimulating factor, inducing haemopoietic differentiation). TPOFLK induced a higher CD133(+) HSPC proliferation (up to 60-fold more, at week 8) and maintained a higher frequency of the primitive colony-forming cells than K36EG. Quantitative polymerase chain reaction analysis revealed opposite expression patterns for specific HOX genes in expanding cord blood CD133(+) HSPC. After 8 weeks in liquid culture, TPOFLK increased the expression of HOX B3, B4 and A9 (associated with uncommitted HSPC) and reduced the expression of HOX B8 and A10 (expressed in committed myeloid cells) when compared to K36EG. These results suggest that TPOFLK induces CD133(+) HSPC proliferation, self-renewal and maintenance, up-regulation of HOX B3, B4 and A9 and down-regulation of HOX B8 and A10 gene expression.

Activated oncogenic pathways and therapeutic targets in extranodal nasal-type NK/T cell lymphoma revealed by gene expression profiling

We performed comprehensive genome-wide gene expression profiling (GEP) of extranodal nasal-type natural killer/T-cell lymphoma (NKTL) using formalin-fixed, paraffin-embedded tissue (n = 9) and NK cell lines (n = 5) in comparison with normal NK cells, with the objective of understanding the oncogenic pathways involved in the pathogenesis of NKTL and to identify potential therapeutic targets. Pathway and network analysis of genes differentially expressed between NKTL and normal NK cells revealed significant enrichment for cell cycle-related genes and pathways, such as PLK1, CDK1, and Aurora-A. Furthermore, our results demonstrated a pro-proliferative and anti-apoptotic phenotype in NKTL characterized by activation of Myc and nuclear factor kappa B (NF-kappa B), and deregulation of p53. In corroboration with GEP findings, a significant percentage of NKTLs (n = 33) overexpressed c-Myc (45.4%), p53 (87.9%), and NF-kappa B p50 (67.7%) on immunohistochemistry using a tissue microarray containing 33 NKTL samples. Notably, overexpression of survivin was observed in 97% of cases. Based on our findings, we propose a model of NKTL pathogenesis where deregulation of p53 together with activation of Myc and NF-kappa B, possibly driven by EBV LMP-1, results in the cumulative up-regulation of survivin. Down-regulation of survivin with Terameprocol (EM-1421, a survivin inhibitor) results in reduced cell viability and increased apoptosis in tumour cells, suggesting that targeting survivin may be a potential novel therapeutic strategy in NKTL. Copyright (C) 2011 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

An international standardization programme towards the application of gene expression profiling in routine leukaemia diagnostics: the Microarray Innovations in LEukemia study prephase

Gene expression profiling has the potential to enhance current methods for the diagnosis of haematological malignancies. Here, we present data on 204 analyses from an international standardization programme that was conducted in 11 laboratories as a prephase to the Microarray Innovations in LEukemia (MILE) study. Each laboratory prepared two cell line samples, together with three replicate leukaemia patient lysates in two distinct stages: (i) a 5-d course of protocol training, and (ii) independent proficiency testing. Unsupervised, supervised, and r(2) correlation analyses demonstrated that microarray analysis can be performed with remarkably high intra-laboratory reproducibility and with comparable quality and reliability.

Oxidised, glycated LDL selectively influences tissue inhibitor of metalloproteinase-3 gene expression and protein production in human retinal capillary pericytes

Aims/hypothesis: Matrix metalloproteinases (MMPs) and their natural inhibitors, tissue inhibitor of metalloproteinases (TIMPs), regulate important biological processes including the homeostasis of the extracellular matrix, proteolysis of cell surface proteins, proteinase zymogen activation, angiogenesis and inflammation. Studies have shown that their balance is altered in retinal microvascular tissues in diabetes. Since LDLs modified by oxidation/glycation are implicated in the pathogenesis of diabetic vascular complications, we examined the effects of modified LDL on the gene expression and protein production of MMPs and TIMPs in retinal pericytes. Methods: Quiescent human retinal pericytes were exposed to native LDL (N-LDL), glycated LDL (G-LDL) and heavily oxidised and glycated LDL (HOG-LDL) for 24 h. We studied the expression of the genes encoding MMPs and TIMPs mRNAs by analysis of microarray data and quantitative PCR, and protein levels by immunoblotting and ELISA. Results: Microarray analysis showed that MMP1, MMP2, MMP11, MMP14 and MMP25 and TIMP1, TIMP2, TIMP3 and TIMP4 were expressed in pericytes. Of these, only TIMP3 mRNA showed altered regulation, being expressed at significantly lower levels in response to HOG- vs N-LDL. Quantitative PCR and immunoblotting of cell/matrix proteins confirmed the reduction in TIMP3 mRNA and protein in response to HOG-LDL. In contrast to cellular TIMP3 protein, analysis of secreted TIMP1, TIMP2, MMP1 and collagenase activity indicated no changes in their production in response to modified LDL. Combined treatment with N- and HOG-LDL restored TIMP3 mRNA expression to a level comparable with that after N-LDL alone. Conclusions/interpretation: Among the genes encoding for MMPs and TIMPs expressed in retinal pericytes, TIMP3 is uniquely regulated by HOG-LDL. Reduced TIMP3 expression might contribute to microvascular abnormalities in diabetic retinopathy. © 2007 Springer-Verlag.

Effects of oxidized and glycated LDL on gene expression in human retinal capillary pericytes

Modified (oxidized and/or glycated) low-density lipoproteins (LDLs) have been implicated in retinal pericyte loss, one of the major pathologic features of early-stage diabetic retinopathy. To delineate underlying molecular mechanisms, the present study was designed to explore the global effects of modified LDL on pericyte gene expression.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Gene expression in normal Urothelium depends on location within the bladder:A possible link to bladder carcinogenesis

Objectives: Clinical studies have shown that more than 70% of primary bladder tumours arise in the area around the ureteric orifice. In this study a genomic approach was taken to explore the molecular mechanisms that may influence this phenomenon.

Methods: RNA was isolated from each individual normal ureteric orifice and the dome biopsy from 33 male patients. Equal amounts of the pooled ureteric orifice and dome mRNAs were labelled with Cy3 and Cy5, respectively before hybridising to the gene chip (UniGEM 2.0, Incyte Genomics Inc., Wilmington, Delaware, USA). Results: Significant changes (more than a twofold difference) in gene expression were observed in 3.1% (312) of the 10,176 gene array: 211 genes upregulated and 101 downregulated. Analysis of Cdc25B, TK1, PKM, and PDGFra with RT-PCR supported the reliability of the microarray result. Seladin-1 was the most upregulated gene in the ureteric orifice: 8.3-fold on the microarray and 11.4-fold by real time PCR.

Conclusions: Overall, this study suggests significant altered gene expression between these two anatomically distinct areas of the normal human bladder. Of particular note is Seladin-1, whose significance in cancer is yet to be clarified. Further studies of the genes discovered by this work will help clarify which of these differences influence primary bladder carcinogenesis. (c) 2006 European Association of Urology. Published by Elsevier B.V. All rights reserved.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

QUADrATiC: scalable gene expression connectivity mapping for repurposing FDA-approved therapeutics

Background: Gene expression connectivity mapping has proven to be a powerful and flexible tool for research. Its application has been shown in a broad range of research topics, most commonly as a means of identifying potential small molecule compounds, which may be further investigated as candidates for repurposing to treat diseases. The public release of voluminous data from the Library of Integrated Cellular Signatures (LINCS) programme further enhanced the utilities and potentials of gene expression connectivity mapping in biomedicine. Results: We describe QUADrATiC (http://go.qub.ac.uk/QUADrATiC), a user-friendly tool for the exploration of gene expression connectivity on the subset of the LINCS data set corresponding to FDA-approved small molecule compounds. It enables the identification of compounds for repurposing therapeutic potentials. The software is designed to cope with the increased volume of data over existing tools, by taking advantage of multicore computing architectures to provide a scalable solution, which may be installed and operated on a range of computers, from laptops to servers. This scalability is provided by the use of the modern concurrent programming paradigm provided by the Akka framework. The QUADrATiC Graphical User Interface (GUI) has been developed using advanced Javascript frameworks, providing novel visualization capabilities for further analysis of connections. There is also a web services interface, allowing integration with other programs or scripts.Conclusions: QUADrATiC has been shown to provide an improvement over existing connectivity map software, in terms of scope (based on the LINCS data set), applicability (using FDA-approved compounds), usability and speed. It offers potential to biological researchers to analyze transcriptional data and generate potential therapeutics for focussed study in the lab. QUADrATiC represents a step change in the process of investigating gene expression connectivity and provides more biologically-relevant results than previous alternative solutions.

Immune-derived PD-L1 gene expression defines a subgroup of stage II/III colorectal cancer patients with favorable prognosis that may be harmed by adjuvant chemotherapy

A recent phase 2 study of metastatic colorectal carcinoma (CRC) patients showed that mismatch repair gene status was predictive of clinical response to PD-1-targeting immune checkpoint blockade. Further examination revealed strong correlation between PD-L1 protein expression and microsatellite instability (MSI) in stage IV CRC, suggesting that the amount of PD-L1 protein expression could identify late stage patients who may benefit from immunotherapy. To assess whether the clinical associations between PD-L1 gene expression and MSI identified in metastatic CRC are also present in stage II/III CRC, we used in silico analysis to elucidate the cell types expressing the PD-L1 gene. We found a significant association of PD-L1 gene expression with MSI in early stage CRC (P < 0.001) and show that unlike in non-CRC tumors, PD-L1 is derived predominantly from the immune infiltrate. We demonstrate that PD-L1 gene expression has positive prognostic value in the adjuvant disease setting (PD-L1low v PD-L1high HR = 9.09; CI, 2.11-39.10). PD-L1 gene expression had predictive value, as patients with high PD-L1 expression appear to be harmed by standard-of-care treatment (HR = 4.95; CI,1.10-22.35). Building on the promising results from the metastatic CRC PD-1-targeting trial, we provide compelling evidence that PD-L1high/MSI/immunehigh stage II/III CRC patients should not receive standard chemotherapy. This conclusion supports the rationale to clinically evaluate this patient subgroup for PD-1 blockade treatment.

Association Between Results of a Gene Expression Signature Assay and Recurrence-Free Interval in Patients With Stage II Colon Cancer in Cancer and Leukemia Group B 9581 (Alliance)

PURPOSE: Conventional staging methods are inadequate to identify patients with stage II colon cancer (CC) who are at high risk of recurrence after surgery with curative intent. ColDx is a gene expression, microarray-based assay shown to be independently prognostic for recurrence-free interval (RFI) and overall survival in CC. The objective of this study was to further validate ColDx using formalin-fixed, paraffin-embedded specimens collected as part of the Alliance phase III trial, C9581.

PATIENTS AND METHODS: C9581 evaluated edrecolomab versus observation in patients with stage II CC and reported no survival benefit. Under an initial case-cohort sampling design, a randomly selected subcohort (RS) comprised 514 patients from 901 eligible patients with available tissue. Forty-nine additional patients with recurrence events were included in the analysis. Final analysis comprised 393 patients: 360 RS (58 events) and 33 non-RS events. Risk status was determined for each patient by ColDx. The Self-Prentice method was used to test the association between the resulting ColDx risk score and RFI adjusting for standard prognostic variables.

RESULTS: Fifty-five percent of patients (216 of 393) were classified as high risk. After adjustment for prognostic variables that included mismatch repair (MMR) deficiency, ColDx high-risk patients exhibited significantly worse RFI (multivariable hazard ratio, 2.13; 95% CI, 1.3 to 3.5; P < .01). Age and MMR status were marginally significant. RFI at 5 years for patients classified as high risk was 82% (95% CI, 79% to 85%), compared with 91% (95% CI, 89% to 93%) for patients classified as low risk.

CONCLUSION: ColDx is associated with RFI in the C9581 subsample in the presence of other prognostic factors, including MMR deficiency. ColDx could be incorporated with the traditional clinical markers of risk to refine patient prognosis.