Abstract
Although testicular germ cell tumor (TGCT) overall is highly curable, patients may experience late effects after treatment. An increased understanding of the mechanisms behind the development of TGCT may pave the way for better outcome for patients. To elucidate molecular changes prior to TGCT diagnosis we sequenced small RNAs in serum from 69 patients who were later diagnosed with TGCT and 111 matched controls. The deep RNA profiles, with on average 18 million sequences per sample, comprised of nine classes of RNA, including microRNA. We found that circulating RNA signals differed significantly between cases and controls regardless of time to diagnosis. Different levels of TSIX related to X-chromosome inactivation and TEX101 involved in spermatozoa production are among the interesting findings. The RNA signals differed between seminoma and non-seminoma TGCT subtypes, with seminoma cases showing lower levels of RNAs and non-seminoma cases showing higher levels of RNAs, compared with controls. The differentially expressed RNAs were typically associated with cancer related pathways. Our results indicate that circulating RNA profiles change during TGCT development according to histology and may be useful for early detection of this tumor type.
Introduction
Testicular germ cell tumor (TGCT) is diagnosed in 1% of men worldwide and is the most common malignancy in males between 20 and 39 years of age. The incidence rates are still rising (), and some of the highest incidences are found in Northern Europe (). The etiology is largely unknown, although genetic components and conditions during pregnancy seem to play a role (). The susceptibility to TGCT is shown to have a strong familial link, with a fourfold increased risk for fathers and eightfold for brothers (, ). The polygenic nature of TGCT has been recognized and more than 50 susceptibility genes have been identified (–). The susceptibility loci contain genes linked to germ cell development and sex determination, as well as genes related to tumor growth/suppression. Cryptorchidism is a well-known risk factor, and around 10% of males with TGCT have a history of cryptorchidism (). The main histologic subtypes of TGCT are seminoma, non-seminoma and mixed, making up around 50, 30, and 20%, respectively ().
TGCT is a highly curable disease since the introduction of cisplatin-based chemotherapy with a 5-year survival rate of above 95% (, ). However, there is an increased risk of long-term side effects including secondary non-germ cell (GC) cancers, cardiovascular disease, and hypogonadism (–).
The current guidelines presented by the European Association of Urology () for the diagnosis of TGCT include several techniques, including ultrasound imaging for initial diagnosis and serum markers for specific subtype diagnosis and prognosis (–). The serum markers used are alpha-fetoprotein, hCG and LDH, however alpha-fetoprotein is only seen in yolk-sac tumor and teratomas whilst hCG is expressed by trophoblasts only (). Whilst LDH is a less specific serum marker, the serum level is usually proportional to tumor volume and levels may be increased in advanced TGCT (). The presence of histological markers is also used for diagnosis. For germ cell neoplasia in situ (GCNIS), the markers include OCT3/4, PLAP, CD117, and SALL 4 (). Furthermore, it has been shown that the metabolic biomarkers leptin and resisting may predict cancer mortality due to their function of inducing pro-tumorigenic environment that further promotes tumor initiation, angiogenesis, and metastasis ().
RNAs, such as microRNA (miRNA), piwi-interacting RNA (piRNA), and long non-coding RNA (lncRNA), regulate gene expression on transcriptional and post-transcriptional levels and have been shown to be present in serum (). The differential composition of RNAs in circulation can help determine diagnosis as well as the developmental stage of the tumor (). Circulating miRNAs have been identified as both prognostic and diagnostic markers in biliary tract cancer (BTC), which has led to earlier diagnosis and less invasive procedure compared with previously used techniques (). Exosomal miRNAs have also been observed in colorectal cancer (CRC). Specifically, miR-150-5p and miR-99b-5p were found to be downregulated in CRC patients compared to healthy patients (). Differential levels of circulating RNA were also observed in pre-diagnostic serum samples in lung cancer patients, with an overall dynamic trend when advancing clinically ().
In TGCT patients, high expression of miRNAs belonging to the clusters miR-302/367 and miR-371-373 has been found in serum (–). Small RNA sequencing performed on TGCT tissue samples revealed miRNAs profiles to differ between normal and TGCT tissues, as well as between histological subtypes. A genome wide downregulation or loss of piRNAs was observed in TGCT, through mechanisms such as hypermethylation in CpG islands on genes associated with piRNAs (, –).
Circulating RNAs in TGCT patients have allowed for a less invasive and a more specific histology at diagnosis (). However, pre-diagnostic circulating RNA changes in TGCT have yet to be studied.
The aim of this study was to investigate the role of RNA in the development of TGCT in serum samples collected in a 10-years period prior to diagnosis. Furthermore, we investigated how differences in circulating RNAs were related to histologies and time periods before diagnosis. Differences in RNAs profiles were then used in functional enrichment analysis, and RNAs with stable patterns were investigated further.
Materials and Methods
Study Design and Participants
All samples included were retrieved from the Janus Serum Bank (JSB). JSB is a population-based cancer research biobank containing serum samples from 318 628 individuals collected from 1972 to 2004 (, ). The TGCT cases were identified by linking the JSB to the Cancer Registry of Norway using the Norwegian individually unique national identity numbers. The prediagnostic serum samples were donated up to 10 years before the diagnosis. We drew 111 cancer-free Janus participants for comparison of RNA levels with the cancer cases. The control subjects had to be alive and free from cancer, except for non-melanoma skin cancer, at the time of their matched case's diagnosis, and up to 10 years after blood collection. Controls were matched on age, time of blood collection and blood donor group (dependent on their county of residence, year of collection and method of storage), and as a result of these matching criteria, the average age at sample donation was 35 years (standard deviation of 6.5/6.7, respectively) for both cases and controls.
Previous studies have shown the feasibility of using long-term archived serum samples for RNA analyses and the variability of the RNA levels (, , ). The methods and analyses of the TGCT serum samples were similar to studies by Umu et al. (). The donors have given broad consent for the use of the samples in cancer research. The study was approved by the Norwegian regional committee for medical and health research ethics (REC no: 24 846, 2012/1590).
Laboratory Processing of Serum to RNA Profiles
We extracted RNAs from 400 μl serum using phenol-chloroform phase separation and the miRNeasy Serum/Plasma kit (Cat. no 1071073, Qiagen) using a QIAcube (Qiagen). NEBNext® Small RNA Library Prep Set for Illumina (Cat. No E7300, New England Biolabs Inc.) was used for small RNA-seq library preparation. RNA molecules from 17 to 47 nt in length were selected. We sequenced 12 samples per lane of a HiSeq 2500 (Illumina). Additional information is available in our previous study ().
Bioinformatics Analyses
Total number of reads generated was 3.5 billion with an average sampling depth of 18.4 million raw reads. AdapterRemoval v2.1.7 () was used to trim for adapters. We then mapped the collapsed reads to human genome version hg38 with Bowtie2 v2.2.9, with 10 alignments per read being allowed. Our annotation set consist of miRBase(v.22) () for miRNAs, pirBAse for piRNAs () and GENCODE () for other RNAs. IsomiR and tRF profiles were obtained through SeqBuster () and MINTmap, respectively (42). In the analyses, we included RNAs with at least five reads in more than 20% of the samples. Details are available in our previous study ().
The optmatch R package (github.com/markmfredrickson/optmatch) allowed us to find optimally matched sets of controls (Supplementary Table 1) for each analysis. The analyses were matched on age, histology when appropriate and technical artifacts. The technical artifacts accounted for differences in pre-analytical treatment and storage time (, ), termed blood donor groups. The DESeq2 R package (v1.18.1) (43) was used for the differential expression analyses using the default generalized linear model with a negative binomial distribution. We consider differences in expression outside the range{−1,1} log2 fold change (log2fc) as biologically meaningful and are counted in the summary bar charts. Venn diagrams, produced with the R package (VennDiagram, v. 1.6.20), represent common and unique differentially expressed RNA trends among the categories.
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed using kegga function from the limma R package. The inputs included mRNA, miRNA and isomiR targets extracted from miRDB (v5.0) predictions (44) with a cut off score of >60. With the mRNA and miRNA targets we also performed a pathway visualization (pathview, v.1.22.3) of selected pathway seen in the KEGG enrichment analysis (45, 46).
To exploit the full statistical power of the dataset when identifying RNAs that differ in cases and controls regardless of other factors, we first compared all cases and controls. To identify histology (determined by ICD-O (3rd revision) codes) specific RNA signals, we compared the RNA levels of seminoma with those of the matched controls, and similarly for non-seminomas and matched controls. To identify RNA signals according to the time between blood draw and diagnosis we divided prediagnostic time into four discrete time intervals, 0–2, 2–5, 5–8, and 8–10 years. We selected these intervals to optimize resolution on time to diagnosis while still having sufficient statistical power. To make the time windows comparable with respect to statistical power, the cut points were chosen to secure the same number of cases and controls, and similar proportions of histologies.
Results
Reduced Levels of RNAs in TGCT Cases
To identify differentially expressed RNAs in all TGCT cases (n = 79) against matched controls (n = 111), we analyzed 4231 RNAs that passed our inclusion criteria. Of these, we identified 818 RNAs that were differentially expressed with a p-adjusted value ≤0.05, and 88 of these had a log2 fold change (log2fc) outside the range {−1,1} (Figure 1A). The majority of these RNAs (82) had reduced levels in cases compared to controls, with an average log2fc = −1.37 and p-adjusted = 0.00032. The RNAs showing reduced levels were primarily isomiRs and mRNA fragments, whereas the few elevated RNAs consisted of lncRNAs, mRNA fragments, and a tRF (Figure 1B).
Figure 1
RNA Signals Differ With Histological Subtype
We compared RNA levels in seminoma and non-seminoma cases with control samples. As seen for all cases combined, seminomas exhibited a similar pattern with reduced level for the majority of the differentially expressed RNAs.
In total, 112 RNA signals were different in seminomas compared to controls, 102 with lower RNA levels and 10 with elevated levels (Figure 2A). Similar to all TGCT cases, the majority of the downregulated RNAs, with a mean log2fc values of −1.35, consisted of mRNA fragments and isomiRs, whereas the upregulated RNAs consisted of lncRNAs, mRNA fragments, piRNAs, and tRFs (Figure 2B).
Figure 2
In contrast to TGCT and seminomas, where the differences compared to controls were generally large, non-seminomas had 63 RNAs with reduced levels and 46 RNAs with elevated levels (Figure 2C), with mean log2fc values of −0.131 and 1.34, respectively. The majority of the RNAs with reduced levels consisted of mRNA fragments and isomiRs, as was also observed with all TGCT, a result driven by the seminomas. However, in non-seminomas the elevated RNAs showed higher proportions of isomiRs, and mRNA fragments, piRNAs and miRNAs as the second most abundant types (Figure 2D).
RNA Signals Are Not Associated With Time to Diagnosis
To investigate if there is any association between RNA signals and time to diagnosis, four different time frames during the prediagnostic period were examined. The RNA signals in all time frames exhibited similar patterns to the overall prediagnostic TGCT patterns, with the majority of significantly differentially expressed RNAs having reduced levels. This was strongest in the 2–5 years time frame with 282 RNA levels being reduced and 64 RNAs with elevated levels. The 0–2 years time frame showed the least amount of differentially expressed RNAs, with 93 RNAs showing reduced levels and 18 with elevated levels (Figure 3A). Fragments containing mRNAs constitute the major proportion of the significantly changed RNAs in all time frames. However, in the 2–5 years time frame, a higher number of piRNAs (76 out of 244) had different levels, and furthermore, in the 8–10 years, the higher proportion was isomiRs (64 out of 155) (Figure 3B).
Figure 3
Specific RNA Signals Are Independent of Histology and Time to Diagnosis
A number of RNAs are specific or common to histologies and time frames as illustrated by the venn diagrams in Figure 4. A total of 83 out of 88 of the RNAs identified were observed in both seminoma separately and all TGCT profiles. Furthermore, seminomas and non-seminomas shared 35 RNAs, and the majority of these were either mRNA fragments or isomiRs (Figures 4A,B).
Figure 4
An investigation into the common RNAs between time frames showed that 18 RNAs are consistently differentially expressed at all time frames, and that these RNAs were piRNAs and mRNA fragments. RNAs from the 2–5 years time frame, which has the highest number of RNAs above our cut-off, consisted of primarily piRNAs (Figures 4C,D).
mRNA Fragments Have Consistently Reduced Levels in TGCT Cases Compared to Controls
To investigate the stability of TGCT RNAs over time, we looked at nine RNAs with the lowest adjusted p-values which appeared in all four time frames and their associated log2fc. Primarily, the significantly reduced RNAs were mRNA fragments which showed a relatively stable log2fc over time. BHLHE41 showed some dynamic changes across the time frames with a log2fc range between −1.91 to −3.45, however, this could still be considered stable due to the overall log2fc remaining negative. Of all the mRNAs, a slight overall increase in log2fc was observed over time, with a mean log2fc of −2.91 at 0–2 years and −2.52 at 8–10 years. However, the lowest average log2fc appears at 2–5 years with a mean of −3.08, followed by a mean of −3.06 at 5–8 years. This showed that during the extended time frame of 0–8 years, there was a stability in the RNA levels, with a slight increase at 8–10 years (Figure 5).
Figure 5
Common Cancer Related Pathways Were Enriched in Specific Time Periods and Histologies
Enrichment analysis showed multiple pathways enriched for significant miRNA targets and mRNA fragments. Of the significant pathways, only one [Axon guidance [mean p-adjusted = 0.0094]] appeared in all seven subgroups, and two appeared in six of the seven sub-analyses [Ras signaling pathway [mean p-adj = 0.10], MAPK signaling pathway [mean p-adj = 0.043]]. In total there are 19 pathways that span at least three of the time frames, with the remaining six only appearing in two of the four. Non-seminoma and 0–2 years time frame had the lowest number of significant enriched pathways with only three pathways in each.
Looking at each sub analysis, enrichment of miRNAs and mRNA fragments among the 88 significant RNAs in all TGCTs revealed several cancer related pathways, including mTOR, MAPK and ErbB2 (Figure 6). In seminomas the expressed miRNAs and mRNA fragments also showed the cancer related pathways mTOR, MAPK and AMPK as well as the proteoglycans in cancer pathway. Similar pathways were present in non-seminomas as well, with miRNAs and mRNA fragments showing cancer related pathways, with mTOR, MAPK, ERBb2, and PI3K-Akt showing highest significance (Figure 6). Finally, all four time frames showed significant pathways, including the mTOR signaling pathway which was one of the top 10 more significant pathways throughout, with the MAPK and Erb2 pathways also showing up once more within the time frames (Figure 6). Visualization of the pathway “Pathway in Cancer,” selected since it summarizes the most relevant cancer paths, showed us which of the associated genes were present throughout the time frames and across histologies (Figures 7, 8). These pathview plots show a notable difference between seminoma and non-seminoma. For example the KITLG-IGFR/EGFR-Grb2-Sos-Raf-Ras-ERK path is specific to seminomas. Of particular interest to us was the presence of the PTEN in the pathways, known from TGCT-genome-wide association studies, which were significantly associated with all four time periods and both histologies.
Figure 6
Figure 7
Figure 8
Discussion
We found the RNA profiles in pre-diagnostic serum of TGCT patients to be different from that of controls, and observed that the difference was relatively stable across the pre-diagnostic time period of 10 years. Although previous studies have identified TGCT specific RNA profiles in patients (–), this has not been described for the pre-diagnostic period, as seen in our study. In contrast to our findings, highly dynamic pre-diagnostic RNA signals have been observed in lung cancer and breast cancer samples (, 47, 48). The reason for this contrast may partly be due to the localization of the tumor within the testis. Typically, the developed GCNIS is found above the basement membrane in seminiferous tubules, which are located in the intratesticular environment and therefore separate from the bloodstream via the blood testis barrier (BTB) (49, 50). The BTB's primary function is to prevent both antibodies and T lymphocytes from affecting the testis (51). Furthermore, angiogenesis in breast and lung cancer is well-documented, showing a closer physical interaction between these tumors and circulatory system than what applies to TGCT (52, 53). GCNIS may also be less affected by the immune system (54, 55) as the BTB provides adequate protection against immune response during the early stages of the malignant growth. Another potential implication of the presence of the BTB function is a delay in serum signals. Changes in the GCNIS embedded in the seminiferous tubules do not necessarily appear in blood unless the tumor causes disruption of the BTB (56, 57). Further studies are needed to show how the interference of TGCT with the BTB affects the levels of circulating RNAs.
In this study we found that differentially expressed levels of mRNAs were related to both testis development and known cancer pathways, as well as the lncRNA, TSIX which exhibited increased levels in TGCT patients. The role of TSIX in X-chromosome inactivation (XCI) has been widely reported (58–60), acting as the antisense repressor of XCI. The expression of TSIX is limited to pluripotent cells and the testis. A zinc finger protein, GLI, binds to the 5′ ends of TSIX to reduce TSIX expression and block the initiation of XCI. Physiologically, XCI is observed during spermatogenesis, but can also be expressed in TGCTs through supernumerary X chromosome constitution (61). XCI has been observed in TGCT of different histogenesis and may be a potential biomarker for some types of TGCT (62).
Alongside increased levels of TSIX, we also observed a positive log2fc with RAB21, an oncogene involved in mediating endocytosis and connected to the tumorigenesis and vesicle transport mechanisms of the other Rab proteins (63). The silencing of RAB21 has also been utilized to induce apoptosis in glioma cells and can also be used to significantly inhibit cell growth (64). The increased presence of RAB21 in the serum of patients after diagnosis of TGCT has not been observed before, however, previous studies have elucidated the function of a similar RAS family associated protein, RAB12, in rat testis development, especially noting the high expression in Sertoli cells (SC) (65). RAB GTPases have been identified as mediators for vesicle trafficking in cancer, including both RAB21 and RAB12 (63).
Investigating of the different histology patterns showed that despite the increased positive log2fc in non-seminomas, there are still RNA signals, such as mRNA TEX101, that are independent of the histology. TEX101 is involved in the production of spermatozoa, and studies using animal models have shown that males with a disrupted TEX101 gene often produce spermatozoa that are unable to fertilize, despite looking normal (66). Previous studies have found links between TGCT and changes in fertility (67). These findings include men with TGCT having fewer children than average, as well as a lower proportion of those children being male. Other studies have found that abnormal semen characteristics can be observed in men who later develop TGCT, indicating that these two aspects are aetiologically linked (68).
As hypothesized previously, the functions of the mRNAs identified as stable across all time frames could give some indication as to how GCNIS develops into tumors. Fragments of the mRNA for AP000295.1 had lower levels in cases compared to controls. AP000295.1 mRNA expression has been previously noted in the Human Protein Atlas (HPA) project, and HPA shows RNA-seq tissue data with high protein-coding transcripts per million (pTPM) counts of AP000295.1 in the blood constituents, granulocytes and monocytes (69). AP000295.1 belongs to a network of genes associated with kinase binding and type I interferon binding, AP000295.1 is also a paralog of the gene IFNAR2 with 100% match between target and query genes as seen using Ensembl (release 99) (70). Of note here is that IFNAR2 has been detected on the surface of pre-sertoli cells (pSC) (71). Previous studies have demonstrated the importance of Sertoli cells in the BTB formation and homeostasis (72), thus potentially explaining the presence of an RNA signal in serum despite the impermeable BTB. Through the reduction in expression of IFNAR2 in pSC's surface, inadequate formation and maintenance of the BTB could occur, causing a breakdown of this filter between the testis and the blood, thus explaining higher presence of signals of tumor development in serum.
Circulating miRNAs have been identified as both prognostic and diagnostic markers in biliary tract cancers (BTC), and this has led to earlier diagnosis and less invasive procedure of BTC than previously used techniques (). The circulating miRNAs used for BTC diagnosis were reported to be stable in serum and showed significantly different miRNA expression in patients with BTC and controls, and there also may be ethnic differences. Circulating biomarkers has also been identified in colorectal cancer (CRC), specifically the exosomal miRNAs, miR-150-5p, and miR-99b-5p were found to be downregulated in CRC patients compared to healthy individuals (, 73). This demonstrates the potential for their later use in clinical settings as a non-invasive diagnostic technique.
Pathways enrichment analysis showed overall many associated cancer related pathways across the 10-years period prior to TGCT diagnosis. The time frames 2–5 and 5–8 years prior to diagnosis showed the highest numbers of enriched pathways, many of which were cancer related. A possible explanation for the absence of enriched pathways in the 0–2 year time interval, could be an overall increase of all RNA types in serum as the tumor develops, this would effectively mask cancer pathways from being detected due to the large amounts of noise from aberrant RNA expression.
The presence of PTEN and KRAS genes in our pathway maps (Figures 7, 8), shows that these genes have been significantly differentially expressed in the prediagnostic cases over controls and can be found in all four time frames and both histologies. PTEN is a tumor suppressor gene that has been previously found to play an important role in the progression of intratubular germ cell neoplasias into the more mature invasive germ cell tumors. PTEN is often highly expressed in normal germ cells, but it is absent in the most common forms of TGCT (74).
Genome-wide association studies for TGCT have consistently shown that KITLG (ligand for the receptor tyrosine kinase) is implicated in TGCT susceptibility and with high effect sizes (75, 76). In these studies, the marker genotypes in KITLG were associated with both seminoma and non-seminoma TGCT without indication that genotype associations differed between the two subtypes. Our data thus suggest that epigenetics might play a more important role in the heterogeneity between seminoma and non-seminoma than hitherto recognized. More specifically, our results fit nicely with the TGCT histology specific DNA methylation patterns for KIT, shown by Shen et al. (77). Shen et al. (77) showed histology specific DNA methylation patterns for KIT and we have now shown differential RNA expression of this pathway to be specific to seminomas. Together this indicates histology-specific epigenetic mechanisms at play in the KIT–TGCT association.
A strength of this study is the large and robust RNA dataset. To our knowledge this is the largest study of the RNA profiles from TGCT patients prior to diagnosis also including harmonized confounders collected from health surveys (). The time interval throughout which the samples were collected is also a strength of this study. With up to 10 years time before diagnosis we were able to observe TGCT development and elucidate the body's response.
The sample size is still a limitation, and stratified analyses on histology and time to diagnosis could benefit from an increase in statistical power. Serum samples and survey data were collected over a period of more than 30 years, and lifestyle factors such as smoking and BMI have changed within this time period as well as carcinogens in the environment, such as organochlorine pesticides (78–80). Cases and controls were matched based on year of donation/storage time and sampling protocol, to minimize the possible bias of temporal/secular trends of exposure variables that could have biased the study.
Primarily, the new insights into TGCT carcinogenesis helps us better understand the development of the disease, which is thought to be initiated in utero. Stable serum RNAs in the pre-diagnostic period have the potential to be biomarkers for earlier detection. Early diagnosis of TGCT would lead to less use of cisplatin treatment, and thereby reduce long-term adverse effects, such as risk of second primary cancer and cardiovascular disease. Hazard ratios for all secondary cancers after a single cisplatin-based chemotherapy cycle were significantly lower than hazard ratios after two or more cycles (81).
In conclusion, RNA levels associated with cancer related pathways were different between individuals who developed TGCT, compared to matched controls. There is some loss of pathway signals closer to diagnosis, however, there is inherent reduced levels of serum RNAs in cases compared to controls for most of the 10 years pre-diagnostic follow-up time. The presence of these stable RNA signals may help in the identification of biomarkers give insights to molecular mechanisms driving TGCT development.
Statements
Data availability statement
The datasets generated for this article are not readily available because of the principles and conditions set out in articles 6 (1) (e) and 9 (2) (j) of the General Data Protection Regulation (GDPR). National legal basis as per the Regulations on population-based health surveys and ethical approval from the Norwegian Regional Committee for Medical and Health Research Ethics (REC) is also required. Requests to access the datasets should be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Norwegian regional committee for medical and health research ethics. The patients/participants provided their written informed consent to participate in this study.
Author contributions
TH and TR designed the study. JB and SU performed the analyses of the data. JB, TH, and TR drafted the manuscript. All authors discussed the results, contributed to the writing, and approved the final manuscript.
Funding
This work was supported by the Research Council of Norway's program Human Biobanks and Health Data (229621/H10 and 248791/H10) and internal funds of OsloMet—Oslo Metropolitan University and Cancer Registry of Norway.
Acknowledgments
We would like to acknowledge Cecilie Bucher-Johannessen, Marianne Lauritzen, Kari Furu, and Magnus Leithaug for performing laboratory and coordination tasks. We acknowledge the Norwegian Institute of Public Health for access to survey data in this study. The sequencing service was provided by the Norwegian Sequencing Center (http://www.sequencing.uio.no), a national technology platform hosted by Oslo University Hospital and the University of Oslo supported by the Research Council of Norway and the Southeastern Regional Health Authority. This manuscript has been released as a pre-print at medRxiv (82).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.574977/full#supplementary-material
References
1.
TrabertBChenJDevesaSSBrayFMcGlynnKA. International patterns and trends in testicular cancer incidence, overall and by histologic subtype, 1973-2007. Andrology. (2015) 3:4–12. 10.1111/andr.293
2.
ZnaorALortet-TieulentJJemalABrayF. International variations and trends in testicular cancer incidence and mortality. Eur Urol. (2014) 65:1095–106. 10.1016/j.eururo.2013.11.004
3.
Elzinga-TinkeJEDohleGRLooijengaLH. Etiology and early pathogenesis of malignant testicular germ cell tumors: towards possibilities for preinvasive diagnosis. Asian J Androl. (2015) 17:381–93. 10.4103/1008-682X.148079
4.
HemminkiKLiX. Familial risk in testicular cancer as a clue to a heritable and environmental aetiology. Br J Cancer. (2004) 90:1765–70. 10.1038/sj.bjc.6601714
5.
DongCLönnstedtIHemminkiK. Familial testicular cancer and second primary cancers in testicular cancer patients by histological type. Eur J Cancer. (2001) 37:1878–85. 10.1016/S0959-8049(01)00172-1
6.
KristiansenWKarlssonRRoungeTBWhitingtonTAndreassenBKMagnussonPKet al. Two new loci and gene sets related to sex determination and cancer progression are associated with susceptibility to testicular germ cell tumor [Internet]. Hum Mol Genet. (2015) 24:4138–46. 10.1093/hmg/ddv129
7.
LitchfieldKLovedayCLevyMDudakiaDRapleyENsengimanaJet al. Large-scale sequencing of Testicular Germ Cell Tumour (TGCT) cases excludes major TGCT predisposition gene. Eur Urol. (2018) 73:828–831. 10.1016/j.eururo.2018.01.021
8.
DasMKKleppaLHaugenTB. Functions of genes related to testicular germ cell tumour development. Andrology. (2019) 7:527–35. 10.1111/andr.12663
9.
WangZMcGlynnKARajpert-De MeytsEBishopDTChungCCDalgaardMDet al. Meta-analysis of five genome-wide association studies identifies multiple new loci associated with testicular germ cell tumor. Nat Genet. (2017) 49:1141–7. 10.1038/ng.3879
10.
FergusonLAgoulnikAI. Testicular cancer and cryptorchidism. Front Endocrinol. (2013) 4:32. 10.3389/fendo.2013.00032
11.
OosterhuisJWWolter OosterhuisJLooijengaLHJ. Testicular germ-cell tumours in a broader perspective. Nat Rev Cancer. (2005) 5:210–22. 10.1038/nrc1568
12.
TramaAFoschiRLarrañagaNSantMFuentes-RaspallRSerrainoDet al. Survival of male genital cancers (prostate, testis and penis) in Europe 1999-2007: results from the EUROCARE-5 study. Eur J Cancer. (2015) 51:2206–16. 10.1016/j.ejca.2015.07.027
13.
MeytsER-DMcGlynnKAOkamotoKJewettMASBokemeyerC. Testicular germ cell tumours. Lancet. (2016) 387:1762–74. 10.1016/S0140-6736(15)00991-5
14.
HaugnesHSBoslGJBoerHGietemaJABrydøyMOldenburgJet al. Long-term and late effects of germ cell testicular cancer treatment and implications for follow-up. J Clin Oncol. (2012) 30:3752–63. 10.1200/JCO.2012.43.4431
15.
RaphaelMJWeiXKarimSRobinsonAGBedardPLBoothCM. Neurotoxicity among survivors of testicular cancer: a population-based study. Clin Oncol. (2019) 31:653–8. 10.1016/j.clon.2019.04.008
16.
Bucher-JohannessenCPageCMHaugenTBWojewodzicMWFossåSDGrotmolTet al. Cisplatin treatment of testicular cancer patients introduces long-term changes in the epigenome. Clin Epigenet. (2019) 11:179. 10.1186/s13148-019-0764-4
17.
LagunaMPAlbersPAlgabaFBokemeyerCBoormansBLFischerSet al. European Association of Urology Guidelines. Edition. Arnhem: European Association of Urology Guidelines Office (2020).
18.
RichieJPBirnholzJGarnickMB. Ultrasonography as a diagnostic adjunct for the evaluation of masses in the scrotum. Surg Gynecol Obstet. (1982) 154:695–8.
19.
Germà-LluchJ. Clinical pattern and therapeutic results achieved in. patients with germ-cell tumours of the testis: the experience of the spanish germ-cell Cancer Group (GG). Eur Urol. (2002). 42:553–63. 10.1016/S0302-2838(02)00439-6
20.
MurrayMJHuddartRAColemanN. The present and future of serum diagnostic tests for testicular germ cell tumours. Nat Rev Urol. (2016) 13:715–25. 10.1038/nrurol.2016.170
21.
WuDShenNLinXChenX. Prepubertal testicular tumors in China: a 10-year experience with 67 cases. Pediatr Surg Int. (2018) 34:1339–43. 10.1007/s00383-018-4366-6
22.
AkinyemijuTMooreJXJuddSEPisuMGoodmanMHowardVJet al. Pre-diagnostic biomarkers of metabolic dysregulation and cancer mortality. Oncotarget. (2018) 9:16099–109. 10.18632/oncotarget.24559
23.
UmuSULangsethHBucher-JohannessenCFrommBKellerAMeeseEet al. A comprehensive profile of circulating RNAs in human serum. RNA Biol. (2018) 15:242–50. 10.1080/15476286.2017.1403003
24.
Fernandez-MercadoMManterolaL. The circulating transcriptome as a source of non-invasive cancer biomarkers: concepts and controversies of non-coding and coding RNA in body fluids. J Cell. (2015) 19:2307–23. 10.1111/jcmm.12625
25.
LetelierPRiquelmeIHernándezAHGuzmánNFaríasJGRoaJC. Circulating microRNAs as biomarkers in biliary tract cancers. Int J Mol Sci. (2016) 17:791. 10.3390/ijms17050791
26.
ZhaoYJSongXNiuLTangYSongXXieL. Circulating exosomal miR-150-5p and miR-99b-5p as diagnostic biomarkers for colorectal cancer. Front Oncol. (2019) 9:1129. 10.3389/fonc.2019.01129
27.
UmuSULangsethHKellerAMeeseEHellandÅLyleRet al. A 10-year prediagnostic follow-up study shows that serum RNA signals are highly dynamic in lung carcinogenesis. Mol Oncol. (2019) 14:235–47. 10.1002/1878-0261.12620
28.
PalmerRDMurrayMJSainiHKvan DongenSAbreu-GoodgerCMuralidharBet al. Malignant germ cell tumors display common microRNA profiles resulting in global changes in expression of messenger RNA targets. Cancer Res. (2010) 70:2911–23. 10.1158/0008-5472.CAN-09-3301
29.
MurrayMJHalsallDJHookCEWilliamsDMNicholsonJCColemanN. Identification of microRNAs From the miR-371~373 and miR-302 clusters as potential serum biomarkers of malignant germ cell tumors. Am J Clin Pathol. (2011) 135:119–25. 10.1309/AJCPOE11KEYZCJHT
30.
LooijengaLHJGillisAJMStoopHJHersmusROosterhuisJW. Chromosomes and expression in human testicular germ-cell tumors: insight into their cell of origin and pathogenesis. Ann N Y Acad Sci. (2007) 1120:187–214. 10.1196/annals.1411.000
31.
RoungeTBLauritzenMLangsethHEnerlyELyleRGislefossRE. microRNA biomarker discovery and high-throughput DNA sequencing are possible using long-term archived serum samples [Internet]. Cancer Epidemiol Biomarkers Prev. (2015) 24:1381–7. 10.1158/1055-9965.EPI-15-0289
32.
FerreiraHJHeynHGarcia del MuroXVidalALarribaSMuñozCet al. Epigenetic loss of the PIWI/piRNA machinery in human testicular tumorigenesis. Epigenetics. (2014) 9:113–8. 10.4161/epi.27237
33.
RoungeTBFuruKSkotheimRIHaugenTBGrotmolTEnerlyE. Profiling of the small RNA populations in human testicular germ cell tumors shows global loss of piRNAs. Mol Cancer. (2015) 14:153. 10.1186/s12943-015-0411-4
34.
LangsethHGislefossREMartinsenJIDillnerJUrsinG. cohort profile: the janus serum bank cohort in norway. Int J Epidemiol. (2017) 46:403–4g. 10.1093/ije/dyw027
35.
HjerkindKVGislefossRETretliSNystadWBjørgeTEngelandAet al. Cohort profile update: the janus serum bank cohort in norway. Int J Epidemiol. (2017) 46:1101–2f. 10.1093/ije/dyw302
36.
RoungeTBUmuSUKellerAMeeseEUrsinGTretliSet al. Circulating small non-coding RNAs associated with age, sex, smoking, body mass and physical activity. Sci Rep. (2018) 8:17650. 10.1101/247155
37.
SchubertMLindgreenSOrlandoL. Adapter removal v2: rapid adapter trimming, identification, and read merging. BMC Res Notes. (2016) 9:88. 10.1186/s13104-016-1900-2
38.
KozomaraAGriffiths-JonesS. miRBase: annotating high confidence microRNAs using deep sequencing data [Internet]. Nucleic Acids Res. (2014) 42:D68–73. 10.1093/nar/gkt1181
39.
ZhangPSiXSkogerbøG.WangJCuiDLiYet al. piRBase: a web resource assisting piRNA functional study. Database. (2014) 2014:bau110. 10.1093/database/bau110
40.
HarrowJFrankishAGonzalezJMTapanariEDiekhansMKokocinskiFet al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. (2012) 22:1760–74. 10.1101/gr.135350.111
41.
PantanoLEstivillXMartíE. SeqBuster, a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells. Nucleic Acids Res. (2010) 38:e34. 10.1093/nar/gkp1127
42.
LoherPTelonisAGRigoutsosI. MINTmap: fast and exhaustive profiling of nuclear and mitochondrial tRNA fragments from short RNA-seq data. Sci Rep. (2017) 7:41184. 10.1038/srep41184
43.
LoveMIHuberWAndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. 10.1186/s13059-014-0550-8
44.
WongNWangX. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. (2015) 43:D146–52. 10.1093/nar/gku1104
45.
LuoWBrouwerC. Pathview: an R/Biocondutor package for pathway-based data integration and visualization. Bioinformatics. (2013) 29:1830–1. 10.1093/bioinformatics/btt285
46.
LuoWPantGBhavnasiYKBlanchardSGBrouwerC. (2017) Pathview web: user friendly pathway visualization and data integration. Nucleic Acids Res.45:W501–8. 10.1093/nar/gkx372
47.
LundEHoldenLBøvelstadHPlancadeSModeNGüntherC-Cet al. A new statistical method for curve group analysis of longitudinal gene expression data illustrated for breast cancer in the NOWAC postgenome cohort as a proof of principle. BMC Med Res Methodol. (2016) 16:28. 10.1186/s12874-016-0129-z
48.
HoldenMHoldenLOlsenKSLundE. Local in Time Statistics for detecting weak gene expression signals in blood – illustrated for prediction of metastases in breast cancer in the NOWAC Post-genome Cohort. Adv Genom Genet. (2017) 7:11–28. 10.2147/agg.s130004
49.
SpillerCMBowlesJ. Germ cell neoplasia in situ: the precursor cell for invasive germ cell tumors of the testis. Int J Biochem Cell Biol. (2017) 86:22–5. 10.1016/j.biocel.2017.03.004
50.
ReuterVE. Origins and molecular biology of testicular germ cell tumors. Mod Pathol. (2005) 18:S51–60. 10.1038/modpathol.3800309
51.
BartJGroenHJMvan der GraafWTAHollemaHHendrikseNHVaalburgWet al. An oncological view on the blood–testis barrier. Lancet Oncol. (2002) 3:357–63. 10.1016/S1470-2045(02)00776-3
52.
HerbstRSOnnASandlerA. Angiogenesis and lung cancer: prognostic and therapeutic implications. J Clin Oncol. (2005) 23:3243–56. 10.1200/JCO.2005.18.853
53.
Longatt FilhoALopesJMSchmittFC. Angiogenesis and breast cancer. J Oncol. (2010). 10.1155/2010/576384
54.
EndoHInoueM. Dormancy in cancer. Cancer Sci. (2019) 110:474–80. 10.1111/cas.13917
55.
WeisSMChereshDA. A wake-up call for hibernating tumour cells [Internet]. Nat Cell Biol. (2013) 15:721–3. 10.1038/ncb2794
56.
ChengCYMrukDD. The blood-testis barrier and its implications for male contraception. Pharmacol Rev. (2012) 64:16–64. 10.1124/pr.110.002790
57.
KalavskaKKucerovaLSchmidtovaSChovanecMMegoM. Cancer stem cell niche and immune-active tumor microenvironment in testicular germ cell tumors. Adv Exp Med Biol. (2020) 1226:111–21. 10.1007/978-3-030-36214-0_9
58.
GayenSMaclaryEButtigiegEHintenMKalantryS. A primary role for the tsix lncrna in maintaining random X-chromosome inactivation. Cell Rep. (2015) 11:1251–65. 10.1016/j.celrep.2015.04.039
59.
TurnerJMAMahadevaiahSKElliottDJGarchonH-JPehrsonJRJaenischRet al. Meiotic sex chromosome inactivation in male mice with targeted disruptions of Xist. J Cell Sci. (2002) 115:4097–105. 10.1242/jcs.00111
60.
Del RosarioBCDel RosarioAMAnselmoAWangPISadreyevRILeeJT. Genetic intersection of tsix and hedgehog signaling during the initiation of X-chromosome inactivation. Dev Cell. (2017) 43:359–71.e6. 10.1016/j.devcel.2017.09.027
61.
LoboJNunesSPGillisAJMBarros-SilvaDMiranda-GonçalvesVvan den BergAet al. XIST-promoter demethylation as tissue biomarker for testicular germ cell tumors and spermatogenesis quality. Cancers. (2019) 11:1385. 10.3390/cancers11091385
62.
LooijengaLHJGillisAJvan GurpRJVerkerkAJOosterhuisJW. X inactivation in human testicular tumors. XIST expression and androgen receptor methylation status. Am J Pathol. (1997) 151:581–90.
63.
TzengH-TWangY-C. Rab-mediated vesicle trafficking in cancer. J Biomed Sci. (2016) 23:70. 10.1186/s12929-016-0287-7
64.
GeJChenQLiuBWangLZhangSJiB. Knockdown of Rab21 inhibits proliferation and induces apoptosis in human glioma cells. Cell Mol Biol Lett. (2017) 22:30. 10.1186/s11658-017-0062-0
65.
IidaHNodaMKanekoTDoiguchiMMoriT. Identification of rab12 as a vesicle-associated small GTPase highly expressed in sertoli cells of rat testis. Mol Reprod Dev. (2005) 71:178–85. 10.1002/mrd.20294
66.
FujiharaYTokuhiroKMuroYKondohGArakiYIkawaMet al. Expression of TEX101, regulated by ACE, is essential for the production of fertile mouse spermatozoa. Proc Natl Acad Sci USA. (2013) 110:8111–6. 10.1073/pnas.1222166110
67.
SkakkebækNERajpert-De MeytsEMainKM. Testicular dysgenesis syndrome: an increasingly common developmental disorder with environmental aspects. Hum Reprod. (2001) 16:972–8. 10.1093/humrep/16.5.972
68.
JacobsenR. Risk of testicular cancer in men with abnormal semen characteristics: cohort study. BMJ. (2000) 321:789–92. 10.1136/bmj.321.7264.789
69.
UhlenMFagerbergLHallstromBMLindskogCOksvoldPMardinogluAet al. Tissue-based map of the human proteome. Science. (2015) 347:1260419. 10.1126/science.1260419
70.
CunninghamFAchuthanPAkanniWAllenJAmodeMRArmeanIMet al. Ensembl. Nucleic Acids Res. (2019). 47:D745–51. 10.1093/nar/gky1113
71.
EdgarRMazorYRinonABlumenthalJGolanYBuzhorEet al. Life Map discoverytm: the embryonic development, stem cells, and regenerative medicine research portal. PLoS ONE. (2013) 8:e66629. 10.1371/journal.pone.0066629
72.
GerberJHeinrichJBrehmR. Blood–testis barrier and sertoli cell function: lessons from SCCx43KO mice. Reproduction. (2016) 151:R15–27. 10.1530/REP-15-0366
73.
FerracinMLupiniLMangoliniANegriniM. Circulating non-coding RNA as biomarkers in colorectal cancer. Adv Exp Med Biol. (2016) 937:171–81. 10.1007/978-3-319-42059-2_9
74.
Di VizioDCitoLBocciaAChieffiPInsabatoLPettinatoGet al. Loss of the tumor suppressor gene PTEN marks the transition from intratubular germ cell neoplasias (ITGCN) to invasive germ cell tumors. Oncogene. (2005) 24:1882–94. 10.1038/sj.onc.1208368
75.
KanetskyPAMitraNVardhanabhutiSLiMVaughnDJLetreroRet al. Common variation in KITLG and at 5q31.3 predisposes to testicular germ cell cancer. Nat Genet. (2009) 41:811–815. 10.1038/ng.393
76.
RapleyEATurnbullCAl OlamaAADermitzakisETLingerRHuddartRAet al. A genome-wide association study of testicular germ cell tumor. Nat Genet. (2009) 41:807–10. 10.1038/ng.394
77.
ShenHShihJHollernDPWangLBowlbyRTickooSKet al. Integrated molecular characterization of testicular germ cell tumors. Cell Rep. (2018) 23:3392–406. 10.1016/j.celrep.2018.05.039
78.
BelpommeDIrigarayPHardellLClappRMontagnierLEpsteinSet al. The multitude and diversity of environmental carcinogens. Environ Res. (2007) 105:414–29. 10.1016/j.envres.2007.07.002
79.
McGlynnKATrabertB. Adolescent and adult risk factors for testicular cancer [Internet]. Nat Rev Urol. (2012) 9:339–49. 10.1038/nrurol.2012.61
80.
IrigarayPNewbyJAClappRHardellLHowardVMontagnierLet al. Lifestyle-related factors and environmental agents causing cancer: an overview. Biomed Pharmacother. (2007) 61:640–58. 10.1016/j.biopha.2007.10.006
81.
HellesnesRKvammenØMyklebustTÅBremnesRMKarlsdottirÁNegaardHFSet al. Continuing increased risk of second cancer in long-term testicular cancer survivors after treatment in the cisplatin era. Int J Cancer. (2020) 147:21–32. 10.1002/ijc.32704
82.
BurtonJUmuSULangsethHGrotmolTGrimsrudTGHaugenTBet al. Serum RNA profiling in the 10-year period prior to diagnosis of testicular germ cell tumour. medRxiv. (2020) 147:21–32. 10.1101/2020.06.03.20121087
Summary
Keywords
RNA profiling, serum, testicular cancer, pre-diagnostic, seminoma, non-seminoma, sequencing, miRNA
Citation
Burton J, Umu SU, Langseth H, Grotmol T, Grimsrud TK, Haugen TB and Rounge TB (2020) Serum RNA Profiling in the 10-Years Period Prior to Diagnosis of Testicular Germ Cell Tumor. Front. Oncol. 10:574977. doi: 10.3389/fonc.2020.574977
Received
01 July 2020
Accepted
18 September 2020
Published
28 October 2020
Volume
10 - 2020
Edited by
Margaret Jane Currie, University of Otago, New Zealand
Reviewed by
Rossano Lattanzio, University of Studies G. d'Annunzio Chieti and Pescara, Italy; Adam J. Dupuy, The University of Iowa, United States
Updates
Copyright
© 2020 Burton, Umu, Langseth, Grotmol, Grimsrud, Haugen and Rounge.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Trine B. Haugen tribha@oslomet.noTrine B. Rounge trine.rounge@kreftregisteret.no
This article was submitted to Cancer Genetics, a section of the journal Frontiers in Oncology
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.