Whole-Transcriptome profiling of formalin-fixed, paraffin-embedded renal cell carcinoma by RNA-seq
© Li et al.; licensee BioMed Central. 2014
Received: 2 July 2014
Accepted: 26 November 2014
Published: 11 December 2014
Formalin-fixed paraffin-embedded (FFPE) tissue samples are routinely archived in the course of patient care and can be linked to clinical outcomes with long-term follow-up. However, FFPE tissues have degraded RNA which poses challenges for analyzing gene expression. Next-generation sequencing (NGS) is rapidly becoming accepted as an effective tool for measuring gene expressions for research and clinical use. However, the feasibility of NGS has not been firmly established when using FFPE tissue.
We optimized strategies for whole transcriptome sequencing (RNA-seq) using FFPE tissue. Ribosomal RNA (rRNA) was successfully depleted by competitive hybridization using the Ribo-zero™ Kit (Epicentre Biotechnologies), and rRNA sequence content was less than one percent for each library. Gene expression measured by FFPE RNA-seq was compared to two different standards: RNA-seq from fresh frozen (FF) tissue and quantitative PCR (qPCR). Both FF and FFPE tumors were sequenced on an Illumina Genome Analyzer IIX with an average of 10 million reads. The distribution of FPKMs (fragments per kilobase of exon per million fragments mapped) and number of detected genes were similar between FFPE and FF. RNA-seq expressions from FF and FFPE samples from the same renal cell carcinoma (RCC) correlated highly (r = 0.919 for tumor 1 and r = 0.954 for tumor 2). On hierarchical cluster analysis, samples clustered by patient identity rather than method of preservation. TaqMan qPCR of 424 RCC-related genes correlated highly with FFPE RNA-seq expressions (r = 0.775 for FFPE tumor 1, r = 0.803 for FFPE tumor 2). Expression fold changes were considered, to assess biologic relevance of gene expressions. Expression fold changes between FFPE tumors (tumor 1/tumor 2) correlated well when comparing qPCR and RNA-seq (r = 0.890). Expression fold changes between tumors from different risk groups (our high risk RCC/The Cancer Genome Atlas, TCGA, low risk RCC) also correlated well when comparing RNA-seq from FF and FFPE tumors (r = 0.887).
FFPE RNA-seq provides reliable genes expression data, comparable to that obtained from fresh frozen tissue. It represents a useful tool for discovery and validation of biomarkers.
KeywordsFormalin-fixed paraffin-embedded (FFPE) qPCR Renal cell carcinoma (RCC) RNA-seq Gene expression
RNA expression profiling may lead to the discovery of molecular markers for disease diagnosis, assessing prognosis, and targeting with drugs. Quantitative (qPCR) has been the “gold standard” for measuring gene expressions [1, 2] due to its high sensitivity and specificity, reproducibility, and large dynamic range [3–5]. However, next-generation sequencing is rapidly becoming accepted  as an effective and more versatile tool for measuring gene expressions for research and clinical use [7–9]. Compared to qPCR, the major advantages of next generation sequencing (NGS) include the ability to analyze a sample’s whole transcriptome in an unbiased way, to discover novel transcripts, and to detect gene fusions, which are common in cancer .
The optimal tissue for RNA-seq is fresh frozen (FF) tissue with high quality RNA. Unfortunately, frozen tumors are not widely available because they are costly to collect and maintain. However, formalin-fixed paraffin-embedded (FFPE) tissue samples are routinely archived in the course of patient care and can often be linked to clinical outcomes with long-term follow-up. Unfortunately, FFPE tissues yield relatively low quantities of degraded RNA. A small number of studies have reported using FFPE tissue for whole transcriptome mRNA expression profiling [11–14]. Here, we characterize the performance of RNA-seq on FFPE renal tumors by comparing results to RNA-seq on FF renal tumors and qPCR, which is considered the “gold standard” for measuring gene expression.
In this study, transcriptome-wide RNA-seq was successfully performed on matching FFPE and FF clear cell renal cell carcinomas (ccRCCs). The expression profiles generated from FFPE and FF tumors correlated well. The tumor RNA was also assessed by qPCR using the OpenArray® NT Cycler system, which we’ve previously validated for use with FFPE tumors . RNA quantities measured by RNA-seq and qPCR also correlated well. Expression fold changes between our tumors and tumors from The Cancer Genome Atlas (TCGA)  correlated well when expressions from FFPE and FF were compared, suggesting that FFPE RNA-seq can provide biologically meaningful information. We establish the feasibility of using RNA-seq with FFPE tissue and recommend its use in future large-scale RNA-seq studies.
Expression levels determined by qPCR
Library preparation and RNA-seq
Average cluster density for the 4 libraries was approximately 800 K/mm2. The average depth of sequencing for each library was 10 million reads. The distribution of FPKMs for all the genes was similar between FFPE and FF (Additional file 2: Figure S2). The number of genes detected by RNA-seq was comparable between FF and FFPE (Additional file 3: Table S1). The distribution of detected genes was also comparable between FF and FFPE (Additional file 4: Figure S3). RNA-seq quality metrics showed that for each library, greater than 70% of the reads were uniquely mapped to the genome with less than 1% being rRNA. The strand specificity was higher than 80%. (Additional file 3: Table S2).
Comparison of FFPE RNA-seq to established RNA profiling strategies
FFPE RNA-seq results contain biologically meaningful information
The past decade has seen accelerated development of genomic and transcriptomic techniques. A wealth of biologic insights has come from data generated from thousands of tumor samples. Both the cost and turn-around time have dropped dramatically, which in turn, has put this powerful technology within reach of most wet-lab scientists. RNA-seq has become the method of choice for transcriptome profiling. In contrast to microarray and qPCR technologies, RNA-seq can identify novel transcripts, examine all RNA species, and identify alternative splices and mutations.
RCC has been studied using RNA-seq [7, 22]. These studies provide transcriptomic profiles for understanding disease development and prognosis. However, all of these studies used FF tumors as the resource of RNA, and FF tumors are relatively scarce. On the other hand, FFPE samples are widely available because they are routinely stored by all pathology departments. Importantly, FFPE samples often come from patients who have had years of clinical follow-up. Therefore, this clinical annotation can provide the phenotype needed to understand the genomic and molecular data that can now be readily generated.
There are several studies of RNA-seq on FFPE tissue [11–13]. However, all these studies have important limitations. One of the studies used poly-A selection to capture mRNA and exclude ribosomal RNA. This results in the loss of all non-mRNAs and because FFPE RNA is highly degraded, poly-A capture provides only a limited sequence from the 3’ end of the mRNA. Therefore, no information is usually available on alternative splicing and many mutations will be missed . Another study used lungs from two patients who died of influenza to perform FFPE RNA-seq. However, the study profiled viral RNA and not human RNA . A third study from Genomic Health Inc. used FFPE breast cancers; however, there was no comparison with RNA-seq performed with FF breast cancer . This study did establish that the expression of a prognostic signature based on a limited number of genes can be determined using FFPE qPCR or FFPE RNA-seq.
We performed whole-transcriptome RNA-seq on FFPE and FF samples. We demonstrated that gene expressions measured by RNA-seq using FFPE tissue correlate well with FF RNA-seq. Our findings are consistent with a recent report using a similar next-generation sequencing platform as ours to compare RNA from paired FF and FFPE tissue from a variety of tumor types to report a correlation coefficient of 0.90 . This study examined RNA-seq in 38 paired FF and FFPE normal and cancer tissues from bladder, prostate, liver, colon and tonsil. In our study, we examined clear cell RCC. In addition to RNA-seq, we also used qPCR as an established standard for expression profiling and showed good correlation with FFPE RNA-seq. By assessing fold change between patient tumors and performing hierarchical cluster analysis, we showed that FFPE RNA-seq contains potentially meaningful biologic information that can be used for basic discoveries as well as patient care.
A limitation of our study is that we only examined two patient tumors. However, RNA-seq outputs covering the entire transcriptome provide ample data for a robust analysis. Furthermore, use of multiple technology platforms and use of the publically available TCGA results strengthen our conclusion. FFPE tumors, routinely stored by medical centers, are associated with a wealth of clinical data and follow-up information. However, to take advantage of these annotated tissues, RNA-seq needs to be applied to large numbers of archival tissue. Our study, along with other recent, independent reports, establishes the feasibility of RNA-seq for these future studies.
We demonstrated that gene expressions measured by RNA-seq using FFPE tissue correlate well with two different gold standards: FF RNA-seq and qPCR.
FFPE RNA-seq also provided biological information: The expression profiles from the matched FF and FFPE tissue clustered together. The fold change between high risk FF tumor and TCGA low risk tumors correlated well with the fold change between high risk FFPE tumor and TCGA low risk tumors.
Our study suggests that FFPE RNA-seq can be used for basic discoveries as well as clinical applications.
Patients and tissues
Cedars-Sinai Medical Center Biobank (Los Angeles, CA) provided 4 renal tumors from 2 matching pairs of FFPE and FF ccRCCs. The study was approved by Cedars-Sinai Medical Center’s Institutional Review Board (IRB). For prospectively collected tumor tissue, patients gave written informed consent for data and tissue collection. For FFPE tissue, which already existed in our biobank, IRB provided an exemption for obtaining informed consent. The clinical information associated with each tumor is listed in Additional file 3: Table S1. The FFPE tissues had been stored for about two years at room temperature prior to use in this study.
Our method for RNA extraction has previously been described [15, 16, 23]. Briefly, for FFPE tissue, RNA from six 10-μm sections were extracted using the MasterPure™ RNA Purification Kit (Epicentre Biotechnologies, Madison, WI) with minor modification: 200 ug Proteinase K was added to the Cell Lysis Solution and used for a 3-hr incubation. For FF tissue, RNA was extracted from 100 mg of snap-frozen tissues using TRIzol®. RNA from both FFPE and FF was treated with 20 units DNase I. Complete digestion of genomic DNA was confirmed by TaqMan qPCR for ACTB. The final RNA concentration and purity were measured using a NanoDrop ND-2000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE). The integrity of RNA was measured with the Eukaryote Total RNA Nano Assay of the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). 200 ng of RNA was loaded for each sample.
RNA-seq Library preparation and sequencing
Ribosomal RNA was depleted from 4 micrograms of total RNA using the Ribo-Zero™ Magnetic Kit (Epicentre Biotechnologies, Madison, WI) following the manufacturer’s procedure. Sequencing libraries for whole transcriptome analysis were prepared using ScriptSeq™ v2 RNA-Seq Library Preparation Kit (Epicentre Biotechnologies, Madison, WI) following the manufacturer’s procedure. After 3′-terminal tagging, the di-tagged cDNA was purified using MinElute PCR Purification Kit (Qiagen, Valencia, CA). Four 6-base index sequences were used to prepare bar-coded libraries for sample multiplexing (ScriptSeq™ Index PCR Primers (Set 1); Epicentre Biotechnologies, Madison, WI). PCR was carried out through 13 cycles to generate the second strand of cDNA, incorporate barcodes, and amplify libraries. The amplified libraries were size selected by a solid phase reversible immobilization, paramagnetic bead-based process (Agencourt AMPure XP system; Beckman Coulter Genomics, Danvers, MA). The beads: PCR product ratio was 1.8:1.
Libraries were quantified by Qubit® dsDNA HS Assay (Life Technologies, Carlsbad, CA). 15ng libraries were visualized on a 2% E-Gel (Life Technologies, Carlsbad, CA). The four RNA-seq Libraries were then multiplexed in a single lane on a paired end flow cell; clonal amplification was performed using an Illumina cBot (Illumina Inc.; San Diego, CA). Sequencing-by-Synthesis chemistry was performed on an Illumina Genome Analyzer IIX using 35 bp single end sequencing chemistry.
BCL files were converted to FASTQ files using BCLconverter. The FASTQ files were aligned against human genome reference version 19 (or hg19), using TopHat2 (http://ccb.jhu.edu/software/tophat/index.shtml) software package. Transcripts were assembled using Cufflink2 (http://cufflinks.cbcb.umd.edu/). The number of sequenced reads that align to a gene of interest was conventionally called a tag count. To make tag count comparable among samples, we transformed the tag count into FPKM/RPKM (Reads per kilobase per million reads mapped) for normalization, which is widely used in RNA-seq analysis. FPKM value was shown in Additional file 3: Table S4.
To get RNA-seq quality metrics, bedtools utility (http://bedtools.readthedocs.org) was used to overlap TopHat2 mappings with mapping quality of at least 20 and defined exons and introns with a minimum overlap fraction of .8. To reconcile differing exon and intron annotations for a single gene due to alternative splicing, a read was defined as intronic if it was (1) within an intron and (2) overlapped with no exons; reads overlapping any exon were counted as exonic. The bedtools utility was used to determine concordance between the strands of reads mapping to exons and the strands of the exons. To determine the presence of rRNA in the samples, reads were mapped to rRNA human sequences (NR_003287/RNA28S5, NR_003286/ RNA18S5, NR_003285/RNA5-8S5, and NR_023379/ RNA5S17) using the BWA utility (http://bio-bwa.sourceforge.net). Reads which mapped to an rRNA sequence with a mapping quality of at least 20 were considered to be rRNA derived.
TaqMan qPCR using the open array platform
Reverse transcription (RT) was performed using the High Capacity cDNA Reverse Transcription Kit (Life Technologies, Grand Island, NY) following the manufacturer’s procedure. Each RT reaction contained 150ng of total RNA, 1 μl of 10× RT buffer, 0.5 μl of 25× dNTP mixture, 1 μl of 10× random reverse primers, 1 μl of 10× gene-specific reverse primers (1μM) and 0.5 μl of MultiScribe RT (50 U/μl). The 10 μl reactions were incubated in a Life Technologies Thermocycler for 10 min at 25°C, 2 hours at 37°C and then held at 4°C. The same primers were used for the pre-amplification and the TaqMan qPCR. The 3’end primers used in the PCR were used for the RT.
Preamplification was performed using the TaqMan® PreAmp Master Mix Kit (Life Technologies, Grand Island, NY). Each reaction included 2.5 μl of 2× TaqMan® PreAmp Master Mix, 1.25 μl of 0.09× pooled TaqMan assays (TaqMan primers and probe) and 1.25 μl of cDNA. The reactions were incubated in a Life Technologies Thermocycler for 10 min at 95°C following by 13 cycles of 95°C for 15 seconds and 60°C for 4 min, and then held at 4°C. 0.09× pooled Taqman assays were prepared by combining equal volumes of each 20× Taqman assay (needed for PCR on the Openarray®, 218 assays for each set). Each cDNA was preamplified on two sets of 218 pooled assays.
The two sets of OpenArray plates were custom made. The preamplified products were diluted 1:10; 10 μl of diluted cDNA was mixed with 10 μl of TaqMan Real-time PCR OpenArray Master Mix (Life Technologies, Grand Island, NY). The pre-amplified cDNA samples were dispensed using the Accufill System onto the corresponding OpenArray plate containing 218 gene assays. Twelve cDNA samples were tested simultaneously per plate, with 36 samples per qPCR run on the Openarray® NT Cycler system. Post-acquisition data processing generated fluorescence amplification for each assay, from which cycle threshold (CT) were computed. CT value was shown in Additional file 3: Table S5.
TCGA RNA-seq data mining and analysis
Clinical information and Level 3 RNA-seq FPKM data were retrieved from the TCGA data portal (https://tcga-data.nci.nih.gov/tcga/). Our tumors were clinically localized tumors that posed high risk for metastasis. To compare these tumors to low risk-RCC using expression fold change, all pT1N0M0, grade 1 or 2 ccRCCs were selected from TCGA (n = 53).
To assess the distribution of FPKMs for all the genes, box plots were generated (Additional file 2: Figure S2). The goal was to identify the most highly expressed genes. The FPKM distribution was different between TCGA data and our data. One of the reasons is that we used different platform for sequencing. TCGA used Hiseq2000. We used GAIIx. Another reason is that we sequenced different cohorts. Based on the distribution of FPKMs, an average FPKM >10 was used to select genes from the TCGA samples. For our samples, an average FPKM >5 was used to select genes. A total of 3,195 overlapping genes common to both datasets were identified and used to calculate expression fold changes between our tumors and TCGA tumors.
Statistics and bioinformatics
Fold changes in gene expression were calculated: for qPCR results, log2 fold change between sample a and b = Δ CT for b - Δ CT for a, and for RNA-seq results, fold change between sample a and b = log2 (FPKM for a/FPKM for b). Complete-linkage Hierarchical clustering was performed using Gene Cluster 3.0 software. The similarity (distance) between gene expression data was defined using Pearson correlation (uncentered). The result was viewed using Java TreeView software.
Availability of supporting data
The TCGA sequence data and clinical data discussed in this publication are available from the TCGA Data Portal .
The supplemental data was available from LabArchive with DOI of 10.6070/H4Q23X6N DOI: 10.6070/H4Q23X6N#doi.
Deep sequencing files have been deposited into Sequence Read Archive under accession number SRP050335.
This study was supported in part by NIH (R01CA133072-01), Winnick Family Foundation, and Lauhlere Family Fund.
- Mackay IM, Arden KE, Nitsche A: Real-time PCR in virology. Nucleic Acids Res. 2002, 30: 1292-1305. 10.1093/nar/30.6.1292.PubMed CentralPubMedView ArticleGoogle Scholar
- van’t Veer LJ, Bernards R: Enabling personalized cancer medicine through analysis of gene-expression patterns. Nature. 2008, 452: 564-570. 10.1038/nature06915.PubMedView ArticleGoogle Scholar
- Arya M, Shergill IS, Williamson M, Gommersall L, Arya N, Patel HR: Basic principles of real-time quantitative PCR. Expert Rev Mol Diagn. 2005, 5: 209-219. 10.1586/1473718.104.22.168.PubMedView ArticleGoogle Scholar
- Wilhelm J, Pingoud A: Real-time polymerase chain reaction. Chembiochem. 2003, 4: 1120-1128. 10.1002/cbic.200300662.PubMedView ArticleGoogle Scholar
- Wong ML, Medrano JF: Real-time PCR for mRNA quantitation. Biotechniques. 2005, 39: 75-85. 10.2144/05391RV01.PubMedView ArticleGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralPubMedView ArticleGoogle Scholar
- Cancer Genome Atlas Research Network: Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013, 499: 43-49. 10.1038/nature12222.View ArticleGoogle Scholar
- Tuch BB, Laborde RR, Xu X, Gu J, Chung CB, Monighetti CK, Stanley SJ, Olsen KD, Kasperbauer JL, Moore EJ, Broomer AJ, Tan R, Brzoska PM, Muller MW, Siddiqui AS, Asmann YW, Sun Y, Kuersten S, Barker MA, De La Vega FM, Smith DI: Tumor transcriptome sequencing reveals allelic expression imbalances associated with copy number alterations. PLoS One. 2010, 5: e9317-10.1371/journal.pone.0009317.PubMed CentralPubMedView ArticleGoogle Scholar
- Guo G, Sun X, Chen C, Wu S, Huang P, Li Z, Dean M, Huang Y, Jia W, Zhou Q, Tang A, Yang Z, Li X, Song P, Zhao X, Ye R, Zhang S, Lin Z, Qi M, Wan S, Xie L, Fan F, Nickerson ML, Zou X, Hu X, Xing L, Lv Z, Mei H, Gao S, Liang C, et al: Whole-genome and whole-exome sequencing of bladder cancer identifies frequent alterations in genes involved in sister chromatid cohesion and segregation. Nat Genet. 2013, 45: 1459-1463. 10.1038/ng.2798.PubMedView ArticleGoogle Scholar
- Maher CA, Kumar-Sinha C, Cao X, Kalyana-Sundaram S, Han B, Jing X, Sam L, Barrette T, Palanisamy N, Chinnaiyan AM: Transcriptome sequencing to detect gene fusions in cancer. Nature. 2009, 458: 97-101. 10.1038/nature07638.PubMed CentralPubMedView ArticleGoogle Scholar
- Sinicropi D, Qu K, Collin F, Crager M, Liu ML, Pelham RJ, Pho M, Dei Rossi A, Jeong J, Scott A, Ambannavar R, Zheng C, Mena R, Esteban J, Stephans J, Morlan J, Baker J: Whole transcriptome RNA-Seq analysis of breast cancer recurrence risk using formalin-fixed paraffin-embedded tumor tissue. PLoS One. 2012, 7: e40092-10.1371/journal.pone.0040092.PubMed CentralPubMedView ArticleGoogle Scholar
- Beck AH, Weng Z, Witten DM, Zhu S, Foley JW, Lacroute P, Smith CL, Tibshirani R, van de Rijn M, Sidow A, West RB: 3’-end sequencing for expression quantification (3SEQ) from archival tumor samples. PLoS One. 2010, 5: e8768-10.1371/journal.pone.0008768.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiao YL, Kash JC, Beres SB, Sheng ZM, Musser JM, Taubenberger JK: High-throughput RNA sequencing of a formalin-fixed, paraffin-embedded autopsy lung tissue sample from the 1918 influenza pandemic. J Pathol. 2013, 229: 535-545. 10.1002/path.4145.PubMed CentralPubMedView ArticleGoogle Scholar
- Hedegaard J, Thorsen K, Lund MK, Hein AM, Hamilton-Dutoit SJ, Vang S, Nordentoft I, Birkenkamp-Demtröder K, Kruhøffer M, Hager H, Knudsen B, Andersen CL, Sørensen KD, Pedersen JS, Ørntoft TF, Dyrskjøt L: Next-generation sequencing of RNA and DNA isolated from paired fresh-frozen and formalin-fixed paraffin-embedded samples of human cancer and normal tissue. PLoS One. 2014, 9 (5): e98187-10.1371/journal.pone.0098187.PubMed CentralPubMedView ArticleGoogle Scholar
- Li P, Grigorenko E, Funari V, Enright E, Zhang H, Kim HL: Evaluation of a high-throughput, microfluidics platform for performing TaqMan™ qPCR using formalin-fixed paraffin-embedded tumors. Bioanalysis. 2013, 5: 1623-1633. 10.4155/bio.13.125.PubMedView ArticleGoogle Scholar
- Glenn ST, Jones CA, Liang P, Kaushik D, Gross KW, Kim HL: Expression profiling of archival renal tumors by quantitative PCR to validate prognostic markers. Biotechniques. 2007, 43: 639-640. 10.2144/000112562. 642-643, 647PubMedView ArticleGoogle Scholar
- Thompson E, Burt AD, Barker CE, Kirby JA, Brain JG: Development of a robust protocol for gene expression analysis using formalin-fixed, paraffin-embedded liver transplant biopsy specimens. J Clin Pathol. 2013, 66: 815-818. 10.1136/jclinpath-2013-201548.PubMed CentralPubMedView ArticleGoogle Scholar
- Takahashi M, Rhodes DR, Furge KA, Kanayama H, Kagawa S, Haab BB, Teh BT: Gene expression profiling of clear cell renal cell carcinoma: gene identification and prognostic classification. Proc Natl Acad Sci U S A. 2001, 98: 9754-9759. 10.1073/pnas.171209998.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao H, Ljungberg B, Grankvist K, Rasmuson T, Tibshirani R, Brooks JD: Gene expression profiling predicts survival in conventional renal cell carcinoma. PLoS Med. 2006, 3: e13-10.1371/journal.pmed.0030013.PubMed CentralPubMedView ArticleGoogle Scholar
- Vasselli JR, Shih JH, Iyengar SR, Maranchie J, Riss J, Worrell R, Torres-Cabala C, Tabios R, Mariotti A, Stearman R, Merino M, Walther MM, Simon R, Klausner RD, Linehan WM: Predicting survival in patients with metastatic kidney cancer by gene-expression profiling in the primary tumor. Proc Natl Acad Sci U S A. 2003, 100: 6958-6963. 10.1073/pnas.1131754100.PubMed CentralPubMedView ArticleGoogle Scholar
- Brannon AR, Reddy A, Seiler M, Arreola A, Moore DT, Pruthi RS, Wallen EM, Nielsen ME, Liu H, Nathanson KL, Ljungberg B, Zhao H, Brooks JD, Ganesan S, Bhanot G, Rathmell WK: Molecular stratification of clear cell renal cell carcinoma by consensus clustering reveals distinct subtypes and survival patterns. Genes Cancer. 2010, 1: 152-163. 10.1177/1947601909359929.PubMed CentralPubMedView ArticleGoogle Scholar
- Sato Y, Yoshizato T, Shiraishi Y, Maekawa S, Okuno Y, Kamura T, Shimamura T, Sato-Otsubo A, Nagae G, Suzuki H, Nagata Y, Yoshida K, Kon A, Suzuki Y, Chiba K, Tanaka H, Niida A, Fujimoto A, Tsunoda T, Morikawa T, Maeda D, Kume H, Sugano S, Fukayama M, Aburatani H, Sanada M, Miyano S, Homma Y, Ogawa S: Integrated molecular analysis of clear-cell renal cell carcinoma. Nat Genet. 2013, 45: 860-867. 10.1038/ng.2699.PubMedView ArticleGoogle Scholar
- Glenn ST, Head KL, Teh BT, Gross KW, Kim HL: Maximizing RNA yield from archival renal tumors and optimizing gene expression analysis. J Biomol Screen. 2010, 15: 80-85. 10.1177/1087057109355059.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.