Skip to main content

Advertisement

Limited utility of qPCR-based detection of tumor-specific circulating mRNAs in whole blood from clear cell renal cell carcinoma patients

Abstract

Background

RNA sequencing data is providing abundant information about the levels of dysregulation of genes in various tumors. These data, as well as data based on older microarray technologies have enabled the identification of many genes which are upregulated in clear cell renal cell carcinoma (ccRCC) compared to matched normal tissue. Here we use RNA sequencing data in order to construct a panel of highly overexpressed genes in ccRCC so as to evaluate their RNA levels in whole blood and determine any diagnostic potential of these levels for renal cell carcinoma patients.

Methods

A bioinformatics analysis with Python was performed using TCGA, GEO and other databases to identify genes which are upregulated in ccRCC while being absent in the blood of healthy individuals. Quantitative Real Time PCR (RT-qPCR) was subsequently used to measure the levels of candidate genes in whole blood (PAX gene) of 16 ccRCC patients versus 11 healthy individuals. PCR results were processed in qBase and GraphPadPrism and statistics was done with Mann-Whitney U test.

Results

While most analyzed genes were either undetectable or did not show any dysregulated expression, two genes, CDK18 and CCND1, were paradoxically downregulated in the blood of ccRCC patients compared to healthy controls. Furthermore, LOX showed a tendency towards upregulation in metastatic ccRCC samples compared to non-metastatic.

Conclusions

This analysis illustrates the difficulty of detecting tumor regulated genes in blood and the possible influence of interference from expression in blood cells even for genes conditionally absent in normal blood. Testing in plasma samples indicated that tumor specific mRNAs were not detectable. While CDK18, CCND1 and LOX mRNAs might carry biomarker potential, this would require validation in an independent, larger patient cohort.

Peer Review reports

Background

In the United States, it is estimated that 65,340 new cases and 14,970 deaths from kidney cancer will occur in 2018 [1]. ccRCC is the most common renal malignancy, accounting for around 80% of the cases [2]. Together with papillary and chromophobe carcinoma, it comprises 2% of all cancers worldwide [3]. Renal cell carcinoma (RCC) incidence increases markedly with age, peaking at 50–70 years, with males being affected twice more frequently than females [4]. The major risk factors for RCC include excess body weight, hypertension and cigarette smoking [5] and associations have also been made with different lifestyle, dietary, occupational and environmental factors [6]. Primary RCC displays no early clinical symptoms as most renal masses remain asymptomatic until the late stages of the disease, with over 50% of all cases of RCC being discovered by chance during imaging studies for other comorbidities [7, 8]. Only 10% of patients have the classical triad of symptoms: hematuria, flank pain and weight loss. Around 25% of RCCs have already metastasized by the time of diagnosis [9]. RCC is mostly unresponsive to conventional chemotherapy and radiation, which is the main reason for treatment failures [10, 11]. The gold standard for the management of renal masses is nephrectomy, in spite of which approximately 30% of patients develop recurrence or metastases [12, 13], which require systemic therapies and are associated with high mortality.

As current prognostic models based on conventional clinicopathological and imaging data have limited accuracy, new biomarkers are needed for early detection, improved diagnostics and the prediction of the clinical outcome of patients with RCC [14,15,16,17]. The ideal biomarker or biomarker panel should have high specificity, sensitivity, and reproducibility. Plasma, serum, and urine have recently gained interest in the field of cancer biomarker discovery. They represent potential sources of valuable biomarkers, containing proteins, DNA, and various RNA species, with blood being especially suitable in terms of kidney disease and low invasiveness. Steady progress in the field is being made, however to date, none of the identified ccRCC biomarkers have been clinically validated [18].

RNA circulating in blood is highly degraded (usually less than 100 bp in length [19]) and even after the introduction of systems that enable the stabilization and storage of whole blood mRNA (e.g. the PAXgene platform) studies tend to be limited to shorter RNA subspecies, or those protected from degradation due to their specific structure or association with proteins or membranous vesicular structures such as exosomes. Analyzing ccRCC biomarkers in urine would be particularly convenient, however this field is far less fruitful compared to blood studies. Similarly as for blood, urine is problematic in terms of RNA detection because of the presence of RNAses, but also because of PCR inhibition [20], which is steering the focus into the analysis of shorter RNA subspecies. Especially in the case of microRNA (miRNA), liquid biopsy has expanded from use in plasma to other bodily fluids in an increasing number of malignancies, making rapid progress since 2008 [21]. Along with miRNA, the use of Circulating Tumor Cells (CTCs), cell free DNA (cfDNA), and more recently circular RNA (circRNA) and long non-coding RNAs (lncRNA), are proving to be far more viable strategies, as for most tumors there is somewhat sporadic progress in detecting tumor-derived mRNA in blood and associating it with cancer prognosis, for example as described here [22,23,24]. A further issue that complicates tumor-derived RNA detection in blood is the uncertainty regarding the exact RNA origin, i.e. whether it is derived from solid tumor or CTCs, and what percentage of source cells are living and actively secreting RNA as opposed to undergoing apoptosis [25].

The first of the five proposed stages in biomarker development is the comparison of tumor with nontumor tissue [26]. Here techniques such as microarrays and more recently RNAseq are used to assess gene expression, while protein expression profiles are based on immunohistochemistry and mass spectroscopy, with the goal of discovering genes displaying dysregulation (usually overexpression in tumor compared to normal tissue). This phase is followed by the development of a clinical assay utilizing blood for non-invasive screening. Of course, the blood levels of selected genes do not have to precicely mirror the expression in tissue, as a result of e.g. the specific rate of mRNA release from cancer tissue into blood. The approach of using mRNA expression of tumor tissue as a starting point and analyzing the levels od respective transcripts in blood by RT-qPCR has been previously used with success, resulting in promising assays deserving of clinical validation. A recent study showed the validation of an RT-PCR assay based on prostate-specific RNA in whole blood from patients with metastatic castration-resistant prostate cancer (mCRPC) [27]. Several databases were consulted to select a panel of genes that were overexpressed in prostate tissue but showed no detection in peripheral blood mononuclear cells (PBMC). This was followed by RT-PCR analysis of blood samples of cancer patients and volunteers, resulting in the establishment of a 5-gene panel that enhances and complements the previously established CTC enumeraration assay. Similarly, in another study focusing on early detection of colorectal cancer [28], meta-analysis of microarray data was used to identify RNAs with highest diffeential expression between cancer tissue and normal blood samples. Subsequent RT-qPCR analysis revealed that blood expression of 3 specific genes shows promising sensitivity and specificity with respect to detection of this cancer.

In this study, TCGA database was used as a starting point to identify genes which are most highly overexpressed in the tissue of ccRCC patients, after which a subset containing genes that according to other databases have no blood expression was evaluated by qPCR in whole blood samples from ccRCC patients and healthy individuals. While RNA transcripts of some of these genes had good detectability in blood, none of the genes were significantly up-regulated in blood from ccRCC patients and two genes paradoxically displayed downregulation.

Methods

Patients and samples

The staging and grading of the tumor samples were done according to the 2002 TNM classification and the Fuhrman grading system [29, 30]. The ccRCC tissue samples were obtained during partial or radical nephrectomy at the University Hospital Charité in Berlin in 2011 and blood samples in the period between 2010 and 2016. Tissue samples were frozen in liquid nitrogen directly after surgical resection and stored at − 80 °C until RNA extraction. They came from tumor and matched normal tissue of 3 male patients without diagnosed metastasis (ages: 47–71; tumor stages: 2 x pT1, and pT3; grading: G1, G2, G3). PAXgene blood samples were obtained from 27 individuals and included 16 ccRCC samples, out of which 10 were non-metastatic (8 male and 2 female patients; median age 70, range 47–84 years; tumor staging: 1x pT1, 2x pT2, 7x pT3; grading: 2x G1, 7x G2, 1x G3) and 6 metastatic: (5 male and 1 female patients; median age 67, range 47–72 years; tumor staging: 6x pT3; grading: 5x G2, 1x G3). In total there were 11 samples without diagnosed cancer, 4 from patients suffering from non-cancer kidney illnesses, and 7 healthy volunteers (7 male and 4 female; median age 47, range 29–80 years).

Bioinformatics analysis

The first stage in gene selection was the analysis of ccRCC expression in TCGA database, followed by the use of GEO and GTEx databases to remove genes present in blood (Fig. 1). Subsequently, in order to evaluate candidate genes with respect to their suitability to serve as blood biomarkers by distinguishing ccRCC vs. normal patients, their expression was first tested by RT-qPCR in ccRCC and normal tissues, and secondly in blood samples of cancer patients versus non-cancer patients and healthy donors. According to the bioinformatics analysis, the higher expression in ccRCC tissue would be expected to be confirmed as compared to normal tissue, and subsequently, when PAXgene blood samples are tested, the higher expression of at least some of the candidate genes would hopefully be determined in PAX blood from cancer patients compared to healthy. In order to obtain RNA seq based expression profiles in ccRCC and compare them to normal tissue as well as blood, the Cancer Genome Atlas database (TCGA, [31]) was used. The TCGA data portal is the largest and most commonly used public resource providing somatic and germline mutation, gene expression, gene methylation and copy number variation (CNV) data sets, amongst others, for several thousands of tumor samples. Data was obtained for 470 ccRCC patients, including 68 samples from matched normal tissue. In cases where multiple samples corresponded to a single patient, average expression values were calculated. Out of 20,533 TCGA genes in total, blood expression data from sources described below was found for 20,466 genes. Ideally candidate genes should not have wide expression domains; so as to provide a measure of kidney specificity for a gene, Tissue-specific Gene Expression and Regulation database (TiGER, [32]) was consulted, which is based on the analysis of the NCBI EST database [33] for 30 human tissues and contains tissue-specific expression profiles for 20,000 UniGenes. From 458 enriched in kidney, genes also expressed in blood, liver, prostate and bladder were deducted, leaving a list of 95 conditionally named ‘kidney specific’ genes.

Fig. 1
figure1

Workflow diagram

In order to obtain blood expression profiles, a comprehensive search for RNA seq expression data from healthy individuals was made in literature and online databases. Gene Expression Omnibus database (GEO, [34]) archives and freely distributes microarray, next-generation sequencing, and other forms of high-throughput functional genomics data. This database was searched by variations of ‘blood [Sample Source] AND Homo sapiens [Organism] AND high throughput sequencing [Platform Technology Type]’ providing a total of seven usable datasets altogether comprising 91 individual blood samples. Further 376 blood samples were obtained from GTEx database [35] and an additional source of one blood sample pooled from five individuals was kindly provided by Dr. Zhao and Dr. Zhang of Pfizer.

In order for expression profiles in important organs or organs related to the urological system to be given some relevance, RNA seq data from normal tissue was also considered in the analysis. From TCGA database, data was obtained for normal liver and bladder (9 and 11 samples respectively) and an analogous GEO search resulted in a recovery of a small number of samples for kidney, liver and bladder. Additional samples for kidney and liver (pooled from multiple donors) were included from RNA seq Atlas [36] (Table 1).

Table 1 Sources of expression profile datasets

Processing of data downloaded from TiGER database, RNA seq expression data, calculation of rpkm values (reads per kilobase million) where necessary, translation of gene names and statistics were done in Python. Rpkm values were calculated according to the formula: raw count × 1,000,000/(gene length x library size). Translation of gene names was done using BioMart [38]. Mann-Whitey U test was used to discern between cancer and matched normal samples from TCGA with statistical significance defined as p < 0.05. In cases of gene expression entries with multiple isoforms, replicate samples, duplicate gene names, absolute highest values were taken, so as not to underestimate the possible presence in blood.

RNA isolation and RT-qPCR analysis

Total tissue RNA (1 μg) was purified using miRNeasy Kit (Qiagen, Hilden, Germany) following homogenization using TissueLyser II (Qiagen). Total RNA from PAXgene blood tubes was purified using the PAXgene Blood miRNA Kit (Qiagen). Total RNA concentration was determined by NanoDrop 1000 Spectrometer (Thermo Fisher Scientific Inc., Wilmington, DE, USA) by measuring the absorbance at 260 nm and RNA purity by measuring A260/280 ratios. The integrity and size distribution of the tissue and blood derived RNA was analyzed using Bioanalyzer (Agilent RNA 6000 Nano Kit). Only samples with RNA integrity number values equal or above 7 were included. RNA samples from normal tissues were pooled together and the same was done with the ones from cancer, producing one normal pool (NN) and one cancer pool (NC). Complementary DNA synthesis was performed using the Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science, Mannheim, Germany) with a mix of random hexamer and anchored-oligo(dT) primers. RNA was also isolated and transcribed from the renal cell carcinoma cell line 786–0 to assess the quality of all newly made cDNA from tissue and PAXgene blood samples. Normalization of the RT-qPCR data was done using the kidney reference gene peptidylproline isomerase A (PPIA) [39].

Primers were designed for SYBR Green using NCBI’s PrimerBlast and Primer3 (see Additional file 1), so as to cover the maximum number of isoforms. The criteria for primer design were: amplicon length 60–150 nt, primer length 18–30 nt, intron spanning (intron length > 1000 nt), GC content 40–60%. For certain genes UPL probes were used in which case primers were automatically suggested with a given probe by the online tool (Universal Probe Library, Roche [40]), and common assays were selected for genes with multiple isoforms.

The relative quantification of transcripts was done on the Light Cycler 480 (Roche) using the QuantiTec SYBR Green PCR Kit (Qiagen) as previously described [37]. In case of UPL probes LightCycler 480 Probes Master Kit (Roche) was used. PCR was done on 96-well plates, with kidney cancer cell line 786–0 and ccRCC tissues as positive controls. PCR conditions were optimized where necessary and the size of PCR products was confirmed by electrophoresis using Bioanalyzer (Agilent DNA 1000 Kit). PCR data were analyzed by qBasePLUS software (Biogazelle NV, Gent, Belgium). With respect to the qBasePLUS processing, the samples were divided in 2 or 3 groups: normal vs. all cancer samples, i.e. cancer and metastatic cancer in a single group, as shown in the table ‘qBasePLUS results: normal vs. all cancer samples’ (see Additional file 2), normal vs. non-metastatic cancer, normal vs. metastatic cancer, and non-metastatic cancer vs. metastatic, as shown in the table ‘qBasePLUS results: non-metastatic cancer vs. metastatic samples’ (see Additional file 3). Results were calculated for 100% PCR efficiency and ´unpaired´ experimental design.

Statistics

Statistical analysis was done with GraphPad Prism 6.07 (GraphPad Software, San Diego, CA, USA) and qBasePLUS, using the Mann-Whitney U-test. P values < 0.05 were considered statistically significant. Graphs were generated in GraphPad Prism using the Mann-Whitney U-test.

Results

Candidate gene selection

In order to obtain a list of genes potentially useful as biomarkers, only genes with supposedly no blood expression, favorable statistical distance between distributions of cancer and normal values and high expression in cancer were taken into account. Regarding blood expression, values below 1 rpkm were considered low enough as to signify possible non expression, with respect to the sensitivity of detection. As a measure of distance of cancer and matched normal tissue distributions, the ratio of 5th percentile of cancer distribution with 95th percentile from normal was taken, and values above 0.5 considered favorable. Another measure of distance was calculated where the score represents the multiplication of probabilities of patients from each distribution falling within the overlap interval (score = Xprob x Yprob). Individual probabilities are calculated as the number of patients whose rpkm values fall within the overlap interval, divided by the total number of patients in the distribution (Xprob = patients within the overlap interval/total number of patients). Score is assigned 0 if the distributions do not overlap, and 1 for identical distributions. In cases when one distribution is within the other, but there are no patients from the larger one which fall into the overlap interval (they are distributed on both sides of it) score is assigned 1 as those genes are not valuable for further analysis. This method of calculating statistical distance is generally stricter then the percentile ratio, with favorable distance represented by values less than 0.3.

For genes of interest, the expression levels in liver, bladder, prostate and kidney in healthy individuals were also taken into account, giving preferential ranking to genes with lower rpkm values. Literature, the Human Protein Atlas [41], and OMIM [42] were consulted to gather information regarding gene function and expression domains of selected genes. Gene functions related to metabolic pathways in kidney or implicated in cancer (especially genes linked to ccRCC and hypoxia-inducible factors HIF1α and HIF2α) together with absence of expression in bone marrow and immune system, low or no expression in most tissues, and enrichment in kidney were considered favorable with respect to gene ranking.

A group of 20 genes were found to strictly fulfill expression criteria (defined as: blood expression GEO sources 95th percentile <1rpkm, GTEx 95th percentile <=1; fold change TCGA cancer median/matched normal tissue median > 1; distribution distance 5th percentile TCGA cancer/95th percentile matched normal tissue > 0.5, TCGA cancer median > 5rpkm) (Table 2, first 20 genes). The first 13 genes have median cancer values above 10. In addition, when it is considered that the rate of release of RNA from ccRCC into blood may be much higher than from normal kidney, as well as the presence of circulating tumor cells, the fold change median cancer/matched normal tissue, as well as the percentile ratio distribution distance measure become less relevant and may be relaxed in terms of gene selection. A similar argument follows for blood expression considering that individual blood sources may not be fully reliable and false outliers may be present. This enables the inclusion of certain genes which do not satisfy all the criteria fully, but may excel in some (last 11 genes in the table).

Table 2 Candidate genes

Many of these genes have previously been implicated in ccRCC, largely in micro-array studies [43,44,45,46,47,48,49,50]. This analysis identified this group of genes as having zero or low RNA blood presence, suggesting their potential use as ccRCC biomarkers in blood.

Analysis of the expression of candidate genes

To obtain an approximate overview of the levels in tissue, the expression was analyzed in cancer versus normal tissue for 15 genes of highest interest and the bioinformatics analysis was confirmed, as all these genes showed increased levels in cancer, notably so CA9 and NDUFA4L2 (Fig. 2). Certain genes were excluded from the analysis because of detection issues (multiple isoforms, etc.).

Fig. 2
figure2

Confirmation of TCGA data by RT-qPCR: Candidate genes were overexpressed in ccRcc compared to normal tissue (all values are above 1). Fold change is calculated as 2exp (Cqnormal-Cqcancer)

Blood testing consisted of two stages: in the first stage 3 PAXgene cancer samples were used (Table 3), and only genes with good (Cq < 33) detectability were selected for the second stage of blood testing with 24 more PAXgene samples (13 cancer and 11 normal/healthy), in order to broadly evaluate their expression. Especially good detectability was shown for following genes: CDK18 (Cq = 27), EGLN3 (Cq = 26), TMEM45A (Cq = 28), CAV2 (Cq = 26). A larger number of genes was undetectable or with extremely high Cq values. The genes with the supposedly highest potential based on bioinformatics and tissue PCR analysis (NDUFA4L2 and CA9) had low detectability with very high Cq values (around 34). Nevertheless, NDUFA4L2 was tested on all 27 samples and was confirmed to be undetectable. In summary, 9 genes were finally selected for the second stage of testing (CAV2a, FABP7, ESM1, NOL3, LOX, CDK18, EGLN3, TMEM45A, CCND1). In the second stage, the expression levels turned out to be similar in cancer versus normal blood for most genes, except for CDK18 and CCND1 which paradoxically turned out to be downregulated in cancer blood (Table 4). Additional testing with 10 plasma samples indicated non-measurable expression. There was no correlation between the expression levels in blood for CDK18, CCND1, and LOX, and patient data such as age, tumor grade and stage.

Table 3 Evaluation of candidate genes by RT-qPCR in tissue and 3 blood samples
Table 4 Evaluation of expression in the second stage of blood testing with 27 samples

The downregulation of CDK18 RNA in cancer blood (metastatic and non-metastatic grouped together) compared to normal was significant with p value = 0.001, whereas CCND1 was downregulated with p = 0.039 (Fig. 3). For both genes, there was no significant difference in levels when non-metastatic and metastatic samples were compared with each other. The results also showed a tendency towards upregulation for LOX when non-metastatic were compared to metastatic cancer samples, with the p value very close to significant (p = 0.058) (Fig. 3).

Fig. 3
figure3

Blood relative mRNA expression of CDK18, CCND1, and LOX based on the qBase exported relative quantity (RQ) values, calculated from Cq values, according to the formula: RQ = 2 (meanCq-Cq); results from qBase (RQ values) were processed in GraphPad Prism in order to generate graphs using the Mann-Whitney U-test. N- normal patient samples; T- tumor patient samples; mT- metastatic. a CDK18 was underexpressed in PAX blood tumor samples compared to normal PAX blood. b There is no significant difference in expression of CDK18 between tumor and metastatic tumor PAX blood samples. c CCND1 was underexpressed in PAX blood tumor samples compared to normal PAX blood. d There is no significant difference in expression of CCND1 between tumor and metastatic tumor PAX blood samples. e There is no significant difference in expression of LOX in PAX blood tumor samples compared to normal PAX blood. f Lox shows a tendency towards upregulation in metastatic compared to non-metastatic tumor PAX blood samples

Discussion

In this study, a gene panel was made comprising the most highly overexpressed genes in ccRCC tissue whose mRNA also had the potential of being absent in the blood of healthy individuals. The first stage in the construction of this panel was the TCGA database- to select a panel of the most overexpressed genes in ccRCC, followed by the GEO and GTEx databases- to deduct from this panel genes showing measurable expression in the blood of healthy individuals. After confirming the tissue overexpression on samples from ccRCC patients in the next step, RT-qPCR analysis was done to evaluate mRNA levels in whole blood of ccRCC patients vs. patients without ccRCC and healthy donors. Measurable genes did not show overexpression in normal blood while two genes exhibited downregulation.

Whole blood analysis of selected genes does not show increased mRNA levels

The genes that had highest potential based on the bioinformatics analysis were NDUFA4L2 and CA9. According to the TCGA, the former has very high median expression in ccRCC tissue (701 rpkm) whereas the latter has the highest overexpression in ccRCC compared to normal tissue (1218). However, we found both to be undetectable in whole blood by qPCR. Several other candidate genes were found to be undetectable in whole blood, while most of the genes that were detectable (EGLN3, CAV2, ESM1, TMEM45A, NOL3, FABP7) did not exhibit significant dysregulation in expression between cancer and healthy PAXgene samples. A plausible path to overcoming this outcome was to examine these genes in plasma, as the mRNA levels (supposedly originating from the expression in blood cells) could significantly drop in healthy samples compared to cancer ones once the blood cells are removed, revealing the effect of the tumor derived RNA. However, after testing 10 plasma samples, our conclusion was that gene expression in the plasma was not measurable. The PAXgene system is used for stabilization and isolation of mRNA and other classes of nucleic acids (such as genomic DNA and miRNA). Blood specimens are collected in tubes containing a stabilization reagent preventing nuclease degradation and transcriptional changes in anticoagulated whole blood, and stabilizing RNA for up to 3 days at room temperature, for the purpose of expression profiling [51]. All handling of RNA was done with special care, and although it is reasonable to assume that for many or all of the candidate genes RNA was degraded by blood RNAses, the RNA integrity of whole PAXgene samples was indeed satisfactory, reflected by their high RIN values. Apart from the issues regarding RNA stability and concerning expression interference from blood cells, possible limitations of this work design may stem from the bioinformatics phase. Acquired GEO datasets, which were used to screen for genes absent from blood (with supposedly no blood expression) may not be 100% reliable; they came from many different sources and were not in a perfect mutual accordance. A separate issue is the cutoff value of < 1 rpkm to signify absence of blood presence for a gene. Most authors somewhat arbitrarily define the expression threshold as 1 rpkm (more generally anywhere between 0.3 rpkm and 1 rpkm), below which the sensitivity of RNA sequencing is insufficient to confirm expression and distinguish it from the background [52, 53]. Our cutoff potentially may have enabled genes with minute expression in blood cells to be included in the wet lab analysis.

CDK18, CCND1, lox

The study also revealed the downregulation of two genes, CDK18 and CCND1, in ccRCC blood compared to healthy samples, as well as a tendency towards upregulation for LOX in metastatic compared to non-metastatic ccRCC. These results may be suitable for additional analysis in a larger patient cohort.

Cyclin-dependent kinase 18 (CDK18, PCTK3, PCTAIRE, PCTAIRE3) belongs to the PCTAIRE protein kinases, which are a subfamily of cdc2-related serine/threonine protein kinases named for a cysteine-for-serine substitution in the PCTAIRE motif conserved in the initially characterized CDK proteins (PCTAIRE sequence instead of the PSTAIRE sequence). They have unique N and C-terminal domains that extend forth from a serine/threonine kinase domain that is highly homologous to cdc2 [54]. PCTAIRE kinase subfamily includes three members, PCTK1/CDK16, PCTK2/CDK17, and PCTK3/CDK18 which are poorly researched. Insights into the activation of CDK18 have recently been obtained- it binds cyclin A2 and cyclin E1 (pulldown experiment with HEK293T cells) and is activated by cyclin A2 and PKA (cAMP-dependent protein kinase) [55]. CDK18 has recently been shown to regulate cell migration and adhesion in HEK293T cells by negatively modulating FAK (focal adhesion kinase) activity and reorganizing actin and associated skeletal/adhesion proteins such as cofilin, and has also been implicated in vesicular transport via interaction with Sec23Ap [56]. Overexpression of CDK18 also led to the formation of filopodia during the early stages of cell adhesion in HeLa cells [57]. Interestingly, it has also been recently found to play a role in replication stress and positively regulate genome stability, by associating with RAD proteins [58]. Lastly, PCTAIRE-3 as well as PCTAIRE-2 have been implicated in Alzheimer’s disease [59, 60]. CDK18 was induced by CTS-1 (Chimeric tumor suppressor-1, p53-derived synthetic tumor suppressor) and mediated growth arrest and death in glioma cells [61]. Apart from its activation by cyclin A2, in the same study CDK18 was shown to phosphorylate retinoblastoma tumorsupressor protein (Rb) in vitro [55]. Although PCTAIRE1 has been found to be upregulated in many cancers, so far there is no such data for CDK18.

Cyclin D1 (CCND1) regulates CDK4 or CDK6, whose activity is necessary for G1/S transition of the cell cycle. CCND1 is more frequently dysregulated in human cancers and therefore more studied than cyclin D2 or D3. Its overexpression leads to aberrant CDK activation resulting in rapid growth and division and is correlated with tumor stage, increased metastasis and poor prognosis in various cancers [62]. It is also involved in processes such as DNA repair and the control of mitochondrial activity and cell migration; it may assume CDK-independent functions as well [63]. CCND1 was investigated by microarray and TMA in ccRCC, and found to be upregulated and a potential therapeutic target [64]. In another study, CCND1 was found to be a useful immunohistochemical marker to discriminate between chromophobe renal cell carcinoma and renal oncocytoma [65].

Lysyl oxidase (LOX) performs covalent cross linking in elastin and collagen by oxidizing lysine residues, and therefore is important for the integrity of extracellular matrix [66]. It has both intracellular and extracellular functions and is involved in a number of pathological processes that affect the connective tissue [67]. It is upregulated in many cancers and involved in tumor progression, although it has been reported to function as a tumor-suppressor as well. Its concrete roles in various aspects of tumorigenesis have been reviewed recently [68]. LOX is a HIF target [69] and in ccRCC, LOX has been shown to be strongly overexpressed compared to normal tissue; it is one of the genes postranscriptionally regulated by miR-141-3p and miR-145-5p; and has prognostic relevance for the overall survival of ccRCC patients [70]. In ccRCC cell cultures, it has been found to function in a positive-regulative loop with HIF-1α, and to influence ccRCC progression by modifying cellular adhesion, migration, and the rigidity of the collagen matrix [71].

Conclusions

In summary, with the goal of finding potential blood-based biomarkers for ccRCC, this study investigated the blood presence of genes highly overexpressed in ccRCC tissue and compared their blood mRNA levels between healthy and ccRCC patient samples. The overexpression in tissue was not reflected in the increase in the levels of mRNAs circulating in the blood of ccRCC patients. The analysis revealed the transcripts of CDK18 and CCND1 as underexpressed in the blood of ccRCC patients, and LOX as showing a tendency towards upregulation in metastatic ccRCC compared to non-metastatic. Further analysis of the selected gene panel by using a larger patient cohort may prove useful.

Availability of data and materials

All the data are available from the corresponding author upon request.

Abbreviations

CCND1:

Cyclin d1

ccRCC:

Clear cell renal cell carcinoma

CDK18:

Cyclin-dependent kinase 18

CEA:

Carcinoembryonic antigen

cfDNA:

Circulating cell-free dna

circRNA:

Circular rna

CNV:

Copy number variation

CTCs:

Circulating tumor cells

GEO:

Gene expression omnibus database

HCC:

Hepatocellular carcinoma

HIF:

Hypoxia inducible factor

lncRNA:

Long non-coding rnas

LOX:

Lysyl oxidase

miRNAs:

Micrornas

PSA:

Prostate-specific antigen

RCC:

Renal cell carcinoma

RT-qPCR:

Quantitative real time pcr

TCGA:

Cancer Genome Atlas database

TiGER:

Tissue-specific Gene Expression and Regulation database

References

  1. 1.

    Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68:7–30.

  2. 2.

    Srigley JR, Delahunt B, Eble JN, Egevad L, Epstein JI, Grignon D, et al. The International Society of Urological Pathology (ISUP) Vancouver classification of renal neoplasia. Am J Surg Pathol. 2013;37:1469–89.

  3. 3.

    Rydzanicz M, Wrzesiński T, Bluyssen HA, Wesoły J. Genomics and epigenomics of clear cell renal cell carcinoma: recent developments and potential applications. Cancer Lett. 2013;341:111–26.

  4. 4.

    Martel CL, Lara PN. Renal cell carcinoma: current status and future directions. Crit Rev Oncol Hematol. 2003;45(2):177–90.

  5. 5.

    McLaughlin J, Lipworth L, Tarone R, Blot W. Cancer Epidemiology and Prevention. Oxford Univ. Press; 2006. p. 1087–100.

  6. 6.

    Chow WH, Dong LM, Devesa SS. Epidemiology and risk factors for kidney cancer. Nat Rev Urol. 2010;7:245–57.

  7. 7.

    Gospodarowicz MK, Miller D, Groome PA, Greene FL, Logan PA, Sobin LH. The process for continuous improvement of the TNM classification. Cancer. 2004;100(1):1–5.

  8. 8.

    Powles T, Staehler M, Ljungberg B, Bensalah K, Canfield SE, Dabestani S, et al. European Association of Urology guidelines for clear cell renal cancers that are resistant to vascular endothelial growth factor receptor-targeted therapy. Eur Urol. 2016;70(5):705–6.

  9. 9.

    Jonasch E, Gao J, Rathmell WK. Renal cell carcinoma. BMJ. 2014;349:4797.

  10. 10.

    Ljungberg B, Bensalah K, Bex A, Canfield S, Dabestani S, Hofmann F, et al. Guidelines on renal cell carcinoma: 2014 update. Eur Assoc Urol. 2015;67:913–24.

  11. 11.

    Singer EA, Gupta GN, Srinivasan R. Update on targeted therapies for clear cell renal cell carcinoma. Curr Opin Oncol. 2011;23(3):283–9.

  12. 12.

    MacLennan S, et al. Systematic review of oncological outcomes following surgical management of localised renal cancer. Eur Urol. 2012;61:972–93.

  13. 13.

    Janzen NK, Kim HL, Figlin RA, Belldegrun AS. Surveillance after radical or partial nephrectomy for localized renal cell carcinoma and management of recurrent disease. Urol Clin North Am. 2003;30:843–52.

  14. 14.

    Sun M, Shariat SF, Cheng C, et al. Prognostic factors and predictive models in renal cell carcinoma: a contemporary review. Eur Urol. 2011;60:644–61.

  15. 15.

    Atzpodien J, Royston P, Wandert T, et al. Metastatic renal carcinoma comprehensive prognostic system. Br J Cancer. 2003;88:348–53.

  16. 16.

    Kim HL, Seligson D, Liu X, et al. Using tumor markers to predict the survival of patients with metastatic renal cell carcinoma. J Urol. 2005;173:1496–501.

  17. 17.

    Adashek JJ, Salgia MM, Posadas EM, Figlin RA, Gong J. Role of Biomarkers in Prediction of Response to Therapeutics in Metastatic Renal-Cell Carcinoma. Clin Genitourin Cancer. 2019;17(3):e454–e460. Review.

  18. 18.

    Majer W, Kluzek K, Bluyssen H, Wesoły J. Potential approaches and recent advances in biomarker discovery in clear-cell renal cell carcinoma. J Cancer. 2015;6(11):1105–13.

  19. 19.

    Reddi KK, Holland JF. Elevated serum ribonuclease in patients with pancreatic cancer. Proc Natl Acad Sci U S A. 1976;73:2308–10.

  20. 20.

    Huggett JF, Novak T, Garson JA, Green C, Morris-Jones SD, Miller RF, et al. Differential susceptibility of PCR reactions to inhibitors: an important and unrecognised phenomenon. BMC Res Notes. 2008;1:70.

  21. 21.

    Connolly ID, et al. The “liquid biopsy”: the role of circulating DNA and RNA in central nervous system tumors. Curr Neurol Neurosci Rep. 2016;16(3):25.

  22. 22.

    Yörüker EE, Holdenrieder S, Gezer U. Blood-based biomarkers for diagnosis, prognosis and treatment of colorectal cancer. Clin Chim Acta. 2016;455:26–32.

  23. 23.

    Zane K, et al. Emerging blood-based biomarkers for detection of gastric cancer. World J Gastroenterol. 2015;21(41):11636–53.

  24. 24.

    Nandagopal L, Sonpavde G. Circulating biomarkers in bladder cancer. Bladder Cancer. 2016;2(4):369–79.

  25. 25.

    Stroun M, Lyautey J, Lederrey C, et al. About the possible origin and mechanism of circulating DNA apoptosis and active DNA release. Clin Chim Acta. 2001;313(1–2):139–42.

  26. 26.

    Pepe MS, Etzioni R, Feng Z, Potter JD, Thompson ML, Thornquist M, et al. Phases of biomarker development for early detection of cancer. J Natl Cancer Inst. 2001;93(14):1054–61.

  27. 27.

    Danila DC, Anand A, Schultz N, et al. Analytic and clinical validation of a prostate cancer-enhanced messenger RNA detection assay in whole blood as a prognostic biomarker for survival. Eur Urol. 2014;65(6):1191–7.

  28. 28.

    Rodia MT, et al. Systematic large-scale meta-analysis identifies a panel of two mRNAs as blood biomarkers for colorectal cancer detection. Oncotarget. 2016;7(21):30295–306.

  29. 29.

    Fuhrman SA, Lasky LC, Limas C. Prognostic significance of morphologic parameters in renal cell carcinoma. Am J Surg Pathol. 1982;6:655–63.

  30. 30.

    Sobin LH, Wittekind C. TNM classification of malignant tumours. New York: Wiley-Liss; 2002. p. 193–5.

  31. 31.

    Cancer Genome Atlas Research Network. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013;499:43–9.

  32. 32.

    Liu X, Yu X, Zack DJ, Zhu H, Qian J. TIGER: a database for tissue-specific gene expression and regulation. BMC Bioinformatics. 2008;9:271.

  33. 33.

    Boguski MS, Lowe TM, Tolstoshev CM. dbEST--database for “expressed sequence tags”. Nat Genet. 1993;4:332–3.

  34. 34.

    Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.

  35. 35.

    Lonsdale J, et al. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45:580–5.

  36. 36.

    Krupp M, Marquardt JU, Sahin U, Galle PR, Castle J, Teufel A. RNA-Seq atlas- a reference database for gene expression profiling in normal tissue by next-generation sequencing. Bioinformatics. 2012;28(8):1184–5.

  37. 37.

    Zhao S, et al. Comparison of stranded and non-stranded RNA-seq transcriptome profiling and investigation of gene overlap. BMC Genomics. 2015;16:675.

  38. 38.

    BioMart. http://www.ensembl.org. Accessed 20 Apr 2016.

  39. 39.

    Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, Kristiansen G, et al. In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC Mol Biol. 2007;8:47.

  40. 40.

    Universal Probe Library. www.universalprobelibrary.com. Accessed 5 May 2016.

  41. 41.

    Uhlen M, et al. A human protein atlas for normal and cancer tissues based on antibody proteomics. Mol Cell Proteomics. 2005;4(12):1920–32.

  42. 42.

    Hamosh A, et al. Online Mendelian inheritance in man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res. 2002;30(1):52–5.

  43. 43.

    Meyer HA, Tolle A, Jung M, Fritzsche FR, Haendler B, Kristiansen I, et al. Identification of stanniocalcin 2 as prognostic marker in renal cell carcinoma. Eur Urol. 2009;55:669–78.

  44. 44.

    Thibodeau BJ, Fulton M, Fortier LE, Geddes TJ, Pruetz BL, Ahmed S, et al. Characterization of clear cell renal cell carcinoma by gene expression profiling. Urol Oncol. 2016;34(4):168.

  45. 45.

    Feng JY, Diao XW, Fan MQ, Wang PX, Xiao Y, Zhong X, et al. Screening of feature genes of the renal cell carcinoma with DNA microarray. Eur Rev Med Pharmacol Sci. 2013;17(22):2994–3001.

  46. 46.

    Fu L, Minton DR, Zhang T, Nanus DM, Gudas LJ. Genome-wide profiling of track kidneys shows similarity to the human ccRCC transcriptome. Mol Cancer Res. 2015;13(5):870–8.

  47. 47.

    Schrödter S, Braun M, Syring I, Klümper N, Deng M, Schmidt D, et al. Identification of the dopamine transporter SLC6A3 as a biomarker for patients with renal cell carcinoma. Mol Cancer. 2016;15:10.

  48. 48.

    Minton DR, Fu L, Mongan NP, Shevchuk MM, Nanus DM, Gudas LJ. Role of NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4-like 2 in clear cell renal cell carcinoma. Clin Cancer Res. 2016;22(11):2791–801.

  49. 49.

    Klatte T, Kabbinavar F, Zomordian N, Belldegrun AS, Pantuck AJ. Prospective evaluation of carbonic anhydrase IX (CAIX) as a molecular marker in metastatic renal cell carcinoma: interim results. J Urol. 2007;177:164.

  50. 50.

    Bruick RK, McKnight SL. A conserved family of prolyl-4-hydroxylases that modify HIF. Science. 2001;294:1337–40.

  51. 51.

    Rainen L, Oelmueller U, Jurgensen S, et al. Stabilization of mRNA expression in whole blood samples. Clin Chem. 2002;48:1883–90.

  52. 52.

    Ramsköld D, Wang ET, Burge CB, Sandberg R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009;5(12):e1000598.

  53. 53.

    Magee R, Loher P, Londin E, Rigoutsos I. Threshold-seq: a tool for determining the threshold in short RNA-seq datasets. Bioinformatics. 2017;33(13):2034–6.

  54. 54.

    Meyerson M, Enders GH, Wu CL, et al. A family of human cdc2-related protein kinases. EMBO J. 1992;11:2909–17.

  55. 55.

    Matsuda S, Kominato K, Koide-Yoshida S, et al. PCTAIRE kinase 3/cyclin-dependent kinase 18 is activated through association with cyclin a and/or phosphorylation by protein kinase a. J Biol Chem. 2014;289(26):18387–400.

  56. 56.

    Palmer KJ, Konkel JE, Stephens DJ. PCTAIRE protein kinases interact directly with the COPII complex and modulate secretory cargo transport. J Cell Sci. 2005;118:3839–47.

  57. 57.

    Matsuda S, Kawamoto K, Miyamoto K, Tsuji A, Yuasa K. PCTK3/CDK18 regulates cell migration and adhesion by negatively modulating FAK activity. Sci Rep. 2017;7:45545.

  58. 58.

    Barone G, Staples CJ, Ganesh A, Patterson KW, Bryne DP, Myers KN, et al. Human CDK18 promotes replication stress signaling and genome stability. Nucleic Acids Res. 2016;44(18):8772–85.

  59. 59.

    Herskovits AZ, Davies P. The regulation of tau phosphorylation by PCTAIRE 3: implications for the pathogenesis of Alzheimer's disease. Neurobiol Dis. 2006;23(2):398–408.

  60. 60.

    Chaput D, Kirouac L, Stevens SM Jr, Padmanabhan J. Potential role of PCTAIRE-2, PCTAIRE-3 and P-histone H4 in amyloid precursor protein-dependent Alzheimer pathology. Oncotarget. 2016;7(8):8481–97.

  61. 61.

    Naumann U, Huang H, Wolburg H, Wischhusen J, Weit S, Ohgaki H, et al. PCTAIRE3: a putative mediator of growth arrest and death induced by CTS-1, a dominant-positive p53-derived synthetic tumor suppressor, in human malignant glioma cells. Cancer Gene Ther. 2006;5:469–78.

  62. 62.

    Qie S, Diehl JA. Cyclin D1, cancer progression, and opportunities in cancer treatment. J Mol Med (Berl). 2016;94(12):1313–26.

  63. 63.

    Zwijsen RM, Wientjens E, Klompmaker R, van der Sman J, Bernards R, Michalides RJ. CDK-independent activation of estrogen receptor by cyclin D1. Cell. 1997;88(3):405–15.

  64. 64.

    Karim S, Al-Maghrabi JA, Farsi HM, Al-Sayyad AJ, Schulten HJ, Buhmeida A, et al. Cyclin D1 as a therapeutic target of renal cell carcinoma- a combined transcriptomics, tissue microarray and molecular docking study from the Kingdom of Saudi Arabia. BMC Cancer. 2016;16(Suppl 2):741.

  65. 65.

    Zhao W, Tian B, Wu C, Peng Y, Wang H, Gu WL, et al. DOG1, cyclin D1, CK7, CD117 and vimentin are useful immunohistochemical markers in distinguishing chromophobe renal cell carcinoma from clear cell renal cell carcinoma and renal oncocytoma. Pathol Res Pract. 2015;211(4):303–7.

  66. 66.

    Kumari S, Panda TK, Pradhan T. Lysyl oxidase: its diversity in health and diseases. Indian J Clin Biochem. 2017;32(2):134–41.

  67. 67.

    Zheng W, et al. Low extracellular lysyl oxidase expression is associated with poor prognosis in patients with prostate cancer. Oncol Lett. 2016;12(5):3161–6.

  68. 68.

    Wang TH, Hsia SM, Shieh TM. Lysyl oxidase and the tumor microenvironment. Int J Mol Sci. 2017;8(1):62.

  69. 69.

    Tanaka T, et al. Anthracycline inhibits recruitment of hypoxia-inducible transcription factors and suppresses tumor cell migration and cardiac Angiogenic response in the host. J Biol Chem. 2012;287(42):34866–82.

  70. 70.

    Liep J, Kilic E, Meyer HA, Busch J, Jung K, Rabien A. Cooperative effect of miR-141-3p and miR-145-5p in the regulation of targets in clear cell renal cell carcinoma. PLoS One. 2016;11(6):e0157801.

  71. 71.

    Di Stefano V, Torsello B, Bianchi C, Cifola I, Mangano E, Bovo G, et al. Major action of endogenous Lysyl oxidase in clear cell renal cell carcinoma progression and collagen stiffness revealed by primary cell cultures. Am J Pathol. 2016;186(9):2473–85.

Download references

Acknowledgements

We are grateful to Silke Rabenhorst for excellent technical assistance.

Funding

SS, KJ, KS and AR were supported by funds from the Urological Research Foundation (Stiftung Urologische Forschung Berlin), which was not involved in any of the scientific or practical parts of the study.

Author information

Study design: KMS, KJ, AR, SS; bioinformatics: CH, SS; acquisition of data: JB, MJ, SS; interpretation of data: KMS, KJ, AR, SS; manuscript preparation: SS; substantial revision: JB, KJ, KMS, AR. We confirm that all authors read and approved the final manuscript.

Correspondence to Sinisa Simonovic.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the ethics committee of Charite - Universitätsmedizin Berlin by number EA1/134/12. The study was conducted in compliance with the declaration of Helsinki and written informed consent was obtained.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Simonovic, S., Hinze, C., Schmidt-Ott, K.M. et al. Limited utility of qPCR-based detection of tumor-specific circulating mRNAs in whole blood from clear cell renal cell carcinoma patients. BMC Urol 20, 7 (2020). https://doi.org/10.1186/s12894-019-0542-9

Download citation

Keywords

  • Renal cell carcinoma
  • CDK18
  • CCND1
  • Biomarker
  • Blood
  • PAXgene