Cross platform analysis of methylation, miRNA and stem cell gene expression data in germ cell tumors highlights characteristic differences by tumor histology
BMC Cancer volume 15, Article number: 769 (2015)
Alterations in methylation patterns, miRNA expression, and stem cell protein expression occur in germ cell tumors (GCTs). Our goal is to integrate molecular data across platforms to identify molecular signatures in the three main histologic subtypes of Type I and Type II GCTs (yolk sac tumor (YST), germinoma, and teratoma).
We included 39 GCTs and 7 paired adjacent tissue samples in the current analysis. Molecular data available for analysis include DNA methylation data (Illumina GoldenGate Cancer Methylation Panel I), miRNA expression (NanoString nCounter miRNA platform), and stem cell factor expression (SABiosciences Human Embryonic Stem Cell Array). We evaluated the cross platform correlations of the data features using the Maximum Information Coefficient (MIC).
In analyses of individual datasets, differences were observed by tumor histology. Germinomas had higher expression of transcription factors maintaining stemness, while YSTs had higher expression of cytokines, endoderm and endothelial markers. We also observed differences in miRNA expression, with miR-371-5p, miR-122, miR-302a, miR-302d, and miR-373 showing elevated expression in one or more histologic subtypes. Using the MIC, we identified correlations across the data features, including six major hubs with higher expression in YST (LEFTY1, LEFTY2, miR302b, miR302a, miR 126, and miR 122) compared with other GCT.
While prognosis for GCTs is overall favorable, many patients experience resistance to chemotherapy, relapse and/or long term adverse health effects following treatment. Targeted therapies, based on integrated analyses of molecular tumor data such as that presented here, may provide a way to secure high cure rates while reducing unintended health consequences.
Germ cell tumors (GCTs) include germinomas, comprised of testicular seminomas and ovarian dysgerminomas, and nonseminomas, comprised of yolk sac tumors (YSTs), teratomas and embryonal carcinoma . While GCTs are heterogeneous, they are grouped together due to a presumed common stem cell of origin, the primordial germ cell (PGC). Oosterhuis and Looijenga have proposed classification of GCTs into five distinct entities based on cell of origin, histology, genomic imprinting status, age at and location of clinical presentation, and chromosomal constitution . Type I GCTs are those found predominantly in infants and young children, often manifesting in the first four years of life and always before puberty. Type II GCTs are most commonly found in the testis of adolescent males and young men following puberty, but are also found in the ovaries of adolescent and young adult women and the midline/brain of children and adolescents. Pathologic evidence confirms that GCTs are the neoplastic counterpart of the PGC , and several lines of evidence indicate that GCTs, including adult testicular GCT (TGCT), begin in utero . Thus, alterations in normal embryonic development are likely to be etiologically relevant to GCTs. Of particular interest are the processes the PGCs undergo during normal development, including segregation from the somatic cells, migration to the gonads, complete epigenetic reprogramming, reacquisition of pluripotency and sex determination .
Aberrant DNA methylation has been implicated in cancer etiology, and may be especially relevant in GCTs due to the extensive epigenetic reprogramming that occurs in the germ line and early embryo during normal development . Adult TGCTs have been studied most thoroughly in the context of DNA methylation, and thus a majority of our knowledge regarding methylation is limited to these tumors. Interestingly, methylation patterns in GCT differ by histologic subtype in both adults and children [6–16]. In general, methylation increases with tumor differentiation: the lowest levels of methylation occur in the embryonal carcinomas and the highest in the teratomas [6, 7, 10, 11, 13, 15–18]. Understanding methylation patterns in GCTs, overall and by histologic type, may identify the developmental stage at which the tumor arose. This knowledge in turn may identify the at-risk period when external exposures are most harmful.
MicroRNAs (miRNAs) are small endogenous noncoding RNAs that regulate gene function in a manner specific to cell type and developmental stage [19–23]. Differential miRNA expression is associated with human cancers [24–28], including GCTs in children and adults [29–35]. These studies have reported higher expression of miRNAs in the miR-371–73 and the miR-302 clusters and lower expression of let-7 in Type I and Type II GCTs compared to normal samples [29–37]. Alterations in the serum levels of the miR371–3 and miR-302/367 MiRNAs also show promise as a diagnostic and follow-up tool for TGCT patients , highlighting the potential translational impact of molecular evaluation.
Knowledge of stem cell biology is directly relevant to mechanisms of GCT tumor initiation, maintenance and metastasis, since reacquisition of pluripotency is a key step in early germ cell development . Typically, expression of stem cell markers (e.g., OCT3/4, STELLAR, NANOG, LIN28) is induced following demethylation of early stage germ cells [6, 17] and is turned off following entry to meiosis [40–42]. Expression of pluripotency markers past the appropriate developmental stage is a hypothesized explanation for tumorigenesis in germ cells . Notably, studies of adult TGCT have shown aberrant expression of stem cell markers in intratubular germ cell neoplasia (IGCNU), the precursor of TGCT, and in undifferentiated histologic subtypes of GCTs (seminomas and embryonal carcinomas) [43, 44]. Stem cell markers are also expressed in early germ cells in females [45–47] and have been detected in ovarian dysgerminomas . Marker expression past the appropriate developmental stage is correlated with genetic variation, including mutation in c-KIT  and its ligand (KITLG) , and DNA methylation . Given that pediatric GCTs likely originate from a germ cell at an earlier stage of development than adult TGCTs , stem cell marker expression may be particularly relevant in pediatric tumors.
As described above, previous studies have described variation in methylation, miRNA, and mRNA expression in Type I and Type II GCTs, including studies that have evaluated the interaction between miRNA and mRNA expression [29, 35]. To further explore relevant molecular interactions, we used an integrated approach to understand differences in promoter methylation, miRNA expression, stem cell gene expression, and genotype data by tumor characteristics in a series of GCTs. We evaluated correlations between data based on the assumption that these processes are linked and co-regulated (for example, epigenetic changes within promoter regions and expression of cognate miRNA species determine the level of mRNA). We also find differences in miRNA expression and stem cell gene expression by tumor histology.
Type I and Type II GCT samples from males and females were obtained from the Cooperative Human Tissue Network (Columbus, OH). Tumors were resected at initial diagnosis and snap frozen at −70 °C. Pathology reports were also provided. Data were available for tumor histology (YST, teratoma, germinoma, or mixed/other), tumor location (gonadal or extragonadal), sex, and age at diagnosis (< 10 years, ≥ 10 years). The age categories were chosen based on tumor histology. The majority of the tumors diagnosed between the ages of 4 and 10 were of similar histology to the Type I tumors while tumors diagnosed after age 10 included histologic subtypes typically included in the Type II category. Normal adjacent tissue was also available for seven of the tumors in our case series.
This analysis used existing data with no personal identifiers; therefore, the study was deemed exempt from review (category #4) by the Institutional Review Board of the University of Minnesota.
DNA and RNA extraction
Genomic DNA was isolated from GCT tissue and paired normal adjacent tissue (when available) using either the TRIzol® extraction method (Invitrogen Life Technologies, California) or a QIAamp DNA Mini Kit (Qiagen Sciences, Maryland) according to the manufacturer’s recommended protocol. DNA yield was quantified using 1 μl DNA on a NanoDrop™ spectrophotometer (Thermo Scientific, Maryland). Extracted DNA was stored at −80 °C until further analysis.
Total RNA was extracted from fresh frozen tissue using the TRIzol® extraction method (Invitrogen Life Technologies, California) according to the manufacturer's protocol. Following extraction, RNA was cleaned using the RNeasy Mini Kit (Qiagen, Maryland) according to the manufacturer's recommended protocol. RNA yield was then quantified using 2 μl on a NanoDrop™ spectrophotometer (Thermo Scientific, Maryland). Extracted RNA was stored at -80O C until further analysis.
DNA methylation was measured as previously described . Briefly, prior to methylation analysis, 1 μg genomic DNA was treated with sodium bisulfite to convert unmethylated cytosines to uracil using the EZ DNA Methylation Kit (Zymo Research, Orange, CA) according to manufacturer’s protocol. DNA methylation at 1505 CpG loci in 807 cancer-related genes was evaluated using the GoldenGate Cancer Methylation Panel I (Illumina, Inc.) in the University of Minnesota Genomics Center following the manufacturer’s protocol as described . Replicates were included, including four duplicates that were included on both arrays and five duplicates that were included within one array.
Methylation was calculated as the variable β, which is the ratio of the fluorescent signal from the methylated allele to the sum of the fluorescent signals of both methylated and unmethylated alleles . GenomeStudio software (Illumina, Inc) was used to calculate the average methylation values (β) from the ~30 replicate methylation measurements for each CpG locus. We used raw average β values without normalization. GenomeStudio software was also used to assess data quality for each CpG loci. We omitted all CpG loci where ≥ 25 % of the samples had a detection p-value > 0.05 (N = 16, 1 %). X-linked CpG loci (N = 84) were also removed, resulting in 1,405 loci for analysis.
Expression of 800 miRNAs was measured using the NanoString nCounter miRNA Expression Assay kit (NanoString Technologies, Seattle, WA). The nCounter detects total counts of miRNA through hybridization with fluorescently labeled bar coded probes to the miRNAs of interest followed by scanning and counting to quantify expression . Total RNA samples were analyzed following the manufacturer’s instructions for the Human v2 miRNA Expression Assay Kit (NanoString Technologies, Seattle, WA).
Stem cell factor expression
Real-time quantitative PCR gene expression profiling was performed for 84 pathway-specific genes using the human Embryonic Stem Cells RT2 Profiler PCR Array according to the manufacturer’s protocol (SABiosciences, Frederick, MD). Briefly, the RT2 First Strand Kit (SABioscience, Frederick, MD) was used to synthesize cDNA from 1 μg purified RNA. cDNA was obtained with a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA), and quantitative PCR (qPCR) was performed with a StepOnePlus Real-Time PCR system (Applied Biosystems, Foster City, CA). Expression of all genes was normalized to average expression of five endogenous housekeeping genes (E2M, HPRT1, RPL13A, GAPDH, and ACTB).
Genotype data were generated for four SNPs identified in GWAS of adult TGCT as previously described . Briefly, PCR amplification and sequencing were performed for SNPs in three genes: SPRY4 (rs4324715), BAK1 (rs210138) and DMRT1 (rs755383). The KITLG (rs4474514) SNP was detected using a made-to-order TaqMan® SNP Genotyping Assay from Applied Biosystems Inc (catalog# 4351379, assay# C_26154778_10). Primers and conditions for all assays are available upon request.
We used NanoStriDE for normalization and differential expression analysis of the miRNA data . Positive control normalization was conducted by creating a normalization factor for each sample using the 6 positive assay controls on the array and negative control normalization was conducted as an upper quantile approach as recommended by the manufacturer (NanoString Technologies, Seattle, WA). Discrete count data were compared across demographic and tumor characteristics using a negative binomial distribution as described by Anders and Huber . We included a Benjamini-Hochberg correction for multiple comparisons .
Stem cell factor expression analysis
Raw gene expression values were normalized to endogenous housekeeping genes prior to statistical analyses. Unsupervised hierarchical clustering was conducted using the matrix visualization and analysis platform Gene-E with the city block metric and average linkage . Fold change of gene expression was determined using the 2(−ΔΔCt) method, and compared YST (n = 9) to germinomas (n = 8).
Cross platform analysis
Patient age and sex, tumor location (ovary, testis, and extragonadal), and histology (normal adjacent, teratoma, dysgerminoma, YST and mixed) data were combined with molecular data across a common set of 40 samples into a single two-dimensional matrix. Categorical phenotype and genotype values were arranged in a logical manner (e.g., low severity to high severity for tumor histology) and assigned a numerical code that could be used for correlation analysis. These data were combined with the β values from the methylation analysis, the miRNA counts, and the normalized ΔCt values from the SABiosciences Embryonic Stem Cell array. Additionally, genotypes for the four molecular markers listed above were encoded based on the Kimura matrix  idea that transitions are more likely than transversions and hence grouped together where possible.: −4.0 = t/t; −3.0 = c/t; −2.0 = c/c; −0.5 = g/t; −0.2 = c/g; 0.2 = a/t; 0.5 = a/c; 2.0 = a/a; 3.0 = a/g; 4.0 = g/g.
Two types of correlations were explored using the vector of values for each matrix row (i.e., molecular probe or phenotypic label) against the corresponding vector of values for each of the other rows: standard linear correlations using the Pearson correlation coefficient and more complex nonlinear correlations using the Maximal Information Coefficient (MIC) . MIC analysis was performed using the R implementation of the MINE software package (http://www.exploredata.net/). This resulted in two values for correlation (linear, nonlinear) for every pairwise combination of molecular probes and/or phenotypic variables. This very large result matrix was filtered to retain only those pairs that represented a comparison across platforms (or between a phenotypic variable and a platform) and whose Pearson correlation or MIC values exceeded a threshold of 0.75. This threshold was chosen as it provided a compromise between very dense graph connectivity at lower thresholds and sparse connectivity at high cutoffs.
The resulting pairwise correlations were visualized as a network using Gephi software (https://gephi.org/). In order to visualize the correlation networks, the molecular probes and phenotypic variables served as nodes in the network graph, and an edge was drawn between any nodes that had either a linear or nonlinear correlation that exceeded 0.75. The edge was labeled with the larger of either the Pearson R-value or the MIC. These data were loaded into Gephi, and the Force Atlas 2 layout was executed, running until the nodes were far separated in apparent equilibrium, and then the Fruchterman Reingold layout was selected. All pairwise correlations were grouped by single linkage clustering to create networks. Network hubs with four or more neighbors were identified.
The hub nodes and their nearest neighbors were concatenated into a list, and the standard gene symbols extracted. These gene symbols were uploaded to Ingenuity Pathway Analysis software (IPA, http://www.ingenuity.com/products/ipa), and investigated for upstream activators. To guard against the possibility that functional categories were identified by chance due to the a priori bias in the initial set of profiling molecules, we randomly selected three additional sets of the same size from the same initial profiling molecules, and submitted these random gene lists to IPA.
The entire analysis process was repeated using subsets of the data. In the entire data set of 40 samples with complete genomic data, it was noted via manual analysis of scatterplots of primary hub genes that differences in the 8 YSTs relative to the other samples dominated the signal. Following the assumption that the YSTs were overshadowing any signal from other histology types, we tested a subset consisting of the 32 non-YST samples.
Characteristics of the study samples
Tumor specimens from 38 cases of Type I and Type II GCT ranging in age from 0 to 21 years were included in this analysis, including 9 YSTs, 18 teratomas, 7 dysgerminomas, and 4 with mixed histology (Table 1). In addition, one dysgerminoma from a 45 year old woman (Type II GCT) was included because it had similar methylation, miRNA and stem cell factor expression values as the adolescent dysgerminoma samples. The YSTs were evenly distributed among boys and girls while the majority of cases with a germinoma or teratoma were female. Information on race/ethnicity was not available for the cases.
The four cases with mixed histology all had a teratoma component of the tumor. The three tumors in male cases also had YST. One mixed tumor in an adolescent male and the tumor in an adolescent female also had components of embryonal carcinoma and choriocarcinoma. Normal adjacent DNA was available for seven cases. The seven normal adjacent tissue samples included four ovary or fallopian tube tissues, one testis tissue, one adjacent lymph node and one thymus tissue. The testis sample was from a one year old case and is unlikely to contain IGCNU. Results from methylation and genotyping data have been previously published [16, 52]. and are therefore not included in this report.
In a comparison across all histologic subtypes, we observed significant differences in miRNA expression for five miRNA species (Fig. 1). Of these, all had low expression in normal adjacent tissue. One miRNA was elevated only in dysgerminomas (hsa-miR-371–5p). The other four miRNAs were elevated in YST with varying expression differences in other histologic subtypes. We did not observe differential expression for any miRNA species when we compared gonadal vs. extragonadal tumors, tumors from males vs. females, or tumors from children diagnosed prior to age 10 years vs. those diagnosed at or after 10 years of age (data not shown).
Stem cell factor expression
Unsupervised hierarchical clustering based on tumor histology highlights differences in the expression of stem cell genes by tumor histology (Fig. 2a), with the YST and dysgerminoma separating into distinct clusters, while the teratomas (mature and immature) and normal adjacent tissue clustered together. Mixed GCTs were interspersed throughout the groups. To better understand the expression differences in YST and dysgerminoma, we compared relative expression of each gene in these two groups. Of the 84 genes included in the array, 40 had a statistically significant up- or down-regulation in YST compared with dysgerminoma (p < 0.05). When these genes were categorized by function in pathways regulating initiation and maintenance of cellular stemness, we noted distinct patterns based on tumor histology (Fig. 2b). For example, transcription factors related to stemness had higher expression in dysgerminomas, while endoderm, trophoblast, and mesoderm markers had higher expression in YST.
We observed 760 cross-platform correlations that exceeding our threshold of 0.75 for either the Pearson correlation coefficient or the MIC. When these data correlations were visualized, six network hubs were identified (Fig. 3), including the miRNAs, miR-122, −126, −302a and -302b and the stem cell genes, LEFTY1 and LEFTY2. These changes were influenced largely by differences in the expression in YST vs. tumors of other histologic subtypes. In order to detect differences that may have been overshadowed by the strong influence of the YSTs, we also repeated the analysis including only the dysgerminomas and teratomas. In this analysis, the protein T (brachyury) was identified as a network hub linking a large number of molecular features (data not shown).
Using IPA, we identified many predicted upstream activators that were enriched in the list of genes with large cross platform correlations (N = 45 genes). After comparison with the p-values from the randomly selected gene lists, there were three molecules (TP73, decitabine, and tretinoin) with statistically significant p-values after Benjamini-Hochberg correction, suggesting that these molecules are promising upstream activators.
In this analysis, we used a novel method to integrate molecular data across platforms to gain biological insight into the function of GCTs. These analyses highlighted several network hubs, including miRNA clusters and stem cell genes that distinguish YST from normal germ cell samples and other GCTs. Importantly, our approach confirmed several previously reported alterations in embryonic stem cell specific miRNAs. Finally, our analysis of stem cell factor expression highlights the altered expression of multiple stem cell genes in GCTs, with distinct patterns in YST and dysgerminoma. Collectively, these findings suggest that ectopic, aberrant expression of stem cell genes may underlie the unusual and defining capacity of self-renewal in the face of wide differentiation into cells with characteristics of tissues derived from any of the three germ layers observed in GCT, and it is possible that higher levels of stem cell gene expression correlate with tumor progression and prognosis of GCT and other tumors.
Our findings are consistent with previous studies of Type I and Type II GCTs that have identified an overexpression of the miR-371–73 and miR-302 clusters in GCTs compared with normal samples [29–35]. The miRNA-302 and miRNA-371–373 clusters are highly plausible candidate miRNAs in GCTs given their roles as regulators of embryonic stem cell pluripotency markers [59, 60] which has been discussed extensively in previous miRNA expression studies of GCTs [33, 35]. Pinpointing the relevant targets of miRNA can be daunting given the large number of target proteins for each miRNA. According to the online database for miRNA target prediction and functional annotations, miRDB [61, 62] the number of predicted targets for miRNAs in the miR-371–373 and miR-302 families range from 469 to 529; interestingly, LATS2 is one of the highest ranking targets on the list for all miR-302 family members as well as miR-372 and miR-373 (Target Score > 98) . Functional studies have demonstrated that miRNA-372 and miRNA-373 act as oncogenes in TGCT through interactions with the p53 pathway, in particular through regulation of LATS2 . Bioinformatic algorithms indicated that miRNA expression in these clusters is associated with downregulation of mRNA expression in other pathways with biological relevance to GCT [33, 35]. Data also suggest that the miRNA 302 family can be used to reprogram cancer cells into pluripotent cells with an ES cell-like phenotype . Collectively, these data suggest that these miRNA clusters are in large part responsible for regulating the stem cell phenotype of GCTs.
Four of the six hubs identified in the cross-platform analysis were stem cell related, and highlight the importance of the transforming growth factor (TGF)-β superfamily members in GCT development. Specifically, the miR-302 cluster regulates the Nodal inhibitors LEFTY1 and LEFTY2 . This interaction plays an important role in germ layer specification by promoting the formation of the mesendodermal lineage while suppressing neuroectorderm formation . Knockdown of either protein in mice results in altered differentiation, with Lefty1 knockdown leading to increased differentiation potential and Lefty2 knockdown leading to increased immature neuroepithelium . This pathway also plays other important roles in germ cell development, including regulating meiosis in the male germ cells [67, 68]. Finally, the nodal signaling pathway regulates the bone morphogenic protein pluripotency pathway , which was previously shown to play an important role in the development of GCT . The importance of miR-122 and miR-126 in GCT is less clear and will require further study. Interestingly, miR-122 was also overexpressed in YST compared with germinoma and normal gonad tissue in a previous study of miRNA in GCT . Notably, none of the methylation differences that were so striking in our comparison of YST to other GCTs  were identified as being hubs in the cross-platform analysis, suggesting that these changes may be consequence of the altered expression of stem cell genes and miRNAs rather than drivers of the oncogenic process. Given that the associations between the four SNPs evaluated did not differ by tumor histology [52, 71, 72], it was not as surprising that none of these were identified as hubs in the comparison of tumors by histology.
We chose to use platforms that were highly enriched for cancer and stem cell genes due to the higher information yield when compared to an unbiased search; however, this limited our ability to conduct an unbiased search for over-represented gene categories in the list of features with high cross-platform correlations. We were able to evaluate upstream regulators of these highly correlated features, and we identified three highly enriched activators (TP73, Tretinoin, and Decitabine) in the set of hub-connected genes in the entire dataset. Given the known importance of retinoic acid in germ cell development , it is not surprising that Tretinoin, a topical retinoid, would be identified as a potential regulator in GCTs. The DNA hypomethylating agent Decitabine (5-aza-2’-deoxycytidine or 5-aza) is also intriguing as a potential therapy given relevant data in the literature. Decitabine is effective for treatment of hematologic malignancies [74, 75] and has also been evaluated in studies of solid tumors with varying success rates . Preclinical data suggesting that hypomethylating agents may re-sensitize cells to platinum based chemotherapy [77–79] are of particular relevance for GCT; however, early phase clinical trials in epithelial ovarian cancer have provided mixed results of the combination of decitabine and carboplatinum [80–82]. In a previous study of embryonal carcinoma, high expression of the pluripotency-associated DNA methyltransferase 3B (DNMT3B) was associated with sensitivity to 5-Aza [77, 83]. The elevated expression of DNMT3B in YST observed here and in a previous DNA methylation study  suggests that this may also be a relevant alternative therapy for YST, which often have poorer outcomes than other histologic subtypes .
Most cancers are heterogeneous and multifactorial . This complexity stems from molecular events on the level of gene integrity, epigenetic modification and transcription, stability of mRNA transcripts (e.g., by miRNAs), translation, protein activation (e.g., by phosphorylation), and cellular interactions, including within the tumor microenvironment. Because of this complexity, it is important to evaluate the joint effects of these alterations. There is not one clear method for evaluating these complex interactions on large datasets. In this analysis, we chose to utilize the MIC as a tool to evaluate the correlations in the dataset in addition to simple linear Pearson correlation. MIC is able to detect linear relationships and correlations among variables that are not strictly linear (such as quadratic or oscillatory associations). Care should be taken when using the MIC with small datasets (e.g., N < 30), as it may identify spurious correlations in this extreme . Also, despite a recent debate  as to which information-theory based measure provides the highest statistical power (maximal information coefficient or a related measure, mutual information), both measures clearly identify nonlinear trends that are missed by Pearson correlation. Additionally, our analysis included a small number of samples which could have limited our power to detect relevant associations. Finally, our analysis focused on a selected list of genes with a priori significance in cancer and stem cells. It is possible that a genome-wide, agnostic approach would identify additional relevant characteristics of GCTs in additional pathways.
Our analysis suggests that the stem cell phenotype of GCTs is a defining characteristic of GCTs, especially in YSTs. While prognosis for GCTs is overall favorable, subgroups of patients experience resistance to chemotherapy and/or high rates of relapse [88–90]. In addition, the nonspecific and highly cytotoxic chemotherapy used can cause many adverse health effects including cardiovascular disease, hearing loss, and second cancers [91–95]. Targeted therapies, based on integrated analyses of molecular tumor data such as that presented here, may provide a way to improve cure rates in the subgroup of patients who fail to respond to current therapies and may also provide an opportunity to reduce the unintended health consequences associated with current chemotherapeutic agents.
Germ cell tumor
Yolk sac tumor
Maximum information coefficient
Primordial germ cell
testicular germ cell tumor
Ingenuity Pathway Analysis
Transforming growth factor
DNA methyltransferase 3B
Cushing B, Perlman EJ, Marina NM, Castleberry RP. Germ Cell Tumors. In: Pizzo P, Poplack D, editors. Principles and practice of pediatric oncology. 5th ed. Philadelphia: Lippincott, Williams and Wilkins; 2006.
Oosterhuis JW, Looijenga LH. Testicular germ-cell tumours in a broader perspective. Nat Rev Cancer. 2005;5(3):210–22.
Moller H. Decreased testicular cancer risk in men born in wartime. J Natl Cancer Inst. 1989;81(21):1668–9.
Wylie C. Germ cells. Cell. 1999;96(2):165–74.
Smallwood SA, Kelsey G. De novo DNA methylation: a germ cell perspective. Trends Genet. 2012;28(1):33–42.
Gillis AJ, Stoop H, Biermann K, van Gurp RJ, Swartzman E, Cribbes S, et al. Expression and interdependencies of pluripotency factors LIN28, OCT3/4, NANOG and SOX2 in human testicular germ cells and tumours of the testis. Int J Androl. 2011;34(4 Pt 2):e160–174.
Lind GE, Skotheim RI, Lothe RA. The epigenome of testicular germ cell tumors. Apmis. 2007;115(10):1147–60.
Honecker F, Oosterhuis JW, Mayer F, Hartmann JT, Bokemeyer C, Looijenga LH. New insights into the pathology and molecular biology of human germ cell tumors. World J Urol. 2004;22(1):15–24.
McIntyre A, Summersgill B, Lu YJ, Missiaglia E, Kitazawa S, Oosterhuis JW, et al. Genomic copy number and expression patterns in testicular germ cell tumours. Br J Cancer. 2007;97(12):1707–12.
Netto GJ, Nakai Y, Nakayama M, Jadallah S, Toubaji A, Nonomura N, et al. Global DNA hypomethylation in intratubular germ cell neoplasia and seminoma, but not in nonseminomatous male germ cell tumors. Mod Pathol. 2008;21(11):1337–44.
Smiraglia DJ, Szymanska J, Kraggerud SM, Lothe RA, Peltomaki P, Plass C. Distinct epigenetic phenotypes in seminomatous and nonseminomatous testicular germ cell tumors. Oncogene. 2002;21(24):3909–16.
Koul S, Houldsworth J, Mansukhani MM, Donadio A, McKiernan JM, Reuter VE, et al. Characteristic promoter hypermethylation signatures in male germ cell tumors. Mol Cancer. 2002;1:8.
Looijenga LH, Gillis AJ, Stoop H, Biermann K, Oosterhuis JW. Dissecting the molecular pathways of (testicular) germ cell tumour pathogenesis; from initiation to treatment-resistance. Int J Androl. 2011;34(4 Pt 2):e234–251.
Kremenskoy M, Kremenska Y, Ohgane J, Hattori N, Tanaka S, Hashizume K, et al. Genome-wide analysis of DNA methylation status of CpG islands in embryoid bodies, teratomas, and fetuses. Biochem Biophys Res Commun. 2003;311(4):884–90.
Jeyapalan JN, Noor DA, Lee SH, Tan CL, Appleby VA, Kilday JP, et al. Methylator phenotype of malignant germ cell tumours in children identifies strong candidates for chemotherapy resistance. Br J Cancer. 2011;105(4):575–85.
Amatruda JF, Ross JA, Christensen B, Fustino NJ, Chen KS, Hooten AJ, et al. DNA methylation analysis reveals distinct methylation signatures in pediatric germ cell tumors. BMC Cancer. 2013;13:313.
Wermann H, Stoop H, Gillis AJ, Honecker F, van Gurp RJ, Ammerpohl O, et al. Global DNA methylation in fetal human germ cells and germ cell tumours: association with differentiation and cisplatin resistance. J Pathol. 2010;221(4):433–42.
Godmann M, Lambrot R, Kimmins S. The dynamic epigenetic program in male germ cells: Its role in spermatogenesis, testis cancer, and its response to the environment. Microsc Res Tech. 2009;72(8):603–19.
Wyman SK, Knouf EC, Parkin RK, Fritz BR, Lin DW, Dennis LM, et al. Post-transcriptional generation of miRNA variants by multiple nucleotidyl transferases contributes to miRNA transcriptome complexity. Genome Res. 2011;21(9):1450–61.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
Holley CL, Topkara VK. An introduction to small non-coding RNAs: miRNA and snoRNA. Cardiovasc Drugs Ther. 2011;25(2):151–9.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Djuranovic S, Nahvi A, Green R. A parsimonious model for gene regulation by miRNAs. Science. 2011;331(6017):550–3.
Lynam-Lennon N, Maher SG, Reynolds JV. The roles of microRNA in cancer and apoptosis. Biol Rev Camb Philos Soc. 2009;84(1):55–71.
Looijenga LH, Gillis AJ, Stoop H, Hersmus R, Oosterhuis JW. Relevance of microRNAs in normal and malignant development, including human testicular germ cell tumours. Int J Androl. 2007;30(4):304–14. discussion 314–305.
Li X, Chen J, Hu X, Huang Y, Li Z, Zhou L, et al. Comparative mRNA and microRNA expression profiling of three genitourinary cancers reveals common hallmarks and cancer-specific molecular events. PLoS One. 2011;6(7):e22570.
Novotny GW, Sonne SB, Nielsen JE, Jonstrup SP, Hansen MA, Skakkebaek NE, et al. Translational repression of E2F1 mRNA in carcinoma in situ and normal testis correlates with expression of the miR-17-92 cluster. Cell Death Differ. 2007;14(4):879–82.
Tellez CS, Juri DE, Do K, Bernauer AM, Thomas CL, Damiani LA, et al. EMT and stem cell-like properties associated with miR-205 and miR-200 epigenetic silencing are early manifestations during carcinogen-induced transformation of human lung epithelial cells. Cancer Res. 2011;71(8):3087–97.
Gillis AJ, Stoop HJ, Hersmus R, Oosterhuis JW, Sun Y, Chen C, et al. High-throughput microRNAome analysis in human germ cell tumours. J Pathol. 2007;213(3):319–28.
Krausz C, Looijenga LH. Genetic aspects of testicular germ cell tumors. Cell Cycle. 2008;7(22):3519–24.
van de Geijn GJ, Hersmus R, Looijenga LH. Recent developments in testicular germ cell tumor research. Birth Defects Res C Embryo Today. 2009;87(1):96–113.
Voorhoeve PM, le Sage C, Schrier M, Gillis AJ, Stoop H, Nagel R, et al. A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors. Cell. 2006;124(6):1169–81.
Murray MJ, Saini HK, van Dongen S, Palmer RD, Muralidhar B, Pett MR, et al. The two most common histological subtypes of malignant germ cell tumour are distinguished by global microRNA profiles, associated with differential transcription factor expression. Mol Cancer. 2010;9:290.
Fustino N, Rakheja D, Ateek CS, Neumann JC, Amatruda JF. Bone morphogenetic protein signalling activity distinguishes histological subsets of paediatric germ cell tumours. Int J Androl. 2011;34(4 Pt 2):e218–233.
Palmer RD, Murray MJ, Saini HK, van Dongen S, Abreu-Goodger C, Muralidhar B, et al. Malignant germ cell tumors display common microRNA profiles resulting in global changes in expression of messenger RNA targets. Cancer Res. 2010;70(7):2911–23.
Viswanathan SR, Powers JT, Einhorn W, Hoshida Y, Ng TL, Toffanin S, et al. Lin28 promotes transformation and is associated with advanced human malignancies. Nat Genet. 2009;41(7):843–8.
West JA, Viswanathan SR, Yabuuchi A, Cunniff K, Takeuchi A, Park IH, et al. A role for Lin28 in primordial germ-cell development and germ-cell malignancy. Nature. 2009;460(7257):909–13.
Gillis AJ, Rijlaarsdam MA, Eini R, Dorssers LC, Biermann K, Murray MJ, et al. Targeted serum miRNA (TSmiR) test for diagnosis and follow-up of (testicular) germ cell cancer patients: a proof of principle. Mol Oncol. 2013;7(6):1083–92.
Yabuta Y, Kurimoto K, Ohinata Y, Seki Y, Saitou M. Gene expression dynamics during germline specification in mice identified by quantitative single-cell gene expression profiling. Biol Reprod. 2006;75(5):705–16.
Pesce M, Wang X, Wolgemuth DJ, Scholer H. Differential expression of the Oct-4 transcription factor during mouse germ cell differentiation. Mech Dev. 1998;71(1–2):89–98.
Western PS, van den Bergen JA, Miles DC, Sinclair AH. Male fetal germ cell differentiation involves complex repression of the regulatory network controlling pluripotency. FASEB J. 2010;24(8):3026–35.
Clark AT, Rodriguez RT, Bodnar MS, Abeyta MJ, Cedars MI, Turek PJ, et al. Human STELLAR, NANOG, and GDF3 genes are expressed in pluripotent cells and map to chromosome 12p13, a hotspot for teratocarcinoma. Stem Cells. 2004;22(2):169–79.
Hart AH, Hartley L, Parker K, Ibrahim M, Looijenga LH, Pauchnik M, et al. The pluripotency homeobox gene NANOG is expressed in human germ cell tumors. Cancer. 2005;104(10):2092–8.
Rijlaarsdam MA, van Herk HA, Gillis AJ, Stoop H, Jenster G, Martens J, et al. Specific detection of OCT3/4 isoform A/B/B1 expression in solid (germ cell) tumours and cell lines: confirmation of OCT3/4 specificity for germ cell tumours. Br J Cancer. 2011;105(6):854–63.
Zuccotti M, Merico V, Sacchi L, Bellone M, Brink TC, Stefanelli M, et al. Oct-4 regulates the expression of Stella and Foxj2 at the Nanog locus: implications for the developmental competence of mouse oocytes. Hum Reprod. 2009;24(9):2225–37.
Yeom YI, Fuhrmann G, Ovitt CE, Brehm A, Ohbo K, Gross M, et al. Germline regulatory element of Oct-4 specific for the totipotent cycle of embryonal cells. Development. 1996;122(3):881–94.
Rosner MH, Vigano MA, Ozato K, Timmons PM, Poirier F, Rigby PW, et al. A POU-domain transcription factor in early stem cells and germ cells of the mammalian embryo. Nature. 1990;345(6277):686–92.
Hoei-Hansen CE, Kraggerud SM, Abeler VM, Kaern J, Rajpert-De Meyts E, Lothe RA. Ovarian dysgerminomas are characterised by frequent KIT mutations and abundant expression of pluripotency markers. Mol Cancer. 2007;6:12.
Stoop H, Honecker F, van de Geijn GJ, Gillis AJ, Cools MC, de Boer M, et al. Stem cell factor as a novel diagnostic marker for early malignant germ cells. J Pathol. 2008;216(1):43–54.
Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, et al. High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006;16(3):383–93.
Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008;26(3):317–25.
Poynter JN, Hooten AJ, Frazier AL, Ross JA. Associations between variants in KITLG, SPRY4, BAK1, and DMRT1 and pediatric germ cell tumors. Genes Chromosomes Cancer. 2012;51(3):266–71.
Brumbaugh CD, Kim HJ, Giovacchini M, Pourmand N. NanoStriDE: Normalization and Differential Expression Analysis of NanoString nCounter Data. BMC Bioinformatics. 2011;12(1):479.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20.
Reshef DN, Reshef YA, Finucane HK, Grossman SR, McVean G, Turnbaugh PJ, et al. Detecting novel associations in large data sets. Science. 2011;334(6062):1518–24.
Barroso-delJesus A, Romero-Lopez C, Lucena-Aguilar G, Melen GJ, Sanchez L, Ligero G, et al. Embryonic stem cell-specific miR302-367 cluster: human gene structure and functional characterization of its core promoter. Mol Cell Biol. 2008;28(21):6609–19.
Suh MR, Lee Y, Kim JY, Kim SK, Moon SH, Lee JY, et al. Human embryonic stem cells express a unique set of microRNAs. Dev Biol. 2004;270(2):488–98.
Wang X, El Naqa IM. Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics. 2008;24(3):325–32.
Wang X. miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA. 2008;14(6):1012–7.
Lin SL, Chang DC, Chang-Lin S, Lin CH, Wu DT, Chen DT, et al. Mir-302 reprograms human skin cancer cells into a pluripotent ES-cell-like state. RNA. 2008;14(10):2115–24.
Barroso-delJesus A, Lucena-Aguilar G, Sanchez L, Ligero G, Gutierrez-Aranda I, Menendez P. The Nodal inhibitor Lefty is negatively modulated by the microRNA miR-302 in human embryonic stem cells. FASEB J. 2011;25(5):1497–508.
Rosa A, Spagnoli FM, Brivanlou AH. The miR-430/427/302 family controls mesendodermal fate specification via species-specific target selection. Dev Cell. 2009;16(4):517–27.
Kim DK, Cha Y, Ahn HJ, Kim G, Park KS. Lefty1 and lefty2 control the balance between self-renewal and pluripotent differentiation of mouse embryonic stem cells. Stem Cells Dev. 2014;23(5):457–66.
Souquet B, Tourpin S, Messiaen S, Moison D, Habert R, Livera G. Nodal signaling regulates the entry into meiosis in fetal germ cells. Endocrinology. 2012;153(5):2466–73.
Wu Q, Kanata K, Saba R, Deng CX, Hamada H, Saga Y. Nodal/activin signaling promotes male germ cell fate and suppresses female programming in somatic cells. Development. 2013;140(2):291–300.
Galvin KE, Travis ED, Yee D, Magnuson T, Vivian JL. Nodal signaling regulates the bone morphogenic protein pluripotency pathway in mouse embryonic stem cells. J Biol Chem. 2010;285(26):19747–56.
Neumann JC, Chandler GL, Damoulis VA, Fustino NJ, Lillard K, Looijenga L, et al. Mutation in the type IB bone morphogenetic protein receptor Alk6b impairs germ-cell differentiation and causes germ-cell tumors in zebrafish. Proc Natl Acad Sci U S A. 2011;108(32):13153–8.
Kanetsky PA, Mitra N, Vardhanabhuti S, Li M, Vaughn DJ, Letrero R, et al. Common variation in KITLG and at 5q31.3 predisposes to testicular germ cell cancer. Nat Genet. 2009;41(7):811–5.
Rapley EA, Turnbull C, Al Olama AA, Dermitzakis ET, Linger R, Huddart RA, et al. A genome-wide association study of testicular germ cell tumor. Nat Genet. 2009;41(7):807–10.
Bowles J, Knight D, Smith C, Wilhelm D, Richman J, Mamiya S, et al. Retinoid signaling determines germ cell fate in mice. Science. 2006;312(5773):596–600.
Issa JP, Garcia-Manero G, Giles FJ, Mannari R, Thomas D, Faderl S, et al. Phase 1 study of low-dose prolonged exposure schedules of the hypomethylating agent 5-aza-2'-deoxycytidine (decitabine) in hematopoietic malignancies. Blood. 2004;103(5):1635–40.
Joeckel TE, Lubbert M. Clinical results with the DNA hypomethylating agent 5-aza-2'-deoxycytidine (decitabine) in patients with myelodysplastic syndromes: an update. Semin Hematol. 2012;49(4):330–41.
Nie J, Liu L, Li X, Han W. Decitabine, a new star in epigenetic therapy: the clinical application and biological mechanism in solid tumors. Cancer Lett. 2014;354(1):12–20.
Beyrouthy MJ, Garner KM, Hever MP, Freemantle SJ, Eastman A, Dmitrovsky E, et al. High DNA methyltransferase 3B expression mediates 5-aza-deoxycytidine hypersensitivity in testicular germ cell tumors. Cancer Res. 2009;69(24):9360–6.
Balch C, Yan P, Craft T, Young S, Skalnik DG, Huang TH, et al. Antimitogenic and chemosensitizing effects of the methylation inhibitor zebularine in ovarian cancer. Mol Cancer Ther. 2005;4(10):1505–14.
Zeller C, Dai W, Steele NL, Siddiq A, Walley AJ, Wilhelm-Benartzi CS, et al. Candidate DNA methylation drivers of acquired cisplatin resistance in ovarian cancer identified by methylome and expression profiling. Oncogene. 2012;31(42):4567–76.
Fang F, Balch C, Schilder J, Breen T, Zhang S, Shen C, et al. A phase 1 and pharmacodynamic study of decitabine in combination with carboplatin in patients with recurrent, platinum-resistant, epithelial ovarian cancer. Cancer. 2010;116(17):4043–53.
Glasspool RM, Brown R, Gore ME, Rustin GJ, McNeish IA, Wilson RH, et al. A randomised, phase II trial of the DNA-hypomethylating agent 5-aza-2'-deoxycytidine (decitabine) in combination with carboplatin vs carboplatin alone in patients with recurrent, partially platinum-sensitive ovarian cancer. Br J Cancer. 2014;110(8):1923–9.
Matei D, Fang F, Shen C, Schilder J, Arnold A, Zeng Y, et al. Epigenetic resensitization to platinum in ovarian cancer. Cancer Res. 2012;72(9):2197–205.
Biswal BK, Beyrouthy MJ, Hever-Jardine MP, Armstrong D, Tomlinson CR, Christensen BC, et al. Acute hypersensitivity of pluripotent testicular cancer-derived embryonal carcinoma to low-dose 5-aza deoxycytidine is associated with global DNA Damage-associated p53 activation, anti-pluripotency and DNA demethylation. PLoS One. 2012;7(12):e53003.
De Backer A, Madern GC, Pieters R, Haentjens P, Hakvoort-Cammel FG, Oosterhuis JW, et al. Influence of tumor site and histology on long-term survival in 193 children with extracranial germ cell tumors. Eur J Pediatr Surg. 2008;18(1):1–6.
Potter JD. Morphogens, morphostats, microarchitecture and malignancy. Nat Rev Cancer. 2007;7(6):464–74.
de Siqueira SS, Takahashi DY, Nakata A, Fujita A. A comparative study of statistical methods used to identify dependencies between gene expression signals. Brief Bioinform. 2013;15(6):906–18.
Kinney JB, Atwal GS. Equitability, mutual information, and the maximal information coefficient. Proc Natl Acad Sci U S A. 2014;111(9):3354–9.
Baranzelli MC, Kramar A, Bouffet E, Quintana E, Rubie H, Edan C, et al. Prognostic factors in children with localized malignant nonseminomatous germ cell tumors. J Clin Oncol. 1999;17(4):1212.
International Germ Cell Consensus Classification: a prognostic factor-based staging system for metastatic germ cell cancers. International Germ Cell Cancer Collaborative Group. J Clin Oncol. 1997; 15(2):594–603
Koychev D, Oechsle K, Bokemeyer C, Honecker F. Treatment of patients with relapsed and/or cisplatin-refractory metastatic germ cell tumours: an update. Int J Androl. 2011;34(4 Pt 2):e266–273.
Travis LB, Fossa SD, Schonfeld SJ, McMaster ML, Lynch CF, Storm H, et al. Second cancers among 40,576 testicular cancer patients: focus on long-term survivors. J Natl Cancer Inst. 2005;97(18):1354–65.
Travis LB, Beard C, Allan JM, Dahl AA, Feldman DR, Oldenburg J, et al. Testicular cancer survivorship: research strategies and recommendations. J Natl Cancer Inst. 2010;102(15):1114–30.
Meinardi MT, Gietema JA, van der Graaf WT, van Veldhuisen DJ, Runne MA, Sluiter WJ, et al. Cardiovascular morbidity in long-term survivors of metastatic testicular cancer. J Clin Oncol. 2000;18(8):1725–32.
van den Belt-Dusebout AW, de Wit R, Gietema JA, Horenblas S, Louwman MW, Ribot JG, et al. Treatment-specific risks of second malignancies and cardiovascular disease in 5-year survivors of testicular cancer. J Clin Oncol. 2007;25(28):4370–8.
Fligor BJ, Krasker JD, Villaluna D, Krailo M, Frazier AL. “Accelerated ear-age”: a new measure of chemotherapy-induced ototoxicity. Pediatr Blood Cancer. 2012;59(5):947–9.
The authors would like to acknowledge the University of Minnesota Genomics Center for performing the DNA methylation analysis and the miRNA analysis. Supported by grants from the National Institutes of Health (R03 CA141481 to J.N.P., K05 CA157439 and T32 CA099936 to J.A.R.) and the Children’s Cancer Research Fund, Minneapolis, MN
The authors declare that they have no competing interests.
JNP conceived of the study and participated in its design and coordination, analyzed data and wrote the manuscript. JRBMB analyzed miRNA data and contributed to manuscript writing. KATS conducted cross platform data analysis. AJH performed methylation, miRNA, and genotyping analysis. CL performed stem cell factor analysis. JAR conceived of the study and participated in its design and coordination. JT conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.
About this article
Cite this article
Poynter, J.N., Bestrashniy, J.R.B.M., Silverstein, K.A.T. et al. Cross platform analysis of methylation, miRNA and stem cell gene expression data in germ cell tumors highlights characteristic differences by tumor histology. BMC Cancer 15, 769 (2015). https://doi.org/10.1186/s12885-015-1796-6