Small RNAs in the peripheral blood discriminate metastasized from non-metastasized seminoma
Molecular Cancer volume 13, Article number: 47 (2014)
We aimed to better discriminate metastasized (lymphogen/occult/both combined) from non-metastasized seminoma based on post-transcriptional changes examined in the peripheral blood.
Total RNAs including small RNAs were isolated from the peripheral blood of patients suffering from metastasized testicular tumours (lymphogen, n = 5, clinical stage IIb/c; occult, n = 5, clinical stage I) and non-metastasized patients (n = 5, clinical stage I). Small RNA next generation sequencing (SOLID, Life Technologies) was employed to examine post-transcriptional changes. We searched for small RNAs showing at least 50 reads and a significant ≥ 2-fold difference using peripheral blood small RNAs of non-metastasized tumours as the reference group. Candidate small RNAs were examined in univariate logistic regression analysis and combinations of two small RNAs were further examined using support vector machines.
On average 1.3x107, 1.2x107 and 1.2x107 small RNA reads were detectable in non-metastasized, lymphogen and occult metastasized seminoma, respectively of which 73-76% remained after trimming. From these between 80-82% represented annotated reads and 7.2-7.8% (1.6-1.7x104) were annotated small RNA tags. Of them 137 small RNAs showed > 50 reads and a ≥ two-fold difference to the reference. In univariate analysis we detected 33-35 different small RNAs which significantly discriminated lymphogen/occult/combined metastasized from non-metastasized seminoma and among these different comparisons it were the same small RNAs in 44-79%. Many combinations of two of these small RNAs completely discriminated metastasized from non-metastasized seminoma irrespective of the metastasis subtype.
Metastasized (either lymphogen or occult) seminoma can be completely discriminated from non-metastasized seminoma with a combination of two small RNAs measured in the peripheral blood.
Testicular tumor, as the most common tumor in young men, is associated with a 5 year survival rate close to 100% in early stages. Pure seminoma are the most frequent histological subtype (55%) nowadays and more than 70% of patients present without visible metastasis at primary staging . Gold standard for primary staging is computed tomography (CT) of the chest, abdomen and pelvis to detect metastases. In truly non-metastasized clinical stage I (cS I) patients are cured by orchidectomy alone, but despite modern staging and classification procedures up to 30% of cS I seminoma patients’ bear occult metastasis in primary staging and relapse after orchidectomy alone [2, 3]. Until today no reliable biological parameters exist and clinical parameter are showing a concordance of 65% only in differentiating occult metastasized stages from non-metastasized seminoma . Identification of occult metastasized patients is one of the main goals to prevent toxicity (e.g. cardiovascular and kidney disease, secondary malignancies and decreased fertility) caused by unnecessary adjuvant treatment or diagnostic procedures (additional radiation exposure due to quarterly CT scans) during follow up .
Recently, other authors started to examine whether a certain set of micro RNAs (miRNA) might be suitable for discriminating between seminoma bearing patients and healthy persons [6–8]. Expression of miRNAs in testicular germ cell cancer is also known to be associated with the histologic subtype  as well as cisplatin resistance [9, 10]. Additionally miRNAs are known to be involved in different processes of metastatic spread in other tumours . Among seminoma bearing patients circulating tumour cells are already detected in the peripheral blood . We wondered whether changes in microRNA expression in the peripheral blood might be able to discriminate metastasized (either lymphogen, occult or a combination of both subtypes) from non-metastasized seminoma. We utilized an agnostic approach investigating the whole genome for any kind of small RNA species suitable to discriminate metastatic stage in seminoma employing next generation sequencing (NGS) on peripheral blood samples drawn at the time of the primary tumour’s diagnosis.
Characteristics of seminoma groups
The average age at diagnosis was 39.1 (+/- 7.2) years for non-metastasized, higher (45.2 years, +/- 10.8) for lymphogen metastasized and lower (32.1 years, +/- 5.4) for occult metastasized seminoma. Primary tumor size was comparable between lymphogen and occult metastasized seminoma (37.8 mm and 38.6 mm, respectively), but smaller (23.8 mm) in non-metastasized seminoma (Table 1).
Per 2.5 ml whole blood we isolated about the same amount (mean) of total RNA from seminoma patients without metastasis (7.7 μg, +/-1.9), lymphogen metastasis (4.9 μg, +/-2.0 μg) and occult metastasis (6.8 μg, +/-2.4, Table 1). The RNA integrity number (RIN) ranged between 7.0-8.6 per group. No DNA contamination could be detected in our samples.
Analysis of small RNA next generation sequencing results
The average total number of reads for lymphogen/occult metastasized and non-metastasized seminoma was 1.3x107, 1.2x107 and 1.2x107 and on average 73-76% remained for further analysis after trimming of the reads (Table 2). From these reads 80-84% appeared annotated reads with 7.2-7.8% (1.6-1.7x104) representing annotated small RNAs. Of them 137 small RNAs showed > 50 reads and a two-fold difference to the reference. In univariate analysis we identified 33/34 and 35 small RNA species which significantly discriminated lymphogen/occult metastasized and both metastasis subtypes combined from non-metastasized seminoma, respectively (Table 3). The overlap of small RNAs separating either lymphogen or occult metastasis from non-metastasized seminoma was 47-48% and increased up to 79% when comparing lymphogen and occult small RNAs candidates with the small RNAs based on the combined metastasis subtypes (Figure 1). We finally employed support vector machines which completely separated lymphogen, occult metastasis and the combined metastasis subtypes from non-metastasized seminoma using a combination of two small RNA species only (Table 4). Altogether, 891, 668 and 87 combinations allowed a complete separation of lymphogen (n = 5), occult metastasis (n = 5) and the combined metastasis subtypes (n = 10) from non-metastasized seminoma (n = 5), respectively.
In our study we aimed to better discriminate metastasized (either lymphogen or occult subtypes or both combined) from non-metastasized seminoma based on small RNA changes examined in the peripheral blood. A whole genome screening on small RNA species was performed employing NGS. We demonstrated each metastasis subtype and its combination to be completely discriminated from non-metastasized seminoma using two small RNAs combined.
MicroRNAs are previously described to be involved in different distant (e.g. migration, angiogenesis and colonization of distant organs) and local (e.g. changes in tumor microenvironment) processes in metastatic cascade. Interestingly, the microRNAs identified in our study are involved in local tumour microenvironmental changes only. For instance, members of the Let-7 family are involved in inflammatory processes as a part of the metastatic cascade. The Let-7 family targets oncogenes such as HMGA2 and KRAS and seem to be a key component in a so called epigenetic switch [13–15]. Cancer associated fibroblasts are involved in tumor formation and progression where miR-15 and miR-16 regulate FGF2 and FGFR1 in prostate cancer  and miR-18 in breast cancer .
Recently other authors demonstrated that microRNA 371-73 cluster and microRNA 302 discriminated seminoma bearing patients from healthy persons [6–8]. We did the next step and examined the suitability of these microRNA to predict the metastasis status, these microRNA 371-73 cluster and miRNA 302 was detectable in our analysis, but appeared not significantly associated with the metastasis status as it is true for protein tumour markers (e.g. HCG) as well.
The complete separation of metastasis from non-metastasized seminoma using a combination of two small RNA species points to the significant diagnostic potential of these biological marker which is in line with previous examinations on the transcriptional level showing the superiority of molecular marker over epidemiological or clinical-histological parameter [18, 19]. After validation of these results on a larger group, these findings could help for therapy decision making in favour of adjuvant chemotherapy or surveillance based on miRNA expression changes predicting metastasis spread, even if metastasis spreads are not visible using a CT scan.
Interesting, the discrimination in our study occurred irrespective of the metastasis subtype and with the same microRNAs in up to 79%. This was expected, since we recently demonstrated lymphogen and occult metastasized seminoma to be indistinguishable on the transcriptional and post-transcriptional level (revised version submitted).
Noteworthy, these changes were examined in the peripheral blood of seminoma patients taken at the time of the primary tumour’s diagnosis. It can be hypothesized that our measurements are linked to seminoma tumour cells circulating in the peripheral blood.
Our study has certain weaknesses such as the low number of cases examined (total n = 15). However, the molecular biological methodology applied provides a deepness, which under financial considerations does not allow large scale studies. Hence, this study provides hints towards certain small RNA species comprising a significant diagnostic potential for prediction of seminoma metastasis based on examinations in the peripheral blood, but certainly these candidate small RNA species have to be examined on a larger independent group for validation purposes. The validation should also include haematogenous metastasized seminoma (clinical Stage III), not included in this study, to prove the usability in this tumor stage as well. Furthermore, we did not adjust the p-values for multiple comparisons, because of the explorative character of this study and the low sample size. Based on experiences with previous work we expect about 10-20% of our candidates being false positives/negatives . However, this study represents a screening approach. About 100 most promising candidate genes will be considered for validation in a following study using another group and another methodology (qRT-PCR).
Metastasized (either lymphogen or occult) seminoma can be completely discriminated from non-metastasized seminoma with a combination of two small RNAs measured in the peripheral blood.
Non-metastasized seminoma (n = 5) received no adjuvant treatment and were free of relapse/progress for at least two years of follow up. Occult metastasized patients (n = 5) presented without visible metastasis at primary staging, received no adjuvant treatment, and developed retroperitoneal tumour progress during the follow up of 12 month. For patients with detectable metastasis at primary staging (n = 5) we focused on clinical stage IIb and IIc, to include lymphogenic metastatic spread only and to provide a high level of diagnostic accuracy (avoiding doubtful lymph nodes). Lymphogen metastasis and non-metastasized patients were matched with occult metastasized patients considering demographic and histological parameters where applicable as well as quality criteria for isolated RNA (Table 1). All patients included in this study were treated between 2008 and 2010 in one testis cancer centre. These 15 patients were selected out of about 300 seminoma patients, according to the criteria mentioned above.
All testicular tumours (n = 15) were examined by an experienced pathologist for histological and TNM classification (Table 1) and pure seminoma was diagnosed in all cases.
Blood sample taking
Whole blood samples (2.5 ml) were taken intraoperative from cubital vein directly into the PAXgene Blood RNA system (BD Diagnostics, PreAnalytiX GmbH, Hombrechtikon, Switzerland). The tube was gently inverted (10 times), settled at room temperature overnight and stored at -20° until use. The Ethical Review Committee of the Medical Association Hamburg approved the study and all human samples were obtained with written informed consent.
After thawing the PAXGene tubes, and washing and centrifugation of the samples the cells became lysed (Proteinase K) followed by adding Lysis/Binding Solution taken from the mirVana Kit (Life Technologies, Darmstadt, Germany). With the mirVana kit total RNA including small RNA species was isolated by combining a Phenol-Chloroform RNA precipitation with further processing using a silica membrane. After several washing procedures DNA residuals became digested on the membrane (RNAse free DNAse Set, Qiagen, Hilden, Germany), washed, RNA was eluted in a collection tube and frozen at -80°C. Quality and quantity of isolated total RNA were measured spectrophotometrically (NanoDrop, PeqLab Biotechnology, Erlangen, Germany). RNA integrity was assessed by the 2100 Agilent Bioanalyser (Life Science Group, Penzberg, Germany) and DNA contamination was controlled by conventional PCR using actin primer.
For analysis, we used only RNA specimens with a ratio of A260/A280 ≥ 2.0 (Nanodrop) and RNA integrity number (RIN) between 6.0-9.3 for small RNA Next Generation Sequencing (IMGM Laboratories, Martinsried, Germany/CeGat, Tübingen, Germany). For RNA samples with RIN below 7.0 or questionable 28/18S bands (indications of RNA degradation) we performed additional checks to exclude the presence of RNAses and true RNA-degradation.
Small RNA next generation sequencing and data analysis
We performed a genome wide small RNA sequencing using the SOLiD5500xl Next Generation Sequencing Technology (Life Technologies, Penzberg, Germany). In brief, total RNA was purified (PureLink microRNA isolation Kit), enriched small RNAs were ligated to SOLiD adaptors, reverse transcribed (SOLiD RT primer and ArrayScript RT), cDNA was purified (MinElute PCR purification Kit, Qiagen), a cut-off size of 60-80 nt was selected (Novex pre-cast gel products, invitrogen), cDNA was in-gel amplified and samples were barcoded using SOLiD 3′Barcode primer at the same time. Amplified cDNA was purified (PureLink PCR Micro Kit, Life Technologies) and used in emulsion PCR (SOLiD EZ Bead). Thereafter, the emulsion was broken to recover enriched beads and the so-called di-base probes were used by the SOLiD system in the sequencing-by-ligation procedure. Beside the SOLiD5500xl inherent software (LifeScope) used for image and signal processing, software CLC Genomics Workbench 5.1 (CLC bio) was used for clustering, counting, and annotation of all generated 75 bp reads. After discarding reads without adaptor sequence and being shorter than 15 bp (trimming) reads were assigned to known microRNAs (Sanger miRBase release 18, http://www.mirbase.org/) and known non-coding RNAs (ensembl database Homo_sapiens.GRCh37.67.ncrna.fa, http://www.ensembl.org). Small RNAs showing a significant and ≥ 2-fold differences in gene expression among groups and at least 50 reads were further analysed on their ability to separate both groups using univariate logistic regression analysis  and combinations of two small RNAs were examined employing support vector machines (linear Kermel).
Copy deoxyribonucleic acids
- CT scan:
micro ribonucleic acid
Next generation sequencing
Quantitative real time polymerase chain reaction
RNA integrity number
Tumor nodus metastasis classification.
Ruf CG, Isbarn H, Wagner W, Fisch M, Matthies C, Dieckmann KP: Changes in epidemiologic features of testicular germ cell cancer: age at diagnosis and relative frequency of seminoma are constantly and significantly increasing(). Urol Oncol. 2014, 32 (1): 33-e1-6. doi:10.1016/j.urolonc.2012.12.002. Epub 2013 Feb 6
Krege S, Beyer J, Souchon R, Albers P, Albrecht W, Algaba F, Bamberg M, Bodrogi I, Bokemeyer C, Cavallin-Stahl E, Classen J, Clemm C, Cohn-Cedermark G, Culine S, Daugaard G, De Mulder PH, De Santis M, de Wit M, de Wit R, Derigs HG, Dieckmann KP, Dieing A, Droz JP, Fenner M, Fizazi K, Flechon A, Fosså SD, del Muro XG, Gauler T, Geczi L: European consensus conference on diagnosis and treatment of germ cell cancer: a report of the second meeting of the European Germ Cell Cancer Consensus group (EGCCCG): part I. Eur Urol. 2008, 53: 478-496. 10.1016/j.eururo.2007.12.024
Albers P, Albrecht W, Algaba F, Bokemeyer C, Cohn-Cedermark G, Fizazi K, Horwich A, Laguna MP: [EAU guidelines on testicular cancer: 2011 update. European Association of Urology]. Actas Urol Esp. 2012, 36: 127-145.
Ruf CG, Khalili-Harbi N, Sachs S, Isbarn H, Wagner W, Matthies C, Meineke V, Fisch M, Chun FK, Abend M: The search for biomarkers of metastatic seminoma. J Urol. 2013, 190 (3): 1046-1051. doi:10.1016/j.juro.2013.04.022. Epub 2013 Apr 10
Kollmannsberger C, Tyldesley S, Moore C, Chi KN, Murray N, Daneshmand S, Black P, Duncan G, Hayes-Lattin B, Nichols C: Evolution in management of testicular seminoma: population-based outcomes with selective utilization of active therapies. Ann Oncol. 2011, 22: 808-814. 10.1093/annonc/mdq466
Dieckmann KP, Spiekermann M, Balks T, Flor I, Loning T, Bullerdiek J, Belge G: MicroRNAs miR-371-3 in serum as diagnostic tools in the management of testicular germ cell tumours. Br J Cancer. 2012, 107: 1754-1760. 10.1038/bjc.2012.469
Gillis AJ, Stoop HJ, Hersmus R, Oosterhuis JW, Sun Y, Chen C, Guenther S, Sherlock J, Veltman I, Baeten J, van der Spek PJ, de Alarcon P, Looijenga LH: High-throughput microRNAome analysis in human germ cell tumours. J Pathol. 2007, 213: 319-328. 10.1002/path.2230
Palmer RD, Murray MJ, Saini HK, van Dongen S, Abreu-Goodger C, Muralidhar B, Pett MR, Thornton CM, Nicholson JC, Enright AJ, Coleman N, : Malignant germ cell tumors display common microRNA profiles resulting in global changes in expression of messenger RNA targets. Cancer Res. 2010, 70: 2911-2923. 10.1158/0008-5472.CAN-09-3301
Port M, Glaesener S, Ruf C, Riecke A, Bokemeyer C, Meineke V, Honecker F, Abend M: Micro-RNA expression in cisplatin resistant germ cell tumor cell lines. Mol Cancer. 2011, 10: 52- 10.1186/1476-4598-10-52
Liu X, Yu J, Jiang L, Wang A, Shi F, Ye H, Zhou X: MicroRNA-222 regulates cell invasion by targeting matrix metalloproteinase 1 (MMP1) and manganese superoxide dismutase 2 (SOD2) in tongue squamous cell carcinoma cell lines. Cancer Genomics Proteomics. 2009, 6: 131-139.
Profumo V, Gandellini P: MicroRNAs: cobblestones on the road to cancer metastasis. Crit Rev Oncog. 2013, 18: 341-355. 10.1615/CritRevOncog.2013007182
Ruf CG, Nastaly P, Becker P, Isbarn H, Honecker F, Pantel K, Riethdorf S, Hoeppner D, Fisch M, Wagner W, Ahyai S: Circulating tumor cells can be detected in patients with testicular germ cell tumors. J Urol. 2013, 189: e289-
Mayr C, Hemann MT, Bartel DP: Disrupting the pairing between let-7 and Hmga2 enhances oncogenic transformation. Science. 2007, 315: 1576-1579. 10.1126/science.1137999
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014
Iliopoulos D, Hirsch HA, Struhl K: An epigenetic switch involving NF-kappaB, Lin28, Let-7 MicroRNA, and IL6 links inflammation to cell transformation. Cell. 2009, 139: 693-706. 10.1016/j.cell.2009.10.014
Musumeci M, Coppola V, Addario A, Patrizii M, Maugeri-Sacca M, Memeo L, Colarossi C, Francescangeli F, Biffoni M, Collura D, Giacobbe A, D'Urso L, Falchi M, Venneri MA, Muto G, De Maria R, Bonci D: Control of tumor and microenvironment cross-talk by miR-15a and miR-16 in prostate cancer. Oncogene. 2011, 30: 4231-4242. 10.1038/onc.2011.140
Yu Z, Willmarth NE, Zhou J, Katiyar S, Wang M, Liu Y, McCue PA, Quong AA, Lisanti MP, Pestell RG: microRNA 17/20 inhibits cellular invasion and tumor metastasis in breast cancer by heterotypic signaling. Proc Natl Acad Sci U S A. 2010, 107: 8231-8236. 10.1073/pnas.1002080107
Ruf CG, Linbecker M, Port M, Riecke A, Schmelz HU, Wagner W, Meineke V, Abend M: Predicting metastasized seminoma using gene expression. BJU Int. 2012, 110: E14-E20. 10.1111/j.1464-410X.2011.10778.x
Port M, Wang Y, Schmelz HU, Pottek T, Meineke V, Ruf C, Abend M: A gene signature of primary tumor identifies metastasized seminoma. Urol Oncol. 2011, 29: 764-773. 10.1016/j.urolonc.2009.08.008
Port M, Seidl C, Ruf CG, Riecke A, Meineke V, Abend M: Reliable and sample saving gene expression analysis approach for diagnostic tool development. Health Phys. 2012, 103: 159-168.
Baggerly KA, Deng L, Morris JS, Aldaz CM: Overdispersed logistic regression for SAGE: modelling multiple groups and covariates. BMC Bioinform. 2004, 5: 144-10.1186/1471-2105-5-144.
We would like to thank S. Senf for his skilful technical assistance. The German Ministry of Defence supported this work.
The German Ministry of Defence supported this work. All authors declare that they have no competing interests.
CR has developed the study design, acquisition and interpretation of data and drafted the manuscript, DD has substantially contributed in the acquisition of data, MP and H-US have substantially contributed in the study design, interpretation of data and in revising the manuscript, WW and CM have made substantial contributions to acquisition of data and in revising the manuscript, BM-M has performed the statisitcal analysis of data and revised the manuscript, VM has made substantial contributions to technical procedures, data analysis and revised the manuscript, MA has developed the study design, supervised all technical/laboratory procedures and quality checks, analysis and interpretation of data and drafted the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ruf, C.G., Dinger, D., Port, M. et al. Small RNAs in the peripheral blood discriminate metastasized from non-metastasized seminoma. Mol Cancer 13, 47 (2014). https://doi.org/10.1186/1476-4598-13-47