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.

IGFBP7 promoter methylation and gene expression analysis in prostate cancer

PURPOSE: IGFBP7 belongs to a family of insulin-like growth factor-1 regulatory binding proteins. IGFBP7 hypermethylation is associated with its down-regulation in various carcinomas. In prostate cancer IGFBP7 down-regulation has been widely reported but to our knowledge the mechanisms behind this event are unknown. We performed a denaturing high performance liquid chromatography screening and validation strategy to profile the methylation status of IGFBP7 in prostate cancer.

MATERIALS AND METHODS: We combined denaturing high performance liquid chromatography and bisulfite sequencing to examine IGFBP7 methylation in a panel of prostate cancer cell lines. Quantitative methylation specific polymerase chain reaction was used to determine methylation levels in prostate tissue specimens of primary prostate cancer, histologically benign prostate adjacent to tumor, high grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. IGFBP7 gene expression was measured by quantitative methylation specific polymerase chain reaction in cell lines and tissue specimens.

RESULTS: IGFBP7 was methylated in the 4 prostate cancer cell lines DU145, LNCaP, PC-3 and 22RV1. Quantitative methylation specific polymerase chain reaction analysis revealed that promoter methylation was associated with decreased IGFBP7 expression. Quantitative methylation specific polymerase chain reaction showed that IGFBP7 methylation was more frequently detected in prostate cancer (60% (31/52)) and high grade prostatic intraepithelial neoplasia (40% (6/15)) samples compared to histologically benign prostate adjacent to tumor (10%) and benign prostatic hyperplasia (0%) samples.

CONCLUSIONS: To our knowledge this is the first report of aberrant IGFBP7 promoter hypermethylation and concurrent IGFBP7 gene silencing in prostate cancer cell lines. Results demonstrate that CpG methylation of IGFBP7 may represent a novel biomarker of prostate cancer and pre-invasive neoplasms. Thus, future examination of IGFBP7 methylation and expression in a larger patient cohort, including bodily fluids, is justified to further evaluate its role in a diagnostic and prognostic setting.

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.

Insights into the multistep transformation of MGUS to myeloma using microarray expression analysis.

To define specific pathways important in the multistep transformation process of normal plasma cells (PCs) to monoclonal gammopathy of uncertain significance (MGUS) and multiple myeloma (MM), we have applied microarray analysis to PCs from 5 healthy donors (N), 7 patients with MGUS, and 24 patients with newly diagnosed MM. Unsupervised hierarchical clustering using 125 genes with a large variation across all samples defined 2 groups: N and MGUS/MM. Supervised analysis identified 263 genes differentially expressed between N and MGUS and 380 genes differentially expressed between N and MM, 197 of which were also differentially regulated between N and MGUS. Only 74 genes were differentially expressed between MGUS and MM samples, indicating that the differences between MGUS and MM are smaller than those between N and MM or N and MGUS. Differentially expressed genes included oncogenes/tumor-suppressor genes (LAF4, RB1, and disabled homolog 2), cell-signaling genes (RAS family members, B-cell signaling and NF-kappaB genes), DNA-binding and transcription-factor genes (XBP1, zinc finger proteins, forkhead box, and ring finger proteins), and developmental genes (WNT and SHH pathways). Understanding the molecular pathogenesis of MM by gene expression profiling has demonstrated sequential genetic changes from N to malignant PCs and highlighted important pathways involved in the transformation of MGUS to MM.

The expression of pain in infants and toddlers:developmental changes in facial action

Age-related changes in the facial expression of pain during the first 18 months of life have important implications for our understanding of pain and pain assessment. We examined facial reactions video recorded during routine immunization injections in 75 infants stratified into 2-, 4-, 6-, 12-, and 18-month age groups. Two facial coding systems differing in the amount of detail extracted were applied to the records. In addition, parents completed a brief questionnaire that assessed child temperament and provided background information. Parents' efforts to soothe the children also were described. While there were consistencies in facial displays over the age groups, there also were differences on both measures of facial activity, indicating systematic variation in the nature and severity of distress. The least pain was expressed by the 4-month age group. Temperament was not related to the degree of pain expressed. Systematic variations in parental soothing behaviour indicated accommodation to the age of the child. Reasons for the differing patterns of facial activity are examined, with attention paid to the development of inhibitory mechanisms and the role of negative emotions such as anger and anxiety.

Neonatal facial and cry responses to invasive and non-invasive procedures

Evaluation of pain in neonates is difficult due to their limited means of communication. The aim was to determine whether behavioural reactions of cry and facial activity provoked by an invasive procedure could be discriminated from responses to non-invasive tactile events. Thirty-six healthy full-term infants (mean age 2.2 h) received 3 procedures in counterbalanced order: intramuscular injection, application of triple dye to the umbilical stub, and rubbing thigh with alcohol. Significant effects of procedure were found for total face activity and latency to face movement. A cluster of facial actions comprised of brow bulging, eyes squeezed shut, deepening of the naso-labial furrow and open mouth was associated most frequently with the invasive procedure. Comparisons between the 2 non-invasive procedures showed more facial activity to thigh swabbing and least to application of triple dye to the umbilical cord. Acoustic analysis of cry showed statistically significant differences across procedures only for latency to cry and cry duration for the group as a whole. However, babies who cried to two procedures showed higher pitch and greater intensity to the injection. There were no significant differences in melody, dysphonation, or jitter. Methodological difficulties for investigators in this area were examined, including criteria for the selection of cries for analysis, and the logical and statistical challenges of contrasting cries induced by different conditions when some babies do not always cry. It was concluded that facial expression, in combination with short latency to onset of cry and long duration of first cry cycle typifies reaction to acute invasive procedures.

Pain expression in neonates:facial action and cry

Pain expression in neonates instigated by heel-lance for blood sampling purposes was systematically described using measures of facial expression and cry and compared across sleep/waking states and sex. From gate-control theory it was hypothesized that pain behavior would vary with the ongoing functional state of the infant, rather than solely reflecting tissue insult. Awake-alert but inactive infants responded with the most facial activity, consistent with current views that infants in this state are most receptive to environmental stimulation. Infants in quiet sleep showed the least facial reaction and the longest latency to cry. Fundamental frequency of cry was not related to sleep/waking state. This suggested that findings from the cry literature on qualities of pain cry as a reflection of nervous system 'stress', in unwell newborns, do not generalize directly to healthy infants as a function of state. Sex differences were apparent in speed of response, with boys showing shorter time to cry and to display facial action following heel-lance. The findings of facial action variation across sleep/waking state were interpreted as indicating that the biological and behavioral context of pain events affects behavioral expression, even at the earliest time developmentally, before the opportunity for learned response patterns occurs. Issues raised by the study include the importance of using measurement techniques which are independent of preconceived categories of affective response.

Network signatures based on gene pair expression ratios improve classification and the analysis of muscle-invasive urothelial cancer

Urothelial cancer (UC) is highly recurrent and can progress from non-invasive (NMIUC) to a more aggressive muscle-invasive (MIUC) subtype that invades the muscle tissue layer of the bladder. We present a proof of principle study that network-based features of gene pairs can be used to improve classifier performance and the functional analysis of urothelial cancer gene expression data. In the first step of our procedure each individual sample of a UC gene expression dataset is inflated by gene pair expression ratios that are defined based on a given network structure. In the second step an elastic net feature selection procedure for network-based signatures is applied to discriminate between NMIUC and MIUC samples. We performed a repeated random subsampling cross validation in three independent datasets. The network signatures were characterized by a functional enrichment analysis and studied for the enrichment of known cancer genes. We observed that the network-based gene signatures from meta collections of proteinprotein interaction (PPI) databases such as CPDB and the PPI databases HPRD and BioGrid improved the classification performance compared to single gene based signatures. The network based signatures that were derived from PPI databases showed a prominent enrichment of cancer genes (e.g., TP53, TRIM27 and HNRNPA2Bl). We provide a novel integrative approach for large-scale gene expression analysis for the identification and development of novel diagnostical targets in bladder cancer. Further, our method allowed to link cancer gene associations to network-based expression signatures that are not observed in gene-based expression signatures.

A global expression-based analysis of the consequences of the t(4;14) translocation in myeloma

Purpose: Our purpose in this report was to define genes and pathways dysregulated as a consequence of the t(4;14) in myeloma, and to gain insight into the downstream functional effects that may explain the different prognosis of this subgroup.Experimental Design: Fibroblast growth factor receptor 3 (FGFR3) overexpression, the presence of immunoglobulin heavy chain-multiple myeloma SET domain (IgH-MMSET) fusion products and the identification of t(4;14) breakpoints were determined in a series of myeloma cases. Differentially expressed genes were identified between cases with (n = 55) and without (n = 24) a t(4;14) by using global gene expression analysis.Results: Cases with a t(4;14) have a distinct expression pattern compared with other cases of myeloma. A total of 127 genes were identified as being differentially expressed including MMSET and cyclin D2, which have been previously reported as being associated with this translocation. Other important functional classes of genes include cell signaling, apoptosis and related genes, oncogenes, chromatin structure, and DNA repair genes. Interestingly, 25% of myeloma cases lacking evidence of this translocation had up-regulation of the MMSET transcript to the same level as cases with a translocation.Conclusions: t(4;14) cases form a distinct subgroup of myeloma cases with a unique gene signature that may account for their poor prognosis. A number of non-t(4;14) cases also express MMSET consistent with this gene playing a role in myeloma pathogenesis.

Individual mRNA expression profiles reveal the effects of specific microRNAs

BACKGROUND: MicroRNAs (miRNAs) are oligoribonucleotides with an important role in regulation of gene expression at the level of translation. Despite imperfect target complementarity, they can also significantly reduce mRNA levels. The validity of miRNA target gene predictions is difficult to assess at the protein level. We sought, therefore, to determine whether a general lowering of predicted target gene mRNA expression by endogenous miRNAs was detectable within microarray gene expression profiles. RESULTS: The target gene sets predicted for each miRNA were mapped onto known gene expression data from a range of tissues. Whether considering mean absolute target gene expression, rank sum tests or 'ranked ratios', many miRNAs with significantly reduced target gene expression corresponded to those known to be expressed in the cognate tissue. Expression levels of miRNAs with reduced target mRNA levels were higher than those of miRNAs with no detectable effect on mRNA expression. Analysis of microarray data gathered after artificial perturbation of expression of a specific miRNA confirmed the predicted increase or decrease in influence of the altered miRNA upon mRNA levels. Strongest associations were observed with targets predicted by TargetScan. CONCLUSION: We have demonstrated that the effect of a miRNA on its target mRNAs' levels can be measured within a single gene expression profile. This emphasizes the extent of this mode of regulation in vivo and confirms that many of the predicted miRNA-mRNA interactions are correct. The success of this approach has revealed the vast potential for extracting information about miRNA function from gene expression profiles.

Expression of the phosphonoalanine-degradative gene cluster from Variovorax sp Pal2 is induced by growth on phosphonoalanine and phosphonopyruvate

The phosphonopyruvate hydrolase (PalA) found in Variovorax sp., Pal2, is a novel carbon-phosphorus bond cleavage enzyme, which is expressed even in the presence of high levels of phosphate, thus permitting phosphonopyruvate to be used as the sole carbon and energy source. Analysis of the regions adjacent to the palA gene revealed the presence of the five structural genes that constitute the 2-amino-3-phosphonopropionic acid (phosphonoalanine)-degradative operon. Reverse transcriptase-PCR (RT-PCR) experiments demonstrated that all five genes in the operon are transcribed as a single mRNA and that their transcription is induced by phosphonoalanine or phosphonopyruvate. Transcriptional fusions of the regulatory region of the phosphonoalanine degradative operon with the gfp gene were constructed. Expression analysis indicated that the presence of a LysR-type regulator (encoded by the palR gene) is essential for the transcription of the structural genes of the operon. Similar gene clusters were found in the sequenced genomes of six bacterial species from the Alpha-, Beta- and Gammaproteobacteria, and analysis of metagenomic libraries revealed that sequences related to palA are widely spread in the marine environment.

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.

Protein expression profiling during chick retinal maturation: a proteomics-based approach

Background: The underlying pathways that drive retinal neurogenesis and synaptogenesis are still relatively poorly understood. Protein expression analysis can provide direct insight into these complex developmental processes. The aim of this study was therefore to employ proteomic analysis to study the developing chick retina throughout embryonic (E) development commencing at day 12 through 13, 17, 19 and post-hatch (P) 1 and 33 days.

Results: 2D proteomic and mass spectrometric analysis detected an average of 1514 spots per gel with 15 spots demonstrating either modulation or constitutive expression identified via MS. Proteins identified included alpha and beta-tubulin, alpha enolase, B-creatine kinase, gamma-actin, platelet-activating factor (PAF), PREDICTED: similar to TGF-beta interacting protein 1, capping protein (actin filament muscle Z line), nucleophosmin 1 (NPM1), dimethylarginine dimethylaminohydrolase, triosphoaphate isomerase, DJ1, stathmin, fatty acid binding protein 7 (FABP7/B-FABP), beta-synuclein and enhancer of rudimentary homologue.

Conclusion: This study builds upon previous proteomic investigations of retinal development and represents the addition of a unique data set to those previously reported. Based on reported bioactivity some of the identified proteins are most likely to be important to normal retinal development in the chick. Continued analysis of the dynamic protein populations present at the early stages and throughout retinal development will increase our understanding of the molecular events underpinning retinogenesis.