|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Kompetenzzentrum fuer Drug Design und Target Monitoring, Luebeck, Germany.
2 Institut fuer Molekulare Medizin,3
Klinik und Poliklinik für Urologie, and4
Institut fuer Medizinische Biometrie und Statistik, UK-S-H, Campus Luebeck, Luebeck, Germany.
aAddress correspondence to this author at: Institut fuer Molekulare Medizin, Universitaetsklinikum Schleswig-Holstein, Campus Luebeck, Ratzeburger Allee 160, 23538 Luebeck, Germany. Fax 49-451-500-2729; e-mail warnecke{at}imm.uni-luebeck.de.
| Abstract |
|---|
|
|
|---|
Methods: RNA isolated from urinary cell pellet, cell-depleted fraction, and whole urine was quantified by reverse transcription quantitative–PCR. The donor-dependent influence of urine background on individual steps of the standardized procedure was analyzed using an external RNA standard. Using a test set of samples from 61 patients with bladder cancer and 37 healthy donors, we compared 4 putative RNA tumor markers identified in whole urine with 5 established, tissue-derived RNA tumor markers for the detection of bladder cancer.
Results: Of the markers analyzed by this system, the RNA ratio of v-ets erythroblastosis virus E26 oncogene homolog 2 (avian; ETS2) to urokinase plasminogen activator (uPA) enabled the most specific (100%) and sensitive (75.4%) detection of bladder cancer from whole urine, with an area under the curve of 0.929 (95% CI 0.882–0.976).
Conclusions: The described methodology for RNA marker detection in urine appears to be clinically applicable. The ratio of ETS2 mRNA to uPA mRNA in urine is a potential marker for bladder cancer.
| Introduction |
|---|
|
|
|---|
To avoid differences in processing, the use of body fluids without any pretreatment (e.g., centrifugation) is most desirable for a noninvasive diagnostic test. This application is challenging, since body fluids from different donors differ with respect to the amount, origin, and integrity of cells and nucleic acids. The cells may contain RNA that is fragmented by necrotic and apoptotic processes that also give rise to cell-free RNA in urine(4). Cell-free RNA is resistant to the activity of urinary RNases because of packaging into apoptotic bodies. In addition, cell death–independent processes have been described for the generation of extracellular RNA(5). To avoid exclusion of patients, however, a diagnostic routine application based on body fluids ideally is not dependent on the integrity of the RNA.
Despite these obstacles, the detection of RNA tumor markers in urine has been reported as an emerging tool for noninvasive tumor diagnosis(6). One of the most frequently studied markers for bladder cancer is the mRNA of human telomerase reverse transcriptase (hTERT),2 because its concentration correlates with telomerase activity, which is absent in most human somatic cells but detectable in 85% of human tumors(7). RNA extracts from tumor material(8)(10), bladder washings(11)(12), and exfoliated cells(13) have been used for RT-qPCR–based hTERT mRNA quantification for the detection of bladder cancer, resulting in reported diagnostic sensitivities of 50% to 100%. When urine sediment was used as sample material in a larger prospective study, however, hTERT mRNA expression was often found to be near the detection limit, resulting in a sensitivity of 55%(14).
In this study, we aimed to develop a RT-qPCR–based test of whole urine for robust and sensitive diagnosis of bladder cancer in a clinical setting. Toward this goal, we 1st analyzed factors influencing the RNA isolation and quantification procedure. We then applied a standardized process to compare RNA tumor markers originating from our screening experiments with established markers for bladder cancer.
| Materials and Methods |
|---|
|
|
|---|
To investigate selected tumor markers, we used whole urine from 98 donors. For a control group, we chose healthy donors who best matched the age and sex distribution of the bladder cancer group. Detailed clinical information is shown in Table 1 in the Data Supplement that accompanies the online version of this article athttp://www.clinchem.org/content/vol53/issue12. An overview of the experimental procedures is given in Fig. 1
.
|
preparation of the external rna standard (rnaluc)
We synthesized truncated luciferase RNA (RNALUC) by in vitro runoff transcription using T7 RNA-polymerase. The T7-transcript was prepared from a purified EcoRV-restricted pet24a(+) vector (Promega) fragment containing the Photinus pyralis luciferase gene. The reaction mixture (50 µL) for T7 transcription contained 1x reaction buffer (MBI Fermentas), 2 mmol/L each nucleotide triphosphate, 40 units RNase inhibitor (Ambion), 30 units T7 RNA-Polymerase (MBI Fermentas), and 0.8 µg DNA template. After transcription of RNALUC (1570 nucleotides), the DNA template was digested using DNase I (Ambion). The reaction mixture was gel-filtrated through a NAP-5 column (Amersham), extracted by chloroform-phenol, and precipitated with ethanol. Integrity and size of RNALUC were controlled by denaturing agarose gel electrophoresis. RNA concentration was determined by absorbance at 260 nm (1 A260 = 40 mg RNA/L) of a dilution series of the purified RNALUC in vitro transcript using NanoDrop ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies). We used mean values to calculate RNALUC copy numbers.
processing of urine samples
Spontaneously voided urine was adjusted to a final concentration of 3 mol/L guanidinium thiocyanate (GTC), 0.025 mol/L sodium acetate, and 0.25% N-lauroylsarcosine. For the monovette-based sample collection, we transferred 3.54 g GTC powder into a 10-mL urine monovette (Sarstedt) and adjusted the stamp to the 7-mL mark of the monovette. Further handling of the monovette was performed according to manufacturers instructions. After the GTC powder was dissolved in the collected urine, the sample was transported to the analytic laboratory. Adjustment to 0.025 mol/L sodium acetate and 0.25% N-lauroylsarcosine and addition of 1 mol/L HEPES (pH 7) to a final volume of 10 mL were carried out before storage at –80 °C.
preparation of urinary cell pellet
The cellular fraction was obtained by centrifugation of 3 mL voided urine at 400g for 5 min. We discarded the supernatant and dissolved the pellet in a special lysis buffer(15). The lysate was frozen in liquid nitrogen and stored at –80 °C until RNA isolation.
preparation of cell-depleted urine
Cell-depleted urine was obtained by passing voided urine through a 5-µm filter (Sartorius). The urine filtrate was treated as described for whole urine.
standard procedure for the isolation and reverse transcription of total rna
Before RNA isolation, 107 copies of RNALUC were added. Total RNA was isolated using the RNeasy Midi Kit (Qiagen) according to the manufacturers instructions, except that special lysis buffer was used instead of RLT buffer. The procedure included an on-column digestion of genomic DNA with DNase I. RNA was eluted twice with 160 µL nuclease-free H2O and lyophilized. The RNA pellets were resolved in 20 µL nuclease-free H2O. We used 10 µL RNA extract for cDNA synthesis and non-reverse transcription reaction. Additionally, 107 copies of RNALUC were directly reverse-transcribed to verify applied copies. A detailed description of the cDNA synthesis protocol is provided in the online Data Supplement.
quantitative pcr and data analysis
TaqMan-based quantification of cDNA was always run in triplicate, and the non–reverse transcription reaction was run in duplicate. The protocol is described in detail in the online Data Supplement, together with primer sequences and amplicon characteristics (see Table 2 in the online Data Supplement). We performed data analysis with SDS 2.1 software (Applied Biosystems). The threshold cycle values of amplified targets were transformed into absolute RNA copy numbers using standard curves, allowing an absolute quantification of target RNA copy numbers.
assay characterization
RT-qPCR.
RNALUC (107 copies) were directly reverse-transcribed, and we quantified RNALUC cDNA by use of qPCR in triplicate. We performed independent qPCR assay runs using a standard curve for every run to determine the interrun variation of RNALUC quantification. To determine the donor-dependent influence of the urine background on the qPCRs, 5 x 105 copies of RNALUC were added to RNA extracts from different donors. Copy numbers of RNALUC were determined by the standard RT-qPCR procedure.
Interrun imprecision of the complete RNA isolation and quantification procedure.
We added 107 copies of in vitro–transcribed RNALUC to 8 aliquots of spontaneously voided urine of a single donor after adjustment to 3 mol/L GTC. For each sample, we determined copy numbers of RNALUC, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein large P0 (RPLP0), and ubiquitin C (UBC) by independent isolation and quantification procedures.
Intradonor variability.
For daytime profiles, spontaneously voided urine was collected in the morning, at midday, and in the evening. We added 107 copies of RNALUC to whole urine after adjustment to 3 mol/L GTC. RNA was isolated in duplicate and quantified in triplicate via RT-qPCR. Copy numbers of RNALUC, GAPDH, RPLP0, and UBC were determined by standardized RT-qPCR procedure.
Interdonor variability.
For the investigation of donor-dependent influence of urine background on RNA recovery, we added 107 copies of RNALUC to whole urine, cell-depleted urine, or the cellular fraction before RNA isolation and after adjustment to 3 mol/L GTC. RNA was isolated in duplicate, and we determined copy numbers of RNALUC by the standard procedure for RNA isolation and quantification. In addition, we determined copy numbers of endogenous GAPDH.
statistical analysis and software
The statistical analysis of RNA marker ratios is described in detail in the online Data Supplement. The diagnostic power of selected markers was analyzed by ROC curves. For the comparison of 2 datasets (comparison of RNALUC yields and housekeeping gene ratios), we used the nonparametric paired Wilcoxon test. For all statistical tests, 2-sided P values
0.05 were considered statistically significant.
Statistical analyses were performed using SAS statistical software (version 9.1, SAS Institute), SPSS statistical software (version 13.0, SPSS GmbH Software), R language and environment for statistical computing (version 2.3.1, R Foundation for Statistical Computing), and STATA/SE statistical software (version 9.1, Stata).
| Results |
|---|
|
|
|---|
assay characterization
The assay characteristics are summarized in Tables 1
and 2
.
|
|
interrun imprecision
Interrun imprecision (CV) for the isolation and quantification process of RNALUC from 8 aliquots of whole urine was 0.16 (Table 1
). Endogenous RNAs for GAPDH, RPLP0, and UBC were detected with a comparable imprecision (CVGAPDH 0.15; CVRPLP0 0.18; CVUBC 0.16).
interdonor variability
We next determined the donor-dependent influence of urine background on assay imprecision (Table 2
). The overall RNALUC yield of the RNA isolation and quantification process was quantified for different fractions (cell-depleted, cellular fraction, and whole urine) prepared from urine of different donors (n = 21). The yield of RNALUC was best for the cellular fraction and lowest for cell-depleted urine (see Fig. 1 in the online Data Supplement). RNALUC yields for the cellular fraction (mean 27%) were significantly higher (Wilcoxon test P = 0.002) than for whole urine (mean 16%). The highest variability is introduced by the RNA isolation step, since downstream steps turned out to be robust (CVRT-qPCR 0.14; Table 2
). Accordingly, for all subsequent experiments, RNALUC-normalized copy numbers were determined to allow a comparison of urine samples from different donors.
intradonor variability
The amount of endogenous housekeeping genes such as GAPDH RNA in urine varied by 4 orders of magnitude (Table 2
). To allow a comparison between individual samples, we calculated the ratios of housekeeping gene RNAs. Systematic differences in RNA composition were analyzed by investigation of housekeeping gene ratios at 3 times of the day (n = 10).
Ratios were calculated from RNA copy numbers of GAPDH, RPLP0, and UBC. The intradonor variability for GAPDH:RPLP0 (Fig. 2A
) was smaller (mean CV 0.25; range 0.09–0.61) than for GAPDH:UBC (mean CV 0.41; range 0.07–0.71). Pairwise comparison using nonparametric Wilcoxon test provides no evidence for a significant influence of time of day on the GAPDH:UBC ratio (P = 0.13). Only the comparison of the GAPDH:RPLP0 ratio from morning and midday showed a significant difference (P = 0.05). To obtain conclusive data independently of patients compliance with instructions, we used only spontaneously voided urine (midday) instead of morning urine for further experiments.
|
rna composition in urine fractions
The GAPDH:UBC ratios in total differed significantly from the cellular fraction (Wilcoxon test, P = 0.01; n = 21) in contrast to the GAPDH:RPLP0 RNA ratios (Wilcoxon test, P = 0.07; n = 21). The cell-depleted fraction differed significantly from whole urine and the cellular fraction for the 2 RNA ratios (Wilcoxon test, P <0.001 for GAPDH:UBC and GAPDH:RPLP0, n = 21).
We next determined the amount of cell-free RNA that contributes to the detected RNA in whole urine (Fig. 2B
). Quantification of housekeeping gene RNA in the cell-depleted urine fractions (n = 21) revealed high amounts of GAPDH RNA (mediancell-depleted = 1.06 x 106 copy numbers/mL; medianwhole urine = 1.89 x 106 copy numbers/mL) and UBC RNA (mediancell-depleted = 5.12 x 105 copy numbers/mL; medianwhole urine = 5.58 x 106 copy numbers/mL). In some samples, the copy numbers of GAPDH RNA detected in whole urine consisted almost exclusively of cell-free RNA (see Fig. 2B
, samples 2, 11, 15, and 21). The cell-free RNA portion of the whole urine fraction was higher for GAPDH RNA (mean 46.9%) compared with UBC RNA (mean 16.9%).
These results provided a rationale for the reinvestigation of different cell-based tumor markers for their applicability using whole urine. Their diagnostic performance was analyzed together with markers derived from RT-qPCR–based screening of RNA isolated from whole urine (data not shown).
investigation of rna tumor markers in whole urine
RNA was isolated from urine of 37 healthy donors and 61 patients with bladder cancer. A set of 8 different RNA tumor markers (Table 3
) containing putative as well as established bladder cancer markers and the housekeeping gene GAPDH were quantified as specified in Fig. 1
, right panel.
|
Multivariate analysis of individual marker ratios was performed to test the ability to separate the 98 donors into those with cancer and those without. We selected 3 important marker ratios using classification trees. In a logistic regression model with backward selection including the selected variables and all interactions, only the ratio of v-ets erythroblastosis virus E26 oncogene homolog 2 (avian; ETS2) to urokinase plasminogen activator (uPA) allowed a statistically significant separation. To confirm the result of the logistic regression model, we performed 2 nonparametric diagonal linear discriminant analyses, one with the 3 important marker ratios and another with ETS2:uPA. The classification results of the diagonal linear discriminant analyses were not different (P >0.05). Therefore, the RNA ratio of ETS2:uPA was shown to be the only independent marker for the detection of bladder cancer. The ETS2:uPA ratios in whole urine of patients with bladder cancer and healthy controls are presented in Fig. 3A
. ROC curves were calculated that represent the diagnostic power of the RNA marker combinations. ROC analysis (Fig. 3B
) of the study population by use of the ETS2:uPA ratio revealed an area under the curve (AUC) of 0.929 (95% CI 0.882–0.976), indicating the strong diagnostic power of the test. Setting the specificity at 100% (cutoff value 0.96), a sensitivity of 75.4% was achieved. For the group of low-grade tumors, we determined a sensitivity of 53.9%. Sensitivity can be enhanced to 79.9% for low-grade tumors and 89.1% for the group of all tumors, if the specificity is reduced to 89.2%. As demonstrated by the ROC analysis in Fig. 3B
, the use of GAPDH-normalized ETS2 RNA and GAPDH-normalized uPA RNA as individual markers did not turn out to be informative. The AUC values for the other marker combinations tested are listed in Table 3 in the online Data Supplement.
|
| Discussion |
|---|
|
|
|---|
The initial step of sampling, i.e., the transfer of urine into the monovette containing GTC powder, does not require any laboratory equipment. It is therefore simple to implement into hospital routine, thereby avoiding increased processing time. Inactivation of RNA-degrading enzymes occurs rapidly by the dissolved GTC. The urine samples were thus stabilized immediately, transported, and frozen within 6 h after voiding. This is an important step toward the development of a standard method for the sampling of urine, comparable to EDTA monovettes for the collection of blood samples. To obtain a high informative rate, the RNA isolates have to be concentrated. We have recently demonstrated the feasibility of ethanol precipitation for this purpose(15). To allow maximum standardization, this step was substituted by lyophilization.
The use of the "1st void of the day" (morning) urine is quite common in diagnostic studies because of the higher content of nucleic acids. Investigation of intraday variance indicated that the RNA ratio of GAPDH:RPLP0 showed a significant difference between morning and midday urine. This RNA ratio–dependent intraday variation indicates an advantage of spontaneously voided urine (midday), since the results are independent of the patients compliance with instructions.
Using whole urine, differences in yield of the complete RNA isolation and quantification process ranged from 0.7% to 47.9%. Accordingly, RNA tumor markers with low abundance have to be excluded from this application. This holds particularly for the analysis of hTERT RNA, in which our analysis indicated insufficient abundance for reliable quantification (data not shown). This is in line with results from the prospective study of Weikert et al.(14), in which hTERT mRNA expression in exfoliated cells was often found to be near the detection limit. In contrast to hTERT mRNA, ETS2 and uPA RNA copy numbers/mL were sufficiently high in whole urine (medianETS2 = 4.64 x 105 copy numbers/mL; medianuPA= 4.84 x 105 copy numbers/mL).
The emerging interest in using RNA-tumor marker detection in whole urine for the diagnosis of other tumor entities(16)(17) accentuates the need to investigate the variability and source of RNA species in whole urine. Our experiments indicate that RNA tumor markers derived from gene expression analysis of cells or tumor material (e.g., hTERT, UPK1A, HTATIP2) cannot be transferred unconditionally to RT-qPCR–based analysis of whole urine. Based on our data, one major reason is the presence of cell-free RNA in significant amounts. When GAPDH is used for normalization, the high concentrations of cell-free GAPDH RNA, in addition to the well-known overexpression of GAPDH RNA in tumor tissue, mask the diagnostic power of cellular markers, as demonstrated here for uPA/GAPDH (AUC 0.359; see Fig. 3B
). In contrast to extracellular DNA in urine(18)(19), the presence of specific extracellular RNA has not been reported. In addition to apoptosis and necrosis, cell death–independent generation of extracellular RNA has been demonstrated(5).
RT-qPCR array–based screening experiments (data not shown) using whole urine and the subsequent evaluation of putative tumor markers revealed the ratio of ETS2:uPA to be a suitable marker for the detection of bladder cancer. When the specificity was set at 100%, a sensitivity of 75.4% was achieved. This result outperforms the sensitivity (28%–76%) of urinary cytology(20) as the standard noninvasive method and most of the molecular bladder cancer markers when the specificity is set at 100%(21)(22). Stratification by tumor grade (Fig. 3A
) indicated a higher sensitivity for high-grade tumors (81.3%, grade 2–3), an observation that has been described for other urine-based RNA markers(14)(23)(24).
ETS2 is a member of the ETS family of transcription factors that regulate a variety of biological processes(25). ETS2 RNA and protein concentrations were found to be increased in breast cancer samples compared with normal tissue(26). Buggy et al.(26) showed a significant correlation of protein expression between ETS2 and uPA, which contains ETS2 binding sites in its promoter region. The role of the uPA system in tumorigenesis has been intensively studied(27). In bladder cancer, urinary protein concentrations of uPA [e.g., (28)] as well as the uPA mRNA content in tissue [e.g.,(29)], have been suitable for tumor diagnostics. The experimental data of our study population indicate a higher ETS2 RNA concentration compared with uPA in the case of bladder cancer, resulting in an increased ETS2:uPA RNA ratio.
In conclusion, we describe a fully standardized process for the isolation and quantification of RNA tumor markers from total urine. The new tumor marker ratio of ETS2:uPA is promising for the detection of bladder cancer with high specificity and sensitivity in a clinical setting.
| Acknowledgments |
|---|
Financial disclosures: None declared.
Acknowledgments: We thank Drs. C. Höppner, M. Horn, T. Menke, S. Thomas, and J. Träder for collecting samples, A. Krüger and B. Thode for technical assistance, and G. Sczakiel for critically reading the manuscript.
| Footnotes |
|---|
2 Human genes: hTERT, human telomerase reverse transcriptase; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; RPLP0, ribosomal protein, large, P0; UBC, ubiquitin C; ETS2, v-ets erythroblastosis virus E26 oncogene homolog 2 (avian); uPA, urokinase plasminogen activator; HTATIP2, HIV-1 Tat interactive protein 2, 30kDa; UPK1A, uroplakin 1A. ![]()
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |