Introduction

Intracranial germ cell tumors (iGCT) constitute a heterogeneous group of brain tumors affecting children, adolescents, and young adults. Histologically, they comprise pure germinoma (GER) and non-germinomatous tumors (NGGCT), such as embryonal carcinoma (EC), choriocarcinoma (CHC), yolk sac tumor (YST) or teratoma (Murray et al. 2014). Localized in brain midline structures, especially in the pineal and suprasellar regions, they may present as bifocal iGCT at both sites simultaneously (Calaminus et al. 1994).

One-fifth of iGCT are detectable by elevated tumor markers in serum or cerebrospinal fluid (CSF) classifying these tumors as marker-positive iGCT (Calaminus et al. 1997). In all other iGCT diagnosis relies on invasive biopsy, which although nowadays an overall safe procedure still carries a risk (Schulz et al. 2021). Moreover, biopsies may not represent the entire tumor as the histologic composition of iGCT may be heterogeneous (Takami et al. 2019). Platin-based chemotherapy, surgery of residual lesions, and adjuvant radiotherapy are the treatment modalities of choice in iGCT in Europe. Yet, 10% of GER and 30% of NGGCT relapse within 5 years after intense multimodality therapy (Calaminus et al. 2017, 2013). As discrimination from other intracranial tumors and early detection are essential for immediate treatment initiation in iGCT, less invasive methods are warranted that facilitate reliable diagnosis.

MicroRNAs (miR/miRNA) are small non-coding RNAs that have emerged as non-invasive biomarkers. They are involved in the regulation of gene expression and are frequently dysregulated in cancers. Among these, the miR-371~373 and miR-302/367-clusters are overexpressed in malignant extracranial GCT regardless of the histologic subtype, tumor site, or patient age (Palmer et al. 2010; Gillis et al. 2007). We and others have previously analyzed the miR-371~373-cluster and miR-367-3p in serum of adults suffering from testicular cancer and were able to discriminate GCT patients from healthy subjects with high sensitivity (85–98%) and specificity (94–99%) (Syring et al. 2015; Gillis et al. 2013; Dieckmann et al. 2019). Testicular and intracranial GCT showed identical histologic patterns and similar genomic profiles (Takami et al. 2019; Schneider et al. 2006; Ichimura et al. 2016). Therefore, analysis of the GCT-specific miR-371~373 and miR-302/367-clusters in serum and/or CSF of iGCT patients seems a promising approach that is underlined by a pioneering publication on two iGCT patients (Murray et al. 2016).

In a global miRNA expression analysis of 12 tumor specimens from primary pediatric iGCT, GER-specific upregulation of miR-142-5p and miR-146a-5p was documented in contrast to overexpression of miR-335-5p and miR-654-3p in NGGCT (Wang et al. 2010). Yet, it is still unknown whether detection of these miRNAs is discriminatory in serum of iGCT patients versus healthy controls.

Here we present a comprehensive miRNA analysis in a case series of iGCT patients exploring 12 different microRNAs (miR-371a-3p, miR-372-3p, miR-373-3p, miR-302a-3p, miR-302b-3p, miR-302c-3p, miR-302d-3p, miR-367-3p, miR-142-5p, miR-146a-5p, miR-335-5p and miR-654-3p) in both serum and CSF at time of diagnosis.

Patients and methods

Patient enrollment

The study was approved by the local ethics committee of the medical faculty of the Rheinische Friedrich-Wilhelms-Universität of Bonn, Germany (No. 048/18). Adult patients or parents of minors provided full informed consent. Serum and/or CSF of eight iGCT patients was collected at diagnosis for miRNA analysis in eight different oncology centers in Germany between September 2018 and January 2020. From these eight patients, serum (n = 6) and/or CSF (n = 4) was obtained prior to treatment. In two iGCT patients, both serum and CSF was available at diagnosis. Diagnosis of NGGCT (n = 3) was based on elevated α-fetoprotein (AFP, > 25 ng/ml) and/or ß-human chorionic gonadotropin (HCG, > 50 IU/l) in serum and/or CSF. In marker-negative cases (n = 5) biopsy was performed and revealed pure GER. Mean patient age was 17.5 years (range 10.3–33.2). Patient characteristics are summarized in Table 1. Serum from healthy donors (n = 3) and one leukemia patient as well as CSF from children suffering from other malignant diseases (acute lymphoblastic leukemia, n = 5; other malignant brain tumors, n = 3) and from one case of pituitary inflammation was included as controls.

Table 1 Clinical data of eight patients with intracranial germ cell tumors (iGCT)

MicroRNA analysis

RNA isolation and microRNA quantification was performed based on a published protocol, which had already been implemented in miRNA analysis in GCT (Murray et al. 2016; Bell et al. 2017) and adopted in this study to achieve higher miRNA concentrations. Briefly, prior to RNA extraction MS2 carrier RNA (Roche, Welwyn Garden City, USA) was added to serum and CSF samples to increase RNA yield. To enhance RNA quantity for further microRNA analysis, RNA isolation was performed twice per sample and eluted in 20 µl each. Non-human spike-in cel-miR-39-3p was used to control for RNA extraction efficiency. Two samples presented an increased ΔCtmiR23a-miR-451 between 9 and 11 indicating some degree of hemolysis. We proceeded to further analyze these two cases as miR-30b-5p was in the range of the other patient and control samples. TaqMan® microRNA reverse transcription kit (Applied Biosystems, USA) was applied for multiplex reverse transcription of 5 µl RNA followed by a multiplex pre-amplification step with 12.5 µl cDNA using the TaqMan® PreAmp Master Mix (Applied Biosystems, USA). Afterwards, the final product of 50 µl was diluted 1:1 (50 µl + 50 µl) with RNAse-free water. TaqMan® Fast Advanced Master mix (Applied Biosystems, USA) was employed for the final singleplex qPCR of GCT-specific microRNAs being performed in duplicates and by adding 9 µl of pre-amplified cDNA to a final reaction volume of 20 µl per sample. Only samples yielding a maximum fluorescence value of at least 30 and exhibiting sigmoid-like shaped curves were included in the analysis. The mean of Ct values was normalized against the reference genes miR-30b-5p for serum and CSF as well as miR-204-5p for CSF (ΔCt). In addition, patient ΔCt values were standardized against ΔCt values of control samples and expressed as \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\). In this manuscript miRNAs targeting the 5′-untranslated region are labeled with “-5p” (miR-142-5p, miR-146a-5p, miR-335-5p, miR-30b-5p, miR-204-5p), miRNAs without additional notation bind to the 3′-untranslated region.

Statistical analysis

The Mann–Whitney U Test was used to determine differences between two groups, e.g., iGCT patients and controls. p values < 0.05 were considered statistically significant.

Results

Evaluation of reference genes miR-30b-5p and miR-204-5p

We compared Ct values of miR-30b-5p (serum + CSF) and miR-204-5p (CSF) for evaluation of reference gene expression (Fig. 1). Across all initial serum samples in patients (n = 6) and controls (n = 4) the median of miR-30b-5p Ct values was 15.48 (range 12.68–21.41) indicating robust reference gene expression in serum. In CSF, the median miR-30b-5p Ct value was 22.74 (range 14.02–26.41) of patients (n = 4) and controls (n = 9), while Ct values of the alternative reference gene miR-204-5p were significantly lower with a median of 18.14 and a narrower range (16.55–24.83; p = 0.007). Consequently, miR-204-5p was selected for normalization of target CSF miRNA.

Fig. 1
figure 1

Comparison of miRNA reference gene expression of miR-30b-5p and miR-204-5p in CSF samples. Ct values of miRNA reference genes miR-30b-5p and miR-204-5p of initial patient (n = 4) and control (n = 9) CSF samples are presented as median and range

Serum levels of miR-371a, miR-372, miR-367, miR-302a and miR-302d identify patients with intracranial germ cell tumors

In serum, Ct values were normalized to the reference gene miR-30b-5p (ΔCt). Among the 12 tested microRNAs, miR-371a (p = 0.0159), miR-372 (p = 0.0095), miR-367 (p = 0.0190), miR-302a (p = 0.0381) as well as miR-302d (p = 0.0159) differed significantly between iGCT patients and control subjects (Table 2, Fig. 2a). In contrast, ΔCt values of miR-373 (p = 0.0667) and miR-302b (p = 0.0952) known to be overexpressed in serum of patients with extracranial GCT failed to discriminate between patients suffering from iGCT and controls. MiR-302c-qPCR repetitively failed to present sigmoidal curves but exhibited rather linear slopes with a fluorescence level < 25 in both patients and controls. Consequently, analysis of miR-302c was entirely excluded from further analysis. MiR-142-5p, miR-146a-5p, miR-335-5p, and miR-654 differentially expressed in tissue of iGCT patients showed appropriate amplification curves yet without any significant difference between serum of patients and controls.

Table 2 Median ΔCt and range of 11 different microRNAs in serum of controls and patients suffering from intracranial germ cell tumors (iGCT)
Fig. 2
figure 2

Analysis of different serum miRNAs in patients diagnosed with iGCT and controls. a Mean of CT duplicates was normalized against the mean of CT values of the reference gene miR-30b-5p and expressed as a grey point (patient sample) or a clear triangle (control sample). Horizontal black lines indicate the median values of patient (n = 5–6) or control samples (n = 3–4). MiRNAs being statistically different overexpressed in iGCT patients compared to controls are marked with an asterisk. b Mean of CT duplicates was normalized against both miR-30b-5p and normalized CT values of controls and expressed as \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) level for every patient suffering from non-germinomatous (black bars) or germinomatous (grey bars) iGCT

For further analysis, patient ΔCt values were standardized to the median ΔCt of controls and expressed as relative miRNA level (\(2^{{ - \Delta \Delta C_{{\text{t}}} }}\); Fig. 2b). Abundant relative levels of miR-371a and miR-302d were found across all patient samples. The median relative level was 1333 for miR-371a (range: 102–5615) and 283 for miR-302d (94–2548). In addition, miR-372 and miR-367 displayed pronounced expression with median \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) level of 68 (11–116) and 19 (5–180), respectively. In miR-373 an increase of relative miRNA expression was observed in 4/6 iGCT patients with a median of 243 (115–395). With exception of NGGCT#1, miR-302b was markedly expressed in 5/6 patients with a median \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) level of 7 (2–73). Relative expression levels > 2 were also detected in 3/6 patients for miR-302a.

In keeping with results from testicular GCT on the 371~373-- and 302/367-clusters, \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) levels of each serum miRNA did not significantly differ between patients with non-germinomatous and germinomatous iGCT (each p > 0.5). Of note, in two biomarker-negative patients (GER#1, GER#2) without elevated AFP and/or ß-HCG in serum and/or CSF at diagnosis, miR-371a and miR-302d were markedly overexpressed and thus informative. In the other biomarker-negative patient GER#3, miR-371a and miR-302d yielded only Ct values of low fluorescence intensity. All other \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) levels of miRNAs in this patient from the 371~373- and 302/367-clusters were strikingly elevated ranging from 75–471 and thus allowing diagnosis of iGCT.

While expression of miR-142-5p and miR-146a-5p was enhanced in one patient (NGGCT#2) suffering from a NGGCT comprising CHC and GER histology, no patient presented elevated \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) levels of miR-335 and miR-654. As miR-142-5p, miR-146a-5p, miR-335 and miR-654 harbored no additional value on discrimination of iGCT from controls as already shown in Fig. 1, further miRNA analysis in CSF focused on the 371~373- and 302/367-clusters.

In CSF miR-371a, miR-372, miR-367 and miR-302d also allowed for discrimination of patients with intracranial germ cell tumors from controls

In CSF, the housekeeper miR-204-5p and most of the analyzed miRNAs presented with adequate maximum fluorescence levels. Only in miR-371a analysis of 5/9 control samples maximum fluorescence was exceedingly low and although considered negative had to be excluded from the analysis as no valid CT value could be attributed. The remaining four control and four patient samples were evaluable and ΔCt values of miR-371a were discriminatory (p = 0.0286). Discriminatory pattern of the other miRNAs in CSF was similar to serum as miR-372 (p = 0.0028), miR-367 (p = 0.0167) and miR-302d (p = 0.0061) in iGCT patients differed from controls. In contrast to serum, expression of miR-302b in CSF also separated patients from controls (p = 0.0040), while miR-373 did not differentiate between these groups (p = 0.2601; Table 3, Fig. 3a). In CSF in spite of adequate maximum fluorescence levels, miR-302c and miR-302a did not retrieve sigmoidal-like curves for both patients and controls with the established primer set resulting in the exclusion of these miRNAs from analysis.

Table 3 Median ΔCt and range of 7 different microRNAs in cerebrospinal fluid of controls and patients suffering from intracranial germ cell tumors (iGCT)
Fig. 3
figure 3

Analysis of different CSF miRNAs in patients diagnosed with iGCT and controls. a Mean of Ct duplicates was normalized against the mean of Ct values of the reference gene miR-204-5p and expressed as a grey point (patient sample) or a clear triangle (control sample). Horizontal black lines indicate the median values of patient (n = 3–4) or control samples (n = 4–9). MiRNAs being statistically different overexpressed in iGCT patients compared to controls are marked with an asterisk. b Mean of Ct duplicates was normalized against both miR-204-5p and normalized Ct values of controls and expressed as \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) level for every patient suffering from non-germinomatous (black bars) or germinomatous (grey bars) iGCT

Upon calculation of CSF miRNA levels in iGCT patients relative to controls employing the \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) method (Fig. 3b), \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) levels were highest in miR-371a across all patient samples with a median of 10,524 (range 203–27,600). Furthermore, every patient displayed very high relative levels of miR-372 (median, range: 61, 11–171), miR-367 (108, 19–273), miR-302b (9, 4–113) and miR-302d (103, 10–339). All germinoma patients but not the patient suffering from NGGCT exhibited overexpression of miR-373 with median relative levels of 3 (2–64). Of note, the three patients negative for the biomarkers AFP and ß-HCG from whom CSF was obtained at diagnosis (GER#1, GER#4, GER#5) exhibited pronounced expression of miR-372 and miR-302d with relative levels > 10 each.

Discussion

Previous studies have shown that members of the miR-371 ~ 373 and miR-302/367 clusters are frequently and aberrantly expressed in extracranial GCT. Therefore, analysis of these GCT-specific miRNAs in iGCT may be an innovative and promising diagnostic tool. Here, we present the to date largest study on miRNAs as liquid biomarkers in iGCT delineating GCT-specific miRNAs in serum and CSF.

In our study miR-371a, miR-372, miR-367, miR-302a and miR-302d in serum discriminated between iGCT of different histology and controls comprising mostly healthy donors. This is in line with our (Syring et al. 2015) and others’ miRNA analyses in patients suffering from testicular GCT including a total of > 1000 serum samples (Dieckmann et al. 2019, 2012; Agthoven and Looijenga 2017). In iGCT, relative miR-371a levels in two smaller case series were similar to our study but formerly only detectable in the serum of 2/5 patients (Murray et al. 2016, 2020), maybe explained by the adaption of the miRNA quantification protocol in our study.

The analysis of a panel of miRNAs may enhance the diagnostic accuracy. Thus, in our study serum miR-302d also distinguished iGCT patients from healthy controls with a comparable expression level to miR-371a. Although miR-302d is detectable in GCT tissues of both pediatric and adult patients and is overexpressed in YST compared to germinoma (Murray et al. 2010), data on its usefulness as liquid biomarker have been lacking. Serum miR-302d has not been investigated in testicular cancer patients focusing on the miR-371 ~ 373- cluster (Gillis et al. 2013; Dieckmann et al. 2019, 2012; Agthoven and Looijenga 2017) or has been excluded from detailed analysis in extracranial GCT as its levels ranged only within two Ct values between patients and non-tumor control (Murray et al. 2016). Here we document profound expression of both miR-371a and miR-302d in all iGCT sera including serum of biomarker-negative patients. Thus, not only miR-371a but also miR-302d added value in the diagnosis of iGCT especially in biomarker-negative patients, but did not outperform miR-371a as recently been noticed in a multi-institutional pooled miRNA analysis in testicular GCT (Piao et al. 2021).

MiRNA analysis may be influenced by various pre- and post-analytical variables and normalization strategies (Myklebust et al. 2019). The latter requires an endogenous stably expressed reference gene miRNA. MiR-30b-5p fulfilled this purpose in several analyses on testicular GCT in serum (Dieckmann et al. 2019; Agthoven and Looijenga 2017; Terbuch et al. 2018; Rosas Plaza et al. 2019) and was also employed in two studies on intracranial GCT both in serum as well as in CSF (Murray et al. 2016, 2020). However, miR-30b-5p was barely detectable in another analysis in 2/4 CSF samples from healthy donors by RT-qPCR (Lobo et al. 2019). In our study in a total of 13 different CSF samples, Ct values of miR-30b-5p markedly varied within a range of twelve Ct cycles underscoring the need for a more homogeneously expressed reference gene miRNA in CSF. Previous RNAseq-based miRNA examination of 129 samples of 12 different body fluids from healthy donors identified abundant miR-204-5p expression in 11 CSF samples, thus ranking among the top five miRNAs in CSF in this study (Godoy et al. 2018). Here, we confirm the suitability of miR-204-5p for normalization purposes in miRNA analysis of the CSF by demonstrating its higher expression and lower variability compared to the broadly used reference gene miR-30b-5p.

Of note, under these optimized conditions in our study miR-371a, miR-372, miR-367 and miR-302d in CSF presented as strong indicators of iGCT disease. In addition, miR-302b was beneficial in CSF to differentiate iGCT from other diseases, such as low-grade gliomas or pituitary inflammation. In another study analyzing miRNA expression in CSF of two iGCT patients and utilizing miR-30b-5p as housekeeper, only miR-371a and miR-372 discriminated patients from those with a benign brain tumor with miR-302d being not included in the analysis (Murray et al. 2016). In addition, miR-371a analysis in CSF seems to discriminate iGCT from Langerhans cell histiocytosis (LCH) as shown in a recent case report (Murray et al. 2020). Since LCH often presents with similar radiological changes as iGCT such as thickening of the pituitary gland, miRNA analysis may add value to the routine CSF assessment in the differential diagnosis of iGCT.

While healthy donors should be preferred as comparator group for normalization, sampling of CSF remains challenging, since lumbar puncture is mainly performed in symptomatic individuals. This issue was addressed in one study by normalization of CSF-levels against ΔCt values of serum rather than CSF controls (Murray et al. 2016). This resulted in comparatively high miR-371a expression up to \(2^{{ - \Delta \Delta C_{{\text{t}}} }}\) levels of 100,000 in the CSF of the two iGCT patients (Murray et al. 2016). Similarly, normalization of CSF values in our iGCT patients against the serum ΔCt values of our predominantly healthy controls resulted in extraordinary high relative miR-371a CSF levels (3500–470,000) in spite of raw CSF Ct values of miR-371a in our iGCT cohort left unchanged. Since miRNA distribution between CSF and serum differed in comparative microRNA analysis of healthy donors by RNAseq (Godoy et al. 2018) and qPCR (Weber et al. 2010), usage of serum ΔCt control values for normalization of CSF miRNA data remains a matter of debate and selection of appropriate controls is critical. Five of our nine CSF control samples presented an exceedingly low maximum fluorescence in the miR-371a-qPCR although fluorescence in qPCR of the housekeeper and the other miRNAs was sufficient. The amount of miR-371a in these samples seemed to be so low that miR-371a Ct values could not be determined and the samples had to be excluded from the analysis. In the other four control samples valid Ct values could be attributed. Thus, sufficient numbers of CSF control samples for miR-371a analysis are also essential.

In our cohort of iGCT patients we documented a similar GCT-specific miRNA expression pattern between CSF and serum indicating the passage of these miRNAs across the blood–brain barrier. Similarly, in testicular cancer circulating miR-371a is specifically derived from GCT tissue and its tissue levels correlate with levels of the tumor surrounding hydrocele and levels of serum being indicative for miRNA transfer through the blood-testis barrier (Belge et al. 2020; Dieckmann et al. 2016). However, low total RNA yields and different miRNA proportions of total RNA can affect miRNA analysis in CSF. We employed a well-established, robust and quality-controlled pipeline for miRNA analysis including pre-amplification of samples (Murray et al. 2016) and increased total RNA yield by performing RNA isolation twice. Nevertheless, miR-371a analysis in CSF was impaired by low amount of miR-371a in CSF control samples and qPCR of miR-302a and miR-302c repetitively revealed non-sigmoidal but linear amplification curves in both patients and controls. As analysis of GCT-specific miRNAs in serum showed a similar pattern to miRNA assessment in CSF, serum may be preferred to CSF as it yields more RNA and can be obtained by a less invasive and easier access in both patients and healthy controls.

In iGCT, an important clinical question is the distinction between pure germinoma and NGGCT as the latter require a more intensive treatment due to inferior outcomes (Calaminus et al. 2017, 2013). Known GCT-specific miRNAs such as miR-371~ 373 do not permit discrimination between these two iGCT subtypes in our and two other recent studies (Murray et al. 2016, 2020). Thus, we analyzed miR-142-5p and miR-146a-5p in serum of iGCT patients, which were highly expressed in germinoma compared to non-germinoma tissue in a global miRNA analysis of 12 iGCT (Wang et al. 2010). Of note, as a key regulator of anti-tumor immunity, miR-142-5p inhibited PD-L1 expression on tumor cells in vitro enhancing tumor infiltration by CD4+- and CD8+-lymphocytes (Jia et al. 2017). In germinoma, detection of miR-142-5p may, therefore, reflect the immunohistochemically documented high lymphocytic infiltration in the tumor microenvironment in the absence of PD-L1 (Zapka et al. 2018). In our study overexpression of miR-142-5p and miR-146a-5p was detectable in the serum of 1/3 patients who suffered from an iGCT with partly germinoma histology. As miR-142-5p and miR-146a-5p may also be highly expressed in the peripheral blood of healthy subjects due to their regulating role in the immune response, these miRNAs may only discriminate between germinoma and NGGCT in tissue but not in serum.

Perspectively, serum miRNAs of the C19MC cluster may prove beneficial to identify NGGCT, as these miRNAs were indicative for pure embryonal carcinoma histology not only in tissue but also in the serum of patients suffering from testicular GCT (Flor et al. 2016). In a large miRNA tissue analysis of 103 testicular GCT, overexpressed miRNAs from the C19MC cluster were also indicative for NGGCT (Qin et al. 2020).

In the largest data set so far on miRNAs as liquid biomarkers in iGCT, we underlined the utility of miR-371a, miR-372, miR-367 and miR-302d in patient serum for diagnosis of iGCT, a result particularly noteworthy in germinoma patients negative for other biomarkers, such as AFP and ß-HCG. Furthermore, we confirmed the suitability of these miRNAs as liquid biomarkers in both serum and CSF of iGCT patients. Using our miRNA protocol, diagnosis of iGCT in serum samples was a feasible and valid approach, particularly as serum can be readily obtained from patients and healthy controls by a less invasive procedure. This approach complements miRNA analysis in CSF, which may be complicated due to low RNA yield. Our study has some limitations. First, only a small number of patients could be included due to the rarity of iGCT. Second, the multi-institutional approach may have resulted in different pre-analytical sample handling including different time intervals until centrifugation and cryopreservation of probes, which could have affected miRNA quality.

Verification of miR-371 ~ 373 and miR-302d for iGCT diagnosis in larger cohorts of iGCT patients and control subjects suffering from brain tumors of different histologies may facilitate the miRNA analysis into routine clinical practice with the potential to replace neurosurgical interventions by liquid biopsy. Future iGCT studies should also be geared to investigate other miRNAs in serum and/or CSF that enable for discrimination of germinoma from NGGCT. Furthermore, miRNA analysis during therapy may allow to monitor response to therapy as well as detect minimal residual disease.