Inhibition of the galactosyltransferase C1GALT1 reduces osteosarcoma cell proliferation by interfering with ERK signaling and cell cycle progression.

Novel therapeutic strategies are urgently required for osteosarcoma, given the early age at onset and persistently high mortality rate. Modern transcriptomics techniques can identify differentially expressed genes (DEGs) that may serve as biomarkers and therapeutic targets, so we screened for DEGs in osteosarcoma. We found that osteosarcoma cases could be divided into fair and poor survival groups based on gene expression profiles. Among the genes upregulated in the poor survival group, siRNA-mediated knockdown of the glycosylation-related gene C1GALT1 suppressed osteosarcoma cell proliferation in culture. Gene expression, phosphorylation, and glycome array analyses also demonstrated that C1GALT1 is required to maintain ERK signaling and cell cycle progression. Moreover, the C1GALT1 inhibitor itraconazole suppressed osteosarcoma cell proliferation in culture, while doxycycline-induced shRNA-mediated knockdown reduced xenograft osteosarcoma growth in mice. Elevated C1GALT1 expression is a potential early predictor of poor prognosis, while pharmacological inhibition may be a feasible treatment strategy for osteosarcoma.

Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Myeloid/natural killer (NK) cell precursor acute leukemia as a distinct leukemia type.

Myeloid/natural killer (NK) cell precursor acute leukemia (MNKPL) has been described on the basis of its unique immunophenotype and clinical phenotype. However, there is no consensus on the characteristics for identifying this disease type because of its rarity and lack of defined distinctive molecular characteristics. In this study, multiomics analysis revealed that MNKPL is distinct from acute myeloid leukemia, T cell acute lymphoblastic leukemia, and mixed-phenotype acute leukemia (MPAL), and NOTCH1 and RUNX3 activation and BCL11B down-regulation are hallmarks of MNKPL. Although NK cells have been classically considered to be lymphoid lineage-derived, the results of our single-cell analysis using MNKPL cells suggest that NK cells and myeloid cells share common progenitor cells. Treatment outcomes for MNKPL are unsatisfactory, even when hematopoietic cell transplantation is performed. Multiomics analysis and in vitro drug sensitivity assays revealed increased sensitivity to l-asparaginase and reduced levels of asparagine synthetase (ASNS), supporting the clinically observed effectiveness of l-asparaginase.

Heterozygous loss-of-function SMC3 variants are associated with variable and incompletely penetrant growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 13 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated a milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, some instead having intriguing symptomatologies with rational biological links to SMC3 including bone marrow failure, acute myeloid leukemia, and Coats retinal vasculopathy. Analyses of transcriptomic and epigenetic data suggest that SMC3 pLoF variants reduce SMC3 expression but do not result in a blood DNA methylation signature clustering with that of CdLS, and that the global transcriptional signature of SMC3 loss is model-dependent. Our finding of substantial population-scale LoF intolerance in concert with variable penetrance in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multi-layered genomic data paired with careful phenotyping.

Myelodysplasia after clonal hematopoiesis with APOBEC3-mediated CYBB inactivation in retroviral gene therapy for X-CGD.

Stem cell gene therapy using the MFGS-gp91phox retroviral vector was performed on a 27-year-old patient with X-linked chronic granulomatous disease (X-CGD) in 2014. The patient's refractory infections were resolved, whereas the oxidase-positive neutrophils disappeared within 6 months. Thirty-two months after gene therapy, the patient developed myelodysplastic syndrome (MDS), and vector integration into the MECOM locus was identified in blast cells. The vector integration into MECOM was detectable in most myeloid cells at 12 months after gene therapy. However, the patient exhibited normal hematopoiesis until the onset of MDS, suggesting that MECOM transactivation contributed to clonal hematopoiesis, and the blast transformation likely arose after the acquisition of additional genetic lesions. In whole-genome sequencing, the biallelic loss of the WT1 tumor suppressor gene, which occurred immediately before tumorigenesis, was identified as a potential candidate genetic alteration. The provirus CYBB cDNA in the blasts contained 108 G-to-A mutations exclusively in the coding strand, suggesting the occurrence of APOBEC3-mediated hypermutations during the transduction of CD34-positive cells. A hypermutation-mediated loss of oxidase activity may have facilitated the survival and proliferation of the clone with MECOM transactivation. Our data provide valuable insights into the complex mechanisms underlying the development of leukemia in X-CGD gene therapy.

Prevalence of pathogenic variants in cancer-predisposing genes in second cancer after childhood solid cancers.

BACKGROUND: Second malignant neoplasms (SMNs) are one of the most severe late complications after pediatric cancer treatment. However, the effect of genetic variation on SMNs remains unclear. In this study, we revealed germline genetic factors that contribute to the development of SMNs after treatment of pediatric solid tumors. METHODS: We performed whole-exome sequencing in 14 pediatric patients with SMNs, including three brain tumors. RESULTS: Our analysis revealed that five of 14 (35.7%) patients had pathogenic germline variants in cancer-predisposing genes (CPGs), which was significantly higher than in the control cohort (p < 0.01). The identified genes with variants were TP53 (n = 2), DICER1 (n = 1), PMS2 (n = 1), and PTCH1 (n = 1). In terms of the type of subsequent cancer, leukemia and multiple episodes of SMN had an exceptionally high rate of CPG pathogenic variants. None of the patients with germline variants had a family history of SMN development. Mutational signature analysis showed that platinum drugs contributed to the development of SMN in three cases, which suggests the role of platinum agents in SMN development. CONCLUSIONS: We highlight that overlapping effects of genetic background and primary cancer treatment contribute to the development of second cancers after treatment of pediatric solid tumors. A comprehensive analysis of germline and tumor samples may be useful to predict the risk of secondary cancers.

A novel Menin-MLL1 inhibitor, DS-1594a, prevents the progression of acute leukemia with rearranged MLL1 or mutated NPM1.

BACKGROUND: Mixed lineage leukemia 1-rearranged (MLL1-r) acute leukemia patients respond poorly to currently available treatments and there is a need to develop more effective therapies directly disrupting the Menin‒MLL1 complex. Small-molecule-mediated inhibition of the protein‒protein interaction between Menin and MLL1 fusion proteins is a potential therapeutic strategy for patients with MLL1-r or mutated-nucleophosmin 1 (NPM1c) acute leukemia. In this study, we preclinically evaluated the new compound DS-1594a and its salts. METHODS: We evaluated the preclinical efficacy of DS-1594a as well as DS-1594a·HCl (the HCl salt of DS-1594a) and DS-1594a·succinate (the succinic acid salt of DS-1594a, DS-1594b) in vitro and in vivo using acute myeloid leukemia (AML)/acute lymphoblastic leukemia (ALL) models. RESULTS: Our results showed that MLL1-r or NPM1c human leukemic cell lines were selectively and highly sensitive to DS-1594a·HCl, with 50% growth inhibition values < 30 nM. Compared with cytrabine, the standard chemotherapy drug as AML therapy, both DS-1594a·HCl and DS-1594a·succinate mediated the eradication of potential leukemia-initiating cells by enhancing differentiation and reducing serial colony-forming potential in MLL1-r AML cells in vitro. The results were confirmed by flow cytometry, RNA sequencing, RT‒qPCR and chromatin immunoprecipitation sequencing analyses. DS-1594a·HCl and DS-1594a·succinate exhibited significant antitumor efficacy and survival benefit in MOLM-13 cell and patient-derived xenograft models of MLL1-r or NPM1c acute leukemia in vivo. CONCLUSION: We have generated a novel, potent, orally available small-molecule inhibitor of the Menin-MLL1 interaction, DS-1594a. Our results suggest that DS-1594a has medicinal properties distinct from those of cytarabine and that DS-1594a has the potential to be a new anticancer therapy and support oral dosing regimen for clinical studies (NCT04752163).

Identification of the ultrahigh-risk subgroup in neuroblastoma cases through DNA methylation analysis and its treatment exploiting cancer metabolism.

Neuroblastomas require novel therapies that are based on the exploitation of their biological mechanism. To address this need, we analyzed the DNA methylation and expression datasets of neuroblastomas, extracted a candidate gene characterizing the aggressive features, and conducted functional studies. Based on the DNA methylation data, we identified a subgroup of neuroblastoma cases with 11q loss of heterozygosity with extremely poor prognosis. PHGDH, a serine metabolism-related gene, was extracted as a candidate with strong expression and characteristic methylation in this subgroup as well as in cases with MYCN amplification. PHGDH inhibition suppressed neuroblastoma cell proliferation in vitro and in vivo, indicating that the inhibition of serine metabolism by PHGDH inhibitors is a therapeutic alternative for neuroblastoma. Inhibiting the arginine metabolism, which is closely related to serine metabolism using arginine deiminase, had a combination effect both in vitro and in vivo, especially on extracellular arginine-dependent neuroblastoma cells with ASS1 deficiency. Expression and metabolome analyses of post-dose cells confirmed the synergistic effects of treatments targeting serine and arginine indicated that xCT inhibitors that inhibit cystine uptake could be candidates for further combinatorial treatment. Our results highlight the rational therapeutic strategy of targeting serine/arginine metabolism for intractable neuroblastoma.

Multi-omics analysis defines highly refractory RAS burdened immature subgroup of infant acute lymphoblastic leukemia.

KMT2A-rearranged infant acute lymphoblastic leukemia (ALL) represents the most refractory type of childhood leukemia. To uncover the molecular heterogeneity of this disease, we perform RNA sequencing, methylation array analysis, whole exome and targeted deep sequencing on 84 infants with KMT2A-rearranged leukemia. Our multi-omics clustering followed by single-sample and single-cell inference of hematopoietic differentiation establishes five robust integrative clusters (ICs) with different master transcription factors, fusion partners and corresponding stages of B-lymphopoietic and early hemato-endothelial development: IRX-type differentiated (IC1), IRX-type undifferentiated (IC2), HOXA-type MLLT1 (IC3), HOXA-type MLLT3 (IC4), and HOXA-type AFF1 (IC5). Importantly, our deep mutational analysis reveals that the number of RAS pathway mutations predicts prognosis and that the most refractory subgroup of IC2 possesses 100% frequency and the heaviest burden of RAS pathway mutations. Our findings highlight the previously under-appreciated intra- and inter-patient heterogeneity of KMT2A-rearranged infant ALL and provide a rationale for the future development of genomics-guided risk stratification and individualized therapy.

Quantitative assessment of copy number alterations by liquid biopsy for neuroblastoma.

Liquid biopsy, a method of detecting genomic alterations using blood specimens, has recently attracted attention as a noninvasive alternative to surgical tissue biopsy. We attempted quantitative analysis to detect amplification of MYCN (MYCNamp) and loss of heterozygosity at 11q (11qLOH), which are clinical requisites as prognostic factors of neuroblastoma (NB). In this study, cell-free DNA (cfDNA) was extracted from plasma samples from 24 NB patients at diagnosis. Copy numbers of MYCN and NAGK genes were quantitatively analyzed by droplet digital PCR (ddPCR). 11qLOH was also assessed by detecting allelic imbalances of heterozygous single nucleotide polymorphisms in the 11q region. The results obtained were compared to those of specimens from tumor tissues. The correlation coefficient of MYCN copy number of cfDNA and tumor DNA was 0.88 (p < 0.00001). 11qLOH was also accurately detected from cfDNA, except for one case with localized NB. Given the high accuracy of liquid biopsy, to investigate components of cfDNA, the proportion of tumor-derived DNA was estimated by examining the variant allele frequency of tumor-specific mutations in cfDNA. The proportion of tumor-derived DNA in cfDNA was 42.5% (range, 16.9%-55.9%), suggesting sufficient sensitivity of liquid biopsy for NB. In conclusion, MYCN copy number and 11qLOH could be quantitatively analyzed in plasma cfDNA by ddPCR assay. These results suggest that plasma cfDNA can be substituted for tumor DNA and can also be applied for comprehensive genomic profiling analysis.

Novel TENM3-ALK fusion is an alternate mechanism for ALK activation in neuroblastoma.

The identification of molecular events underlying the pathogenesis of neuroblastoma can likely result in improved clinical outcomes for this disease. In this study, a translocation within chromosome 2p and 4q was found to bring about the formation of an in-frame fusion gene that was composed of portions of the teneurin transmembrane protein 3 (TENM3, also known as ODZ3) gene and the anaplastic lymphoma kinase (ALK) gene in tumor cells from patients with neuroblastoma. Expression of the full length TENM3-ALK cDNA in NIH-3T3 cells led to the formation of a fusion protein that: (1) possesses constitutive tyrosine kinase activity, (2) induces strong activation of the downstream targets of extracellular signal-regulated kinase (ERK), protein kinase B (a.k.a. AKT), and signal transducer and activator of transcription 3 (STAT3), (3) provokes oncogenic transformation in NOD.Cg-PrkdcscidIl2rgtm1Sug/ShiJic mice, and (4) possesses sensitivity to ALK inhibitors in vitro and in vivo. Our findings demonstrated that patients with neuroblastoma may express a transforming fusion kinase, which is a promising candidate for a therapeutic target and a diagnostic molecular marker for neuroblastoma. The in-frame 5' partner gene that fuses with ALK has not been reported previously in neuroblastoma. Our data provide novel biological insights into the mechanism of ALK activation due to translocation, with implications for neuroblastoma tumorigenesis, and could be useful as a vital marker for the accurate diagnosis of this type of neuroblastoma.

Genetic features of B-cell lymphoblastic lymphoma with TCF3-PBX1.

BACKGROUND: Lymphoblastic lymphoma (LBL) and acute lymphoblastic leukemia (ALL) are categorized as the same entity under precursor lymphoid neoplasms in the World Health Organization classification. However, compared to B-cell ALL, the molecular genetic makeup of B-cell LBL remains to be understood, mainly due to its rarity. We performed whole exome sequencing (WES) on seven patients with TCF3-PBX1-positive B-cell LBL. METHODS: WES was performed using DNA extracted from tumor specimens and paired blood samples at remission for six patients, and tumor-only analysis was performed for one patient whose remission sample was not available. For one patient, a relapsed sample was also analyzed. RESULTS: KMT2D variants and 6q LOH were found as recurrent alterations. Somatic variants of KMT2D were identified in three of the seven patients. Of note, the two patients with heterozygous nonsense variant of KMT2D were at stage III, without bone marrow infiltration. 6q LOH was also identified in two others, out of the seven patients. The common 6q deleted region of the two patients ranged from 6q12 to 6q16.3. Both patients had bone marrow infiltration. Analysis of recurrent case also revealed that the relapsed clone might be derived from a minor clone of the bone marrow at diagnosis. CONCLUSION: In this study, through WES for seven patients with TCF3-PBX1-positive B-LBL, we identified KMT2D mutations and 6q LOH as recurrent alterations. In order to elucidate the relationship between these recurrent alterations and disease specificity or outcomes, further studies comparing with TCF3-PBX1-positive B-ALL are required.

Genetic features of precursor B-cell phenotype Burkitt leukemia with IGH-MYC rearrangement.

BACKGROUND: An atypical form of Burkitt leukemia/lymphoma (BL), BL with a phenotype of precursor B-cells (preBLL), is listed in the WHO Classification. Recent reports suggested that preBLL and classical BL could be distinguished by the differences in IG-MYC translocation architecture and an additional mutated genes profile. The characteristics of classical BL are IG-MYC by aberrant somatic hypermutation or class switch recombination, and BL-specific gene mutations such as MYC, ID3, and CCND3. Meanwhile, preBLL is characterized by IG-MYC due to aberrant VDJ recombination and mutations in NRAS and KRAS. However, it is not clear whether all preBLL cases can be differentiated. This report investigated the molecular characteristics of an infant preBLL case, with a more advanced stage of maturity than typical preBLL. CASE: The patient showed BL-like morphology with IGH-MYC rearrangement. In the immunophenotyping, CD20 and surface immunoglobulin were negative, whereas other markers were consistent with BL. To evaluate the genetic contribution, we performed whole-exome sequencing. The breakpoint analysis revealed the IG-MYC occurred due to an aberrant VDJ recombination. Meanwhile, additional somatic mutations were detected in FBXO11, one of the mutant genes specific to BL. In the analysis of the specimen in complete remission, mutation in KRAS, frequently mutated in preBLL, was detected with low frequency, suggesting somatic mosaicism. CONCLUSION: The present case showed the characteristics of both typical preBLL and classical BL. Because preBLL includes atypical cases such as the present case, further studies are required to elucidate preBLL features.

NUDT15 variants confer high incidence of second malignancies in children with acute lymphoblastic leukemia.

The effect of genetic variation on second malignant neoplasms (SMNs) remains unclear. First, we identified the pathogenic germline variants in cancer-predisposing genes among 15 children with SMNs after childhood leukemia/lymphoma using whole-exome sequencing. Because the prevalence was low, we focused on the association between SMNs and NUDT15 in primary acute lymphoblastic leukemia (ALL) cases. NUDT15 is one of the 6-mercaptopurine (6-MP) metabolic genes, and its variants are common in East Asian individuals. The prevalence of NUDT15 hypomorphic variants was higher in patients with SMNs (n = 14; 42.9%) than in the general population in the gnomAD database (19.7%; P = .042). In the validation study with a cohort of 438 unselected patients with ALL, the cumulative incidence of SMNs was significantly higher among those with (3.0%; 95% confidence interval [CI], 0.6% to 9.4%) than among those without NUDT15 variants (0.3%; 95% CI, 0.0% to 1.5%; P = .045). The 6-MP dose administered to patients with ALL with a NUDT15 variant was higher than that given to those without SMNs (P = .045). The 6-MP-related mutational signature was observed in SMN specimens after 6-MP exposure. In cells exposed to 6-MP, a higher level of 6-MP induced DNA damage in NUDT15-knockdown induced pluripotent stem cells. Our study indicates that NUDT15 variants may confer a risk of SMNs after treatment with 6-MP in patients with ALL.

Copy number alteration analysis for neuroblastoma using droplet digital polymerase chain reaction.

BACKGROUND: Neuroblastoma (NB) is a malignant tumor derived from the neural crest. MYCN amplification is a well-known adverse molecular prognostic factor for NB. Genome copy number alterations (CNAs) such as chromosome (Chr) 11q deletion, 1p deletion, and 17q gain are associated with a poor prognosis. Fluorescence in situ hybridization (FISH) and Southern blotting analysis are frequently used to detect MYCN amplification. Although comparative genomic hybridization (CGH) and single-nucleotide polymorphism (SNP) chip arrays can easily detect CNAs, these methods are impractical for clinical use due to their cost and run time. Consequently, genome copy number analysis using digital droplet PCR has become widely used to monitor CNAs. METHODS: In this study, we used digital droplet polymerase chain reaction to detect MYCN amplification and Chr 11q CNA, which was used for risk stratification according to the International Neuroblastoma Risk Group classification system. We compared the results with data from SNP chip arrays in seven NB cell lines and eight primary NB samples. RESULTS: Digital droplet PCR assays successfully detected MYCN amplification and 11q CNA. The results were very consistent with those obtained by SNP chip assay. CONCLUSIONS: Digital droplet PCR can be conducted more rapidly than FISH or Southern blotting. Accordingly, it should be useful for on-site clinical applications aimed at detecting CNAs in NB and performing risk stratification promptly after diagnosis.

BICRA, a SWI/SNF Complex Member, Is Associated with BAF-Disorder Related Phenotypes in Humans and Model Organisms.

SWI/SNF-related intellectual disability disorders (SSRIDDs) are rare neurodevelopmental disorders characterized by developmental disability, coarse facial features, and fifth digit/nail hypoplasia that are caused by pathogenic variants in genes that encode for members of the SWI/SNF (or BAF) family of chromatin remodeling complexes. We have identified 12 individuals with rare variants (10 loss-of-function, 2 missense) in the BICRA (BRD4 interacting chromatin remodeling complex-associated protein) gene, also known as GLTSCR1, which encodes a subunit of the non-canonical BAF (ncBAF) complex. These individuals exhibited neurodevelopmental phenotypes that include developmental delay, intellectual disability, autism spectrum disorder, and behavioral abnormalities as well as dysmorphic features. Notably, the majority of individuals lack the fifth digit/nail hypoplasia phenotype, a hallmark of most SSRIDDs. To confirm the role of BICRA in the development of these phenotypes, we performed functional characterization of the zebrafish and Drosophila orthologs of BICRA. In zebrafish, a mutation of bicra that mimics one of the loss-of-function variants leads to craniofacial defects possibly akin to the dysmorphic facial features seen in individuals harboring putatively pathogenic BICRA variants. We further show that Bicra physically binds to other non-canonical ncBAF complex members, including the BRD9/7 ortholog, CG7154, and is the defining member of the ncBAF complex in flies. Like other SWI/SNF complex members, loss of Bicra function in flies acts as a dominant enhancer of position effect variegation but in a more context-specific manner. We conclude that haploinsufficiency of BICRA leads to a unique SSRIDD in humans whose phenotypes overlap with those previously reported.

Comprehensive genetic analysis of pediatric germ cell tumors identifies potential drug targets.

To elucidate the molecular pathogenesis of pediatric germ cell tumors (GCTs), we performed DNA methylation array analysis, whole transcriptome sequencing, targeted capture sequencing, and single-nucleotide polymorphism array analysis using 51 GCT samples (25 female, 26 male), including 6 germinomas, 2 embryonal carcinomas, 4 immature teratomas, 3 mature teratomas, 30 yolk sac tumors, and 6 mixed germ cell tumors. Among the 51 samples, 11 were from infants, 23 were from young children, and 17 were from those aged ≥10 years. Sixteen of the 51 samples developed in the extragonadal regions. Germinomas showed upregulation of pluripotent genes and global hypomethylation. Pluripotent genes were also highly expressed in embryonal carcinomas. These genes may play essential roles in embryonal carcinomas given that their binding sites are hypomethylated. Yolk sac tumors exhibited overexpression of endodermal genes, such as GATA6 and FOXA2, the binding sites of which were hypomethylated. Interestingly, infant yolk sac tumors had different DNA methylation patterns from those observed in older children. Teratomas had higher expression of ectodermal genes, suggesting a tridermal nature. Based on our results, we suggest that KIT, TNFRSF8, and ERBB4 may be suitable targets for the treatment of germinoma, embryonal carcinomas, and yolk sac tumors, respectively.