Abstract
Astroblastomas (ABs) are rare brain tumors of unknown origin. We performed an integrative genetic and epigenetic analysis of AB-like tumors. Here, we show that tumors traceable to neural stem/progenitor cells (radial glia) that emerge during early to later brain development occur in children and young adults, respectively. Tumors with MN1-BEND2 fusion appear to present exclusively in females and exhibit overexpression of genes expressed prior to 25 post-conception weeks (pcw), including genes enriched in early ventricular zone radial glia and ependymal tumors. Other, histologically classic ABs overexpress or harbor mutations of mitogen-activated protein kinase pathway genes, outer and truncated radial glia genes, and genes expressed after 25 pcw, including neuronal and astrocyte markers. Findings support that AB-like tumors arise in the context of epigenetic and genetic changes in neural progenitors. Selective gene fusion, variable imprinting and/or chromosome X-inactivation escape resulting in biallelic overexpression may contribute to female predominance of AB molecular subtypes.
Similar content being viewed by others
Introduction
Astroblastomas (ABs) are uncommon and historically controversial CNS neoplasms characterized by tumor cells surrounding a central capillary forming a flower petal-like structure called the astroblastic pseudorosette (Fig. 1a)1,2,3. AB-like tumors typically occur in the cerebral hemispheres of children and young adults and show considerable female predominance3,4,5. Their clinical behavior ranges from frequently surgically-amenable benign lesions to occasionally malignant. A long-standing question is whether ABs are of astrocytic or ependymal lineage2.
AB-like tumors are comprised of at least two genetic types based on the presence of the pV600E mutation of the B-Raf serine/threonine kinase gene (BRAFV600E)2,4,6 or rearrangement of the meningioma (disrupted in balanced translocation) 1 (MN1) gene4,6. Fusions between MN1 at chromosome region 22q12.1 and BEN domain containing 2 (BEND2) at Xp22.13 are reported in MN1-rearranged tumors6.
Embryonic NSCs include highly proliferative neuroepithelial cells (NEC), which comprise the primitive neural tube and give rise to early neuronal progenitors (neuroblasts). Around five post-conception weeks (pcw) NEC begin to transform into ventricular zone radial glia (vRG), bipolar cells whose apical and basal processes contact the ventricular and pial surfaces of the developing brain, respectively7,8 (Fig. 1b). NEC and vRG exhibit interkinetic nuclear migration (IKNM), which helps maintain the tightly packed cellular structure of the ventricular lining during cell division. Like NEC, vRG can directly give rise to neuronal progenitors through asymmetric division. Further such vRG divisions produce B1 astrocytes (future adult subventricular zone NSCs), followed by ependymal cells9. The latter replace vRG in lining the brain ventricles.
vRG also produce outer radial glia (oRG), whose nuclei reside in the outer subventricular zone (OSVZ), and which appear at about 12 pcw10. Unlike vRG, oRG lack a prominent apical process and demonstrate mitotic somal translocation (MST) in which the cell body migrates basally prior to mitosis, facilitating expansion of the OSVZ and cortical plate11. They generate intermediate progenitor (IP) cells capable of further division termed transiently amplifying progenitors which generate the majority of cortical neurons. Lastly, oRG also give rise to astrocytes and oligodendrocytes.
Nowakowski et al.12 observed that the long basal processes of vRG become truncated in the human fetus at about 15 pcw. Unlike vRG and oRG, truncated radial glia (tRG) extend their basal processes only into the OSVZ where their end-feet typically surround blood capillaries similar to classic astroblastoma cells1,2 (Fig. 1a, b).
DNA methylation in the human brain is dynamic during development. Some genes on the inactive X chromosome in females can escape X inactivation13,14,15,16. These X-inactivation escape (XIE) genes, can acquire promoter hypomethylation in somatic cells resulting in biallelic expression from both X chromosomes, potentially doubling their gene dosage13.
Here we show that AB-like tumors with MN1-BEND2 fusions appear to exclusively occur in females, and nearly always in children. They exhibit overexpression of IGF2-H19 and DLK1-DIO3 region, and other imprinted genes, XIE genes, and genes highly expressed during brain development prior to 25 pcw, including genes enriched in vRG and ependymal tumors. Other AB-like tumors overexpress mitogen-activated protein kinase (MAPK) pathway, phosphatidylinositol-3-kinase (PI3K) pathway and alternate XIE genes, and tend to occur in young adults of both sexes depending on molecular subtype. These tumors highly express oRG and tRG genes, and genes expressed after 25 pcw that are often involved in neuronal or astrocytic differentiation. Our findings strongly support that AB-like tumors arise in the context of epigenetic and genetic changes in neural progenitors during brain development: early ependymal tumors with MN1-BEND2 from early vRG-derived cells, and MAPK pathway tumors demonstrating more classic astroblastoma histomorphology from oRG- or tRG-derived progenitors, respectively. Data also suggest that increased gene fusion, variable imprinting, and XIE may contribute to the female predominance of certain AB-like tumor molecular subtypes. Lastly, we found that like MN1-BEND2 tumors, nascent immature ependymal cells express IGF2 and may represent an important source of this growth factor in the fetal lateral ventricular zone niche.
Results
Tumors with AB histology segregate into distinct subtypes based on transcriptomic and mutational profiles
We performed RNAseq on 24 supratentorial AB-like cases with available material, including six BRAFV600E-mutant tumors2, six MN1-rearranaged tumors4 and twelve tumors without these genetic lesions (Table 1). Three tumors diagnosed as supratentorial ependymoma (STE) and four papillary tumor of the pineal region (PTPR)17 were also studied for comparison. Affymetrix gene microarray (HTA-2) expression analysis was additionally used to study 15 of the above AB-like tumors, one MN1-rearranaged tumor recurrence (C10r), the four PTPR, six histologically classic pleomorphic xanthoastrocytomas (PXA), four fetal brain samples (28–30 pcw) and 14 normal pediatric and adult brain samples.
Dimensionality reduction of the two expression datasets was then performed individually using tSNE (23,937 mappable genes in the Affymetrix data and 24,917 mappable and encodable RNAseq gene transcripts). Next, we identified 1793 genes differentially expressed among the resultant expression clusters that were common to both datasets using analysis of variance (P < 0.05).
Unsupervised hierarchical clustering of these 1793 genes using the RNAseq data yielded seven tumor clusters (Fig. 1c), similar to tSNE (Fig. 1e). Tumors with MN1 rearrangement, including five containing the MN1-BEND2 gene fusion (Table 1), formed a discrete RNA cluster, as did all BRAFV600E-mutant tumors combined with three tumors without BRAFV600E. The latter showed mutations of other MAPK, and frequently PI3K/AKT/mTOR pathway genes (Fig. 2, Supplementary Data 1). We designated this RNA cluster MAPK-A. Three additional tumors without BRAFV600E mutation, but otherwise showing similar mutation profiles as MAPK-A tumors, closely clustered with them and were designated MAPK-B (Fig. 1c, e). Hence forth, we refer to these two groups collectively as MAPK-AB lesions. Three tumors lacking both BRAFV600E mutation or obvious PI3K pathway mutations, designated MAPK-C, also formed a subgroup within the MAPK supercluster (MAPK-ABC).
Three histologically atypical AB-like tumors (ATAB) also formed a transcriptomic cluster (Fig. 1c, e). These occurred in older adults and showed high-grade histology (Table 1). Lastly, three STEs (two containing the ZFTA-RELA fusion), and the four PTPR grouped into respective RNA clusters.
MN1-BEND2 containing tumors, MAPK supercluster astroblastomas, and pleomorphic xanthoastrocytomas (PXAs) are genetically distinct entities
Hierarchical clustering of the Affymetrix HTA-2 data using the same 1,793 genes resulted in nine clusters: 1) MN1-rearranged tumors; 2) fetal brain control tissue; 3) PTPR; 4) MAPK-ABC tumors; 5) ATAB; 6) PXA; and 7-9) normal adult and pediatric brain control tissue (Fig. 1d, f). Although BRAFV600E-mutant AB-like tumors loosely group with PXA by genomic DNA methylation4, all AB-like tumors clustered separately from PXA in hierarchical clustering and tSNE analyses of tumor transcriptomes. Notably, gene expression in normal control tissues was markedly different from that of tumors (Fig. 1d). These data indicate that, albeit likely ontologically related2, MAPK-ABC tumors are genetically distinct from PXA.
Neural stem cell and early ependymal genes are highly expressed in MN1-BEND2 tumors
CUX2, SHOX, SOX14, IRX2, PAX1, HOXD10, DLX5, and PRRX2 are homeobox genes relatively overexpressed in MN1-BEND2 tumors (Fig. 1c, d and Supplementary Figs. 1–3). Interestingly, the fetal growth factor gene insulin growth factor 2 (IGF2) and the ependymal fate determining gene FOXJ19, both markers of intracranial ependymoma18, were relatively overexpressed by MN1-BEND2 tumors, as well as STE (Supplementary Figs. 1, 2). MN1-BEND2 tumors did not overexpress the STE-associated and neuronal marker gene ANK219,20. However, MN1-BEND2 tumors did overexpress CELSR1, encoding a planar cell polarity protein necessary for neural tube closure21 and maintenance of intercellular ciliary basal body (centriole) alignment in vRG and ependymal cells22. They also highly expressed SUN1, encoding a protein required for IKNM in NEC and vRG23 and CHD3, a part of the Mi-2/NuRD histone deacetylase complex (Fig. 1c, d).
MN1-BEND2 tumors and PTPR share ependymal gene expression and low expression of oRG, glial and neuronal genes
SOX14, IRX2, CELSR1, KCNJ5 encoding the inward rectifier potassium channel Kir3.4, and TFF3 encoding an epithelial cell mobility protein, were overexpressed by MN1-rearranged tumors and PTPR almost exclusively (Fig. 1c, d and Supplementary Figs. 1–3). Both also strongly expressed FOXJ124,25 and the additional ependymoma-associated genes RFX3, encoding a transcription factor essential for development of motile cilia26, YAP127 and NELL2 (Supplementary Figs. 1–3). Overall, these findings argue for ependymal differentiation of both MN1-BEND2 tumors and PTPR. The latter differ from MN1-BEND2 tumors and many STEs, however, in that they are IGF2 relative under-expressers. This may be related to their 3rd ventricular diencephalic derivation, in contrast to the likely telencephalic origin of MN1-BEND2 tumors and STE. Moreover, the fact that PTPR more highly express the NSC markers PAX6 and OTX2 argues they may be more closely related to a RG stem-like cell28, perhaps derived directly from such a cell and lacking an IP possibly required in the genesis of MN1-BEND2 tumors and STE.
In contrast to MAPK supercluster tumors, MN1-rearranged tumors, PTPR and STE showed relatively low expression of the neuronal genes GRIA2 (glutamate ionotropic receptor AMPA type subunit 2), MAOA (monoamine oxidase A, an XIE gene) and MAPT (microtubule associated protein tau). They also demonstrated low expression of the astrocyte markers OLIG2, GFAP, S100β and ALDH1, the oRG markers IDH1 and TNC10, and other genes highly expressed in MAPK-ABC tumors, e.g., TCF4, LRP1B, CTNND2, NCAN, DISC1, HIF1A, MYCBP2, ABCC5, and ANK2. PTPR and STE however do overexpress some genes also highly expressed in MAPK supercluster tumors, but not MN1-rearranged tumors (Fig. 1c, d, Supplementary Figs. 1–3).
TCF4 is a WNT transcription factor. LRP1B and CTNND2 are WNT pathway regulatory genes29. The latter encodes δ-catenin 2, which affects dendritic spine and synapse formation, and is a marker of oRG and neurons10. DISC1 encodes disrupted in schizophrenia 1, which is involved in neurogenesis and astrogenesis, and is a positive regulator of WNT signaling and negative regulator of PI3K/AKT/mTOR signaling30,31. The relatively high expression of these genes in MAPK-ABC tumors suggests that the WNT pathway is important in their biology.
MN1-BEND2 and MAPK supercluster tumors show characteristic mutational profiles
In addition to indel, stop gain and stop lost mutations, only missense or splice site mutations found to be “deleterious” by SIFT, “possibly or probably damaging” by PolyPhen and/or deleterious by FathmmMKL32 (defined as FathmmMKL scores above 0.5), are considered herein unless otherwise specified (Supplementary Data 1). MN1-rearranged tumors demonstrated multiple recurrently mutated genes, including ABCC1, IGF2, ABCF1, and MSH3 (Fig. 2). ABCC1 encodes the ATPase cassette transporter protein multidrug resistance 1, a stem cell marker and MYCN transcriptional target conferring resistance to chemotherapeutics. An ABCC1-SSU72 fusion was present in case C10, and ABCC1-LUCAL3 and ABCC1-ZHHC14 fusions were found in C3 (Fig. 2, Supplementary Data 2). A heterozygous ABCC1 single nucleotide substitution (SNS), p1088P/L was present in the MN1-rearanged tumor lacking MN1-BEND2 fusion from a male patient (C16). Predicted deleterious ABCC1 mutations were additionally noted in ZFTA-RELA STE and one ATAB in which ZFTA-RELA was detected by a single read (C4). ABCC1 was also relatively overexpressed in MN1-rearanged and ZFTA-RELA tumors (Fig. 1c, d and Supplementary Figs. 1, 2), suggesting they share a common ependymal stem cell lineage. The mismatch repair gene MSH3 was mutant in four of six MN1-rearranged cases.
Mutations in PI3K/AKT/mTOR pathway genes, e.g., PIK3 subunits, PTEN, TSC2 and NF1, were found in most MAPK-AB and ATAB tumors (Fig. 2). Deleterious mutations of the NF1 tumor suppressor gene occurred in three MAPK-AB tumors, the MN1-BEND2-negative MN1-rearranged tumor, and an ATAB. PI3K/AKT/mTOR pathway gene mutations were notably absent in MN1-BEND2 and MAPK-C tumors.
TCF4 was mutant in two MAPK-B tumors and one MAPK-C tumor, and involved in a fusion in a BRAFV600E-mutant MAPK-A tumor (Fig. 2, Supplementary Data 1, 2). ANK2 and ANK3 mutations were found in non-BRAF-mutant MAPK-AB tumors. The BRAFV600E-mutation negative case C14 occurred in the oldest patient in the MAPK-A group. This case showed multiple PI3K pathway alterations, including a TSC2 mutation, PTEN-ANK3 fusion and an out-of-frame RP5-864K19.4-TRIM27 fusion. TRIM27 encodes an E3 ligase that ubiquitinates phosphatidylinositol-4-phosphate 3-kinase catalytic subunit type 2 beta (PIK3C2B)33. Its disruption may thus lead to increased PI3K signaling.
MAPK-C and other non-MN1-BEND2, non-BRAFV600E–harboring lesions demonstrate unique mutations
MAPK-C tumors occurred in males and showed non-fusion mutations in genes involved in actin and adherens junction binding, including KALRN, TRIO, DISC1, and ANK3 (Fig. 2, Supplementary Data 1–3). KALRN and TRIO encode Rho guanine nucleotide exchange factors, which are downstream effectors of noncanonical WNT signaling34,35. KALRN, TRIO, and DISC1 are involved in synapse formation and plasticity, axonal growth and adolescence-associated dendritic spine pruning36,37. Interestingly, fusions involving KALRN or TRIO were only detected in BRAFV600E-mutant MAPK-A and STE tumors from females, e.g., KALRN-SESN3, KALRN-PLXND1 and TRIO-NDUFV3, and TRIO-DNAH5 and TTC21B-KALRN, respectively. DISC1 was mutant in MAPK-C and STE tumors.
PARP8, PTEN, ANK2, ANKFY1, and KIAA1109 non-fusion mutations were unique to non-MN1-BEND2, non-BRAFV600E tumors. KIAA1109 is important in cortical development38 and mutant in glioblastoma39. ANK genes encode ankyrins, which link membrane proteins to the actin cytoskeleton and are involved in dendritic spine growth and synaptic plasticity40.
Mutations in CTTN (cortactin) or CTTNBPNL (cortactin-binding protein 2 N-terminal–like) occurred in non-BRAFV600E MAPK-ABC tumors. Cortactin is an actin binding protein involved in cytoskeletal remodeling and implicated in tumor invasion and WNT-dependent synaptic plasticity41,42. Three tumors contained VEZT mutations. Its protein product vezatin binds the cadherin-catenin complex linking adherens junctions to the actin cytoskeleton43.
The most frequently mutated gene in non-MN1-BEND2, non-BRAFV600E tumors was AHNAK, mutated in six such tumors, including MAPK-C tumors. AHNAK encodes neuroblast differentiation-associated protein, a cadherin-associated actin binding protein implicated in neural cell differentiation, migration, and tumor metastasis44.
DNA methylome analysis reveals a distinct methylation class for MAPK-C tumors
We performed genomic DNA methylation analysis of twelve AB-like tumors not previously analyzed. When these data were combined with that from our prior study4, tSNE analysis revealed that all MN1-BEND2 tumors clustered with the methylation class “MN1 high-grade neuroepithelial tumor (MN1-HGNET)” (Fig. 3).
MAPK-A AB tumors containing the BRAFV600E mutation again loosely grouped with the “pleomorphic xanthoastrocytoma (PXA)” methylation class4 (Fig. 3). BRAF-wildtype MAPK-A tumors C14, C20 and C39, and an additional AB-like tumor not studied by RNAseq (C15), appeared on the periphery of the PXA methylation cluster. The three MAPK-C transcriptional cluster tumors, one MAPK-B tumor (C8) and two tumors not studied by RNAseq were among six cases (C8, C18, C19, C22, C27, and C41) forming a distinct genomic DNA methylation class, designated the Group C methylation class.
Of the remaining non-BRAFV600E cases, including the MN1-rearranged tumor lacking MN1-BEND2 (C16), three grouped with the nonspecific “low-grade glioma/supratentorial pilocytic astrocytoma/ganglioglioma” methylation class6, and C1 containing the ZFTA-RELA fusion, grouped with RELA ependymomas. A STE lacking ZFTA-RELA grouped with posterior fossa type A ependymomas by DNA methylation, and two ATAB grouped with glioblastoma subtypes (Fig. 3).
AB-like tumor histology is characterized by astroblastic pseudorosettes, and oRG/tRG-like classic astroblastoma cytomorphology in Group C tumors
It is not possible to precisely differentiate AB-like tumor genetic subtypes solely by histologic features, as tumor cytology is variable from cuboidal to tapered-columnar to spindle-shaped, often within the same tumor2,4. In general, however, MN1-BEND2 tumor cellular morphology tends to be more cuboidal to stout-columnar and these tumors more often demonstrate vascular sclerosis, and lack eosinophilic granular material, in comparison to MAPK supercluster tumors4 (Fig. 4a [a–c]). MAPK-ABC tumors (Fig. 4a [d–i]) more frequently demonstrate tapering basal processes and inconspicuous or absent apical processes, especially Group C tumors (Fig. 4a [g–i]). The latter contained elongate monopolar cells whose endfeet contacted the brain pial surface or basal lamina of blood vessels mimicking oRG and tRG, respectively, similar to astroblastoma cells originally depicted by Bailey and Bucy1 (Fig. 4b [a, b]). Like MAPK-AB tumors, Group C tumors also demonstrated eosinophilic granular bodies and occasional multinucleated cells2, as previously noted in astroblastomas by Bailey and Bucy1.
AB tumor subtypes show characteristic protein expression
MN1-BEND2 tumors exhibited relative underexpression of GFAP mRNA (Fig. 1c, d, Supplementary Figs. 1–3). They were also predominantly GFAP negative by immunohistochemical staining or showed only patchy GFAP immunoreactivity (Fig. 4c [a–c]), similar to ependymomas. Conversely, MAPK-ABC tumors were diffusely GFAP immunoreactive (Fig. 4c [d–f]). Human vRG are present by 5 pcw and first become GFAP-positive by immunohistochemistry at about 11 pcw45. Thus, similar to MN1-BEND2 tumors, early vRG do not show diffuse GFAP immunoreactivity.
Strikingly, all MN1-BEND2 tumors demonstrated diffuse nuclear MN1 positivity (Fig. 4d [a–c]). Some MAPK-ABC tumors showed probable nonspecific cytoplasmic MN1 immunoreactivity; however, none exhibited diffuse MN1 nuclear immunoreactivity (Fig. 4d [d–f]). MN1 and BEND2 immunoreactivity were also confirmed by western blot in a MN1-BEND2 case with available frozen tumor (C29) (Supplementary Fig. 4). MN1-BEND2 tumors also demonstrated strong cytoplasmic IGF2 immunoreactivity (Fig. 4e [a–c]), which was absent or weak in non-MN1-BEND2 AB-like tumors (Fig. 4e [d–f]).
MN1-BEND2 and MAPK-ABC tumor gene expression overlaps that of NEC and vRG, and tRG and oRG, respectively
We next compared expression of relatively overexpressed or mutated genes in MN1-BEND2 versus MAPK-ABC AB-like tumors to that in the developing human brain utilizing the Allen Brain Atlas Human Developmental Transcriptome database46. Several transcripts selectively expressed or mutated in MN1-BEND2 tumors, e.g., MN1, IGF2, ABCC1, and VEZF1, are more highly expressed in the fetal brain prior to 25 pcw (Figs. 4f, g, 5). Many such transcripts, e.g., NXN, CTTNBP2NL and H19, are highly expressed prior to 9 pcw when cellular proliferation in the developing brain peaks. Indeed, H19 message is enriched in human vRG as determined by single cell RNAseq47. The period prior to 25 pcw corresponded to the highest initial expression of NEC genes within the database, e.g., SOX1, DACH1, PROM1, NOTCH1, NOTCH2, ITGA6 and NRARP; vRG genes, e.g., H19, FOXJ1 and HES1; the tRG gene PALLD, and a few oRG genes, e.g., LIFR and SEZ6L. In contrast, genes more highly expressed or mutated in MAPK-ABC tumors, e.g., GFAP, BRINP3 and TRIM9, and certain oRG genes, e.g., FAM107A and HOPX, are more highly expressed subsequent to 25 pcw when genes related to dendritic spine and synapse development, e.g., CLU, DLG1, KALRN, SLITRK2, become more highly expressed48 (Fig. 5). The latter period corresponds to the end of widespread cortical neurogenesis and the beginning of cerebral gliogenesis. In contrast to genes highly expressed in MN1-BEND2 tumors that tended to cluster with NEC and vRG marker genes, those selectively expressed in MAPK-ABC tumors tended to cluster with markers of vRG and oRG (Fig. 5). These data suggest that MN1-BEND2 tumors are derived from cells closely related to early vRG that maintain high expression of many NEC genes, and that MAPK-ABC tumors are derived from cells related to oRG that retain expression of some vRG markers.
In order to further explore possible neural stem/progenitor cells which AB-like tumors might share ontologic lineage, we compared published neural stem/progenitor cell gene expression sets10,49,50 (Supplementary Data 4) to tumor gene expression using hierarchical clustering (Supplementary Fig. 5) and GSEA (Supplementary Fig. 6) NEC, general RG (gRG), vRG, tRG, and oRG signature genes were significantly upregulated in nearly all tumor types compared to normal pediatric and adult control brain tissues (based on our stringent family-wise error rate P value, p[FWER]) of ≤0.025). The only exceptions were vRG genes in PTPR and NEC genes in PXA (Supplementary Fig. 6A–D, Supplementary Data 5).
NEC genes were significantly enriched in MN1-BEND2 tumors versus control fetal brain, and relatively enriched in MN1-BEND2 tumors compared to PTPR, MAPK-ABC and PXA. The latter approached statistical significance (pFEWR = 0.028) (Supplementary Fig. 6E–H, Supplementary Data 5). MN1-BEND2 tumors were relatively enriched in vRG genes compared to all other tumors, however, this was not statistically significant (Supplementary Fig. 6J–L).
The most significant and consistent GSEA findings for NSC genes were enrichment of PTPRs for tRG genes and PXAs for oRG genes. PTPR were significantly enriched for tRG gene expression compared to controls and all other tumor types, except PXA: to fetal brain (pFEWR < 0.0001), MN1-BEND2 tumors (pFEWR = 0.004), ATAB (pFEWR = 0.001), MAPK-ABC (pFEWR = 0.002) and control pediatric and adult brain (pFEWR < 0.001) (Supplementary Fig. 6 M–P, Supplementary Data 5).
Both MAPK–ABC tumors and PXA were significantly enriched for only tRG and oRG compared to fetal brain (Supplementary Fig. 6Q, Supplementary Data 5). Compared to MN1-BEND2 tumors, MAPK-ABC tumors were significantly enriched only for oRG genes, but showed more frequent mutation of tRG and oRG genes (Supplementary Fig. 7). PXA was enriched for tRG, and more so for oRG compared to MN1-BEND2 and MAPK-ABC tumors (Supplementary Fig. 6P and T, pFEWER < 0.001, Supplementary Data 5). ATAB showed no clear stem/progenitor cell gene enrichment pattern likely attributable to the heterogeneity of this group.
The MN1-BEND2 tumor markers MN1 and IGF2 are expressed by differentiating ependymal cells in the P0–P6 mouse ventricular zone
MN1-immunoreactivity was detected in the nucleus and cytoplasm of nascent immature ependymal cells of the mouse pup ventricular zone (Fig. 6e, g, l–o). These MN1-positive ependymal precursors are in an intermediate stage of differentiation, referred to as the individualized stage, and exhibit centrioles and few to multiple short cilia. This stage occurs between the earliest stage of ependymal differentiation, i.e., the halo stage, involving generation of multiple pro-centrioles, and fully mature ependyma demonstrating mature centrioles and long motile cilia51 (Fig. 6). Cytoplasmic staining for IGF2 was also observed in these cells (Fig. 6b, h–k), suggesting that MN1-BEND2 or other ependymal tumors may originate from similar immature ependymal cells derived from vRG9,52. IP cells basal to the OSVZ also stained for MN1, but not IGF2. Nevertheless, such cells could also theoretically give rise to MN1-BEND2 tumors.
MN1-BEND2 tumors show increased expression of imprinted genes
Unsupervised hierarchical clustering of 139 currently annotated imprinted genes culled from the literature53,54,55 (Supplementary Data 6), revealed overexpression of distinct signature blocks of imprinted genes in MN1-rearranged tumors as compared to MAPK-ABC lesions (Fig. 7a, b and Supplementary Figs. 8, 9). Examples include genes from the IGF2-H19 imprinted locus on chromosome 11p15, e.g., IGF2, H19, and miR483, and DLK1, DIO3, and miR134, which reside within the DLK1-DIO3 imprinted locus at 14q32 and whose misregulation is linked to schizophrenia and bipolar disorder56,57. miR483 was also found to be highly expressed in MN1-BEND2 tumors and RELA STE compared to non-MN1-BEND2 AB-like tumors independently by Nanostring methodology (Supplementary Fig. 10). GNAS, an imprinted NEC gene encoding a G-protein α-subunit activating adenylate cyclase, was also relatively overexpressed in MN1-BEND2 and some other AB-like tumors, but only in female patients (Supplementary Figs. 7, 8). Most of these genes were also upregulated in MN1-rearranged tumors compared to normal and fetal brain (Fig. 7b and Supplementary Fig. 9).
GSEA revealed that imprinted gene transcripts are significantly enriched in MN1-BEND2 compared to MAPK-ABC tumors (ES = |0.43| , pFWER < 0.0001) (Fig. 7d) and all other tumor types and controls (pFWER ≤ 0.0001) (Supplementary Data 7).
Fusions involving imprinted and chromosome X genes occur predominantly in AB-like tumors from female patients
Fusions containing the imprinted gene maternally expressed 8 (MEG8), encoding a small nucleolar RNA, were only present in tumors from females, in four MAPK-AB lesions and four of five MN1-BEND2 tumors (Fig. 4f, Supplementary Fig. 11A, Supplementary Data 2). Most were SNHG23-MEG8 (chromosome 14:14) fusions, with the exception of SNORD114-26-MEG8 (14:14) in two BRAFV600E-mutant cases and RUNX1-MEG8 (21:14) in one MN1-BEND2 case. SNORD114 resides on chromosome 14 within the same imprinted region as MEG8. All were out-of-frame rearrangements demonstrated by multiple reads. These fusions correlate with chromosome 14 instability previously observed in MN1-rearranged tumors4. Fusions containing GNAS (at 20q13) also only occurred in female cases (Supplementary Fig. 11A, Supplementary Data 2).
Fusions involving imprinted genes overall occurred at similar frequency in MN1-BEND2 (2.6/case, n = 5) and MAPK-ABC tumors (2.1/case, n = 15), but were almost twofold more frequent in BRAFV600E cases (4.2/case, n = 6), and approximately tenfold more frequent in females (2.7/case, n = 16) compared to males (0.25/case, n = 4). Mutations in imprinted genes were found at similar frequency in MNI-BEND2, BRAFV600E, and MAPK-ABC tumors (1.0, 0.7, and 1.7/case, respectively), and were most common in males (2.5/case).
Considering all genes in all AB-like tumors, fusions were 4.9-fold more frequent in female versus male patient tumors (Table 2). Fusions involving chromosome X were 5.3-times more common in tumors from females (P = 0.017) and occurred almost entirely in MN1-BEND2 and BRAFV600E-mutant tumors. Fusions involving chromosome 22 were 17-fold more common in female patients (P = 0.0035, Table 2). *Statistically significant p < 0.05.
X-inactivation escape genes are more widely expressed by and mutant in MAPK supercluster AB-like tumors
We also compiled a list of 95 XIE genes and 105 genes exhibiting variable XIE status (Supplementary Data 6)56.
Overall, MAPK-ABC and other non-MN1-BEND2 AB-like tumors showed overexpression of more XIE genes than MN1-BEND2 tumors (Fig. 7c, Supplementary Figs. 12, 13). XIST, an XIE gene encoding a long noncoding RNA that regulates X-inactivation on the inactive X-chromosome58, was relatively overexpressed in MAPK-ABC ABs from female patients compared to other tumors and controls (Supplementary Figs. 1–3, 12, 13). Indeed, GSEA showed that XIE gene expression is significantly enriched in MAPK-ABC tumors versus MN1-rearranged tumors (ES = |0.31|, pFWER = 0.018) (Fig. 7e).
Three BRAFV600E-mutant MAPK tumors additionally contained fusions with XIST as the downstream gene (Supplementary Data 2). Mutations in XIE genes USP9X, NLGN4x, PLXNB3 (involved in axonal guidance and glioma invasion), SYAP1, SRPX2, and ZFX (ZFX-SMARCA1 fusion) occurred in non-MN1-BEND2, predominantly female cases. Some XIE genes were also overexpressed or mutant in MN1-BEND2 tumors, e.g., SHOX and ZFX.
Deleterious SNPs and other non-fusion mutations of XIE genes occurred slightly more frequently in MAPK-ABC tumors (5.0/case) and the BRAFV600E-mutant subgroup (5.3/case) compared to MN1-BEND2 cases (4.0/case), and in females (5.9/case) compared to males (3.0/case) (Supplementary Fig. 11B, Supplementary Data 2). However, fusions involving XIE genes occurred 3.6-fold more frequent in MN1-BEND2 tumors (4.0/case) than MAPK-ABC tumors (1.2/case) and 4.6-times as frequent in females compared to males (2.3 vs. 0.5/case, respectively).
Expression of AB-like tumor genes correlate with gene promoter methylation
In order to confirm that differential expression of AB-like tumor genes can be, in some cases, attributable to differential promoter methylation, we integrated RNAseq expression and DNA methylation data using the mCSEA R package59, which correlates differentially methylated regions (DMRs) with expression of nearby downstream genes. The expression of several MN1-BEND2 overexpressed genes within the IGF2-H19 and DIO3-DLK1 imprinted regions and other MN1-BEND2 and MAPK-ABC overexpressed genes significantly correlated with their promoter methylation status (Supplementary Fig. 14).
Discussion
MN1 coregulates transcription and chromatin remodeling. Its C-terminal truncation is associated with a distinct craniofacial/brain malformation syndrome60,61. MN1 is implicated as a tumor suppressor in meningioma and an oncoprotein in acute myeloid leukemia62,63,64. Alteration of these potential functions could be important in MN1-BEND2 and other MN1-rearranged tumors, e.g., containing MN1-CXXC5 fusions6,65.
BEND2’s function is unknown however BEN domain-containing homologs participate in chromatin remodeling and transcriptional repression66,67. As a putative chromatin remodeling protein, BEND2 dysregulation could potentially influence gene expression. We identified MN1-BEND2 fusion transcripts in five tumors from females, all over-expressing MN1 and BEND2 mRNA. BEND2 overexpression has also been documented in an additional MN1-BEND2 tumor68. An MN1-rearranged tumor from a male patient in our series, lacking MN1-BEND2 (or other MN1 or BEND2 fusions), did not overexpress MN1, but did overexpress BEND2. It clustered with MN1-BEND2 tumors by RNA expression, but not DNA methylation.
Of eight recently reported AB-like tumors demonstrating EWSR1-BEND2 fusions (Supplementary Data 8), one had documented high BEND2 expression69. Most were infratentorial, involved the medulla or spinal cord and some clustered near HGNET-MN1 by tSNE analysis of DNA methylation70. All but one in these anatomic regions occurred in males. Two frontal lobe cases were in females. EWSR1 is located at chromosome 22q12.2 near MN1, thus both translocations involve the same 22q12 general breakpoint region. BEND2 overexpression may thus be significant in the tumorigenesis of MN1-BEND2 and similar tumors. Indeed, BEND2 appears to be normally expressed in the brain only at very low levels (Fig. 5).
All reported tumors with confirmed MN1-BEND2, and available demographic data, occurred in females aged 4–18 years (mean 10.4 yr, n = 13) (Supplementary Data 8). Tumors occurred in the frontal, parietal or occipital lobe, predominantly the parietal lobe, but not the temporal lobe or infratentorial as seen in non-MN1-BEND2 AB-like tumors (Fig. 8a, Table 1, Supplementary Data 8). In contrast, the mean age of presentation of BRAFV600E-mutant AB tumors in our series was 24.7 years (n = 6), and 29.1 years for MAPK-ABC cases (n = 15). It, therefore, appears that AB-like tumors related to earlier developmental stages of NSCs (RG), i.e., MN1-BEND2 tumors, clinically present earlier in life compared to tumors transcriptionally related to later RG, e.g., MAPK-ABC ABs and PTPR17.
The transition of embryonic stem cells to NSCs is accompanied by changes in DNA methylation. Additional methylation changes allow expression of glia-specific genes during the mid-gestational shift from primarily neurogenesis to gliogenesis28,71,72. Further epigenetic changes occur during increased transcriptional activity and synaptic remodeling in childhood through young adulthood73,74. DNA methylation and gene expression patterns in MN1-BEND2 and MAPK-ABC tumors may reflect these developmental dynamics and partially explain their different ages of clinical presentation. Dysregulation or mutation of chromatin remodeling genes, e.g., MN1, BEND2, ATRX, CXXC5, CREBBP, and SMARCAD1; histone modifier genes, e.g., CHD3, KMT2C; and/or possibly X-inactivation effectors could alter DNA methylation and potentially contribute to AB tumorigenesis.
The terms early-, mid-, and late-RG have been used in model systems28,75. Ziller et al.28 defined mid-RG as corresponding to peak neurogenesis and early gliogenesis, and late-RG to mostly gliogenesis. Like vRG, early-RG expressed PAX6 and FOXJ1. Late-RG highly expressed OLIG2 and TCF4, similar to oRG and MAPK-ABC ABs.
MN1-BEND2 tumors exhibit high expression of several NEC and vRG markers and low expression of the oRG markers FAM107A, LIFR, and HOPX (Supplementary Fig. 5A, B, D). Unlike MN1-BEND2 tumors, MAPK-ABC tumors and tRG are enriched for CRYAB12 (Supplementary Fig. 5C). The latter is more highly expressed after 25 pcw, consistent with tRG expressing a later transcriptional program than vRG (Fig. 5).
Our data suggest that MN1-BEND2 tumors ultimately originate from early vRG-derived ependymal progenitors or a residual population of related cells. Such cells may become neoplastic due to epigenetic and/or genetic events in utero or during the dramatic growth of the brain during childhood. The correlation of MN1-BEND2 tumor transcriptomes with early NSC types (NEC and vRG) and MAPK-ABC tumor, PTPR, and PXA transcriptomes with later NSC types (tRG and oRG) suggests that mutagenic events in early brain development may manifest as tumors in both childhood and adulthood, or that some RG or RG-derived progenitors persist beyond adolescence.
NEC and vRG project an apical primary cilium into the lateral ventricle, facilitating IGF2 binding from the CSF76. Unlike adult rodent lateral ventricle NSCs (B1 astrocytes), NSCs in the adult rodent dentate gyrus express IGF2, which is required for their proliferation in an apparent autocrine fashion73,77. This regional difference in IGF2 expression was attributed to differential imprinting73. We observed IGF2 immunoreactivity in immature ependymal cells of the fetal mouse lateral ventricular zone. Thus, IGF2 produced by differentiating ependymal cells may possibly be a critical pro-proliferative signal in the fetal lateral ventricular zone. This seems plausible considering the requirement of the ependymal layer by B1 astrocytes78. Alternatively, this additional IGF2 may be needed by the developing ependyma themselves or other nearby cells.
MAPK supercluster tumors demonstrated high expression, mutation, and fusion of genes reflecting glial and neural differentiation. GSEA supports that these tumors express transcriptomes reflective of tRG and oRG. PXA was more highly enriched for oRG genes compared to tRG genes. The mixed pattern of tRG and oRG gene expression by MAPK-ABC and PXA might be explained by the possibility that gene expression of late RG types and/or their progeny IPs might vary with development over time, and that the published tRG and oRG gene sets represent developmental timeline-dependent snapshots.
MAPK-ABC tumors exhibit overexpression and mutation of genes involved in neurocognitive disorders characteristic of children and young adults. TCF4 is implicated in schizophrenia and Pitt-Hopkins autism-like syndrome79,80. ANK3 has been associated with autism, bipolar disorder, and schizophrenia81. KALRN, TRIO, and DISC1 misexpression or mutation are implicated in attention deficit hyperactivity disorder, autism, bipolar disorder, schizophrenia, and other conditions82,83,84,85. KALRN protein (kalirin) isoforms show increased expression in adolescence and early adulthood when schizophrenia and many MAPK supercluster AB-like tumors occur86,87.
Group C and other non-MN1-BEND2/non-BRAFV600E-mutant AB-like tumors exhibit AHNAK and SYNE1 mutations. Although AHNAK has been characterized as a tumor suppressor88, AHNAK and SYNE1 have also been implicated as cancer genes89. AHNAK is necessary for pseudopod formation, a key step in the epithelial to mesenchymal transition and tumor cell invasion44. SYNE1 encodes spectrin repeat-containing nuclear envelope family member 1, which is involved in nucleolemmal cytoskeletal anchoring23,90.
An intriguing aspect of MN1-BEND2 and MAPK-AB tumors is their striking female predominance of up to 9:1 or greater4. We previously observed chromosome X gains in seven AB-like tumors from female patients and only one in a male patient; X loss was seen in four tumors from females and none from males4. Here, fusions involving chromosome X were over five times more common in female patient tumors, and almost entirely in MN1-BEND2 and BRAFV600E-mutant tumors. Chromosome 22 fusions were even more frequent in tumors from females. Meiotic recombination of chromosome 22 is greater in females than males91. How this relates to non-homologous mitotic recombination is unclear.
High expression of imprinted genes in MN1-BEND2 tumors and XIE genes in MAPK-ABC tumors suggest a potential role for differential DNA methylation in driving their tumorigenesis. BEND2 itself exhibits variable XIE status15. XIE genes can be activated by epigenetic changes in chromatin structure and are frequently mutated in cancer, especially female predominant cancers92,93.
The imprinted homeobox gene DLX5 may escape imprinting and be biallelically expressed in mice and patients with the nearly entirely female disease Rett syndrome94,95. The latter is characterized by mutations in chromosome X gene MECP2 encoding methyl-CpG binding protein 2, which binds methylated DNA and suppresses transcription. We found frequent overexpression of imprinted genes DIO3, IGF2, H19, and others, correlating with decreased methylation of their promoters in MN1-BEND2 tumors. Our data suggest that epigenetic and genetic alterations associated with AB-like tumors possibly involves loss of imprinting, leading to increased expression and sometimes subsequent mutation of affected genes96,97. Loss of imprinting is linked to several cancers and pediatric syndromes, e.g., Beckwith-Wiedemann syndrome involving alteration of the IGF2 imprinting control region and increased risk for Wilms tumor and other neoplasms98. High expression of imprinted genes in MN1-BEND2 tumors could alternatively be due to retention of fetal methylation patterns. Interestingly, DNA methylation and expression of imprinted genes in the brain is not only developmentally time dependent, but is also anatomic region dependent99,100,101, which may possibly relate to anatomic locations in which AB-like tumors occur (Fig. 8a).
We provide multiple lines of evidence that MN1-BEND2 and MAPK-ABC astroblastomas are of radial glia lineage and likely derived from vRG, and oRG and/or tRG or their close progeny, respectively (Fig. 8b, c). This interpretation is limited however by the mostly correlative and non-mechanistic nature of the data. For instance, the mostly adult glioma, glioblastoma, is thought to originate from subventricular zone B1 astrocytes102. We cannot completely rule out that some AB-like tumors may be derived from B1 astrocytes, as the transcriptome of B1 astrocytes is largely undefined.
PTPR has been hypothesized to be derived from the subcomissural organ (SCO) located in the third ventricle. GSEA data strongly suggest that PTPR are closely related to tRG and thus derived from them. tRG resemble an incompletely defined RG-like ventricular zone “ependymoglial” cell known as the tanycyte, which may also project a single cilium into the ventricle and an elongated basal process to a subependymal capillary103. Such cells are notably found in various third ventricle locations including the SCO. Tanycytes may thus represent persistent tRG capable of generating PTPR. Notably, PTPR strongly express PAX6, which is expressed by tRG12. It has also been proposed that ABs may be derived from tanycytes104, potentially consistent with our findings of partial enrichment of tRG genes in MAPK-ABC tumors.
Recently, tumor stem cells expressing the oRG marker PTPRZ1 capable of MST were found in glioblastoma105. This supports the possibility that residual oRG or oRG-derived IP cells persist into adulthood and lead to tumor formation. Invasive properties of these oRG-like tumor stem cells105 may explain why oRG-related MAPK supercluster ABs have a poorer prognosis than MN1-BEND2 tumors4. It may also be that B1 astrocytes or other glioblastoma initiating cells can differentiate to attain oRG-like properties. The relevance of either possibility in adult gliomagenesis remains to be determined.
Because MN1-BEND2 containing AB-like tumors are generally associated with excellent patient survival4, they should not be termed high-grade neuroepithelial tumors6. Because MN1-BEND2 tumors are most likely derived from cells destined to become ependymal cells, they also should not be called astroblastomas, but rather early ependymal tumors with MN1-BEND2 fusion (EET MN1-BEND2). Of practical clinical diagnostic importance, we show that MN1-BEND2 tumors can be differentiated from other AB-like tumors by their strong diffuse MN1 nuclear immunoreactivity and absent diffuse GFAP staining. We recommend that MN1 and GFAP immunohistochemistry, DNA methylation analysis, and confirmation of the MN1-BEND2 fusion (and if negative, MN1-CXXC5 or EWSR1-BEND2 fusions) by FISH, RT-PCR, or sequencing comprise an integrated approach to the diagnosis of EET MN1-BEND2 and other AB-like tumors.
As oRG, and likely tRG, give rise to IPs destined to become neurons or glia, the coinage of “astroblastoma” by Bailey and Cushing106, indicating derivation from an astrocyte progenitor, was amazingly insightful. We propose that the term astroblastoma be reserved for MAPK supercluster AB-like lesions, particularly DNA methylation Group C tumors, which may exhibit classic astroblastoma histomorphology.
Lastly, our data identify the IGF2 pathway and the multidrug transporter ABCC1 as candidate therapeutic targets in EET MN1-BEND2107,108. While BRAF and downstream MEK and/or PI3K/AKT/mTOR pathway inhibitors may be useful in treating MAPK-ABC astroblastomas109,110.
Methods
Tumors
Construction of our clinical cohort of histologically-defined AB-like cases has been described2. The cohort was augmented with three additional AB cases with available material for further study (C39, 43–44) and an additional five tumors with AB-like histological features, for a total of 35 AB-like cases in the current cohort. All control brain tissue was retrospectively obtained from archived material from routine hospital autopsies. Equivalent fetal age control brain tissue was obtained from autopsies of spontaneous abortuses or deceased premature infants at Nationwide Children’s Hospital, Columbus, Ohio. Autopsies were requested by patient’s families or their physicians for medical reasons and performed with the consent of the patient’s next of kin, completely independent of the studies described herein. Informed consent for the use of archival biopsy or autopsy tissue was not obtained, nor required by the institutional review board. No patients or patient families were compensated for these studies. Institutional review board approval from the University of Louisville was obtained for use of control and case material (Approval number 17.0984). MN1-rearrangement was determined by FISH (cases C3, C7, C10, C13, C16, C21, and C295) and confirmed by CNV analysis using the Illumina EPIC BeadChip 850K microarray data from these and two additional samples (C43 and C44).
DNA methylation
Genomic DNA methylation analysis was performed on an additional fourteen tumors (C8, C15-18, C21, C27, C31, C35-36, C41-44) using the Illumina EPIC BeadChip 850K microarray similar to as previously described6,111,112. Genomic DNA was extracted from 4 μm sections of formalin-fixed paraffin-embedded (FFPE) tissues using the AllPrep DNA/RNA FFPE Kit (Qiagen, Hilden, Germany) and quantified using a Qubit High Sensitivity assay (Thermo Fisher Scientific, UK), before bisulfite conversion of up to 250 ng DNA, using the EZ DNA Methylation kit (Zymo Research, Irvine, CA). Bisulphite-converted FFPE DNA was processed with the Infinium FFPE DNA Restore kit (Illumina, UK) and assayed on the Infinium Methylation EPIC BeadChip (Illumina), according to the Infinium HD FFPE Methylation Assay automated protocol (Illumina).
The tSNE plot was generated using the 32,000 most variably methylated probes of the study samples and 2801 German Cancer Center (DKFZ) reference samples as previously described111. The later consists of 2729 tumors (comprising 76 pathological diagnoses) and 72 normal brain tissues.
RNA isolation
Total RNA was extracted from three 30 μm scrolls cut from FFPE tissue blocks containing 90% tumor or more, using the RecoverAll total nucleic acid isolation kit (Ambion/Thermo Fisher Scientific, Waltham, MA) per the manufacturer’s instructions. RNA was then passed through concentration and clean-up columns (Norgen Biotek, Thorold, ON) and quantified using a Nanodrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA).
RNAseq
Capture-based RNAseq was performed using a modified SureSelect strand-specific RNA library prep and target enrichment protocol (Agilent Technologies, Santa Clara, CA). Briefly, 200 ng of FFPE-derived total RNA was heated to 65 °C for 5 min and used for first-strand cDNA synthesis. Second-strand cDNA was synthesized, and the ends repaired. The 3ʹ ends of ds-cDNA were adenylated, ligated with adapters and fragments were PCR-amplified. PCR product (700 ng) was hybridized overnight to SureSelect XT Human All Exon V5 + lncRNA baits. The hybridized fragments were captured on magnetic beads and PCR-amplified using uniquely barcoded index primers. Barcoded libraries were pooled and underwent paired-end 150 bp sequencing on an Illumina HiSeq4000 instrument (Illumina, San Diego, CA).
For gene expression, read alignment was performed with hisat2-2.0.5113 to human genome GRCh38 and aligned files were converted to bam and sorted with samtools v1.3.1114. Alignment QC was generated with RSeQC v2.6.4115. Counts from the subread package v1.5.1116 were used to count the number of reads aligning to Ensembl IDs according to GRCh38.86 annotation.
Fusion detection was performed using JAFFA117 after merging across replicate lanes. Variant calling was performed following GATK’s best practices for variant calling in RNAseq using a 2-pass STAR alignment. Alignments were made with STAR v2.5.2a118 to human reference genome GRCh38. The junctions detected in all samples were merged and used to make a custom STAR genome index for a second pass alignment. Picard v2.4.1 was used to preprocess the alignment files, add read groups, mark duplicates, sort, and index. GATK v3.6119 was used to split reads containing junctions, realign around indels, and recalibrate the base qualities. Samtools v1.3.1 was used to generate mpilup format. BRAF mutations were identified by a visual inspection of the alignment BAM in the Integrated Genomic Viewer120 since they were not detected by VarScan analysis. VarScan v2.4.1121 was used to call and filter variants from mpilup format from the whole genome. Called variants were annotated with Ensembl Variant Effect Predictor (VEP)122. Variants were eliminated if listed in dbSNP alone. Variants with MEDIUM or HIGH impact as designated by VEP were retained. After filtering variants detected by the same set of reads were eliminated. PolyPhen, SIFT, and FathmmMKL predictions were determined using respective versions available July 2018.
cDNA microarrays
Affymetrix GeneChip Human Transcriptome Array 2.0 (HTA-2) microarray analysis was performed according to the manufacturer’s instructions using the GeneChip WT Pico and Sensation Plus sample preparation kits (Thermo Fisher Scientific). In brief, 50 ng of total RNA was converted to single-stranded cDNA with T7 promoter sequence at the 5′ end. A 3′ Adaptor, as a template for double-stranded cDNA synthesis, was added and the single-stranded cDNA was converted to double-stranded cDNA. Six amplification cycles were applied as recommended by the protocol to enrich the samples. The amplified product was converted to cRNA and purified using purification binding beads. Eluted cRNA yield was measured by NanoDrop. cRNA (20 μg) were reverse transcribed to sense-strand cDNA containing dUTP at a fixed ratio relative to dTTP. The cRNA was removed using RNase H and the 2nd-cycle single-stranded cDNA was purified using purification binding beads. The eluted single-stranded cDNA yield was determined using the NanoDrop and 5.5 μg was fragmented by uracil-DNA glycosylase and apurinic/apyrimidinic endonuclease 1 (APE 1). The fragmented ss-cDNA was labeled by terminal deoxynucleotidyl transferase, and then used to prepare the Hybridization Cocktail, injected and hybridized in the array for 16 h at 45 °C in an oven with rotation at 60 rpm. After overnight hybridization, the hybridization cocktail mix was extracted from the array. Each array was washed and stained using Affymetrix GeneChip Fluidics Station 450 according to the user manual (P/N 08-0295) and the arrays were scanned using GeneChip Scanner 7G. Each array image was collected using the GeneChip Command Console, QC, saved as.cel files and saved for further analysis.
RNA expression analysis, dimensionality reduction, and clustering
For the RNAseq data, transcripts with 0 counts > 50% were removed. Data were normalized with Limma-voom R package with quantile normalization123. The Ensembl IDs were mapped to unique gene symbols using the Bioconductor org.Hs.eg.db package (ver: 3.13.0) using median summarization.
For the microarray data, probe-level data were generated by analyzing the raw signal intensities from the raw Affymetrix CEL data files for each probe set using the rma algorithm for background intensity correction and normalization in the oligo package (ver: 1.56). Batch effect related to the sample preparation kit was reduced using the combat function of the sva package (ver: 3.4). The probes were mapped to unique gene symbols using the hta20transcriptcluster.db package (ver: 8.8) with median summarization.
Differential gene expression between groups of samples in RNAseq and Affymetrix datasets was performed using analysis of variance thresholding with a two-sided p-value to 0.05. Expression values of the genes differentially expressed in both datasets were scaled and denoted on a heatmap with the respective samples. Heatmaps in Figs. 1, 5, 8, S1–S3, S7, S8, and S11, S12 were generated by unsupervised hierarchical clustering.
Dimensionality reduction of the expression data with t-Distributed Stochastic Neighbor Embedding (tSNE) was performed to facilitate visualization using the following parameters (dimensions = 2, perplexity = 4, iterations = 5000) using the Rtsne package (ver: 0.15). All analyses were performed using R software (ver: 4.1).
Gene set enrichment analysis (GSEA)
GSEA was performed using GSEA software from the Broad Institute, Cambridge, MA (ver: 4.1)124,125 using the following settings (permutations: 1000; permutation type: gene_set; enrichment statistic: weighted; metric for ranking genes: t Test).
Differential gene promoter methylation and expression correlation
Raw methylation files were read and normalized using the sesame package (ver: 1.10.4) with the default settings. Differentially methylated gene promoters were identified using the mCSEA package (ver: 1.12). Correlation between methylation of gene promoters and expression levels was quantified using the mCSEAIntegrate function. This correlates differentially methylated regions (DMRs) with the expression of nearby downstream genes59.
Nanostring MicroRNA
MicroRNA (miR) profiling of tumor samples was performed using the Nanostring nCounter system (Nanostring Technologies, Seattle, WA)126. One hundred nanograms of total RNA was annealed with multiplexed DNA tags (miR-tags) complementary to specific target miR sequences. Mature miRs were bound to specific miR-tags using a ligase and excess tags were enzymatically removed. Tagged miRs were diluted 1–5, and 5 µl was combined with 5 µl of capture probes and 20 µl of reporter probes in hybridization buffer. Probes were allowed to complex with specific target sequences overnight (16–20 h) at 65 °C. Excess probe was removed using a two-step magnetic bead-based purification (nCounter Prep Station) and target/probe complexes were immobilized onto cartridges. Reporter probe counts were obtained for each sample using the Gen1 nCounter Digital Analyzer, which utilizes a CCD camera-equipped microscope to image the immobilized fluorescent reporters via a high-density scan encompassing 600 fields of view. Raw data output was imported into nSolver for normalization (http://www.nanostring.com/products/nSolver). For analysis, data are the geometric mean normalized and then log10. For heatmap generation, data were scaled and truncated as follows. For scaling, the data are mean centered and divided by the SD of the miR. Truncation was done at a value of 4.0 to adjust for rare very large outlier values in the matrix.
Immunohistochemistry
Immunohistochemistry was done using 4 µ thick sections of FPPE tissue. Immunohistochemical stains for GFAP clone EP672Y (Cell Marque 258R-16, lot 32653, 1:25), MN1 (Proteintech 24697-1-AP, lot 00021048, 1:40) and IGF2 (Abcam; ab9574, lots GR31975-61, and GR31975-64, 1:500) were performed on the Leica Bond III using Leica Epitope Retrieval 1 (Leica; AR9961) for GFAP, and Epitope Retrieval 2 (Leica; AR9640) for MN1 and IGF2. Retrieval was performed for 20 min, followed by 15 min primary antibody incubation, and Bond Polymer Refine Detection (Leica; DS9800). Hematoxylin was used as a counterstain. Positive IGF2 staining was defined as dark granular perinuclear cytoplasmic staining.
Western blotting
Protein lysates of frozen tumor tissue were prepared as previously described127, adding to the lysis buffer 1 mM sodium orthovanadate and 5 mM sodium fluoride. Fifteen micrograms of protein were resolved on 10% bis-tris gels (Invitrogen) and transferred to PVDF membranes. Primary antibodies, BEND2 (Invitrogen PA5-31747, lot UA2709054B, 1:500), IGF2 (Abcam ab262713, lot GR3295940-1, 1:500), MAP3K5 (Millipore MABC632, lot VP1812286, 1:1000), ABCC1 (Abcam ab24102, lot GR3247402-6, 1:50), FOXJ1 (Novus NBP1-87928, lot H119212, 1:500), MN1 (Invitrogen PA5-38666, lot UD2755192, 1:500), TCF4 (Proteintech 22337-1-AP, lot 00050018, 1:500), and β-actin (Sigma A2228, lot 085M4754V, 1:10,000), were incubated for 90 min at room temperature, followed by a 60 min incubation with HRP-conjugated secondary antibodies (Cell Signaling, 1:2000). Uncropped and unprocessed scans of the most important blots are supplied in the Source Data file.
Animals
Animal studies were performed in accordance with the guidelines of the European Community and French Ministry of Agriculture and were approved by the “Direction départementale de la protection des populations - Paris (authorization 9343201702211706561). Mice used in this study are Centrin2-GFP (CB6-Tg(CAG-EGFP/CETN2)3- 4Jgg/J; The Jackson Laboratory). Animals were group-housed with free access to food and water in controlled temperature conditions (room temperature controlled at 21–22 °C, humidity between 40 and 50%), and exposed to a conventional 12-h light/dark cycle. Experiments were performed on embryos at embryonic (E) days E16, and neonate pups at postnatal day 0–10 (P0–10) of both sexes. All procedures were approved by the Ethical Committee CEEA-005 Charles Darwin (authorization 9343-201702211706561) and conducted in accordance with EU Directive 2010/63/EU. All efforts were made to reduce animal suffering and minimize the number of animals, in compliance with all relevant ethical regulations for animal testing and research.
Immunofluorescence
Animal brains were dissected out in cold PBS from embryos or newborn animals euthanized by decapitation. Lateral ventricle whole-mounts were prepared as previously described9. Whole-mount brain sections were fixed for 30 min in 4% PFA, incubated for 1 h in blocking solution (1X PBS with 0.1% Triton X-100 and 10% fetal bovine serum) at room temperature and were incubated overnight at 4 °C in primary antibodies diluted in blocking solution. The primary antibodies used targeted polyglutamylated tubulin (Adipogen AG-20B-0020-C100, 1:1000), MN1 (Proteintech 24697-1-AP, lot 00021048,1:100), or IGF-2 (Abcam ab9574, 1:100). The following day, they were stained with species-specific AlexaFluor fluorophore-conjugated secondary antibodies for 2 h at room temperature (1:400, Thermo Fischer Scientific or Jackson ImmunoResearch Labs). Finally, the whole-mounts were mounted with Fluoromount-G mounting medium (Southern Biotech, 0100-01).
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
The raw and processed RNAseq, Affymetrix, MethylationEpic BeadChip, and Nanostring microRNA data generated in this study are deposited at the public repository Gene Expression Omnibus (GEO) under the accession numbers GSE165351, GSE165813, GSE166569, and GSE196697, respectively. The remaining data are available within the Article, Source Data File, Supplementary Data, and Supplementary Information. Source data are provided with this paper.
Change history
31 January 2023
A Correction to this paper has been published: https://doi.org/10.1038/s41467-023-36255-z
References
Bailey, P. & Bucy, P. C. Astroblastomas of the brain. Acta Psychiatr. Scand. 5, 439–461 (1930).
Lehman, N. L. et al. Morphological and molecular features of astroblastoma, including BRAFV600E mutations, suggest an ontological relationship to other cortical-based gliomas of children and young adults. Neuro Oncol. 19, 31–42 (2017).
Aldape K. D., Rosenblum M. K. Astroblastoma. In: WHO Classification of Tumours of the Central Nervous System (eds Louis, D. N., Ohgaki, H., Wiestler, O. D. & Cavenee, W. K.). International Agency for Research on Cancer (2016).
Lehman, N. L. et al. Genomic analysis demonstrates that histologically-defined astroblastomas are molecularly heterogeneous and that tumors with MN1 rearrangement exhibit the most favorable prognosis. Acta Neuropathol. Commun. 7, 42 (2019).
Mhatre, R. et al. MN1 rearrangement in astroblastoma: study of eight cases and review of literature. Brain Tumor Pathol. 36, 112–120 (2019).
Sturm, D. et al. New brain tumor entities emerge from molecular classification of CNS-PNETs. Cell 164, 1060–1072 (2016).
Beattie, R. & Hippenmeyer, S. Mechanisms of radial glia progenitor cell lineage progression. FEBS Lett. 591, 3993–4008 (2017).
Silbereis, J. C., Pochareddy, S., Zhu, Y., Li, M. & Sestan, N. The cellular and molecular landscapes of the developing human central nervous system. Neuron 89, 248–268 (2016).
Ortiz-Alvarez, G. et al. Adult neural stem cells and multiciliated ependymal cells share a common lineage regulated by the geminin family members. Neuron 102, 159–172 e157 (2019).
Pollen, A. A. et al. Molecular identity of human outer radial glia during cortical development. Cell 163, 55–67 (2015).
Hansen, D. V., Lui, J. H., Parker, P. R. & Kriegstein, A. R. Neurogenic radial glia in the outer subventricular zone of human neocortex. Nature 464, 554–561 (2010).
Nowakowski, T. J., Pollen, A. A., Sandoval-Espinosa, C. & Kriegstein, A. R. Transformation of the radial glia scaffold demarcates two stages of human cerebral cortex development. Neuron 91, 1219–1227 (2016).
Carrel, L. & Willard, H. F. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature 434, 400–404 (2005).
Hassan H. & Hegyi H. Pervasive chromatin remodeling at X-inactivation escape genes in schizophrenic males. bioRxiv, 300624 (2018).
Tukiainen, T. et al. Landscape of X chromosome inactivation across human tissues. Nature 550, 244–248 (2017).
Cotton, A. M. et al. Analysis of expressed SNPs identifies variable extents of expression from the human inactive X chromosome. Genome Biol. 14, R122 (2013).
Heim, S. et al. Papillary tumor of the pineal region: a distinct molecular entity. Brain Pathol. 26, 199–205 (2016).
Korshunov, A. et al. Gene expression patterns in ependymomas correlate with tumor location, grade, and patient age. Am. J. Pathol. 163, 1721–1727 (2003).
Parker, M. et al. ZFTA-RELA fusions drive oncogenic NF-kappaB signalling in ependymoma. Nature 506, 451–455 (2014).
Taylor, M. D. et al. Radial glia cells are candidate stem cells of ependymoma. Cancer Cell 8, 323–335 (2005).
Robinson, A. et al. Mutations in the planar cell polarity genes CELSR1 and SCRIB are associated with the severe neural tube defect craniorachischisis. Hum. Mutat. 33, 440–447 (2012).
Boutin, C. et al. A dual role for planar cell polarity genes in ciliated cells. Proc. Natl Acad. Sci. USA 111, E3129–E3138 (2014).
Zhang, X. et al. SUN1/2 and Syne/Nesprin-1/2 complexes connect centrosome to the nucleus during neurogenesis and neuronal migration in mice. Neuron 64, 173–187 (2009).
Jacquet, B. V. et al. FoxJ1-dependent gene expression is required for differentiation of radial glia into ependymal cells and a subset of astrocytes in the postnatal brain. Development 136, 4021–4031 (2009).
Coy, S. et al. Nuclear CRX and FOXJ1 expression differentiates non-germ cell pineal region tumors and supports the ependymal differentiation of papillary tumor of the pineal region. Am. J. Surg. Pathol. 41, 1410–1421 (2017).
El Zein, L. et al. RFX3 governs growth and beating efficiency of motile cilia in mouse and controls the expression of genes involved in human ciliopathies. J. Cell Sci. 122, 3180–3189 (2009).
Pajtler, K. W. et al. Molecular classification of ependymal tumors across all CNS compartments, histopathological grades, and age groups. Cancer Cell 27, 728–743 (2015).
Ziller, M. J. et al. Dissecting neural differentiation regulatory networks through epigenetic footprinting. Nature 518, 355–359 (2015).
Zilberberg, A., Yaniv, A. & Gazit, A. The low density lipoprotein receptor-1, LRP1, interacts with the human frizzled-1 (HFz1) and down-regulates the canonical Wnt signaling pathway. J. Biol. Chem. 279, 17535–17542 (2004).
Ye, F. et al. DISC1 Regulates neurogenesis via modulating kinetochore attachment of Ndel1/Nde1 during mitosis. Neuron 96, 1204 (2017).
Mao, Y. et al. Disrupted in schizophrenia 1 regulates neuronal progenitor proliferation via modulation of GSK3beta/beta-catenin signaling. Cell 136, 1017–1031 (2009).
Shihab, H. A. et al. An integrative approach to predicting the functional effects of non-coding and coding sequence variation. Bioinformatics 31, 1536–1543 (2015).
Cai, X. et al. Tripartite motif-containing protein 27 negatively regulates CD4 T cells by ubiquitinating and inhibiting the class II PI3K-C2beta. Proc. Natl Acad. Sci. USA 108, 20072–20077 (2011).
Hall, A. & Lalli, G. Rho and Ras GTPases in axon growth, guidance, and branching. Cold Spring Harb. Perspect. Biol. 2, a001818 (2010).
Penzes, P. et al. Rapid induction of dendritic spine morphogenesis by trans-synaptic ephrinB-EphB receptor activation of the Rho-GEF kalirin. Neuron 37, 263–274 (2003).
Afroz, S., Parato, J., Shen, H. & Smith, S. S. Synaptic pruning in the female hippocampus is triggered at puberty by extrasynaptic GABAA receptors on dendritic spines. Elife 5, e15106 (2016).
Briancon-Marjollet, A. et al. Trio mediates netrin-1-induced Rac1 activation in axon outgrowth and guidance. Mol. Cell Biol. 28, 2314–2323 (2008).
Gueneau, L. et al. KIAA1109 Variants are associated with a severe disorder of brain development and arthrogryposis. Am. J. Hum. Genet. 102, 116–132 (2018).
Asif, S., Fatima, R., Krc, R., Bennett, J. & Raza, S. Comparative proteogenomic characterization of glioblastoma. CNS Oncol. 8, CNS37 (2019).
Leussis, M. P., Madison, J. M. & Petryshen, T. L. Ankyrin 3: genetic association with bipolar disorder and relevance to disease pathophysiology. Biol. Mood Anxiety Disord. 2, 18 (2012).
Alicea, D., Perez, M., Maldonado, C., Dominicci-Cotto, C. & Marie, B. Cortactin is a regulator of activity-dependent synaptic plasticity controlled by wingless. J. Neurosci. 37, 2203–2215 (2017).
Beghein, E., Devriese, D., Van Hoey, E. & Gettemans, J. Cortactin and fascin-1 regulate extracellular vesicle release by controlling endosomal trafficking or invadopodia formation and function. Sci. Rep. 8, 15606 (2018).
Kussel-Andermann, P. et al. Vezatin, a novel transmembrane protein, bridges myosin VIIA to the cadherin-catenins complex. EMBO J. 19, 6020–6029 (2000).
Shankar, J. et al. Pseudopodial actin dynamics control epithelial-mesenchymal transition in metastatic cancer cells. Cancer Res. 70, 3780–3790 (2010).
Middeldorp, J. et al. GFAPdelta in radial glia and subventricular zone progenitors in the developing human cortex. Development 137, 313–321 (2010).
Allen Institute for Brain Science. Allen Human Brain Atlas. BrainSpan: Atlas of the Developing Human Brain—Developmental Transcriptome https://www.brainspan.org/rnaseq/search/index.html (accessed 7 July 2021) (2010).
Johnson, M. B. et al. Single-cell analysis reveals transcriptional heterogeneity of neural progenitors in human cortex. Nat. Neurosci. 18, 637–646 (2015).
Kang, H. J. et al. Spatio-temporal transcriptome of the human brain. Nature 478, 483–489 (2011).
Thomsen, E. R. et al. Fixed single-cell transcriptomic characterization of human radial glial diversity. Nat. Methods 13, 87–93 (2016).
Liu, J. et al. The primate-specific gene TMEM14B marks outer radial glia cells and promotes cortical expansion and folding. Cell Stem Cell 21, 635–649 e638 (2017).
Al Jord, A. et al. Centriole amplification by mother and daughter centrioles differs in multiciliated cells. Nature 516, 104–107 (2014).
Spassky, N. et al. Adult ependymal cells are postmitotic and are derived from radial glial cells during embryogenesis. J. Neurosci. 25, 10–18 (2005).
Morison, I. M., Paton, C. J. & Cleverley, S. D. The imprinted gene and parent-of-origin effect database. Nucleic Acids Res. 29, 275–276 (2001).
Morison I. M. Catalogue of parent of origin effects. http://igc.otago.ac.nz; accessed 15 Jan 2020 (2001).
Jirtle R. L. Imprinted genes: by species. https://www.geneimprint.com/site/genes-by-species; accessed 16 Aug 2021 (1997).
Gardiner, E. et al. Imprinted DLK1-DIO3 region of 14q32 defines a schizophrenia-associated miRNA signature in peripheral blood mononuclear cells. Mol. Psychiatry 17, 827–840 (2012).
Rong, H. et al. MicroRNA-134 plasma levels before and after treatment for bipolar mania. J. Psychiatr. Res. 45, 92–95 (2011).
Ross, M. T. et al. The DNA sequence of the human X chromosome. Nature 434, 325–337 (2005).
Martorell-Marugan, J., Gonzalez-Rumayor, V. & Carmona-Saez, P. mCSEA: detecting subtle differentially methylated regions. Bioinformatics 35, 3257–3262 (2019).
Miyake N. et al. Gain-of-function MN1 truncation variants cause a recognizable syndrome with craniofacial and brain abnormalities. Am. J. Hum. Genet. [Epub ahead of print] (2019).
Riedel, S. S. et al. SWI/SNF dysregulation through a prion-like domain causes AML. Blood 134, 2512–2512 (2019).
Heuser, M. et al. MN1 overexpression induces acute myeloid leukemia in mice and predicts ATRA resistance in patients with AML. Blood 110, 1639–1647 (2007).
Lekanne Deprez, R. H. et al. Cloning and characterization of MN1, a gene from chromosome 22q11, which is disrupted by a balanced translocation in a meningioma. Oncogene 10, 1521–1528 (1995).
Riedel, S. S. et al. Intrinsically disordered Meningioma-1 stabilizes the BAF complex to cause AML. Mol. Cell 81, 2332–2348 e2339 (2021).
Lake, J. A. et al. Targeted fusion analysis can aid in the classification and treatment of pediatric glioma, ependymoma, and glioneuronal tumors. Pediatr. Blood Cancer 67, e28028 (2020).
Dai, Q. et al. The BEN domain is a novel sequence-specific DNA-binding domain conserved in neural transcriptional repressors. Genes Dev. 27, 602–614 (2013).
Sathyan, K. M., Shen, Z., Tripathi, V., Prasanth, K. V. & Prasanth, S. G. A BEN-domain-containing protein associates with heterochromatin and represses transcription. J. Cell Sci. 124, 3149–3163 (2011).
Burford, A. et al. The ten-year evolutionary trajectory of a highly recurrent paediatric high grade neuroepithelial tumour with MN1:BEND2 fusion. Sci. Rep. 8, 1032 (2018).
Yamasaki K. et al. Spinal cord astroblastoma with an EWSR1-BEND2 fusion classified as a high-grade neuroepithelial tumour with MN1 alteration. Neuropathol. Appl. Neurobiol. 46, 190–193 (2020).
Lucas, C. G. et al. EWSR1-BEND2 fusion defines an epigenetically distinct subtype of astroblastoma. Acta Neuropathol. 143, 109–113 (2022).
Takizawa, T. et al. DNA methylation is a critical cell-intrinsic determinant of astrocyte differentiation in the fetal brain. Dev. Cell 1, 749–758 (2001).
Wang, Z., Tang, B., He, Y. & Jin, P. DNA methylation dynamics in neurogenesis. Epigenomics 8, 401–414 (2016).
Ferron, S. R. et al. Differential genomic imprinting regulates paracrine and autocrine roles of IGF2 in mouse adult neurogenesis. Nat. Commun. 6, 8265 (2015).
Arain, M. et al. Maturation of the adolescent brain. Neuropsychiatr. Dis. Treat. 9, 449–461 (2013).
Elkabetz, Y. et al. Human ES cell-derived neural rosettes reveal a functionally distinct early neural stem cell stage. Genes Dev. 22, 152–165 (2008).
Lehtinen, M. K. et al. The cerebrospinal fluid provides a proliferative niche for neural progenitor cells. Neuron 69, 893–905 (2011).
Bracko, O. et al. Gene expression profiling of neural stem cells and their neuronal progeny reveals IGF2 as a regulator of adult hippocampal neurogenesis. J. Neurosci. 32, 3376–3387 (2012).
Lim D. A. & Alvarez-Buylla A. The adult ventricular-subventricular zone (V-SVZ) and olfactory bulb (OB) neurogenesis. Cold Spring Harb. Perspect. Biol. 8, a018820 (2016).
Navarrete, K. et al. TCF4 (e2-2; ITF2): a schizophrenia-associated gene with pleiotropic effects on human disease. Am. J. Med. Genet. B Neuropsychiatr. Genet. 162B, 1–16 (2013).
Whalen, S. et al. Novel comprehensive diagnostic strategy in Pitt-Hopkins syndrome: clinical score and further delineation of the TCF4 mutational spectrum. Hum. Mutat. 33, 64–72 (2012).
Schizophrenia Psychiatric Genome-Wide Association Study Consortium. Genome-wide association study identifies five new schizophrenia loci. Nat. Genet. 43, 969–976 (2011).
Lesch, K. P. et al. Molecular genetics of adult ADHD: converging evidence from genome-wide association and extended pedigree linkage studies. J. Neural Transm. 115, 1573–1585 (2008).
Russell, T. A. et al. A schizophrenia-linked KALRN coding variant alters neuron morphology, protein function, and transcript stability. Biol. Psychiatry 83, 499–508 (2018).
Paskus, J. D., Herring, B. E. & Roche, K. W. Kalirin and Trio: RhoGEFs in synaptic transmission, plasticity, and complex brain disorders. Trends Neurosci. 43, 505–518 (2020).
de Bartolomeis, A., Latte, G., Tomasetti, C. & Iasevoli, F. Glutamatergic postsynaptic density protein dysfunctions in synaptic plasticity and dendritic spines morphology: relevance to schizophrenia and other behavioral disorders pathophysiology, and implications for novel therapeutic approaches. Mol. Neurobiol. 49, 484–511 (2014).
McPherson, C. E., Eipper, B. A. & Mains, R. E. Genomic organization and differential expression of Kalirin isoforms. Gene 284, 41–51 (2002).
Remmers, C., Sweet, R. A. & Penzes, P. Abnormal kalirin signaling in neuropsychiatric disorders. Brain Res. Bull. 103, 29–38 (2014).
Lee, I. H. et al. Ahnak functions as a tumor suppressor via modulation of TGFbeta/Smad signaling pathway. Oncogene 33, 4675–4684 (2014).
Cheng, F. et al. A gene gravity model for the evolution of cancer genomes: a study of 3,000 cancer genomes across 9 cancer types. PLoS Comput. Biol. 11, e1004497 (2015).
Fujita, I. et al. Endfoot regeneration restricts radial glial state and prevents translocation into the outer subventricular zone in early mammalian brain development. Nat. Cell Biol. 22, 26–37 (2020).
Tapper, W. J., Ke, X., Morton, N. E. & Collins, A. Recombination, interference and sequence: comparison of chromosomes 21 and 22. Ann. Hum. Genet. 66, 75–86 (2002).
Balaton, B. P. & Brown, C. J. Escape artists of the X chromosome. Trends Genet. 32, 348–359 (2016).
Jager, N. et al. Hypermutation of the inactive X chromosome is a frequent event in cancer. Cell 155, 567–581 (2013).
Horike, S., Cai, S., Miyano, M., Cheng, J. F. & Kohwi-Shigematsu, T. Loss of silent-chromatin looping and impaired imprinting of DLX5 in Rett syndrome. Nat. Genet. 37, 31–40 (2005).
Kimura, M. I. et al. Dlx5, the mouse homologue of the human-imprinted DLX5 gene, is biallelically expressed in the mouse brain. J. Hum. Genet. 49, 273–277 (2004).
Thorvaldsen, J. L., Duran, K. L. & Bartolomei, M. S. Deletion of the H19 differentially methylated domain results in loss of imprinted expression of H19 and Igf2. Genes Dev. 12, 3693–3702 (1998).
Leick, M. B., Shoff, C. J., Wang, E. C., Congress, J. L. & Gallicano, G. I. Loss of imprinting of IGF2 and the epigenetic progenitor model of cancer. Am. J. Stem Cells 1, 59–74 (2012).
Mussa, A. et al. Cancer risk in Beckwith-Wiedemann syndrome: a systematic review and meta-analysis outlining a novel (epi)genotype specific histotype targeted screening protocol. J. Pediatr. 176, 142–149 e141 (2016).
Ho-Shing, O. & Dulac, C. Influences of genomic imprinting on brain function and behavior. Curr. Opin. Behav. Sci. 25, 66–76 (2019).
Wilkinson, L. S., Davies, W. & Isles, A. R. Genomic imprinting effects on brain development and function. Nat. Rev. Neurosci. 8, 832–843 (2007).
Pidsley, R. et al. DNA methylation at the Igf2/H19 imprinting control region is associated with cerebellum mass in outbred mice. Mol. Brain 5, 42 (2012).
Lee, J. H. et al. Human glioblastoma arises from subventricular zone cells with low-level driver mutations. Nature 560, 243–247 (2018).
Reichenbach A., Bringmann A. Comparative anatomy of glial cells in mammals. In: Evolutionary Neuroscience (Second Edition) (ed. Kaas, J. H.). Academic Press (2020).
Rubinstein, L. J. & Herman, M. M. The astroblastoma and its possible cytogenic relationship to the tanycyte. An electron microscopic, immunohistochemical, tissue- and organ-culture study. Acta Neuropathol. 78, 472–483 (1989).
Bhaduri, A. et al. Outer radial glia-like cancer stem cells contribute to heterogeneity of glioblastoma. Cell Stem Cell 26, 48–63 e46 (2020).
Bailey P., Cushing H. A Classification of the Tumors of the Glioma Group on a Histogenetic Basis with a Correlated Study of Prognosis. JB Lippincott Company (1926).
Silbermann, K., Stefan, S. M., Elshawadfy, R., Namasivayam, V. & Wiese, M. Identification of thienopyrimidine scaffold as an inhibitor of the ABC transport protein ABCC1 (MRP1) and related transporters using a combined virtual screening approach. J. Med. Chem. 62, 4383–4400 (2019).
Gupta, P. et al. GSK1904529A, A potent IGF-IR inhibitor, reverses MRP1-mediated multidrug resistance. J. Cell Biochem. 118, 3260–3267 (2017).
Dienstmann, R., Rodon, J., Serra, V. & Tabernero, J. Picking the point of inhibition: a comparative review of PI3K/AKT/mTOR pathway inhibitors. Mol. Cancer Ther. 13, 1021–1031 (2014).
Usubalieva, A. et al. Primary meningeal pleomorphic xanthoastrocytoma with anaplastic features: a report of 2 cases, One with BRAF(V600E) mutation and clinical response to the BRAF inhibitor dabrafenib. J. Neuropathol. Exp. Neurol. 74, 960–969 (2015).
Capper, D. et al. DNA methylation-based classification of central nervous system tumours. Nature 555, 469–474 (2018).
Moran, S., Arribas, C. & Esteller, M. Validation of a DNA methylation microarray for 850,000 CpG sites of the human genome enriched in enhancer sequences. Epigenomics 8, 389–399 (2016).
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667 (2016).
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Wang, L., Wang, S. & Li, W. RSeQC: quality control of RNA-seq experiments. Bioinformatics 28, 2184–2185 (2012).
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).
Davidson, N. M., Majewski, I. J. & Oshlack, A. JAFFA: High sensitivity transcriptome-focused fusion gene detection. Genome Med. 7, 43 (2015).
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Koboldt, D. C. et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 22, 568–576 (2012).
Thorvaldsdottir, H., Robinson, J. T. & Mesirov, J. P. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief. Bioinform. 14, 178–192 (2013).
McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
McLaren, W. et al. The Ensembl variant effect predictor. Genome Biol. 17, 122 (2016).
Law, C. W., Chen, Y., Shi, W. & Smyth, G. K. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 15, R29 (2014).
Mootha, V. K. et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 34, 267–273 (2003).
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Geiss, G. K. et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat. Biotechnol. 26, 317–325 (2008).
Lehman, N. L. et al. Aurora A is differentially expressed in gliomas, is associated with patient survival in glioblastoma and is a potential chemotherapeutic target in gliomas. Cell Cycle 11, 489–502 (2012).
Stylianopoulou, F., Herbert, J., Soares, M. B. & Efstratiadis, A. Expression of the insulin-like growth factor II gene in the choroid plexus and the leptomeninges of the adult rat central nervous system. Proc. Natl Acad. Sci. USA 85, 141–145 (1988).
Acknowledgements
We thank: Elizabeth Rushing, University Hospital of Zurich, Switzerland; Mohanpal Dulai, Trillium Health Partners, Mississauga, Ontario; Michael Nowacki, Norton Hospital, Louisville KY; Roger McLendon, Duke University, Durham, NC; Matthew Schniederjan, Emory University, Atlanta, GA; Marta Couce, Case Western Reserve University, Cleveland, OH; Maria-Magdalena Georgescu, Louisiana State University Health Shreveport, Shreveport, LA; Cheryl Palmer, University of Utah, Salt Lake City, UT; Mousa Al Abbadi, University of Jordan, Amman, Jordan; and Werner Paulus, University of Münster, Münster, Germany for previously contributing cases, Eric Lipp, Duke University, Durham, NC, Paul Fisher, Stanford University, Palo Alto, CA, Kristin Huntoon, Ohio State University, Columbus, OH for clinical follow-up; James Van Brocklyn, Xiaohui Li, Sarah Warner, Pearlly Yan, Huabao Wang, Sushma Nagaraj and Kayla O’Donnell for technical assistance; Leslie Schier for artwork (Fig. 1a); and David Symer and Keiko Akagi, Ohio State University, Columbus, OH for assistance with RNAseq pilot data. Supported in part by NIH Grants K08NS45077 (NLL), RO1NS081125 (NLL), P30CA016058 (Ohio State University Comprehensive Cancer Center Genomics Shared Resource Core), F32 CA224962 (AMM), UL1TR002243 (Vanderbilt Institute for Clinical and Translational Research), 2018 Burroughs Wellcome Fund Physician-Scientist Institutional Award 1018894 (AMM), Vanderbilt-Ingram Cancer Center Voland Neuro-Oncology Award (AMM), and Vanderbilt Institute for Clinical and Translational Research Award, VR52534 (AMM). All relevant funding awarded to each author is described in these acknowledgements. Figures 1b and 8a–c were created with BioRender.com
Author information
Authors and Affiliations
Contributions
Conceptualization, N.L.L.; Methodology, N.L.L., N.S., A.M.M., S.C.; Validation, N.L.L, A.M.M.; Formal analysis, N.L.L., A.M.M., N.S., A.W., C.T., J.P.M., C.T.Z., M.S.; Investigation, N.L.L., N.S., K.J.A., A.M.M., A.W., P.F., Z.A., T.L., C.T., K.H.M., M.S., J.C., and R.X.C.; Writing—Original draft, N.L.L.; Writing—Review & editing, N.L.L, N.S., K.J.A., A.M.M., B.J.W., M.S., T.D.L., K.H.M.; Visualization, N.L.L., A.M.M., N.S., C.T.Z., J.P.M., V.G., M.S.; Funding acquisition, N.L.L. and A.M.M.; Resources, N.L.L., N.S., C.R.P., M.G., J.R.P, A.P., B.C.M., A.M.M., T.P., E.M.H., and R.L.J.; Data curation, N.L.L., A.M.M., A.H., A.U., T.D.L., M.S., C.T.Z., K.J.A., and B.J.W.; Supervision and project administration, N.L.L.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lehman, N.L., Spassky, N., Sak, M. et al. Astroblastomas exhibit radial glia stem cell lineages and differential expression of imprinted and X-inactivation escape genes. Nat Commun 13, 2083 (2022). https://doi.org/10.1038/s41467-022-29302-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-29302-8
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.