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.

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.

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.

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.

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.

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.

Repurposing proscillaridin A in combination with decitabine against embryonal rhabdomyosarcoma RD cells.

PURPOSE: Embryonal rhabdomyosarcoma (eRMS) is the most common type of rhabdomyosarcoma in children. eRMS is characterized by malignant skeletal muscle cells driven by hyperactivation of several oncogenic pathways including the MYC pathway. Targeting MYC in cancer has been extremely challenging. Recently, we have demonstrated that the heart failure drug, proscillaridin A, produced anticancer effects with specificity toward MYC expressing leukemia cells. We also reported that decitabine, a hypomethylating drug, synergizes with proscillaridin A in colon cancer cells. Here, we investigated whether proscillaridin A exhibits epigenetic and anticancer activity against eRMS RD cells, overexpressing MYC oncogene, and its combination with decitabine. METHODS: We investigated the anticancer effects of proscillaridin A in eRMS RD cells in vitro. In response to drug treatment, we measured growth inhibition, cell cycle arrest, loss of clonogenicity and self-renewal capacity. We further evaluated the impact of proscillaridin A on MYC expression and its downstream transcriptomic effects by RNA sequencing. Then, we measured protein expression of epigenetic regulators and their associated chromatin post-translational modifications in response to drug treatment. Chromatin immunoprecipitation sequencing data sets were coupled with transcriptomic results to pinpoint the impact of proscillaridin A on gene pathways associated with specific chromatin modifications. Lastly, we evaluated the effect of the combination of proscillaridin A and the DNA demethylating drug decitabine on eRMS RD cell growth and clonogenic potential. RESULTS: Clinically relevant concentration of proscillaridin A (5 nM) produced growth inhibition, cell cycle arrest and loss of clonogenicity in eRMS RD cells. Proscillaridin A produced a significant downregulation of MYC protein expression and inhibition of oncogenic transcriptional programs controlled by MYC, involved in cell replication. Interestingly, significant reduction in total histone 3 acetylation and on specific lysine residues (lysine 9, 14, 18, and 27 on histone 3) was associated with significant protein downregulation of a series of lysine acetyltransferases (KAT3A, KAT3B, KAT2A, KAT2B, and KAT5). In addition, proscillaridin A produced synergistic growth inhibition and loss of clonogenicity when combined with the approved DNA demethylating drug decitabine. CONCLUSION: Proscillaridin A produces anticancer and epigenetic effects in the low nanomolar range and its combination with decitabine warrants further investigation for the treatment of eRMS.

Invasive growth associated with cold-inducible RNA-binding protein expression drives recurrence of surgically resected brain metastases.

BACKGROUND: Sixty percent of surgically resected brain metastases (BrM) recur within 1 year. These recurrences have long been thought to result from the dispersion of cancer cells during surgery. We tested the alternative hypothesis that invasion of cancer cells into the adjacent brain plays a significant role in local recurrence and shortened overall survival. METHODS: We determined the invasion pattern of 164 surgically resected BrM and correlated with local recurrence and overall survival. We performed single-cell RNA sequencing (scRNAseq) of >15,000 cells from BrM and adjacent brain tissue. Validation of targets was performed with a novel cohort of BrM patient-derived xenografts (PDX) and patient tissues. RESULTS: We demonstrate that invasion of metastatic cancer cells into the adjacent brain is associated with local recurrence and shortened overall survival. scRNAseq of paired tumor and adjacent brain samples confirmed the existence of invasive cancer cells in the tumor-adjacent brain. Analysis of these cells identified cold-inducible RNA-binding protein (CIRBP) overexpression in invasive cancer cells compared to cancer cells located within the metastases. Applying PDX models that recapitulate the invasion pattern observed in patients, we show that CIRBP is overexpressed in highly invasive BrM and is required for efficient invasive growth in the brain. CONCLUSIONS: These data demonstrate peritumoral invasion as a driver of treatment failure in BrM that is functionally mediated by CIRBP. These findings improve our understanding of the biology underlying postoperative treatment failure and lay the groundwork for rational clinical trial development based upon invasion pattern in surgically resected BrM.

Single-cell analysis of childhood leukemia reveals a link between developmental states and ribosomal protein expression as a source of intra-individual heterogeneity.

Childhood acute lymphoblastic leukemia (cALL) is the most common pediatric cancer. It is characterized by bone marrow lymphoid precursors that acquire genetic alterations, resulting in disrupted maturation and uncontrollable proliferation. More than a dozen molecular subtypes of variable severity can be used to classify cALL cases. Modern therapy protocols currently cure 85-90% of cases, but other patients are refractory or will relapse and eventually succumb to their disease. To better understand intratumor heterogeneity in cALL patients, we investigated the nature and extent of transcriptional heterogeneity at the cellular level by sequencing the transcriptomes of 39,375 individual cells in eight patients (six B-ALL and two T-ALL) and three healthy pediatric controls. We observed intra-individual transcriptional clusters in five out of the eight patients. Using pseudotime maturation trajectories of healthy B and T cells, we obtained the predicted developmental state of each leukemia cell and observed distribution shifts within patients. We showed that the predicted developmental states of these cancer cells are inversely correlated with ribosomal protein expression levels, which could be a common contributor to intra-individual heterogeneity in cALL patients.

Heart failure drug proscillaridin A targets MYC overexpressing leukemia through global loss of lysine acetylation.

BACKGROUND: Cardiac glycosides are approved for the treatment of heart failure as Na+/K+ pump inhibitors. Their repurposing in oncology is currently investigated in preclinical and clinical studies. However, the identification of a specific cancer type defined by a molecular signature to design targeted clinical trials with cardiac glycosides remains to be characterized. Here, we demonstrate that cardiac glycoside proscillaridin A specifically targets MYC overexpressing leukemia cells and leukemia stem cells by causing MYC degradation, epigenetic reprogramming and leukemia differentiation through loss of lysine acetylation. METHODS: Proscillaridin A anticancer activity was investigated against a panel of human leukemia and solid tumor cell lines with different MYC expression levels, overexpression in vitro systems and leukemia stem cells. RNA-sequencing and differentiation studies were used to characterize transcriptional and phenotypic changes. Drug-induced epigenetic changes were studied by chromatin post-translational modification analysis, expression of chromatin regulators, chromatin immunoprecipitation, and mass-spectrometry. RESULTS: At a clinically relevant dose, proscillaridin A rapidly altered MYC protein half-life causing MYC degradation and growth inhibition. Transcriptomic profile of leukemic cells after treatment showed a downregulation of genes involved in MYC pathways, cell replication and an upregulation of hematopoietic differentiation genes. Functional studies confirmed cell cycle inhibition and the onset of leukemia differentiation even after drug removal. Proscillaridin A induced a significant loss of lysine acetylation in histone H3 (at lysine 9, 14, 18 and 27) and in non-histone proteins such as MYC itself, MYC target proteins, and a series of histone acetylation regulators. Global loss of acetylation correlated with the rapid downregulation of histone acetyltransferases. Importantly, proscillaridin A demonstrated anticancer activity against lymphoid and myeloid stem cell populations characterized by MYC overexpression. CONCLUSION: Overall, these results strongly support the repurposing of proscillaridin A in MYC overexpressing leukemia.

GenPipes: an open-source framework for distributed and scalable genomic analyses.

BACKGROUND: With the decreasing cost of sequencing and the rapid developments in genomics technologies and protocols, the need for validated bioinformatics software that enables efficient large-scale data processing is growing. FINDINGS: Here we present GenPipes, a flexible Python-based framework that facilitates the development and deployment of multi-step workflows optimized for high-performance computing clusters and the cloud. GenPipes already implements 12 validated and scalable pipelines for various genomics applications, including RNA sequencing, chromatin immunoprecipitation sequencing, DNA sequencing, methylation sequencing, Hi-C, capture Hi-C, metagenomics, and Pacific Biosciences long-read assembly. The software is available under a GPLv3 open source license and is continuously updated to follow recent advances in genomics and bioinformatics. The framework has already been configured on several servers, and a Docker image is also available to facilitate additional installations. CONCLUSIONS: GenPipes offers genomics researchers a simple method to analyze different types of data, customizable to their needs and resources, as well as the flexibility to create their own workflows.

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.

Molecular Signature of CAID Syndrome: Noncanonical Roles of SGO1 in Regulation of TGF-ß Signaling and Epigenomics.

BACKGROUND & AIMS: A generalized human pacemaking syndrome, chronic atrial and intestinal dysrhythmia (CAID) (OMIM 616201), is caused by a homozygous SGO1 mutation (K23E), leading to chronic intestinal pseudo-obstruction and arrhythmias. Because CAID patients do not show phenotypes consistent with perturbation of known roles of SGO1, we hypothesized that noncanonical roles of SGO1 drive the clinical manifestations observed. METHODS: To identify a molecular signature for CAID syndrome, we achieved unbiased screens in cell lines and gut tissues from CAID patients vs wild-type controls. We performed RNA sequencing along with stable isotope labeling with amino acids in cell culture. In addition, we determined the genome-wide DNA methylation and chromatin accessibility signatures using reduced representative bisulfite sequencing and assay for transposase-accessible chromatin with high-throughput sequencing. Functional studies included patch-clamp, quantitation of transforming growth factor-ß (TGF-ß) signaling, and immunohistochemistry in CAID patient gut biopsy specimens. RESULTS: Proteome and transcriptome studies converge on cell-cycle regulation, cardiac conduction, and smooth muscle regulation as drivers of CAID syndrome. Specifically, the inward rectifier current, an important regulator of cellular function, was disrupted. Immunohistochemistry confirmed overexpression of Budding Uninhibited By Benzimidazoles 1 (BUB1) in patients, implicating the TGF-ß pathway in CAID pathogenesis. Canonical TGF-ß signaling was up-regulated and uncoupled from noncanonical signaling in CAID patients. Reduced representative bisulfite sequencing and assay for transposase-accessible chromatin with high-throughput sequencing experiments showed significant changes of chromatin states in CAID, pointing to epigenetic regulation as a possible pathologic mechanism. CONCLUSIONS: Our findings point to impaired inward rectifier potassium current, dysregulation of canonical TGF-ß signaling, and epigenetic regulation as potential drivers of intestinal and cardiac manifestations of CAID syndrome. Transcript profiling and genomics data are as follows: repository URL: https://www.ncbi.nlm.nih.gov/geo; SuperSeries GSE110612 was composed of the following subseries: GSE110309, GSE110576, and GSE110601.

Very long intergenic non-coding RNA transcripts and expression profiles are associated to specific childhood acute lymphoblastic leukemia subtypes.

Very long intergenic non-coding RNAs (vlincRNAs) are a novel class of long transcripts (~50 kb to 1 Mb) with cell type- or cancer-specific expression. We report the discovery and characterization of 256 vlincRNAs from a cohort of 64 primary childhood pre-B and pre-T acute lymphoblastic leukemia (cALL) samples, of which 61% are novel and specifically expressed in cALL. Validation was performed in 35 pre-B and pre-T cALL primary samples. We show that their expression is cALL immunophenotype and molecular subtype-specific and correlated with epigenetic modifications on their promoters, much like protein-coding genes. While the biological functions of these vlincRNAs are still unknown, our results suggest they could play a role in cALL etiology or progression.

Genome wide mapping of ETV6 binding sites in pre-B leukemic cells.

Genetic alterations in the transcriptional repressor ETV6 are associated with hematological malignancies. Notably, the t(12;21) translocation leading to an ETV6-AML1 fusion gene is the most common genetic alteration found in childhood acute lymphoblastic leukemia. Moreover, most of these patients also lack ETV6 expression, suggesting a tumor suppressor function. To gain insights on ETV6 DNA-binding specificity and genome wide transcriptional regulation capacities, we performed chromatin immunoprecipitation experiments coupled to deep sequencing in a t(12;21)-positive pre-B leukemic cell line. This strategy led to the identification of ETV6-bound regions that were further associated to gene expression. ETV6 binding is mostly cell type-specific as only few regions are shared with other blood cell subtypes. Peaks localization and motif enrichment analyses revealed that this unique binding profile could be associated with the ETV6-AML1 fusion protein specific to the t(12;21) background. This study underscores the complexity of ETV6 binding and uncovers ETV6 transcriptional network in pre-B leukemia cells bearing the recurrent t(12;21) translocation.

LncRNAs downregulated in childhood acute lymphoblastic leukemia modulate apoptosis, cell migration, and DNA damage response.

Childhood acute lymphoblastic leukemia (cALL) accounts for 25% of pediatric cancers and is one of the leading causes of disease-related death in children. Although long non-coding RNAs (lncRNAs) have been implicated in cALL etiology, progression, and treatment response, little is known about their exact functional role. We had previously sequenced the whole transcriptome of 56 cALL patients and identified lncRNA transcripts specifically silenced in tumoral cells. Here we investigated the impact of restoring the expression of three of these (RP11-624C23.1, RP11-203E8, and RP11-446E9) in leukemic cell lines had dramatic impact on cancer hallmark cellular phenotypes such as apoptosis, cell proliferation and migration, and DNA damage response. Interestingly, both RP11-624C23.1 and RP11-203E8 had very similar impacts on DNA damage response, specifically displaying lower γ-H2A.X and higher apoptosis levels than control cells in response to genotoxic stress. These results indicate that silencing RP11-624C23.1 or RP11-203E8 could provide a selective advantage to leukemic cells by increasing resistance to genotoxic stress, possibly by modulating the DDR pathway. Since genotoxic agents are fundamental parts of antineoplastic treatment, further investigation of the mechanisms these lncRNAs impact would provide novel and interesting avenues for overcoming treatment resistance.