Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing.

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Intron dynamics reveal principles of gene regulation during the maternal-to-zygotic transition.

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here, we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in nonmodel organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Cell-level metadata are indispensable for documenting single-cell sequencing datasets.

Single-cell RNA sequencing (scRNA-seq) provides an unprecedented view of cellular diversity of biological systems. However, across the thousands of publications and datasets generated using this technology, we estimate that only a minority (<25%) of studies provide cell-level metadata information containing identified cell types and related findings of the published dataset. Metadata omission hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates. We encourage investigators, reviewers, journals, and data repositories to improve their standards and ensure proper documentation of these valuable datasets.

RNA-binding proteins regulate aldosterone homeostasis in human steroidogenic cells.

Angiotensin II (AngII) stimulates adrenocortical cells to produce aldosterone, a master regulator of blood pressure. Despite extensive characterization of the transcriptional and enzymatic control of adrenocortical steroidogenesis, there are still major gaps in the precise regulation of AII-induced gene expression kinetics. Specifically, we do not know the regulatory contribution of RNA-binding proteins (RBPs) and RNA decay, which can control the timing of stimulus-induced gene expression. To investigate this question, we performed a high-resolution RNA-seq time course of the AngII stimulation response and 4-thiouridine pulse labeling in a steroidogenic human cell line (H295R). We identified twelve temporally distinct gene expression responses that contained mRNA encoding proteins known to be important for various steps of aldosterone production, such as cAMP signaling components and steroidogenic enzymes. AngII response kinetics for many of these mRNAs revealed a coordinated increase in both synthesis and decay. These findings were validated in primary human adrenocortical cells stimulated ex vivo with AngII. Using a candidate screen, we identified a subset of RNA-binding protein and RNA decay factors that activate or repress AngII-stimulated aldosterone production. Among the repressors of aldosterone were BTG2, which promotes deadenylation and global RNA decay. BTG2 was induced in response to AngII stimulation and promoted the repression of mRNAs encoding pro-steroidogenic factors indicating the existence of an incoherent feedforward loop controlling aldosterone homeostasis. These data support a model in which coordinated increases in transcription and decay facilitate the major transcriptomic changes required to implement a pro-steroidogenic expression program that actively resolved to prevent aldosterone overproduction.

Fatal COVID-19 and Non-COVID-19 Acute Respiratory Distress Syndrome Is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State without Fibrosis.

Acute respiratory distress syndrome (ARDS) due to coronavirus disease 2019 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; recovery requires epithelial regeneration. During physiological regeneration in mice, type 2 AECs (AEC2s) proliferate, exit the cell cycle, transiently assume a transitional state, then differentiate into type 1 AECs (AEC1s); in humans, persistence of the transitional state is associated with pulmonary fibrosis. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiological regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or single-cell RNA sequencing revealed that transitional cells in mouse models of physiological regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis.

Simultaneous measurement of biochemical phenotypes and gene expression in single cells.

Methods to measure heterogeneity among cells are rapidly transforming our understanding of biology but are currently limited to molecular abundance measurements. We developed an approach to simultaneously measure biochemical activities and mRNA abundance in single cells to understand the heterogeneity of DNA repair activities across thousands of human lymphocytes, identifying known and novel cell-type-specific DNA repair phenotypes.

Recovery and analysis of transcriptome subsets from pooled single-cell RNA-seq libraries.

Single-cell RNA sequencing (scRNA-seq) methods generate sparse gene expression profiles for thousands of single cells in a single experiment. The information in these profiles is sufficient to classify cell types by distinct expression patterns but the high complexity of scRNA-seq libraries often prevents full characterization of transcriptomes from individual cells. To extract more focused gene expression information from scRNA-seq libraries, we developed a strategy to physically recover the DNA molecules comprising transcriptome subsets, enabling deeper interrogation of the isolated molecules by another round of DNA sequencing. We applied the method in cell-centric and gene-centric modes to isolate cDNA fragments from scRNA-seq libraries. First, we resampled the transcriptomes of rare, single megakaryocytes from a complex mixture of lymphocytes and analyzed them in a second round of DNA sequencing, yielding up to 20-fold greater sequencing depth per cell and increasing the number of genes detected per cell from a median of 1313 to 2002. We similarly isolated mRNAs from targeted T cells to improve the reconstruction of their VDJ-rearranged immune receptor mRNAs. Second, we isolated CD3D mRNA fragments expressed across cells in a scRNA-seq library prepared from a clonal T cell line, increasing the number of cells with detected CD3D expression from 59.7% to 100%. Transcriptome resampling is a general approach to recover targeted gene expression information from single-cell RNA sequencing libraries that enhances the utility of these costly experiments, and may be applicable to the targeted recovery of molecules from other single-cell assays.

Dynamic temperature-sensitive A-to-I RNA editing in the brain of a heterothermic mammal during hibernation.

RNA editing diversifies genomically encoded information to expand the complexity of the transcriptome. In ectothermic organisms, including Drosophila and Cephalopoda, where body temperature mirrors ambient temperature, decreases in environmental temperature lead to increases in A-to-I RNA editing and cause amino acid recoding events that are thought to be adaptive responses to temperature fluctuations. In contrast, endothermic mammals, including humans and mice, typically maintain a constant body temperature despite environmental changes. Here, A-to-I editing primarily targets repeat elements, rarely results in the recoding of amino acids, and plays a critical role in innate immune tolerance. Hibernating ground squirrels provide a unique opportunity to examine RNA editing in a heterothermic mammal whose body temperature varies over 30°C and can be maintained at 5°C for many days during torpor. We profiled the transcriptome in three brain regions at six physiological states to quantify RNA editing and determine whether cold-induced RNA editing modifies the transcriptome as a potential mechanism for neuroprotection at low temperature during hibernation. We identified 5165 A-to-I editing sites in 1205 genes with dynamically increased editing after prolonged cold exposure. The majority (99.6%) of the cold-increased editing sites are outside of previously annotated coding regions, 82.7% lie in SINE-derived repeats, and 12 sites are predicted to recode amino acids. Additionally, A-to-I editing frequencies increase with increasing cold-exposure, demonstrating that ADAR remains active during torpor. Our findings suggest that dynamic A-to-I editing at low body temperature may provide a neuroprotective mechanism to limit aberrant dsRNA accumulation during torpor in the mammalian hibernator.

Redistribution of EC-SOD resolves bleomycin-induced inflammation via increased apoptosis of recruited alveolar macrophages.

A human single nucleotide polymorphism (SNP) in the matrix-binding domain of extracellular superoxide dismutase (EC-SOD), with arginine to glycine substitution at position 213 (R213G), redistributes EC-SOD from the matrix into extracellular fluids. We reported that, following bleomycin (bleo), knockin mice harboring the human R213G SNP (R213G mice) exhibit enhanced resolution of inflammation and protection against fibrosis, compared with wild-type (WT) littermates. In this study, we tested the hypothesis that the EC-SOD R213G SNP promotes resolution via accelerated apoptosis of recruited alveolar macrophage (AM). RNA sequencing and Ingenuity Pathway Analysis 7 d postbleo in recruited AM implicated increased apoptosis and blunted inflammatory responses in the R213G strain exhibiting accelerated resolution. We validated that the percentage of apoptosis was significantly elevated in R213G recruited AM vs. WT at 3 and 7 d postbleo in vivo. Recruited AM numbers were also significantly decreased in R213G mice vs. WT at 3 and 7 d postbleo. ChaC glutathione-specific γ-glutamylcyclotransferase 1 (Chac1), a proapoptotic γ-glutamyl cyclotransferase that depletes glutathione, was increased in the R213G recruited AM. Overexpression of Chac1 in vitro induced apoptosis of macrophages and was blocked by administration of cell-permeable glutathione. In summary, we provide new evidence that redistributed EC-SOD accelerates the resolution of inflammation through redox-regulated mechanisms that increase recruited AM apoptosis.-Allawzi, A., McDermott, I., Delaney, C., Nguyen, K., Banimostafa, L., Trumpie, A., Hernandez-Lagunas, L., Riemondy, K., Gillen, A., Hesselberth, J., El Kasmi, K., Sucharov, C. C., Janssen, W. J., Stenmark, K., Bowler, R., Nozik-Grayck, E. Redistribution of EC-SOD resolves bleomycin-induced inflammation via increased apoptosis of recruited alveolar macrophages.

A basal-enriched microRNA is required for prostate tumorigenesis in a Pten knockout mouse model.

MicroRNAs (miRNAs) play important roles in prostate cancer development. However, it remains unclear how individual miRNAs contribute to the initiation and progression of prostate cancer. Here we show that a basal layer-enriched miRNA is required for prostate tumorigenesis. We identify miR-205 as the most highly expressed miRNA and enriched in the basal cells of the prostate. Although miR-205 is not required for normal prostate development and homeostasis, genetic deletion of miR-205 in a Pten null tumor model significantly compromises tumor progression and does not promote metastasis. In Pten null basal cells, loss of miR-205 attenuates pAkt levels and promotes cellular senescence. Furthermore, although overexpression of miR-205 in prostate cancer cells with luminal phenotypes inhibits cell growth in both human and mouse, miR-205 has a minimal effect on the growth of a normal human prostate cell line. Taken together, we have provided genetic evidence for a requirement of miR-205 in the progression of Pten null-induced prostate cancer.

Genome-wide analysis of Musashi-2 targets reveals novel functions in governing epithelial cell migration.

The Musashi-2 (Msi2) RNA-binding protein maintains stem cell self-renewal and promotes oncogenesis by enhancing cell proliferation in hematopoietic and gastrointestinal tissues. However, it is unclear how Msi2 recognizes and regulates mRNA targets in vivo and whether Msi2 primarily controls cell growth in all cell types. Here we identified Msi2 targets with HITS-CLIP and revealed that Msi2 primarily recognizes mRNA 3'UTRs at sites enriched in multiple copies of UAG motifs in epithelial progenitor cells. RNA-seq and ribosome profiling demonstrated that Msi2 promotes targeted mRNA decay without affecting translation efficiency. Unexpectedly, the most prominent Msi2 targets identified are key regulators that govern cell motility with a high enrichment in focal adhesion and extracellular matrix-receptor interaction, in addition to regulators of cell growth and survival. Loss of Msi2 stimulates epithelial cell migration, increases the number of focal adhesions and also compromises cell growth. These findings provide new insights into the molecular mechanisms of Msi2's recognition and repression of targets and uncover a key function of Msi2 in restricting epithelial cell migration.

Overexpression of pendrin in intercalated cells produces chloride-sensitive hypertension.

Inherited and acquired disorders that enhance the activity of transporters mediating renal tubular Na(+) reabsorption are well established causes of hypertension. It is unclear, however, whether primary activation of an Na(+)-independent chloride transporter in the kidney can also play a pathogenic role in this disease. Here, mice overexpressing the chloride transporter pendrin in intercalated cells of the distal nephron (Tg(B1-hPDS) mice) displayed increased renal absorption of chloride. Compared with normal mice, these transgenic mice exhibited a delayed increase in urinary NaCl and ultimately, developed hypertension when exposed to a high-salt diet. Administering the same sodium intake as NaHCO3 instead of NaCl did not significantly alter BP, indicating that the hypertension in the transgenic mice was chloride-sensitive. Moreover, excessive chloride absorption by pendrin drove parallel absorption of sodium through the epithelial sodium channel ENaC and the sodium-driven chloride/bicarbonate exchanger (Ndcbe), despite an appropriate downregulation of these sodium transporters in response to the expanded vascular volume and hypertension. In summary, chloride transport in the distal nephron can play a primary role in driving NaCl transport in this part of the kidney, and a primary abnormality in renal chloride transport can provoke arterial hypertension. Thus, we conclude that the chloride/bicarbonate exchanger pendrin plays a major role in controlling net NaCl absorption, thereby influencing BP under conditions of high salt intake.

Tumor and immune cell types interact to produce heterogeneous phenotypes of pediatric high-grade glioma.

BACKGROUND: Pediatric high-grade gliomas (PHGG) are aggressive brain tumors with 5-year survival rates ranging from <2% to 20% depending upon subtype. PHGG presents differently from patient to patient and is intratumorally heterogeneous, posing challenges in designing therapies. We hypothesized that heterogeneity occurs because PHGG comprises multiple distinct tumor and immune cell types in varying proportions, each of which may influence tumor characteristics. METHODS: We obtained 19 PHGG samples from our institution's pediatric brain tumor bank. We constructed a comprehensive transcriptomic dataset at the single-cell level using single-cell RNA-Seq (scRNA-Seq), identified known glial and immune cell types, and performed differential gene expression and gene set enrichment analysis. We conducted multi-channel immunofluorescence (IF) staining to confirm the transcriptomic results. RESULTS: Our PHGG samples included 3 principal predicted tumor cell types: astrocytes, oligodendrocyte progenitors (OPCs), and mesenchymal-like cells (Mes). These cell types differed in their gene expression profiles, pathway enrichment, and mesenchymal character. We identified a macrophage population enriched in mesenchymal and inflammatory gene expression as a possible source of mesenchymal tumor characteristics. We found evidence of T-cell exhaustion and suppression. CONCLUSIONS: PHGG comprises multiple distinct proliferating tumor cell types. Microglia-derived macrophages may drive mesenchymal gene expression in PHGG. The predicted Mes tumor cell population likely derives from OPCs. The variable tumor cell populations rely on different oncogenic pathways and are thus likely to vary in their responses to therapy.

Multi-pronged analysis of pediatric low-grade glioma reveals a unique tumor microenvironment associated with BRAF alterations.

Pediatric low-grade gliomas (pLGG) comprise 35% of all brain tumors. Despite favorable survival, patients experience significant morbidity from disease and treatments. A deeper understanding of pLGG biology is essential to identify novel, more effective, and less toxic therapies. We utilized single cell RNA sequencing (scRNA-seq), spatial transcriptomics, and cytokine analyses to characterize and understand tumor and immune cell heterogeneity across pLGG. scRNA-seq revealed tumor and immune cells within the tumor microenvironment (TME). Tumor cell subsets revealed a developmental hierarchy with progenitor and mature cell populations. Immune cells included myeloid and lymphocytic cells. There was a significant difference between the prevalence of two major myeloid subclusters between pilocytic astrocytoma (PA) and ganglioglioma (GG). Bulk and single-cell cytokine analyses evaluated the immune cell signaling cascade with distinct immune phenotypes among tumor samples. KIAA1549-BRAF tumors appeared more immunogenic, secreting higher levels of immune cell activators and chemokines, compared to BRAF V600E tumors. Spatial transcriptomics revealed the differential gene expression of these chemokines and their location within the TME. A multi-pronged analysis of pLGG demonstrated the complexity of the pLGG TME and differences between genetic drivers that may influence their response to immunotherapy. Further investigation of immune cell infiltration and tumor-immune interactions is warranted. Key points: There is a developmental hierarchy in neoplastic population comprising of both progenitor-like and mature cell types in both PA and GG.A more immunogenic, immune activating myeloid population is present in PA compared to GG. Functional analysis and spatial transcriptomics show higher levels of immune mobilizing chemokines in KIAA1549-BRAF fusion PA tumor samples compared to BRAF V600E GG samples. Importance of the Study: While scRNA seq provides information on cellular heterogeneity within the tumor microenvironment (TME), it does not provide a complete picture of how these cells are interacting or where they are located. To expand on this, we used a three-pronged approach to better understand the biology of pediatric low-grade glioma (pLGG). By analyzing scRNA-seq, secreted cytokines and spatial orientation of cells within the TME, we strove to gain a more complete picture of the complex interplay between tumor and immune cells within pLGG. Our data revealed a complex heterogeneity in tumor and immune populations and identified an interesting difference in the immune phenotype among different subtypes.

Spatial transcriptomic analysis delineates epithelial and mesenchymal subpopulations and transition stages in childhood ependymoma.

BACKGROUND: The diverse cellular constituents of childhood brain tumor ependymoma, recently revealed by single cell RNA-sequencing, may underly therapeutic resistance. Here we use spatial transcriptomics to further advance our understanding of the tumor microenvironment, mapping cellular subpopulations to the tumor architecture of ependymoma posterior fossa subgroup A (PFA), the commonest and most deadly childhood ependymoma variant. METHODS: Spatial transcriptomics data from intact PFA sections was deconvoluted to resolve the histological arrangement of neoplastic and non-neoplastic cell types. Key findings were validated using immunohistochemistry, in vitro functional assays and outcome analysis in clinically-annotated PFA bulk transcriptomic data. RESULTS: PFA are comprised of epithelial and mesenchymal histological zones containing a diversity of cellular states, each zone including co-existing and spatially distinct undifferentiated progenitor-like cells; a quiescent mesenchymal zone population, and a second highly mitotic progenitor population that is restricted to hypercellular epithelial zones and that is more abundant in progressive tumors. We show that myeloid cell interaction is the leading cause of mesenchymal transition in PFA, occurring in zones spatially distinct from hypoxia-induced mesenchymal transition, and these distinct EMT-initiating processes were replicated using in vitro models of PFA. CONCLUSIONS: These insights demonstrate the utility of spatial transcriptomics to advance our understanding of ependymoma biology, revealing a clearer picture of the cellular constituents of PFA, their interactions and influence on tumor progression.

Integration of single-nuclei RNA-sequencing, spatial transcriptomics and histochemistry defines the complex microenvironment of NF1-associated plexiform neurofibromas.

Plexiform neurofibroma (PN) is a leading cause of morbidity in children with the genetic condition Neurofibromatosis Type 1 (NF1), often disfiguring or threatening vital structures. During formation of PN, a complex tumor microenvironment (TME) develops, with recruitment of neoplastic and non-neoplastic cell types being critical for growth and progression. Due to the cohesive cellularity of PN, single-cell RNA-sequencing is difficult and may result in a loss of detection of critical cellular subpopulations. To bypass this barrier, we performed single-nuclei RNA-sequencing (snRNA-seq) on 8 frozen PN samples, and integrated this with spatial transcriptomics (ST) in 4 PN samples and immunohistochemistry to provide morphological context to transcriptomic data. SnRNA-seq analysis definitively charted the heterogeneous cellular subpopulations in the PN TME, with the predominant fraction being fibroblast subtypes. PN showed a remarkable amount of inter-sample homogeneity regarding cellular subpopulation proportions despite being resected from a variety of anatomical locations. ST analysis identified distinct cellular subpopulations which were annotated using snRNA-seq data and correlated with histological features. Schwann cell/fibroblast interactions were identified by receptor/ligand interaction analysis demonstrating a high probability of Neurexin 1/Neuroligin 1 (NRXN1/NLGN1) receptor-ligand cross-talk predicted between fibroblasts and non-myelinated Schwann cells (NM-SC) and subtypes, respectively. We observed aberrant expression of NRXN1 and NLGN1 in our PN snRNA-seq data compared to a normal mouse sciatic nerve single-cell RNA-seq dataset. This pathway has never been described in PN and may indicate a clear and direct communication pathway between putative NM-SC cells of origin and surrounding fibroblasts, potentially driving disease progression. SnRNA-seq integrated with spatial transcriptomics advances our understanding of the complex cellular heterogeneity of PN TME and identify potential novel communication pathways that may drive disease progression, a finding that could provide translational therapy options for patients with these devastating tumors of childhood and early adulthood.

Multi-omic approach identifies hypoxic tumor-associated myeloid cells that drive immunobiology of high-risk pediatric ependymoma.

Ependymoma (EPN) is a devastating childhood brain tumor. Single-cell analyses have illustrated the cellular heterogeneity of EPN tumors, identifying multiple neoplastic cell states including a mesenchymal-differentiated subpopulation which characterizes the PFA1 subtype. Here, we characterize the EPN immune environment, in the context of both tumor subtypes and tumor cell subpopulations using single-cell sequencing (scRNAseq, n = 27), deconvolution of bulk tumor gene expression (n = 299), spatial proteomics (n = 54), and single-cell cytokine release assays (n = 12). We identify eight distinct myeloid-derived subpopulations from which a group of cells, termed hypoxia myeloid cells, demonstrate features of myeloid-derived suppressor cells, including IL6/STAT3 pathway activation and wound healing ontologies. In PFA tumors, hypoxia myeloid cells colocalize with mesenchymal-differentiated cells in necrotic and perivascular niches and secrete IL-8, which we hypothesize amplifies the EPN immunosuppressive microenvironment. This myeloid cell-driven immunosuppression will need to be targeted for immunotherapy to be effective in this difficult-to-cure childhood brain tumor.

Significant increase of high-risk chromosome 1q gain and 6q loss at recurrence in posterior fossa group A ependymoma: A multicenter study.

BACKGROUND: Ependymoma (EPN) posterior fossa group A (PFA) has the highest rate of recurrence and the worst prognosis of all EPN molecular groups. At relapse, it is typically incurable even with re-resection and re-irradiation. The biology of recurrent PFA remains largely unknown; however, the increasing use of surgery at first recurrence has now provided access to clinical samples to facilitate a better understanding of this. METHODS: In this large longitudinal international multicenter study, we examined matched samples of primary and recurrent disease from PFA patients to investigate the biology of recurrence. RESULTS: DNA methylome derived copy number variants (CNVs) revealed large-scale chromosome gains and losses at recurrence in PFA. CNV changes were dominated by chromosome 1q gain and/or 6q loss, both previously identified as high-risk factors in PFA, which were present in 23% at presentation but increased to 61% at first recurrence. Multivariate survival analyses of this cohort showed that cases with 1q gain or 6q loss at first recurrence were significantly more likely to recur again. Predisposition to 1q+/6q- CNV changes at recurrence correlated with hypomethylation of heterochromatin-associated DNA at presentation. Cellular and molecular analyses revealed that 1q+/6q- PFA had significantly higher proportions of proliferative neuroepithelial undifferentiated progenitors and decreased differentiated neoplastic subpopulations. CONCLUSIONS: This study provides clinically and preclinically actionable insights into the biology of PFA recurrence. The hypomethylation predisposition signature in PFA is a potential risk-classifier for trial stratification. We show that the cellular heterogeneity of PFAs evolves largely because of genetic evolution of neoplastic cells.

Identification of HDAC4 as a potential therapeutic target and prognostic biomarker for ZFTA-fused ependymomas.

Members of the HDAC family are predictive biomarkers and regulate the tumorigenesis in several cancers. However, the role of these genes in the biology of intracranial ependymomas (EPNs) remains unexplored. Here, an analysis of eighteen HDACs genes in an EPN transcriptomic dataset, revealed significantly higher levels of HDAC4 in supratentorial ZFTA fusion (ST-ZFTA) compared with ST-YAP1 fusion and posterior fossa EPNs, while HDAC7 and SIRT2 were downregulated in ST-ZFTA. HDAC4 was also overexpressed in ST-ZFTA as measured by single-cell RNA-Seq, quantitative real time-polymerase chain reaction, and immunohistochemistry. Survival analyses showed a significantly worse outcome for EPNs with higher HDAC4 and SIRT1 mRNA levels. Ontology enrichment analysis showed an HDAC4-high signature consistent with viral processes while collagen-containing extracellular matrix and cell-cell junction were enriched in those with an HDAC4-low signature. Immune gene analysis demonstrated a correlation between HDAC4 expression and low levels of NK resting cells. Several small molecules compounds targeting HDAC4 and ABCG2, were predicted by in silico analysis to be effective against HDAC4-high ZFTA. Our results provide novel insights into the biology of the HDAC family in intracranial ependymomas and reveal HDAC4 as a prognostic marker and potential therapeutic target in ST-ZFTA.