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.

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.

Acute stress drives global repression through two independent RNA polymerase II stalling events in Saccharomyces.

In multicellular eukaryotes, RNA polymerase (Pol) II pauses transcription ~30-50 bp after initiation. While the budding yeast Saccharomyces has its transcription mechanisms mostly conserved with other eukaryotes, it appears to lack this fundamental promoter-proximal pausing. However, we now report that nearly all yeast genes, including constitutive and inducible genes, manifest two distinct transcriptional stall sites that are brought on by acute environmental signaling (e.g., peroxide stress). Pol II first stalls at the pre-initiation stage before promoter clearance, but after DNA melting and factor acquisition, and may involve inhibited dephosphorylation. The second stall occurs at the +2 nucleosome. It acquires most, but not all, elongation factor interactions. Its regulation may include Bur1/Spt4/5. Our results suggest that a double Pol II stall is a mechanism to downregulate essentially all genes in concert.

In vivo genome-wide binding interactions of mouse and human constitutive androstane receptors reveal novel gene targets.

The constitutive androstane receptor (CAR; NR1I3) is a nuclear receptor orchestrating complex roles in cell and systems biology. Species differences in CAR's effector pathways remain poorly understood, including its role in regulating liver tumor promotion. We developed transgenic mouse models to assess genome-wide binding of mouse and human CAR, following receptor activation in liver with direct ligands and with phenobarbital, an indirect CAR activator. Genomic interaction profiles were integrated with transcriptional and biological pathway analyses. Newly identified CAR target genes included Gdf15 and Foxo3, important regulators of the carcinogenic process. Approximately 1000 genes exhibited differential binding interactions between mouse and human CAR, including the proto-oncogenes, Myc and Ikbke, which demonstrated preferential binding by mouse CAR as well as mouse CAR-selective transcriptional enhancement. The ChIP-exo analyses also identified distinct binding motifs for the respective mouse and human receptors. Together, the results provide new insights into the important roles that CAR contributes as a key modulator of numerous signaling pathways in mammalian organisms, presenting a genomic context that specifies species variation in biological processes under CAR's control, including liver cell proliferation and tumor promotion.

Tail and Kinase Modules Differently Regulate Core Mediator Recruitment and Function In Vivo.

Mediator is a highly conserved transcriptional coactivator organized into four modules, namely Tail, Middle, Head, and Kinase (CKM). Previous work suggests regulatory roles for Tail and CKM, but an integrated model for these activities is lacking. Here, we analyzed the genome-wide distribution of Mediator subunits in wild-type and mutant yeast cells in which RNA polymerase II promoter escape is blocked, allowing detection of transient Mediator forms. We found that although all modules are recruited to upstream activated regions (UAS), assembly of Mediator within the pre-initiation complex is accompanied by the release of CKM. Interestingly, our data show that CKM regulates Mediator-UAS interaction rather than Mediator-promoter association. In addition, although Tail is required for Mediator recruitment to UAS, Tailless Mediator nevertheless interacts with core promoters. Collectively, our data suggest that the essential function of Mediator is mediated by Head and Middle at core promoters, while Tail and CKM play regulatory roles.

The Histone Modification Domain of Paf1 Complex Subunit Rtf1 Directly Stimulates H2B Ubiquitylation through an Interaction with Rad6.

The five-subunit yeast Paf1 complex (Paf1C) regulates all stages of transcription and is critical for the monoubiquitylation of histone H2B (H2Bub), a modification that broadly influences chromatin structure and eukaryotic transcription. Here, we show that the histone modification domain (HMD) of Paf1C subunit Rtf1 directly interacts with the ubiquitin conjugase Rad6 and stimulates H2Bub independently of transcription. We present the crystal structure of the Rtf1 HMD and use site-specific, in vivo crosslinking to identify a conserved Rad6 interaction surface. Utilizing ChIP-exo analysis, we define the localization patterns of the H2Bub machinery at high resolution and demonstrate the importance of Paf1C in targeting the Rtf1 HMD, and thereby H2Bub, to its appropriate genomic locations. Finally, we observe HMD-dependent stimulation of H2Bub in a transcription-free, reconstituted in vitro system. Taken together, our results argue for an active role for Paf1C in promoting H2Bub and ensuring its proper localization in vivo.

Genome-Wide Organization of GATA1 and TAL1 Determined at High Resolution.

Erythroid development and differentiation from multiprogenitor cells into red blood cells requires precise transcriptional regulation. Key erythroid transcription factors, GATA1 and TAL1, cooperate, along with other proteins, to regulate many aspects of this process. How GATA1 and TAL1 are juxtaposed along the DNA and their cognate DNA binding site across the mouse genome remains unclear. We applied high-resolution ChIP-exo (chromatin immunoprecipitation followed by 5'-to-3' exonuclease treatment and then massively parallel DNA sequencing) to GATA1 and TAL1 to study their positional organization across the mouse genome during GATA1-dependent maturation. Two complementary methods, MultiGPS and peak pairing, were used to determine high-confidence binding locations by ChIP-exo. We identified ∼10,000 GATA1 and ∼15,000 TAL1 locations, which were essentially confirmed by ChIP-seq (chromatin immunoprecipitation followed by massively parallel DNA sequencing). Of these, ∼4,000 locations were bound by both GATA1 and TAL1. About three-quarters of them were tightly linked to a partial E-box located 7 or 8 bp upstream of a WGATAA motif. Both TAL1 and GATA1 generated distinct characteristic ChIP-exo peaks around WGATAA motifs that reflect their positional arrangement within a complex. We show that TAL1 and GATA1 form a precisely organized complex at a compound motif consisting of a TG 7 or 8 bp upstream of a WGATAA motif across thousands of genomic locations.

Subnucleosomal structures and nucleosome asymmetry across a genome.

Genes are packaged into nucleosomal arrays, each nucleosome typically having two copies of histones H2A, H2B, H3, and H4. Histones have distinct posttranslational modifications, variant isoforms, and dynamics. Whether each histone copy within a nucleosome has distinct properties, particularly in relation to the direction of transcription, is unknown. Here we use chromatin immunoprecipitation-exonuclease (ChIP-exo) to resolve the organization of individual histones on a genomic scale. We detect widespread subnucleosomal structures in dynamic chromatin, including what appear to be half-nucleosomes consisting of one copy of each histone. We also detect interactions of H3 tails with linker DNA between nucleosomes, which may be negatively regulated by methylation of H3K36. Histone variant H2A.Z is enriched on the promoter-distal half of the +1 nucleosome, whereas H2BK123 ubiquitylation and H3K9 acetylation are enriched on the promoter-proximal half in a transcription-linked manner. Subnucleosome asymmetries might serve as molecular beacons that guide transcription.

A universal RNA polymerase II CTD cycle is orchestrated by complex interplays between kinase, phosphatase, and isomerase enzymes along genes.

Transcription by RNA polymerase II (RNAPII) is coupled to mRNA processing and chromatin modifications via the C-terminal domain (CTD) of its largest subunit, consisting of multiple repeats of the heptapeptide YSPTSPS. Pioneering studies showed that CTD serines are differentially phosphorylated along genes in a prescribed pattern during the transcription cycle. Genome-wide analyses challenged this idea, suggesting that this cycle is not uniform among different genes. Moreover, the respective role of enzymes responsible for CTD modifications remains controversial. Here, we systematically profiled the location of the RNAPII phosphoisoforms in wild-type cells and mutants for most CTD modifying enzymes. Together with results of in vitro assays, these data reveal a complex interplay between the modifying enzymes, and provide evidence that the CTD cycle is uniform across genes. We also identify Ssu72 as the Ser7 phosphatase and show that proline isomerization is a key regulator of CTD dephosphorylation at the end of genes.

Systematic identification of fragile sites via genome-wide location analysis of gamma-H2AX.

Phosphorylation of histone H2AX is an early response to DNA damage in eukaryotes. In Saccharomyces cerevisiae, DNA damage or replication-fork stalling results in phosphorylation of histone H2A yielding gamma-H2A (yeast gamma-H2AX) in a Mec1 (ATR)- and Tel1 (ATM)-dependent manner. Here, we describe the genome-wide location analysis of gamma-H2A as a strategy to identify loci prone to engaging the Mec1 and Tel1 pathways. Notably, gamma-H2A enrichment overlaps with loci prone to replication-fork stalling and is caused by the action of Mec1 and Tel1, indicating that these loci are prone to breakage. Moreover, about half the sites enriched for gamma-H2A map to repressed protein-coding genes, and histone deacetylases are necessary for formation of gamma-H2A at these loci. Finally, our work indicates that high-resolution mapping of gamma-H2AX is a fruitful route to map fragile sites in eukaryotic genomes.

Profiling genome-wide histone modifications and variants by ChIP-chip on tiling microarrays in S. cerevisiae.

Chromatin immunoprecipitation coupled to DNA microarray has become a widely used method to study transcription factors and chromatin structure. Here, we provide a detailed protocol for the localization of the variant histone Htz1 in the S. cerevisiae genome. This protocol can easily be adapted to fit other purposes such as profiling histone modifications.

Nuclear receptor ERR alpha and coactivator PGC-1 beta are effectors of IFN-gamma-induced host defense.

Macrophage activation by the proinflammatory cytokine interferon-gamma (IFN-gamma) is a critical component of the host innate response to bacterial pathogenesis. However, the precise nature of the IFN-gamma-induced activation pathway is not known. Here we show using genome-wide expression and chromatin-binding profiling that IFN-gamma induces the expression of many nuclear genes encoding mitochondrial respiratory chain machinery via activation of the nuclear receptor ERR alpha (estrogen-related receptor alpha, NR3B1). Studies with macrophages lacking ERR alpha demonstrate that it is required for induction of mitochondrial reactive oxygen species (ROS) production and efficient clearance of Listeria monocytogenes (LM) in response to IFN-gamma. As a result, mice lacking ERR alpha are susceptible to LM infection, a phenotype that is localized to bone marrow-derived cells. Furthermore, we found that IFN-gamma-induced activation of ERR alpha depends on coactivator PGC-1 beta (peroxisome proliferator-activated receptor gamma coactivator-1 beta), which appears to be a direct target for the IFN-gamma/STAT-1 signaling cascade. Thus, ERR alpha and PGC-1 beta act together as a key effector of IFN-gamma-induced mitochondrial ROS production and host defense.

Genome-wide computational prediction of transcriptional regulatory modules reveals new insights into human gene expression.

The identification of regulatory regions is one of the most important and challenging problems toward the functional annotation of the human genome. In higher eukaryotes, transcription-factor (TF) binding sites are often organized in clusters called cis-regulatory modules (CRM). While the prediction of individual TF-binding sites is a notoriously difficult problem, CRM prediction has proven to be somewhat more reliable. Starting from a set of predicted binding sites for more than 200 TF families documented in Transfac, we describe an algorithm relying on the principle that CRMs generally contain several phylogenetically conserved binding sites for a few different TFs. The method allows the prediction of more than 118,000 CRMs within the human genome. A subset of these is shown to be bound in vivo by TFs using ChIP-chip. Their analysis reveals, among other things, that CRM density varies widely across the genome, with CRM-rich regions often being located near genes encoding transcription factors involved in development. Predicted CRMs show a surprising enrichment near the 3' end of genes and in regions far from genes. We document the tendency for certain TFs to bind modules located in specific regions with respect to their target genes and identify TFs likely to be involved in tissue-specific regulation. The set of predicted CRMs, which is made available as a public database called PReMod (http://genomequebec.mcgill.ca/PReMod), will help analyze regulatory mechanisms in specific biological systems.

Variant histone H2A.Z is globally localized to the promoters of inactive yeast genes and regulates nucleosome positioning.

H2A.Z is an evolutionary conserved histone variant involved in transcriptional regulation, antisilencing, silencing, and genome stability. The mechanism(s) by which H2A.Z regulates these various biological functions remains poorly defined, in part due to the lack of knowledge regarding its physical location along chromosomes and the bearing it has in regulating chromatin structure. Here we mapped H2A.Z across the yeast genome at an approximately 300-bp resolution, using chromatin immunoprecipitation combined with tiling microarrays. We have identified 4,862 small regions--typically one or two nucleosomes wide--decorated with H2A.Z. Those "Z loci" are predominantly found within specific nucleosomes in the promoter of inactive genes all across the genome. Furthermore, we have shown that H2A.Z can regulate nucleosome positioning at the GAL1 promoter. Within HZAD domains, the regions where H2A.Z shows an antisilencing function, H2A.Z is localized in a wider pattern, suggesting that the variant histone regulates a silencing and transcriptional activation via different mechanisms. Our data suggest that the incorporation of H2A.Z into specific promoter-bound nucleosomes configures chromatin structure to poise genes for transcriptional activation. The relevance of these findings to higher eukaryotes is discussed.

From the Cover: Location analysis of estrogen receptor alpha target promoters reveals that FOXA1 defines a domain of the estrogen response.

Nuclear receptors can activate diverse biological pathways within a target cell in response to their cognate ligands, but how this compartmentalization is achieved at the level of gene regulation is poorly understood. We used a genome-wide analysis of promoter occupancy by the estrogen receptor alpha (ERalpha) in MCF-7 cells to investigate the molecular mechanisms underlying the action of 17beta-estradiol (E2) in controlling the growth of breast cancer cells. We identified 153 promoters bound by ERalpha in the presence of E2. Motif-finding algorithms demonstrated that the estrogen response element (ERE) is the most common motif present in these promoters whereas conventional chromatin immunoprecipitation assays showed E2-modulated recruitment of coactivator AIB1 and RNA polymerase II at these loci. The promoters were linked to known ERalpha targets but also to many genes not directly associated with the estrogenic response, including the transcriptional factor FOXA1, whose expression correlates with the presence of ERalpha in breast tumors. We found that ablation of FOXA1 expression in MCF-7 cells suppressed ERalpha binding to the prototypic TFF1 promoter (which contains a FOXA1 binding site), hindered the induction of TFF1 expression by E2, and prevented hormone-induced reentry into the cell cycle. Taken together, these results define a paradigm for estrogen action in breast cancer cells and suggest that regulation of gene expression by nuclear receptors can be compartmentalized into unique transcriptional domains by means of licensing of their activity to cofactors such as FOXA1.