Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Nuclear mTOR acts as a transcriptional integrator of the androgen signaling pathway in prostate cancer.

Androgen receptor (AR) signaling reprograms cellular metabolism to support prostate cancer (PCa) growth and survival. Another key regulator of cellular metabolism is mTOR, a kinase found in diverse protein complexes and cellular localizations, including the nucleus. However, whether nuclear mTOR plays a role in PCa progression and participates in direct transcriptional cross-talk with the AR is unknown. Here, via the intersection of gene expression, genomic, and metabolic studies, we reveal the existence of a nuclear mTOR-AR transcriptional axis integral to the metabolic rewiring of PCa cells. Androgens reprogram mTOR-chromatin associations in an AR-dependent manner in which activation of mTOR-dependent metabolic gene networks is essential for androgen-induced aerobic glycolysis and mitochondrial respiration. In models of castration-resistant PCa cells, mTOR was capable of transcriptionally regulating metabolic gene programs in the absence of androgens, highlighting a potential novel castration resistance mechanism to sustain cell metabolism even without a functional AR. Remarkably, we demonstrate that increased mTOR nuclear localization is indicative of poor prognosis in patients, with the highest levels detected in castration-resistant PCa tumors and metastases. Identification of a functional mTOR targeted multigene signature robustly discriminates between normal prostate tissues, primary tumors, and hormone refractory metastatic samples but is also predictive of cancer recurrence. This study thus underscores a paradigm shift from AR to nuclear mTOR as being the master transcriptional regulator of metabolism in PCa.

Epigenome-wide analysis across the development span of pediatric acute lymphoblastic leukemia: backtracking to birth.

BACKGROUND: Cancer is the leading cause of disease-related mortality in children. Causes of leukemia, the most common form, are largely unknown. Growing evidence points to an origin in-utero, when global redistribution of DNA methylation occurs driving tissue differentiation. METHODS: Epigenome-wide DNA methylation was profiled in surrogate (blood) and target (bone marrow) tissues at birth, diagnosis, remission and relapse of pediatric pre-B acute lymphoblastic leukemia (pre-B ALL) patients. Double-blinded analyses was performed between prospective cohorts extending from birth to diagnosis and retrospective studies backtracking from clinical disease to birth. Validation was carried out using independent technologies and populations. RESULTS: The imprinted and immuno-modulating VTRNA2-1 was hypermethylated (FDR<0.05) at birth in nested cases relative to controls in all tested populations (totaling 317 cases and 483 controls), including European and Hispanic ancestries. VTRNA2-1 methylation was stable over follow-up years after birth and across surrogate, target and other tissues (n=5,023 tissues; 30 types). When profiled in leukemic tissues from two clinical cohorts (totaling 644 cases), VTRNA2-1 methylation exhibited higher levels at diagnosis relative to controls, it reset back to normal levels at remission, and then re-increased to above control levels at relapse. Hypermethylation was significantly associated with worse pre-B ALL patient survival and with reduced VTRNA2-1 expression (n=2,294 tissues; 26 types), supporting a functional and translational role for VTRNA2-1 methylation. CONCLUSION: This study provides proof-of-concept to detect at birth epigenetic precursors of pediatric pre-B ALL. These alterations were reproducible with different technologies, in three continents and in two ethnicities, and can offer biomarkers for early detection and prognosis as well as actionable targets for therapy.

Characterisation of FLT3 alterations in childhood acute lymphoblastic leukaemia.

BACKGROUND: Alterations of FLT3 are among the most common driver events in acute leukaemia with important clinical implications, since it allows patient classification into prognostic groups and the possibility of personalising therapy thanks to the availability of FLT3 inhibitors. Most of the knowledge on FLT3 implications comes from the study of acute myeloid leukaemia and so far, few studies have been performed in other leukaemias. METHODS: A comprehensive genomic (DNA-seq in 267 patients) and transcriptomic (RNA-seq in 160 patients) analysis of FLT3 in 342 childhood acute lymphoblastic leukaemia (ALL) patients was performed. Mutations were functionally characterised by in vitro experiments. RESULTS: Point mutations (PM) and internal tandem duplications (ITD) were detected in 4.3% and 2.7% of the patients, respectively. A new activating mutation of the TKD, G846D, conferred oncogenic properties and sorafenib resistance. Moreover, a novel alteration involving the circularisation of read-through transcripts (rt-circRNAs) was observed in 10% of the cases. Patients presenting FLT3 alterations exhibited higher levels of the receptor. In addition, patients with ZNF384- and MLL/KMT2A-rearranged ALL, as well as hyperdiploid subtype, overexpressed FLT3. DISCUSSION: Our results suggest that specific ALL subgroups may also benefit from a deeper understanding of the biology of FLT3 alterations and their clinical implications.

CircRNAome of Childhood Acute Lymphoblastic Leukemia: Deciphering Subtype-Specific Expression Profiles and Involvement in TCF3::PBX1 ALL.

Childhood B-cell acute lymphoblastic leukemia (B-ALL) is a heterogeneous disease comprising multiple molecular subgroups with subtype-specific expression profiles. Recently, a new type of ncRNA, termed circular RNA (circRNA), has emerged as a promising biomarker in cancer, but little is known about their role in childhood B-ALL. Here, through RNA-seq analysis in 105 childhood B-ALL patients comprising six genetic subtypes and seven B-cell controls from two independent cohorts we demonstrated that circRNAs properly stratified B-ALL subtypes. By differential expression analysis of each subtype vs. controls, 156 overexpressed and 134 underexpressed circRNAs were identified consistently in at least one subtype, most of them with subtype-specific expression. TCF3::PBX1 subtype was the one with the highest number of unique and overexpressed circRNAs, and the circRNA signature could effectively discriminate new patients with TCF3::PBX1 subtype from others. Our results indicated that NUDT21, an RNA-binding protein (RBP) involved in circRNA biogenesis, may contribute to this circRNA enrichment in TCF3::PBX1 ALL. Further functional characterization using the CRISPR-Cas13d system demonstrated that circBARD1, overexpressed in TCF3::PBX1 patients and regulated by NUDT21, might be involved in leukemogenesis through the activation of p38 via hsa-miR-153-5p. Our results suggest that circRNAs could play a role in the pathogenesis of childhood B-ALL.

Development and validation of five behavioral indices of flood adaptation.

BACKGROUND: In the current context of climate change, climate forecasts for the province of Quebec (Canada) are a lengthening of the thunderstorm season and an increase in episodes of intense precipitations. These changes in the distribution of precipitations could heighten the intensity or frequency of floods, a natural hazard that concerns 80% of Quebec's riverside municipalities. For the health and safety of the at-risk population, it is very important to make sure they have acquired necessary adaptive behaviors against flooding hazard. However, there has been no assessment of these flood adaptation behaviors to date. Thus, the aim of this study was to develop and validate five indices of adaptation to flooding. METHODS: A sample of 1951 adults completed a questionnaire by phone. The questionnaire, specifically developed for this study, measured whether they did or did not adopt the behaviors that are proposed by public health officials to protect themselves against flooding. RESULTS: The results of the item, confirmatory factor, and multiple correspondence analyses contributed to the development of five indices corresponding to the adaptation behaviors to adopt according to the chronology of events: (a) pre-alert preventive behaviors, (b) behaviors to carry out after the alert is issued, (c) behaviors to adopt during a flood not requiring evacuation, (d) behaviors to adopt during a flood requiring evacuation, and (e) post-flood behaviors. The results of this study also showed that people who perceive a risk of flooding in their home in the next 5 years tend to adopt more preventive behaviors and adaptation behaviors than those who perceive little or no risk at all. They also reveal that people who feel more adverse effects on their physical or mental health tend to adopt more adaptive behaviors than those who feel little or no adverse effects on their health. CONCLUSION: Across a series of psychometric analyses, the results showed that these flood adaptation indices could properly measure a vast range of adaptive behaviors according to the chronology of events. Therefore, researchers, public health agencies, and professionals can use them to monitor the evolution of individuals' adaptive behaviors during floods.

Control of embryonic stem cell self-renewal and differentiation via coordinated alternative splicing and translation of YY2.

Translational control of gene expression plays a key role during the early phases of embryonic development. Here we describe a transcriptional regulator of mouse embryonic stem cells (mESCs), Yin-yang 2 (YY2), that is controlled by the translation inhibitors, Eukaryotic initiation factor 4E-binding proteins (4E-BPs). YY2 plays a critical role in regulating mESC functions through control of key pluripotency factors, including Octamer-binding protein 4 (Oct4) and Estrogen-related receptor-ß (Esrrb). Importantly, overexpression of YY2 directs the differentiation of mESCs into cardiovascular lineages. We show that the splicing regulator Polypyrimidine tract-binding protein 1 (PTBP1) promotes the retention of an intron in the 5'-UTR of Yy2 mRNA that confers sensitivity to 4E-BP-mediated translational suppression. Thus, we conclude that YY2 is a major regulator of mESC self-renewal and lineage commitment and document a multilayer regulatory mechanism that controls its expression.

Conserved expression of transposon-derived non-coding transcripts in primate stem cells.

BACKGROUND: A significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs). RESULTS: To identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved. CONCLUSIONS: Together, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.

High-dose folic acid supplementation alters the human sperm methylome and is influenced by the MTHFR C677T polymorphism.

Dietary folate is a major source of methyl groups required for DNA methylation, an epigenetic modification that is actively maintained and remodeled during spermatogenesis. While high-dose folic acid supplementation (up to 10 times the daily recommended dose) has been shown to improve sperm parameters in infertile men, the effects of supplementation on the sperm epigenome are unknown. To assess the impact of 6 months of high-dose folic acid supplementation on the sperm epigenome, we studied 30 men with idiopathic infertility. Blood folate concentrations increased significantly after supplementation with no significant improvements in sperm parameters. Methylation levels of the differentially methylated regions of several imprinted loci (H19, DLK1/GTL2, MEST, SNRPN, PLAGL1, KCNQ1OT1) were normal both before and after supplementation. Reduced representation bisulfite sequencing (RRBS) revealed a significant global loss of methylation across different regions of the sperm genome. The most marked loss of DNA methylation was found in sperm from patients homozygous for the methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism, a common polymorphism in a key enzyme required for folate metabolism. RRBS analysis also showed that most of the differentially methylated tiles were located in DNA repeats, low CpG-density and intergenic regions. Ingenuity Pathway Analysis revealed that methylation of promoter regions was altered in several genes involved in cancer and neurobehavioral disorders including CBFA2T3, PTPN6, COL18A1, ALDH2, UBE4B, ERBB2, GABRB3, CNTNAP4 and NIPA1. Our data reveal alterations of the human sperm epigenome associated with high-dose folic acid supplementation, effects that were exacerbated by a common polymorphism in MTHFR.

Transient DNMT1 suppression reveals hidden heritable marks in the genome.

Genome-wide demethylation and remethylation of DNA during early embryogenesis is essential for development. Imprinted germline differentially methylated domains (gDMDs) established by sex-specific methylation in either male or female germ cells, must escape these dynamic changes and sustain precise inheritance of both methylated and unmethylated parental alleles. To identify other, gDMD-like sequences with the same epigenetic inheritance properties, we used a modified embryonic stem (ES) cell line that emulates the early embryonic demethylation and remethylation waves. Transient DNMT1 suppression revealed gDMD-like sequences requiring continuous DNMT1 activity to sustain a highly methylated state. Remethylation of these sequences was also compromised in vivo in a mouse model of transient DNMT1 loss in the preimplantation embryo. These novel regions, possessing heritable epigenetic features similar to imprinted-gDMDs are required for normal physiological and developmental processes and when disrupted are associated with disorders such as cancer and autism spectrum disorders. This study presents new perspectives on DNA methylation heritability during early embryo development that extend beyond conventional imprinted-gDMDs.

Tissue factor expression provokes escape from tumor dormancy and leads to genomic alterations.

The coagulation system links immediate (hemostatic) and late (inflammatory, angiogenic) tissue responses to injury, a continuum that often is subverted in cancer. Here we provide evidence that tumor dormancy is influenced by tissue factor (TF), the cancer cell-associated initiator of the coagulation system and a signaling receptor. Thus, indolent human glioma cells deficient for TF remain viable but permanently dormant at the injection site for nearly a year, whereas the expression of TF leads to a step-wise transition to latent and overt tumor growth phases, a process that is preceded by recruitment of vascular (CD105(+)) and myeloid (CD11b(+) and F4/80(+)) cells. Importantly, the microenvironment orchestrated by TF expression drives permanent changes in the phenotype, gene-expression profile, DNA copy number, and DNA methylation state of the tumor cells that escape from dormancy. We postulate that procoagulant events in the tissue microenvironment (niche) may affect the fate of occult tumor cells, including their biological and genetic progression to initiate a full-blown malignancy.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Allelic expression mapping across cellular lineages to establish impact of non-coding SNPs.

Most complex disease-associated genetic variants are located in non-coding regions and are therefore thought to be regulatory in nature. Association mapping of differential allelic expression (AE) is a powerful method to identify SNPs with direct cis-regulatory impact (cis-rSNPs). We used AE mapping to identify cis-rSNPs regulating gene expression in 55 and 63 HapMap lymphoblastoid cell lines from a Caucasian and an African population, respectively, 70 fibroblast cell lines, and 188 purified monocyte samples and found 40-60% of these cis-rSNPs to be shared across cell types. We uncover a new class of cis-rSNPs, which disrupt footprint-derived de novo motifs that are predominantly bound by repressive factors and are implicated in disease susceptibility through overlaps with GWAS SNPs. Finally, we provide the proof-of-principle for a new approach for genome-wide functional validation of transcription factor-SNP interactions. By perturbing NFκB action in lymphoblasts, we identified 489 cis-regulated transcripts with altered AE after NFκB perturbation. Altogether, we perform a comprehensive analysis of cis-variation in four cell populations and provide new tools for the identification of functional variants associated to complex diseases.

Identification of divergent placental profiles in clinically distinct pregnancy complications revealed by the transcriptome.

INTRODUCTION: Pregnancy complications, including preeclampsia (PE), preterm birth (PTB), and intra-uterine growth restriction (IUGR) have individually been associated with inflammation but the combined comparative analysis of their placental profiles at the transcriptomic and histological levels is lacking. METHODS: Bulk RNA-sequencing of human placental biopsies from uncomplicated term pregnancies (CTL) and pregnancies complicated with early-onset (EO), and late-onset (LO) PE, as well as PTB and term IUGR were used to characterize individual molecular profiles. We also applied immune-cell-specific cellular deconvolution to address local immune cell compositions and analyzed placental lesions by histology to further characterize these complications. RESULTS: Transcriptome analysis revealed that clinically distinct complications differentiated themselves in unique ways compared to CTLs. Only TMEM136 was commonly modulated. Compared to CTLs, we found that PTB and IUGR were the most distinct, with LOPE being the least distinct. PTB and IUGR revealed differently enhanced inflammatory pathways, where PTB had general inflammatory responses and IUGR had immune cell activation. This inflammation was reflected in the histological profile for PTB only, whereas structural lesions were elevated in all complications. Placental lesions additionally had corresponding enhancement in inflammatory and structural biological processes. We observed that having co-complications, particularly for PTB with or without IUGR, impacted placental transcriptomes. Lastly, cellular deconvolution uncovered shared immune features among the complications. DISCUSSION: Overall, we provide evidence that these pregnancy complications are not only distinct in their clinical manifestations but also in their placental profiles, which could be leveraged to understand their underlying mechanisms and could offer therapeutic targets.

Repurposing disulfiram, an alcohol-abuse drug, in neuroblastoma causes KAT2A downregulation and in vivo activity with a water/oil emulsion.

Neuroblastoma, the most common type of pediatric extracranial solid tumor, causes 10% of childhood cancer deaths. Despite intensive multimodal treatment, the outcomes of high-risk neuroblastoma remain poor. We urgently need to develop new therapies with safe long-term toxicity profiles for rapid testing in clinical trials. Drug repurposing is a promising approach to meet these needs. Here, we investigated disulfiram, a safe and successful chronic alcoholism treatment with known anticancer and epigenetic effects. Disulfiram efficiently induced cell cycle arrest and decreased the viability of six human neuroblastoma cell lines at half-maximal inhibitory concentrations up to 20 times lower than its peak clinical plasma level in patients treated for chronic alcoholism. Disulfiram shifted neuroblastoma transcriptome, decreasing MYCN levels and activating neuronal differentiation. Consistently, disulfiram significantly reduced the protein level of lysine acetyltransferase 2A (KAT2A), drastically reducing acetylation of its target residues on histone H3. To investigate disulfiram's anticancer effects in an in vivo model of high-risk neuroblastoma, we developed a disulfiram-loaded emulsion to deliver the highly liposoluble drug. Treatment with the emulsion significantly delayed neuroblastoma progression in mice. These results identify KAT2A as a novel target of disulfiram, which directly impacts neuroblastoma epigenetics and is a promising candidate for repurposing to treat pediatric neuroblastoma.

Outlier Expression of Isoforms by Targeted or Total RNA Sequencing Identifies Clinically Significant Genomic Variants in Hematolymphoid Tumors.

Recognition of aberrant gene isoforms due to DNA events can impact risk stratification and molecular classification of hematolymphoid tumors. In myelodysplastic syndromes, KMT2A partial tandem duplication (PTD) was one of the top adverse predictors in the International Prognostic Scoring System-Molecular study. In B-cell acute lymphoblastic leukemia (B-ALL), ERG isoforms have been proposed as markers of favorable-risk DUX4 rearrangements, whereas deletion-mediated IKZF1 isoforms are associated with adverse prognosis and have been extended to the high-risk IKZF1plus signature defined by codeletions, including PAX5. In this limited study, outlier expression of isoforms as markers of IKZF1 intragenic or 3' deletions, DUX4 rearrangements, or PAX5 intragenic deletions were 92.3% (48/52), 90% (9/10), or 100% (9/9) sensitive, respectively, and 98.7% (368/373), 100% (35/35), or 97.1% (102/105) specific, respectively, by targeted RNA sequencing, and 84.0% (21/25), 85.7% (6/7), or 81.8% (9/11) sensitive, respectively, and 98.2% (109/111), 98.4% (127/129), or 98.7% (78/79) specific, respectively, by total RNA sequencing. Comprehensive split-read analysis identified expressed DNA breakpoints, cryptic splice sites associated with IKZF1 3' deletions, PTD of IKZF1 exon 5 spanning N159Y in B-ALL with mutated IKZF1 N159Y, and truncated KMT2A-PTD isoforms. Outlier isoforms were also effective targeted RNA markers for PAX5 intragenic amplifications (B-ALL), KMT2A-PTD (myeloid malignant cancers), and rare NOTCH1 intragenic deletions (T-cell acute lymphoblastic leukemia). These findings support the use of outlier isoform analysis as a robust strategy for detecting clinically significant DNA events.

Designing small multiple-target artificial RNAs.

MicroRNAs (miRNAs) are naturally occurring small RNAs that regulate the expression of several genes. MiRNAs' targeting rules are based on sequence complementarity between their mature products and targeted genes' mRNAs. Based on our present understanding of those rules, we developed an algorithm to design artificial miRNAs to target simultaneously a set of predetermined genes. To validate in silico our algorithm, we tested different sets of genes known to be targeted by a single miRNA. The algorithm finds the seed of the corresponding miRNA among the solutions, which also include the seeds of new artificial miRNA sequences potentially capable of targeting these genes as well. We also validated the functionality of some artificial miRNAs designed to target simultaneously members of the E2F family. These artificial miRNAs reproduced the effects of E2Fs inhibition in both normal human fibroblasts and prostate cancer cells where they inhibited cell proliferation and induced cellular senescence. We conclude that the current miRNA targeting rules based on the seed sequence work to design multiple-target artificial miRNAs. This approach may find applications in both research and therapeutics.

Identification of new ETV6 modulators through a high-throughput functional screening.

ETV6 transcriptional activity is critical for proper blood cell development in the bone marrow. Despite the accumulating body of evidence linking ETV6 malfunction to hematological malignancies, its regulatory network remains unclear. To uncover genes that modulate ETV6 repressive transcriptional activity, we performed a specifically designed, unbiased genome-wide shRNA screen in pre-B acute lymphoblastic leukemia cells. Following an extensive validation process, we identified 13 shRNAs inducing overexpression of ETV6 transcriptional target genes. We showed that the silencing of AKIRIN1, COMMD9, DYRK4, JUNB, and SRP72 led to an abrogation of ETV6 repressive activity. We identified critical modulators of the ETV6 function which could participate in cellular transformation through the ETV6 transcriptional network.

Whole-transcriptome analysis in acute lymphoblastic leukemia: a report from the DFCI ALL Consortium Protocol 16-001.

The molecular hallmark of childhood acute lymphoblastic leukemia (ALL) is characterized by recurrent, prognostic genetic alterations, many of which are cryptic by conventional cytogenetics. RNA sequencing (RNA-seq) is a powerful next-generation sequencing technology that can simultaneously identify cryptic gene rearrangements, sequence mutations and gene expression profiles in a single assay. We examined the feasibility and utility of incorporating RNA-seq into a prospective multicenter phase 3 clinical trial for children with newly diagnosed ALL. The Dana-Farber Cancer Institute ALL Consortium Protocol 16-001 enrolled 173 patients with ALL who consented to optional studies and had samples available for RNA-seq. RNA-seq identified at least 1 alteration in 157 patients (91%). Fusion detection was 100% concordant with results obtained from conventional cytogenetic analyses. An additional 56 gene fusions were identified by RNA-seq, many of which confer prognostic or therapeutic significance. Gene expression profiling enabled further molecular classification into the following B-cell ALL (B-ALL) subgroups: high hyperdiploid (n = 36), ETV6-RUNX1/-like (n = 31), TCF3-PBX1 (n = 7), KMT2A-rearranged (KMT2A-R; n = 5), intrachromosomal amplification of chromosome 21 (iAMP21) (n = 1), hypodiploid (n = 1), Philadelphia chromosome (Ph)-positive/Ph-like (n = 16), DUX4-R (n = 11), PAX5 alterations (PAX5 alt; n = 11), PAX5 P80R (n = 1), ZNF384-R (n = 4), NUTM1-R (n = 1), MEF2D-R (n = 1), and others (n = 10). RNA-seq identified 141 nonsynonymous mutations in 93 patients (54%); the most frequent were RAS-MAPK pathway mutations. Among 79 patients with both low-density array and RNA-seq data for the Philadelphia chromosome-like gene signature prediction, results were concordant in 74 patients (94%). In conclusion, RNA-seq identified several clinically relevant genetic alterations not detected by conventional methods, which supports the integration of this technology into front-line pediatric ALL trials. This trial was registered at www.clinicaltrials.gov as #NCT03020030.

Single Cell Transcriptomics of Ependymal Cells Across Age, Region and Species Reveals Cilia-Related and Metal Ion Regulatory Roles as Major Conserved Ependymal Cell Functions.

Ependymal cells are ciliated-epithelial glial cells that develop from radial glia along the surface of the ventricles of the brain and the spinal canal. They play a critical role in cerebrospinal fluid (CSF) homeostasis, brain metabolism, and the clearance of waste from the brain. These cells have been implicated in disease across the lifespan including developmental disorders, cancer, and neurodegenerative disease. Despite this, ependymal cells remain largely understudied. Using single-cell RNA sequencing data extracted from publicly available datasets, we make key findings regarding the remarkable conservation of ependymal cell gene signatures across age, region, and species. Through this unbiased analysis, we have discovered that one of the most overrepresented ependymal cell functions that we observed relates to a critically understudied role in metal ion homeostasis. Our analysis also revealed distinct subtypes and states of ependymal cells across regions and ages of the nervous system. For example, neonatal ependymal cells maintained a gene signature consistent with developmental processes such as determination of left/right symmetry; while adult ventricular ependymal cells, not spinal canal ependymal cells, appeared to express genes involved in regulating cellular transport and inflammation. Together, these findings highlight underappreciated functions of ependymal cells, which will be important to investigate in order to better understand these cells in health and disease.