Mining tissue microarray data to uncover combinations of biomarker expression patterns that improve intermediate staging and grading of clear cell renal cell cancer

PURPOSE: Tumor stage and nuclear grade are the most important prognostic parameters of clear cell renal cell carcinoma (ccRCC). The progression risk of ccRCC remains difficult to predict particularly for tumors with organ-confined stage and intermediate differentiation grade. Elucidating molecular pathways deregulated in ccRCC may point to novel prognostic parameters that facilitate planning of therapeutic approaches. EXPERIMENTAL DESIGN: Using tissue microarrays, expression patterns of 15 different proteins were evaluated in over 800 ccRCC patients to analyze pathways reported to be physiologically controlled by the tumor suppressors von Hippel-Lindau protein and phosphatase and tensin homologue (PTEN). Tumor staging and grading were improved by performing variable selection using Cox regression and a recursive bootstrap elimination scheme. RESULTS: Patients with pT2 and pT3 tumors that were p27 and CAIX positive had a better outcome than those with all remaining marker combinations. A prolonged survival among patients with intermediate grade (grade 2) correlated with both nuclear p27 and cytoplasmic PTEN expression, as well as with inactive, nonphosphorylated ribosomal protein S6. By applying graphical log-linear modeling for over 700 ccRCC for which the molecular parameters were available, only a weak conditional dependence existed between the expression of p27, PTEN, CAIX, and p-S6, suggesting that the dysregulation of several independent pathways are crucial for tumor progression. CONCLUSIONS: The use of recursive bootstrap elimination, as well as graphical log-linear modeling for comprehensive tissue microarray (TMA) data analysis allows the unraveling of complex molecular contexts and may improve predictive evaluations for patients with advanced renal cancer.

Selecting control genes for RT-QPCR using public microarray data

BACKGROUND: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e.g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. RESULTS: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/~vpopovic/research/ CONCLUSION: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

Meta-analyses qualify metzincins and related genes as acute rejection markers in renal transplant patients

Definition of acute renal allograft rejection (AR) markers remains clinically relevant. Features of T-cell-mediated AR are tubulointerstitial and vascular inflammation associated with excessive extracellular matrix (ECM) remodeling, regulated by metzincins, including matrix metalloproteases (MMP). Our study focused on expression of metzincins (METS), and metzincins and related genes (MARGS) in renal allograft biopsies using four independent microarray data sets. Our own cases included normal histology (N, n = 20), borderline changes (BL, n = 4), AR (n = 10) and AR + IF/TA (n = 7). MARGS enriched in all data sets were further examined on mRNA and/or protein level in additional patients. METS and MARGS differentiated AR from BL, AR + IF/TA and N in a principal component analysis. Their expression changes correlated to Banff t- and i-scores. Two AR classifiers, based on METS (including MMP7, TIMP1), or on MARGS were established in our own and validated in the three additional data sets. Thirteen MARGS were significantly enriched in AR patients of all data sets comprising MMP7, -9, TIMP1, -2, thrombospondin2 (THBS2) and fibrillin1. RT-PCR using microdissected glomeruli/tubuli confirmed MMP7, -9 and THBS2 microarray results; immunohistochemistry showed augmentation of MMP2, -9 and TIMP1 in AR. TIMP1 and THBS2 were enriched in AR patient serum. Therefore, differentially expressed METS and MARGS especially TIMP1, MMP7/-9 represent potential molecular AR markers.

A subset of metzincins and related genes constitutes a marker of human solid organ fibrosis

Metzincins and functionally related genes play important roles in extracellular matrix remodeling both in healthy and fibrotic conditions. We recently presented a transcriptomic classifier consisting of 19 metzincins and related genes (MARGS) discriminating biopsies from renal transplant patients with or without interstitial fibrosis/tubular atrophy (IF/TA) by virtue of gene expression measurement (Roedder et al., Am J Transplant 9:517-526, 2009). Here we demonstrate that the same algorithm has diagnostic value in non-transplant solid organ fibrosis. We used publically available microarray datasets of 325 human heart, liver, lung, kidney cortex, and pancreas microarray samples (265 with fibrosis, 60 healthy controls). Expression of nine commonly differentially expressed genes was confirmed by TaqMan low-density arrays (Applied Biosystems, USA) in 50 independent archival tissue specimens with matched histological diagnoses to microarray patients. In separate and in combined, integrated microarray data analyses of five datasets with 325 samples, the previously published MARGS classifier for renal post-transplant IF/TA had a mean AUC of 87% and 82%, respectively. These data demonstrate that the MARGS gene panel classifier not only discriminates IF/TA from normal renal transplant tissue, but also classifies solid organ fibrotic conditions of human pancreas, liver, heart, kidney, and lung tissue samples with high specificity and accuracy, suggesting that the MARGS classifier is a cross-platform, cross-organ classifier of fibrotic conditions of different etiologies when compared to normal tissue.

An integrated nano-scale approach to profile miRNAs in limited clinical samples

Profiling miRNA expression in cells that directly contribute to human disease pathogenesis is likely to aid the discovery of novel drug targets and biomarkers. However, tissue heterogeneity and the limited amount of human diseased tissue available for research purposes present fundamental difficulties that often constrain the scope and potential of such studies. We established a flow cytometry-based method for isolating pure populations of pathogenic T cells from bronchial biopsy samples of asthma patients, and optimized a high-throughput nano-scale qRT-PCR method capable of accurately measuring 96 miRNAs in as little as 100 cells. Comparison of circulating and airway T cells from healthy and asthmatic subjects revealed asthma-associated and tissue-specific miRNA expression patterns. These results establish the feasibility and utility of investigating miRNA expression in small populations of cells involved in asthma pathogenesis, and set a precedent for application of our nano-scale approach in other human diseases. The microarray data from this study (Figure 7) has been submitted to the NCBI Gene Expression Omnibus (GEO; http://ncbi.nlm.nih.gov/geo) under accession no. GSE31030.

Comparison of the gene expression profiles from normal and Fgfrl1 deficient mouse kidneys reveals downstream targets of Fgfrl1 signaling

Fgfrl1 (fibroblast growth factor receptor-like 1) is a transmembrane receptor that is essential for the development of the metanephric kidney. It is expressed in all nascent nephrogenic structures and in the ureteric bud. Fgfrl1 null mice fail to develop the metanephric kidneys. Mutant kidney rudiments show a dramatic reduction of ureteric branching and a lack of mesenchymal-to-epithelial transition. Here, we compared the expression profiles of wildtype and Fgfrl1 mutant kidneys to identify genes that act downstream of Fgfrl1 signaling during the early steps of nephron formation. We detected 56 differentially expressed transcripts with 2-fold or greater reduction, among them many genes involved in Fgf, Wnt, Bmp, Notch, and Six/Eya/Dach signaling. We validated the microarray data by qPCR and whole-mount in situ hybridization and showed the expression pattern of candidate genes in normal kidneys. Some of these genes might play an important role during early nephron formation. Our study should help to define the minimal set of genes that is required to form a functional nephron.

Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

The sequence of addition of IL-3 and IgE-dependent/IgE-independent stimuli regulates RANKL expression in human basophils

Background: Receptor Activator of Nuclear Factor kappaB Ligand (RANKL), a member of the TNF superfamily, contributes to the imbalance of bone resorption and immunoregulation in rheumatoid arthritis. In mice, collagen induced arthritis was exacerbated by IL-3 and anti-IgER antibodies, two mediators activating basophils that are known as effector cells of allergy. Interestingly, our unpublished microarray data revealed that IL-3 induces RANKL mRNA in human basophils. Here we further investigate under which conditions human basophils express surface and/or soluble RANKL. Methods: One part of purified human basophils was co-stimulated with IL-3 and either IgE-dependent or IgE-independent stimuli. The other part of purified basophils was first primed with IL-3 and subsequently triggered with IgE-dependent or IgE-independent stimuli. Expression of surface and soluble RANKL were detected by flow cytometry, ELISA and real-time PCR. Results: By flow cytometry we show that IL-3 induces de novo expression of surface RANKL on human basophils in a time and dose dependent manner. Co-stimulation of basophils with IL-3 and an IgE-dependent stimulus reduces IL-3-induced expression of surface RANKL in a dose dependent manner while IgE-independent stimuli have no effect. In contrast, both IgE-dependent and IgE-independent stimuli enhance expression of surface and soluble RANKL in basophils that were first primed with IL-3 and then triggered. Real-time PCR analysis shows that surface hRANKL1 and soluble hRANKL3 are induced by IL-3 and reduced by co-stimulation with IL-3 and an IgE-dependent stimulus and thus confirms our flow cytometry data. Conclusion: RANKL expression in human basophils is not only dependent on IL-3 and IgE-dependent/IgE-independent stimuli but also on the sequence of their addition to cell culture. Based on our data, we suggest that basophils might have previously unidentified functions in bone resorption or immunoregulation via RANKL.

Loss of tumor suppressor mir-203 mediates overexpression of LIM and SH3 Protein 1 (LASP1) in high-risk prostate cancer thereby increasing cell proliferation and migration

Several studies have linked overexpression of the LIM and SH3 domain protein 1 (LASP1) to progression of breast, colon, liver, and bladder cancer. However, its expression pattern and role in human prostate cancer (PCa) remained largely undefined. Analysis of published microarray data revealed a significant overexpression of LASP1 in PCa metastases compared to parental primary tumors and normal prostate epithelial cells. Subsequent gene-set enrichment analysis comparing LASP1-high and -low PCa identified an association of LASP1 with genes involved in locomotory behavior and chemokine signaling. These bioinformatic predictions were confirmed in vitro as the inducible short hairpin RNA-mediated LASP1 knockdown impaired migration and proliferation in LNCaP prostate cancer cells. By immunohistochemical staining and semi-quantitative image analysis of whole tissue sections we found an enhanced expression of LASP1 in primary PCa and lymph node metastases over benign prostatic hyperplasia. Strong cytosolic and nuclear LASP1 immunoreactivity correlated with PSA progression. Conversely, qRT-PCR analyses for mir-203, which is a known translational suppressor of LASP1 in matched RNA samples revealed an inverse correlation of LASP1 protein and mir-203 expression. Collectively, our results suggest that loss of mir-203 expression and thus uncontrolled LASP1 overexpression might drive progression of PCa.

Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia.

BACKGROUND It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Transcriptional response to cardiac injury in the zebrafish: systematic identification of genes with highly concordant activity across in vivo models

BACKGROUND Zebrafish is a clinically-relevant model of heart regeneration. Unlike mammals, it has a remarkable heart repair capacity after injury, and promises novel translational applications. Amputation and cryoinjury models are key research tools for understanding injury response and regeneration in vivo. An understanding of the transcriptional responses following injury is needed to identify key players of heart tissue repair, as well as potential targets for boosting this property in humans. RESULTS We investigated amputation and cryoinjury in vivo models of heart damage in the zebrafish through unbiased, integrative analyses of independent molecular datasets. To detect genes with potential biological roles, we derived computational prediction models with microarray data from heart amputation experiments. We focused on a top-ranked set of genes highly activated in the early post-injury stage, whose activity was further verified in independent microarray datasets. Next, we performed independent validations of expression responses with qPCR in a cryoinjury model. Across in vivo models, the top candidates showed highly concordant responses at 1 and 3 days post-injury, which highlights the predictive power of our analysis strategies and the possible biological relevance of these genes. Top candidates are significantly involved in cell fate specification and differentiation, and include heart failure markers such as periostin, as well as potential new targets for heart regeneration. For example, ptgis and ca2 were overexpressed, while usp2a, a regulator of the p53 pathway, was down-regulated in our in vivo models. Interestingly, a high activity of ptgis and ca2 has been previously observed in failing hearts from rats and humans. CONCLUSIONS We identified genes with potential critical roles in the response to cardiac damage in the zebrafish. Their transcriptional activities are reproducible in different in vivo models of cardiac injury.

Intestinal current measurement for diagnostic classification of patients with questionable cystic fibrosis: validation and reference data

In questionable cystic fibrosis (CF), mild or monosymptomatic phenotypes frequently cause diagnostic difficulties despite detailed algorithms. CF transmembrane conductance regulator (CFTR)-mediated ion transport can be studied ex vivo in rectal biopsies by intestinal current measurement (ICM).

3D Intervertebral Disc Localization and Segmentation from MR Images by Data-Driven Regression and Classification

In this paper we propose a new fully-automatic method for localizing and segmenting 3D intervertebral discs from MR images, where the two problems are solved in a unified data-driven regression and classification framework. We estimate the output (image displacements for localization, or fg/bg labels for segmentation) of image points by exploiting both training data and geometric constraints simultaneously. The problem is formulated in a unified objective function which is then solved globally and efficiently. We validate our method on MR images of 25 patients. Taking manually labeled data as the ground truth, our method achieves a mean localization error of 1.3 mm, a mean Dice metric of 87%, and a mean surface distance of 1.3 mm. Our method can be applied to other localization and segmentation tasks.