Chromatin accessibility maps of chronic lymphocytic leukaemia identify subtype-specific epigenome signatures and transcription regulatory networks

Non-canonical transcriptional regulation of the poor prognostic factor UGT2B17 in chronic lymphocytic leukemic and normal B cells

BACKGROUND

High expression of the glycosyltransferase UGT2B17 represents an independent adverse prognostic marker in chronic lymphocytic leukemia (CLL). It also constitutes a predictive marker for therapeutic response and a drug resistance mechanism. The key determinants driving expression of the UGT2B17 gene in normal and leukemic B-cells remain undefined. The UGT2B17 transcriptome is complex and is comprised of at least 10 alternative transcripts, identified by previous RNA-sequencing of liver and intestine. We hypothesized that the transcriptional program regulating UGT2B17 in B-lymphocytes is distinct from the canonical expression previously characterized in the liver.

RESULTS

RNA-sequencing and genomics data revealed a specific genomic landscape at the UGT2B17 locus in normal and leukemic B-cells. RNA-sequencing and quantitative PCR data indicated that the UGT2B17 enzyme is solely encoded by alternative transcripts expressed in CLL patient cells and not by the canonical transcript widely expressed in the liver and intestine. Chromatin accessible regions (ATAC-Seq) in CLL cells mapped with alternative promoters and non-coding exons, which may be derived from endogenous retrotransposon elements. By luciferase reporter assays, we identified key cis-regulatory STAT3, RELA and interferon regulatory factor (IRF) binding sequences driving the expression of UGT2B17 in lymphoblastoid and leukemic B-cells. Electrophoretic mobility shift assays and pharmacological inhibition demonstrated key roles for the CLL prosurvival transcription factors STAT3 and NF-κB in the leukemic expression of UGT2B17.

CONCLUSIONS

UGT2B17 expression in B-CLL is driven by key regulators of CLL progression. Our data suggest that a NF-κB/STAT3/IRF/UGT2B17 axis may represent a novel B-cell pathway promoting disease progression and drug resistance.

Ehf and Fezf2 regulate late medullary thymic epithelial cell and thymic tuft cell development

Thymic epithelial cells are indispensable for T cell maturation and selection and the induction of central immune tolerance. The self-peptide repertoire expressed by medullary thymic epithelial cells is in part regulated by the transcriptional regulator Aire (Autoimmune regulator) and the transcription factor Fezf2. Due to the high complexity of mTEC maturation stages (i.e., post-Aire, Krt10+ mTECs, and Dclk1+ Tuft mTECs) and the heterogeneity in their gene expression profiles (i.e., mosaic expression patterns), it has been challenging to identify the additional factors complementing the transcriptional regulation. We aimed to identify the transcriptional regulators involved in the regulation of mTEC development and self-peptide expression in an unbiased and genome-wide manner. We used ATAC footprinting analysis as an indirect approach to identify transcription factors involved in the gene expression regulation in mTECs, which we validated by ChIP sequencing. This study identifies Fezf2 as a regulator of the recently described thymic Tuft cells (i.e., Tuft mTECs). Furthermore, we identify that transcriptional regulators of the ELF, ESE, ERF, and PEA3 subfamily of the ETS transcription factor family and members of the Krüppel-like family of transcription factors play a role in the transcriptional regulation of genes involved in late mTEC development and promiscuous gene expression.

Biological, prognostic, and therapeutic impact of the epigenome in CLL

Chronic lymphocytic leukemia (CLL) is characterized by widespread alterations in the genetic and epigenetic landscapes which seem to underlie the variable clinical manifestations observed in patients. Over the last decade, epigenomic studies have described the whole-genome maps of DNA methylation and chromatin features of CLL and normal B cells, identifying distinct epigenetic mechanisms operating in tumoral cells. DNA methylation analyses have identified that the CLL methylome contains imprints of the cell of origin, as well as of the proliferative history of the tumor cells, with both being strong independent prognostic predictors. Moreover, single-cell analysis revealed a higher degree of DNA methylation noise in CLL cells, which associates with transcriptional plasticity and disease aggressiveness. Integrative analysis of chromatin has uncovered chromatin signatures, as well as regulatory regions specifically active in each CLL subtype or in Richter transformed samples. Unique transcription factor (TF) binding motifs are overrepresented on those regions, suggesting that altered TF networks operate from disease initiation to progression as nongenetic factors mediating the oncogenic transcriptional profiles. Multiomics analysis has identified that response to treatment is modulated by an epigenetic imprint, and that treatments affect chromatin through the activity of particular set of TFs. Additionally, the epigenome is an axis of therapeutic vulnerability in CLL, as it can be targeted by inhibitors of histone modifying enzymes, that have shown promising preclinical results. Altogether, this review aims at summarizing the major findings derived from published literature to distill how altered epigenomic mechanisms contribute to CLL origin, evolution, clinical behavior, and response to treatment.

Chromatin accessibility landscape of relapsed pediatric B-lineage acute lymphoblastic leukemia

For around half of the pediatric B-lineage acute lymphoblastic leukemia (B-ALL) patients, the molecular mechanism of relapse remains unclear. To fill this gap in knowledge, here we characterize the chromatin accessibility landscape in pediatric relapsed B-ALL. We observe rewired accessible chromatin regions (ACRs) associated with transcription dysregulation in leukemia cells as compared with normal B-cell progenitors. We show that over a quarter of the ACRs in B-ALL are in quiescent regions with high heterogeneity among B-ALLs. We identify subtype-specific and allele-imbalanced chromatin accessibility by integrating multi-omics data. By characterizing the differential ACRs between diagnosis and relapse in B-ALL, we identify alterations in chromatin accessibility during drug treatment. Further analysis of ACRs associated with relapse free survival leads to the identification of a subgroup of B-ALL which show early relapse. These data provide an advanced and integrative portrait of the importance of chromatin accessibility alterations in tumorigenesis and drug responses.

TEAD4 antagonizes cellular senescence by remodeling chromatin accessibility at enhancer regions

Dramatic alterations in epigenetic landscapes are known to impact genome accessibility and transcription. Extensive evidence demonstrates that senescent cells undergo significant changes in chromatin structure; however, the mechanisms underlying the crosstalk between epigenetic parameters and gene expression profiles have not been fully elucidated. In the present study, we delineate the genome-wide redistribution of accessible chromatin regions that lead to broad transcriptome effects during senescence. We report that distinct senescence-activated accessibility regions (SAAs) are always distributed in H3K27ac-occupied enhancer regions, where they are responsible for elevated flanking senescence-associated secretory phenotype (SASP) expression and aberrant cellular signaling relevant to SASP secretion. Mechanistically, a single transcription factor, TEAD4, moves away from H3K27ac-labled SAAs to allow for prominent chromatin accessibility reconstruction during senescence. The enhanced SAAs signal driven by TEAD4 suppression subsequently induces a robust increase in the expression of adjacent SASP genes and the secretion of downstream factors, which contribute to the progression of senescence. Our findings illustrate a dynamic landscape of chromatin accessibility following senescence entry, and further reveal an insightful function for TEAD4 in regulating the broad chromatin state that modulates the overall transcriptional program of SASP genes.

Nucleosome repositioning in chronic lymphocytic leukemia

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.

IChrom-Deep: An Attention-Based Deep Learning Model for Identifying Chromatin Interactions

Identification of chromatin interactions is crucial for advancing our knowledge of gene regulation. However, due to the limitations of high-throughput experimental techniques, there is an urgent need to develop computational methods for predicting chromatin interactions. In this study, we propose a novel attention-based deep learning model, termed IChrom-Deep, to identify chromatin interactions using sequence features and genomic features. The experimental results based on the datasets of three cell lines demonstrate that the IChrom-Deep achieves satisfactory performance and is superior to the previous methods. We also investigate the effect of DNA sequence and associated features and genomic features on chromatin interactions, and highlight the applicable scenarios of some features, such as sequence conservation and distance. Moreover, we identify a few genomic features that are extremely important across different cell lines, and IChrom-Deep achieves comparable performance with only these significant genomic features versus using all genomic features. It is believed that IChrom-Deep can serve as a useful tool for future studies that seek to identify chromatin interactions.

Recruitment of transcription factor ETS1 to activated accessible regions promotes the transcriptional program of cilia genes

Defects in cilia genes, which are critical for cilia formation and function, can cause complicated ciliopathy syndromes involving multiple organs and tissues; however, the underlying regulatory mechanisms of the networks of cilia genes in ciliopathies remain enigmatic. Herein, we have uncovered the genome-wide redistribution of accessible chromatin regions and extensive alterations of expression of cilia genes during Ellis-van Creveld syndrome (EVC) ciliopathy pathogenesis. Mechanistically, the distinct EVC ciliopathy-activated accessible regions (CAAs) are shown to positively regulate robust changes in flanking cilia genes, which are a key requirement for cilia transcription in response to developmental signals. Moreover, a single transcription factor, ETS1, can be recruited to CAAs, leading to prominent chromatin accessibility reconstruction in EVC ciliopathy patients. In zebrafish, the collapse of CAAs driven by ets1 suppression subsequently causes defective cilia proteins, resulting in body curvature and pericardial oedema. Our results depict a dynamic landscape of chromatin accessibility in EVC ciliopathy patients, and uncover an insightful role for ETS1 in controlling the global transcriptional program of cilia genes by reprogramming the widespread chromatin state.

Immunosenescence: molecular mechanisms and diseases

Infection susceptibility, poor vaccination efficacy, age-related disease onset, and neoplasms are linked to innate and adaptive immune dysfunction that accompanies aging (known as immunosenescence). During aging, organisms tend to develop a characteristic inflammatory state that expresses high levels of pro-inflammatory markers, termed inflammaging. This chronic inflammation is a typical phenomenon linked to immunosenescence and it is considered the major risk factor for age-related diseases. Thymic involution, naïve/memory cell ratio imbalance, dysregulated metabolism, and epigenetic alterations are striking features of immunosenescence. Disturbed T-cell pools and chronic antigen stimulation mediate premature senescence of immune cells, and senescent immune cells develop a proinflammatory senescence-associated secretory phenotype that exacerbates inflammaging. Although the underlying molecular mechanisms remain to be addressed, it is well documented that senescent T cells and inflammaging might be major driving forces in immunosenescence. Potential counteractive measures will be discussed, including intervention of cellular senescence and metabolic-epigenetic axes to mitigate immunosenescence. In recent years, immunosenescence has attracted increasing attention for its role in tumor development. As a result of the limited participation of elderly patients, the impact of immunosenescence on cancer immunotherapy is unclear. Despite some surprising results from clinical trials and drugs, it is necessary to investigate the role of immunosenescence in cancer and other age-related diseases.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm

Comprehensive cis-regulatory landscapes are essential for accurate enhancer prediction and disease variant mapping. Although cis-regulatory element (CRE) resources exist for most tissues and organs, many rare - yet functionally important - cell types remain overlooked. Despite representing only a small fraction of the heart's cellular biomass, the cardiac conduction system (CCS) unfailingly coordinates every life-sustaining heartbeat. To globally profile the mouse CCS cis-regulatory landscape, we genetically tagged CCS component-specific nuclei for comprehensive assay for transposase-accessible chromatin-sequencing (ATAC-Seq) analysis. Thus, we established a global CCS-enriched CRE database, referred to as CCS-ATAC, as a key resource for studying CCS-wide and component-specific regulatory functions. Using transcription factor (TF) motifs to construct CCS component-specific gene regulatory networks (GRNs), we identified and independently confirmed several specific TF sub-networks. Highlighting the functional importance of CCS-ATAC, we also validated numerous CCS-enriched enhancer elements and suggested gene targets based on CCS single-cell RNA-Seq data. Furthermore, we leveraged CCS-ATAC to improve annotation of existing human variants related to cardiac rhythm and nominated a potential enhancer-target pair that was dysregulated by a specific SNP. Collectively, our results established a CCS-regulatory compendium, identified novel CCS enhancer elements, and illuminated potential functional associations between human genomic variants and CCS component-specific CREs.

Integrative epigenomics in chronic lymphocytic leukaemia: Biological insights and clinical applications

Chronic lymphocytic leukaemia (CLL) is not only characterised by driver genetic alterations but by extensive epigenetic changes. Over the last decade, epigenomic studies have described the DNA methylome, chromatin accessibility, histone modifications and the three-dimensional (3D) genome architecture of CLL. Beyond its regulatory role, the DNA methylome contains imprints of the cellular origin and proliferative history of CLL cells. These two aspects are strong independent prognostic factors. Integrative analyses of chromatin marks have uncovered novel regulatory elements and altered transcription factor networks as non-genetic means mediating gene deregulation in CLL. Additionally, CLL cells display a disease-specific pattern of 3D genome interactions. From the technological perspective, we are currently witnessing a transition from bulk omics to single-cell analyses. This review aims at summarising the major findings from the epigenomics field as well as providing a prospect of the present and future of single-cell analyses in CLL.

Dysregulation of PRMT5 in chronic lymphocytic leukemia promotes progression with high risk of Richter's transformation

Richter's Transformation (RT) is a poorly understood and fatal progression of chronic lymphocytic leukemia (CLL) manifesting histologically as diffuse large B-cell lymphoma. Protein arginine methyltransferase 5 (PRMT5) is implicated in lymphomagenesis, but its role in CLL or RT progression is unknown. We demonstrate herein that tumors uniformly overexpress PRMT5 in patients with progression to RT. Furthermore, mice with B-specific overexpression of hPRMT5 develop a B-lymphoid expansion with increased risk of death, and Eµ-PRMT5/TCL1 double transgenic mice develop a highly aggressive disease with transformation that histologically resembles RT; where large-scale transcriptional profiling identifies oncogenic pathways mediating PRMT5-driven disease progression. Lastly, we report the development of a SAM-competitive PRMT5 inhibitor, PRT382, with exclusive selectivity and optimal in vitro and in vivo activity compared to available PRMT5 inhibitors. Taken together, the discovery that PRMT5 drives oncogenic pathways promoting RT provides a compelling rationale for clinical investigation of PRMT5 inhibitors such as PRT382 in aggressive CLL/RT cases.

Analysis of Data on Fludarabine, Cyclophosphamide, and Rituximab Chemoimmunotherapy for Chronic Lymphocytic Leukemia Shows High Patient Heterogeneity and the Need for More Consideration of Individualized Treatment

Methods

Data from single-cell RNA sequencing (RNA-seq) of CLL patients were obtained from the Gene Expression Omnibus database. The R package was utilized to analyze the data, and the relation of results was predicted via the GeneMANIA website. The information of 7 samples covered three stages: observation stage, pretreatment by CIT with rituximab, fludarabine, and cyclophosphamide (pre-CIT), and post-CIT. The differentially expressed genes (DEGs) were identified, and functional enrichment analyses were performed. B cell subpopulations and pseudotime trajectories analysis was conducted.

Results

A total of 70,659 DEGs were identified. Each patient's DEGs presented their own characteristics, with low similarity. Therefore, it is difficult to identify potential hub genes. Similarly, pathway enrichment analysis showed significant tumor heterogeneity among CLL patients. Analysis of relapsed post-CIT compared to the observation stage suggested that the TP53 pathway should be taken seriously as it is closely related to treatment strategy and patient prognosis.

Conclusions

Tumor heterogeneity may be a more common manifestation of CLL. Individualized treatment should be considered for CLL. TP53 abnormality and its regulatory factors should still be the focus of CLL diagnosis and treatment.

In vivo G-CSF treatment activates the GR-SOCS1 axis to suppress IFN-γ secretion by natural killer cells

Natural killer (NK) cells are lymphocytes that are involved in controlling tumors or microbial infections through the production of interferon gamma (IFN-γ). Granulocyte colony-stimulating factor (G-CSF) inhibits IFN-γ secretion by NK cells, but the mechanism underlying this effect remains unclear. Here, by comparing the multi-omics profiles of human NK cells before and after in vivo G-CSF treatment, we identify a pathway that is activated in response to G-CSF treatment, which suppresses IFN-γ secretion in NK cells. Specifically, glucocorticoid receptors (GRs) activated by G-CSF inhibit secretion of IFN-γ by promoting interactions between SOCS1 promoters and enhancers, as well as increasing the expression of SOCS1. Experiments in mice confirm that G-CSF treatment significantly downregulates IFN-γ secretion and upregulates GR and SOCS1 expression in NK cells. In addition, GR blockade by the antagonist RU486 significantly reverses the effects of G-CSF, demonstrating that GRs upregulate SOCS1 and inhibit the production of IFN-γ by NK cells.

Multi-omics analysis of naïve B cells of patients harboring the C104R mutation in TACI

Common variable immunodeficiency (CVID) is the most prevalent form of symptomatic primary immunodeficiency in humans. The genetic cause of CVID is still unknown in about 70% of cases. Ten percent of CVID patients carry heterozygous mutations in the tumor necrosis factor receptor superfamily member 13B gene (TNFRSF13B), encoding TACI. Mutations in TNFRSF13B alone may not be sufficient for the development of CVID, as 1% of the healthy population carry these mutations. The common hypothesis is that TACI mutations are not fully penetrant and additional factors contribute to the development of CVID. To determine these additional factors, we investigated the perturbations of transcription factor (TF) binding and the transcriptome profiles in unstimulated and CD40L/IL21-stimulated naïve B cells from CVID patients harboring the C104R mutation in TNFRSF13B and compared them to their healthy relatives with the same mutation. In addition, the proteome of stimulated naïve B cells was investigated. For functional validation, intracellular protein concentrations were measured by flow cytometry. Our analysis revealed 8% less accessible chromatin in unstimulated naïve B cells and 25% less accessible chromatin in class-switched memory B cells from affected and unaffected TACI mutation carriers compared to healthy donors. The most enriched TF binding motifs in TACI mutation carriers involved members from the ETS, IRF, and NF-κB TF families. Validation experiments supported dysregulation of the NF-κB and MAPK pathways. In steady state, naïve B cells had increased cell death pathways and reduced cell metabolism pathways, while after stimulation, enhanced immune responses and decreased cell survival were detected. Using a multi-omics approach, our findings provide valuable insights into the impaired biology of naïve B cells from TACI mutation carriers.

Towards precision medicine in lymphoid malignancies

Careful histopathologic examination remains the cornerstone in the diagnosis of the clinically and biologically heterogeneous group of lymphoid malignancies. However, recent advances in genomic and epigenomic characterization using high-throughput technologies have significantly improved our understanding of these tumors. Although no single genomic alteration is completely specific for a lymphoma entity, some alterations are highly recurrent in certain entities and thus can provide complementary diagnostic information when integrated in the hematopathological diagnostic workup. Moreover, other alterations may provide important information regarding the clinical course, that is, prognostic or risk-stratifying markers, or response to treatment, that is, predictive markers, which may allow tailoring of the patient's treatment based on (epi)genetic characteristics. In this review, we will focus on clinically relevant diagnostic, prognostic, and predictive biomarkers identified in more common types of B-cell malignancies, and discuss how diagnostic assays designed for comprehensive molecular profiling may pave the way for the implementation of precision diagnostics/medicine approaches. We will also discuss future directions in this rapidly evolving field, including the application of single-cell sequencing and other omics technologies, to decipher clonal dynamics and evolution in lymphoid malignancies.

Drug-microenvironment perturbations reveal resistance mechanisms and prognostic subgroups in CLL

The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited and systematic analyses are lacking. Using chronic lymphocytic leukaemia (CLL) as a model disease, we investigated the influence of 17 microenvironmental stimuli on 12 drugs in 192 genetically characterised patient samples. Based on microenvironmental response, we identified four subgroups with distinct clinical outcomes beyond known prognostic markers. Response to multiple microenvironmental stimuli was amplified in trisomy 12 samples. Trisomy 12 was associated with a distinct epigenetic signature. Bromodomain inhibition reversed this epigenetic profile and could be used to target microenvironmental signalling in trisomy 12 CLL. We quantified the impact of microenvironmental stimuli on drug response and their dependence on genetic alterations, identifying interleukin 4 (IL4) and Toll‐like receptor (TLR) stimulation as the strongest actuators of drug resistance. IL4 and TLR signalling activity was increased in CLL‐infiltrated lymph nodes compared with healthy samples. High IL4 activity correlated with faster disease progression. The publicly available dataset can facilitate the investigation of cell‐extrinsic mechanisms of drug resistance and disease progression.

Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution

Myofibers are the main components of skeletal muscle, which is the largest tissue in the body. Myofibers are highly adaptive and can be altered under different biological and disease conditions. Therefore, transcriptional and epigenetic studies on myofibers are crucial to discover how chromatin alterations occur in the skeletal muscle under different conditions. However, due to the heterogenous nature of skeletal muscle, studying myofibers in isolation proves to be a challenging task. Single-cell sequencing has permitted the study of the epigenome of isolated myonuclei. While this provides sequencing with high dimensionality, the sequencing depth is lacking, which makes comparisons between different biological conditions difficult. Here, we report the first implementation of single myofiber ATAC-Seq, which allows for the sequencing of an individual myofiber at a depth sufficient for peak calling and for comparative analysis of chromatin accessibility under various physiological and disease conditions. Application of this technique revealed significant differences in chromatin accessibility between resting and regenerating myofibers, as well as between myofibers from a mouse model of Duchenne Muscular Dystrophy (mdx) and wild-type (WT) counterparts. This technique can lead to a wide application in the identification of chromatin regulatory elements and epigenetic mechanisms in muscle fibers during development and in muscle-wasting diseases.

scPower accelerates and optimizes the design of multi-sample single cell transcriptomic studies

Single cell RNA-seq has revolutionized transcriptomics by providing cell type resolution for differential gene expression and expression quantitative trait loci (eQTL) analyses. However, efficient power analysis methods for single cell data and inter-individual comparisons are lacking. Here, we present scPower; a statistical framework for the design and power analysis of multi-sample single cell transcriptomic experiments. We modelled the relationship between sample size, the number of cells per individual, sequencing depth, and the power of detecting differentially expressed genes within cell types. We systematically evaluated these optimal parameter combinations for several single cell profiling platforms, and generated broad recommendations. In general, shallow sequencing of high numbers of cells leads to higher overall power than deep sequencing of fewer cells. The model, including priors, is implemented as an R package and is accessible as a web tool. scPower is a highly customizable tool that experimentalists can use to quickly compare a multitude of experimental designs and optimize for a limited budget.

Core-binding factor leukemia hijacks the T-cell-prone PU.1 antisense promoter

The blood system serves as a key model for cell differentiation and cancer. It is orchestrated by precise spatiotemporal expression of crucial transcription factors. One of the key master regulators in the hematopoietic systems is PU.1. Reduced levels of PU.1 are characteristic for human acute myeloid leukemia (AML) and are known to induce AML in mouse models. Here, we show that transcriptional downregulation of PU.1 is an active process involving an alternative promoter in intron 3 that is induced by RUNX transcription factors driving noncoding antisense transcription. Core-binding factor (CBF) fusions RUNX1-ETO and CBFβ-MYH11 in t(8;21) and inv(16) AML, respectively, activate the PU.1 antisense promoter that results in a shift from sense toward antisense transcription and myeloid differentiation blockade. In patients with CBF-AML, we found that an elevated antisense/sense transcript and promoter accessibility ratio represents a hallmark compared with normal karyotype AML or healthy CD34+ cells. Competitive interaction of an enhancer with the proximal or the antisense promoter forms a binary on/off switch for either myeloid or T-cell development. Leukemic CBF fusions thus use a physiological mechanism used by T cells to decrease sense transcription. Our study is the first example of a sense/antisense promoter competition as a crucial functional switch for gene expression perturbation by oncogenes. Hence, this disease mechanism reveals a previously unknown Achilles heel for future precise therapeutic targeting of oncogene-induced chromatin remodeling.

A Functional Polymorphism in Accessible Chromatin Region Confers Risk of Non-Small Cell Lung Cancer in Chinese Population

Background

The disease-associated non-coding variants identified by genome-wide association studies (GWASs) were enriched in open chromatin regions (OCRs) and implicated in gene regulation. Genetic variants in OCRs thus may exert regulatory functions and contribute to non-small cell lung cancer (NSCLC) susceptibility.

Objective

To fine map potential functional variants in GWAS loci that contribute to NSCLC predisposition using chromatin accessibility and histone modification data and explore their functions by population study and biochemical experimental analyses.

Methods

We mapped the chromatin accessible regions of lung tissues using data of assay for transposase-accessible chromatin using sequencing (ATAC-seq) in The Cancer Genome Atlas (TCGA) and prioritized potential regulatory variants within lung cancer GWAS loci by aligning with histone signatures using data of chromatin immunoprecipitation assays followed by sequencing (ChIP-seq) in the Encyclopedia of DNA Elements (ENCODE). A two-stage case–control study with 1,830 cases and 2,001 controls was conducted to explore the associations between candidate variants and NSCLC risk in Chinese population. Bioinformatic annotations and biochemical experiments were performed to further reveal the potential functions of significant variants.

Results

Sixteen potential functional single-nucleotide polymorphisms (SNPs) were selected as candidates from bioinformatics analyses. Three variants out of the 16 candidate SNPs survived after genotyping in stage 1 case–control study, and only the results of SNP rs13064999 were successfully validated in the analyses of stage 2 case–control study. In combined analyses, rs13064999 was significantly associated with NSCLC risk [additive model; odds ratio (OR) = 1.17; 95%CI, 1.07–1.29; p = 0.001]. Functional annotations indicated its potential enhancer bioactivity, and dual-luciferase reporter assays revealed a significant increase in luciferase activity for the reconstructed plasmid with rs13064999 A allele, when compared to the one with wild-type G allele (p_A549 < 0.001, p_SK-MES-1 = 0.004). Further electrophoretic mobility shift assays (EMSA) and super-shift assays confirmed a stronger affinity of HP1γ for the binding motif containing SNP rs13064999 A allele.

Conclusion

These findings suggested that the functional variant rs13064999, identified by the integration of ATAC-seq and ChIP-seq data, contributes to the susceptibility of NSCLC by affecting HP1γ binding, while the exact biological mechanism awaits further exploration.

Chromatin-based, in cis and in trans regulatory rewiring underpins distinct oncogenic transcriptomes in multiple myeloma

Multiple myeloma is a genetically heterogeneous cancer of the bone marrow plasma cells (PC). Distinct myeloma transcriptome profiles are primarily driven by myeloma initiating events (MIE) and converge into a mutually exclusive overexpression of the CCND1 and CCND2 oncogenes. Here, with reference to their normal counterparts, we find that myeloma PC enhanced chromatin accessibility combined with paired transcriptome profiling can classify MIE-defined genetic subgroups. Across and within different MM genetic subgroups, we ascribe regulation of genes and pathways critical for myeloma biology to unique or shared, developmentally activated or de novo formed candidate enhancers. Such enhancers co-opt recruitment of existing transcription factors, which although not transcriptionally deregulated per se, organise aberrant gene regulatory networks that help identify myeloma cell dependencies with prognostic impact. Finally, we identify and validate the critical super-enhancer that regulates ectopic expression of CCND2 in a subset of patients with MM and in chronic lymphocytic leukemia.

Despite extensive genetic heterogeneity, nearly half of all multiple myeloma (MM) cases are driven by cyclin D2 (CCND2) over-expression. Here the authors dissect the chromatin landscape of MM to provide insights into the transcriptional regulatory landscape driving MM and divergent transcriptomes corresponding to different MM genetic subtypes.

Chronic lymphocytic leukemia (CLL) risk is mediated by multiple enhancer variants within CLL risk loci

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in Western countries. It has a strong genetic basis, showing a ~ 8-fold increased risk of CLL in first-degree relatives. Genome-wide association studies (GWAS) have identified 41 risk variants across 41 loci. However, for a majority of the loci, the functional variants and the mechanisms underlying their causal roles remain undefined. Here, we examined the genetic and epigenetic features associated with 12 index variants, along with any correlated (r2 ≥ 0.5) variants, at the CLL risk loci located outside of gene promoters. Based on publicly available ChIP-seq and chromatin accessibility data as well as our own ChIP-seq data from CLL patients, we identified six candidate functional variants at six loci and at least two candidate functional variants at each of the remaining six loci. The functional variants are predominantly located within enhancers or super-enhancers, including bi-directionally transcribed enhancers, which are often restricted to immune cell types. Furthermore, we found that, at 78% of the functional variants, the alternative alleles altered the transcription factor binding motifs or histone modifications, indicating the involvement of these variants in the change of local chromatin state. Finally, the enhancers carrying functional variants physically interacted with genes enriched in the type I interferon signaling pathway, apoptosis, or TP53 network that are known to play key roles in CLL. These results support the regulatory roles for inherited noncoding variants in the pathogenesis of CLL.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

CLL intraclonal fractions exhibit established and recently acquired patterns of DNA methylation

Intraclonal subpopulations of circulating chronic lymphocytic leukemia (CLL) cells with different proliferative histories and reciprocal surface expression of CXCR4 and CD5 have been observed in the peripheral blood of CLL patients and named proliferative (PF), intermediate (IF), and resting (RF) cellular fractions. Here, we found that these intraclonal circulating fractions share persistent DNA methylation signatures largely associated with the mutation status of the immunoglobulin heavy chain locus (IGHV) and their origins from distinct stages of differentiation of antigen-experienced B cells. Increased leukemic birth rate, however, showed a very limited impact on DNA methylation of circulating CLL fractions independent of IGHV mutation status. Additionally, DNA methylation heterogeneity increased as leukemic cells advanced from PF to RF in the peripheral blood. This frequently co-occurred with heterochromatin hypomethylation and hypermethylation of Polycomb-repressed regions in the PF, suggesting accumulation of longevity-associated epigenetic features in recently born cells. On the other hand, transcriptional differences between paired intraclonal fractions confirmed their proliferative experience and further supported a linear advancement from PF to RF in the peripheral blood. Several of these differentially expressed genes showed unique associations with clinical outcome not evident in the bulk clone, supporting the pathological and therapeutic relevance of studying intraclonal CLL fractions. We conclude that independent methylation and transcriptional landscapes reflect both preexisting cell-of-origin fingerprints and more recently acquired hallmarks associated with the life cycle of circulating CLL cells.

A genome-wide association study identifies susceptibility loci for primary central nervous system lymphoma at 6p25.3 and 3p22.1: a LOC Network study

BACKGROUND

Primary central nervous system lymphoma (PCNSL) is a rare form of extra-nodal non-Hodgkin lymphoma. PCNSL is a distinct subtype of non-Hodgkin lymphoma, with over 95% of tumors belonging to the diffuse large B-cell lymphoma (DLBCL) group. We have conducted a genome-wide association study (GWAS) on immunocompetent patients to address the possibility that common genetic variants influence the risk of developing PCNSL.

METHODS

We performed a meta-analysis of 2 new GWASs of PCNSL totaling 475 cases and 1134 controls of European ancestry. To increase genomic resolution, we imputed >10 million single nucleotide polymorphisms using the 1000 Genomes Project combined with UK10K as reference. In addition we performed a transcription factor binding disruption analysis and investigated the patterns of local chromatin by Capture Hi-C data.

RESULTS

We identified independent risk loci at 3p22.1 (rs41289586, ANO10, P = 2.17 × 10-8) and 6p25.3 near EXOC2 (rs116446171, P = 1.95 x 10-13). In contrast, the lack of an association between rs41289586 and DLBCL suggests distinct germline predisposition to PCNSL and DLBCL. We found looping chromatin interactions between noncoding regions at 6p25.3 (rs11646171) with the IRF4 promoter and at 8q24.21 (rs13254990) with the MYC promoter, both genes with strong relevance to B-cell tumorigenesis.

CONCLUSION

To our knowledge this is the first study providing insight into the genetic predisposition to PCNSL. Our findings represent an important step in defining the contribution of common genetic variation to the risk of developing PCNSL.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Systematic clustering algorithm for chromatin accessibility data and its application to hematopoietic cells

The huge amount of data acquired by high-throughput sequencing requires data reduction for effective analysis. Here we give a clustering algorithm for genome-wide open chromatin data using a new data reduction method. This method regards the genome as a string of 1s and 0s based on a set of peaks and calculates the Hamming distances between the strings. This algorithm with the systematically optimized set of peaks enables us to quantitatively evaluate differences between samples of hematopoietic cells and classify cell types, potentially leading to a better understanding of leukemia pathogenesis.

Chromatin accessibility landscapes of skin cells in systemic sclerosis nominate dendritic cells in disease pathogenesis

Systemic sclerosis (SSc) is a disease at the intersection of autoimmunity and fibrosis. However, the epigenetic regulation and the contributions of diverse cell types to SSc remain unclear. Here we survey, using ATAC-seq, the active DNA regulatory elements of eight types of primary cells in normal skin from healthy controls, as well as clinically affected and unaffected skin from SSc patients. We find that accessible DNA elements in skin-resident dendritic cells (DCs) exhibit the highest enrichment of SSc-associated single-nucleotide polymorphisms (SNPs) and predict the degrees of skin fibrosis in patients. DCs also have the greatest disease-associated changes in chromatin accessibility and the strongest alteration of cell-cell interactions in SSc lesions. Lastly, data from an independent cohort of patients with SSc confirm a significant increase of DCs in lesioned skin. Thus, the DCs epigenome links inherited susceptibility and clinically apparent fibrosis in SSc skin, and can be an important driver of SSc pathogenesis.

Evolution of the Epigenetic Landscape in Childhood B Acute Lymphoblastic Leukemia and Its Role in Drug Resistance

Although B-cell acute lymphoblastic leukemia (B-ALL) is the most common malignancy in children and while highly curable, it remains a leading cause of cancer-related mortality. The outgrowth of tumor subclones carrying mutations in genes responsible for resistance to therapy has led to a Darwinian model of clonal selection. Previous work has indicated that alterations in the epigenome might contribute to clonal selection, yet the extent to which the chromatin state is altered under the selective pressures of therapy is unknown. To address this, we performed chromatin immunoprecipitation, gene expression analysis, and enhanced reduced representation bisulfite sequencing on a cohort of paired diagnosis and relapse samples from individual patients who all but one relapsed within 36 months of initial diagnosis. The chromatin state at diagnosis varied widely among patients, while the majority of peaks remained stable between diagnosis and relapse. Yet a significant fraction was either lost or newly gained, with some patients showing few differences and others showing massive changes of the epigenetic state. Evolution of the epigenome was associated with pathways previously linked to therapy resistance as well as novel candidate pathways through alterations in pyrimidine biosynthesis and downregulation of polycomb repressive complex 2 targets. Three novel, relapse-specific superenhancers were shared by a majority of patients including one associated with S100A8, the top upregulated gene seen at relapse in childhood B-ALL. Overall, our results support a role of the epigenome in clonal evolution and uncover new candidate pathways associated with relapse. SIGNIFICANCE: This study suggests a major role for epigenetic mechanisms in driving clonal evolution in B-ALL and identifies novel pathways associated with drug resistance.

Environmental Epigenetics of Diesel Particulate Matter Toxicogenomics

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by disruptions in social communication and behavioral flexibility. Both genetic and environmental factors contribute to ASD risk. Epidemiologic studies indicate that roadway vehicle exhaust and in utero exposure to diesel particulate matter (DPM) are associated with ASD. Using the Comparative Toxicogenomics Database (CTD), we identified genes connected to DPM exposure and ASD, extracted the known enhancers/promoters of the identified genes, and integrated this with Assay for Transposase Accessible Chromatin (ATAC-seq) data from DPM-exposed human neural progenitor cells. Enhancer/promoter elements with significantly different chromosome accessibility revealed enriched DNA sequence motifs with transcription factor binding sites for EGR1. Variant extraction for linkage disequilibrium blocks of these regions followed by analysis through Genome Wide Association Studies (GWAS) revealed multiple neurological trait associations including exploratory eye movement and brain volume measurement. This approach highlights the effects of pollution on the regulatory regions of genes implicated in ASD by genetic studies, indicating convergence of genetic and environmental factors on molecular networks that contribute to ASD. Integration of publicly available data from the CTD, cell culture exposure studies, and phenotypic genetics synergize extensive evidence of chemical exposures on gene regulation for altered brain development.

Disturbed mitochondrial dynamics in CD8+ TILs reinforce T cell exhaustion

The metabolic challenges present in tumors attenuate the metabolic fitness and antitumor activity of tumor-infiltrating T lymphocytes (TILs). However, it remains unclear whether persistent metabolic insufficiency can imprint permanent T cell dysfunction. We found that TILs accumulated depolarized mitochondria as a result of decreased mitophagy activity and displayed functional, transcriptomic and epigenetic characteristics of terminally exhausted T cells. Mechanistically, reduced mitochondrial fitness in TILs was induced by the coordination of T cell receptor stimulation, microenvironmental stressors and PD-1 signaling. Enforced accumulation of depolarized mitochondria with pharmacological inhibitors induced epigenetic reprogramming toward terminal exhaustion, indicating that mitochondrial deregulation caused T cell exhaustion. Furthermore, supplementation with nicotinamide riboside enhanced T cell mitochondrial fitness and improved responsiveness to anti-PD-1 treatment. Together, our results reveal insights into how mitochondrial dynamics and quality orchestrate T cell antitumor responses and commitment to the exhaustion program.

Repurposing of KLF5 activates a cell cycle signature during the progression from a precursor state to oesophageal adenocarcinoma

Oesophageal adenocarcinoma (OAC) is one of the most common causes of cancer deaths. Barrett's oesophagus (BO) is the only known precancerous precursor to OAC, but our understanding about the molecular events leading to OAC development is limited. Here, we have integrated gene expression and chromatin accessibility profiles of human biopsies and identified a strong cell cycle gene expression signature in OAC compared to BO. Through analysing associated chromatin accessibility changes, we have implicated the transcription factor KLF5 in the transition from BO to OAC. Importantly, we show that KLF5 expression is unchanged during this transition, but instead, KLF5 is redistributed across chromatin to directly regulate cell cycle genes specifically in OAC cells. This new KLF5 target gene programme has potential prognostic significance as high levels correlate with poorer patient survival. Thus, the repurposing of KLF5 for novel regulatory activity in OAC provides new insights into the mechanisms behind disease progression.

Network Approaches for Dissecting the Immune System

The immune system is a complex biological network composed of hierarchically organized genes, proteins, and cellular components that combat external pathogens and monitor the onset of internal disease. To meet and ultimately defeat these challenges, the immune system orchestrates an exquisitely complex interplay of numerous cells, often with highly specialized functions, in a tissue-specific manner. One of the major methodologies of systems immunology is to measure quantitatively the components and interaction levels in the immunologic networks to construct a computational network and predict the response of the components to perturbations. The recent advances in high-throughput sequencing techniques have provided us with a powerful approach to dissecting the complexity of the immune system. Here we summarize the latest progress in integrating omics data and network approaches to construct networks and to infer the underlying signaling and transcriptional landscape, as well as cell-cell communication, in the immune system, with a focus on hematopoiesis, adaptive immunity, and tumor immunology. Understanding the network regulation of immune cells has provided new insights into immune homeostasis and disease, with important therapeutic implications for inflammation, cancer, and other immune-mediated disorders.

Epigenomic technologies for precision oncology

Epigenetic patterns in a cell control the expression of genes and consequently determine the phenotype of a cell. Cancer cells possess altered epigenomes which include aberrant patterns of DNA methylation, histone tail modifications, nucleosome positioning and of the three-dimensional chromatin organization within a nucleus. These altered epigenetic patterns are potential useful biomarkers to detect cancer cells and to classify tumor types. In addition, the cancer epigenome dictates the response of a cancer cell to therapeutic intervention and, therefore its knowledge, will allow to predict response to different therapeutic approaches. Here we review the current state-of-the-art technologies that have been developed to decipher epigenetic patterns on the genomic level and discuss how these methods are potentially useful for precision oncology.

Chromatin accessibility landscape of pediatric T-lymphoblastic leukemia and human T-cell precursors

We aimed at identifying the developmental stage at which leukemic cells of pediatric T-ALLs are arrested and at defining leukemogenic mechanisms based on ATAC-Seq. Chromatin accessibility maps of seven developmental stages of human healthy T cells revealed progressive chromatin condensation during T-cell maturation. Developmental stages were distinguished by 2,823 signature chromatin regions with 95% accuracy. Open chromatin surrounding SAE1 was identified to best distinguish thymic developmental stages suggesting a potential role of SUMOylation in T-cell development. Deconvolution using signature regions revealed that T-ALLs, including those with mature immunophenotypes, resemble the most immature populations, which was confirmed by TF-binding motif profiles. We integrated ATAC-Seq and RNA-Seq and found DAB1, a gene not related to leukemia previously, to be overexpressed, abnormally spliced and hyper-accessible in T-ALLs. DAB1-negative patients formed a distinct subgroup with particularly immature chromatin profiles and hyper-accessible binding sites for SPI1 (PU.1), a TF crucial for normal T-cell maturation. In conclusion, our analyses of chromatin accessibility and TF-binding motifs showed that pediatric T-ALL cells are most similar to immature thymic precursors, indicating an early developmental arrest.

Targeting Nuclear NOTCH2 by Gliotoxin Recovers a Tumor-Suppressor NOTCH3 Activity in CLL

NOTCH signaling represents a promising therapeutic target in chronic lymphocytic leukemia (CLL). We compared the anti-neoplastic effects of the nuclear NOTCH2 inhibitor gliotoxin and the pan-NOTCH γ-secretase inhibitor RO4929097 in primary CLL cells with special emphasis on the individual roles of the different NOTCH receptors. Gliotoxin rapidly induced apoptosis in all CLL cases tested, whereas RO4929097 exerted a variable and delayed effect on CLL cell viability. Gliotoxin-induced apoptosis was associated with inhibition of the NOTCH2/FCER2 (CD23) axis together with concomitant upregulation of the NOTCH3/NR4A1 axis. In contrast, RO4929097 downregulated the NOTCH3/NR4A1 axis and counteracted the spontaneous and gliotoxin-induced apoptosis. On the cell surface, NOTCH3 and CD23 expression were mutually exclusive, suggesting that downregulation of NOTCH2 signaling is a prerequisite for NOTCH3 expression in CLL cells. ATAC-seq confirmed that gliotoxin targeted the canonical NOTCH signaling, as indicated by the loss of chromatin accessibility at the potential NOTCH/CSL site containing the gene regulatory elements. This was accompanied by a gain in accessibility at the NR4A1, NFκB, and ATF3 motifs close to the genes involved in B-cell activation, differentiation, and apoptosis. In summary, these data show that gliotoxin recovers a non-canonical tumor-suppressing NOTCH3 activity, indicating that nuclear NOTCH2 inhibitors might be beneficial compared to pan-NOTCH inhibitors in the treatment of CLL.

Genomic and Epigenomic Alterations in Chronic Lymphocytic Leukemia

Chronic lymphocytic leukemia is a common disease in Western countries and has heterogeneous clinical behavior. The relevance of the genetic basis of the disease has come to the forefront recently, with genome-wide studies that have provided a comprehensive view of structural variants, somatic mutations, and different layers of epigenetic changes. The mutational landscape is characterized by relatively common copy number alterations, a few mutated genes occurring in 10-15% of cases, and a large number of genes mutated in a small number of cases. The epigenomic profile has revealed a marked reprogramming of regulatory regions in tumor cells compared with normal B cells. All of these alterations are differentially distributed in clinical and biological subsets of the disease, indicating that they may underlie the heterogeneous evolution of the disease. These global studies are revealing the molecular complexity of chronic lymphocytic leukemia and provide new perspectives that have helped to understand its pathogenic mechanisms and improve the clinical management of patients.

Clinical implications of chromatin accessibility in human cancers

Assay for transposase-accessible chromatin using sequencing (ATAC-seq) has not yet been widely used in cancer research. Clinical implications of chromatin accessibility assessed by ATAC-seq profiling in human cancers especially in a large patient cohort is largely unknown. In this study, we analyzed ATAC-seq data in 404 cancer patients from the Cancer Genome Atlas, representing the largest cancer patient cohort with ATAC-seq data, and correlated chromatin accessibility with patient demographics, tumor histology, molecular subtypes, and survival. Our results showed that chromatin accessibility varies from chromosome to chromosome, and is different in different genomic regions along the same chromosome. Chromatin accessibility especially on the X chromosome is strongly dependent on patient sex, but not much on patient age or tumor stage. Striking difference in chromatin accessibility is observed between lung adenocarcinoma and lung squamous cell carcinoma, the two most common histological subgroups in lung cancer. Furthermore, chromatin accessibility was different between basal and non-basal breast cancer. Finally, we identified prognostic peaks in the promoter regions that were significantly correlated with survival. In particular, we identified six peaks in the ESR1 gene promoter region in the ATAC-seq profiling and found that the peak about 247 bp away from the transcription start site was significantly associated with better survival. In conclusion, our study provides an alternative mechanism underlying tumor prognosis.

Combining signal and sequence to detect RNA polymerase initiation in ATAC-seq data

The assay for transposase-accessible chromatin followed by sequencing (ATAC-seq) is an inexpensive protocol for measuring open chromatin regions. ATAC-seq is also relatively simple and requires fewer cells than many other high-throughput sequencing protocols. Therefore, it is tractable in numerous settings where other high throughput assays are challenging to impossible. Hence it is important to understand the limits of what can be inferred from ATAC-seq data. In this work, we leverage ATAC-seq to predict the presence of nascent transcription. Nascent transcription assays are the current gold standard for identifying regions of active transcription, including markers for functional transcription factor (TF) binding. We combine mapped short reads from ATAC-seq with the underlying peak sequence, to determine regions of active transcription genome-wide. We show that a hybrid signal/sequence representation classified using recurrent neural networks (RNNs) can identify these regions across different cell types.

Chromatin-Based Classification of Genetically Heterogeneous AMLs into Two Distinct Subtypes with Diverse Stemness Phenotypes

Global investigation of histone marks in acute myeloid leukemia (AML) remains limited. Analyses of 38 AML samples through integrated transcriptional and chromatin mark analysis exposes 2 major subtypes. One subtype is dominated by patients with NPM1 mutations or MLL-fusion genes, shows activation of the regulatory pathways involving HOX-family genes as targets, and displays high self-renewal capacity and stemness. The second subtype is enriched for RUNX1 or spliceosome mutations, suggesting potential interplay between the 2 aberrations, and mainly depends on IRF family regulators. Cellular consequences in prognosis predict a relatively worse outcome for the first subtype. Our integrated profiling establishes a rich resource to probe AML subtypes on the basis of expression and chromatin data.

From reads to insight: a hitchhiker's guide to ATAC-seq data analysis

Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) is widely used in studying chromatin biology, but a comprehensive review of the analysis tools has not been completed yet. Here, we discuss the major steps in ATAC-seq data analysis, including pre-analysis (quality check and alignment), core analysis (peak calling), and advanced analysis (peak differential analysis and annotation, motif enrichment, footprinting, and nucleosome position analysis). We also review the reconstruction of transcriptional regulatory networks with multiomics data and highlight the current challenges of each step. Finally, we describe the potential of single-cell ATAC-seq and highlight the necessity of developing ATAC-seq specific analysis tools to obtain biologically meaningful insights.

Chromatin mapping and single-cell immune profiling define the temporal dynamics of ibrutinib response in CLL

The Bruton tyrosine kinase (BTK) inhibitor ibrutinib provides effective treatment for patients with chronic lymphocytic leukemia (CLL), despite extensive heterogeneity in this disease. To define the underlining regulatory dynamics, we analyze high-resolution time courses of ibrutinib treatment in patients with CLL, combining immune-phenotyping, single-cell transcriptome profiling, and chromatin mapping. We identify a consistent regulatory program starting with a sharp decrease of NF-κB binding in CLL cells, which is followed by reduced activity of lineage-defining transcription factors, erosion of CLL cell identity, and acquisition of a quiescence-like gene signature. We observe patient-to-patient variation in the speed of execution of this program, which we exploit to predict patient-specific dynamics in the response to ibrutinib based on the pre-treatment patient samples. In aggregate, our study describes time-dependent cellular, molecular, and regulatory effects for therapeutic inhibition of B cell receptor signaling in CLL, and it establishes a broadly applicable method for epigenome/transcriptome-based treatment monitoring.

Ibrutinib, a Bruton tyrosine kinase inhibitor, provides effective treatment for chronic lymphocytic leukemia (CLL). Here, the authors describe time-dependent molecular changes to malignant cells and to the immune system in patients undergoing ibrutinib therapy, with can be used for therapy monitoring.

Life-long impairment of glucose homeostasis upon prenatal exposure to psychostimulants

Maternal drug abuse during pregnancy is a rapidly escalating societal problem. Psychostimulants, including amphetamine, cocaine, and methamphetamine, are amongst the illicit drugs most commonly consumed by pregnant women. Neuropharmacology concepts posit that psychostimulants affect monoamine signaling in the nervous system by their affinities to neurotransmitter reuptake and vesicular transporters to heighten neurotransmitter availability extracellularly. Exacerbated dopamine signaling is particularly considered as a key determinant of psychostimulant action. Much less is known about possible adverse effects of these drugs on peripheral organs, and if in utero exposure induces lifelong pathologies. Here, we addressed this question by combining human RNA-seq data with cellular and mouse models of neuroendocrine development. We show that episodic maternal exposure to psychostimulants during pregnancy coincident with the intrauterine specification of pancreatic β cells permanently impairs their ability of insulin production, leading to glucose intolerance in adult female but not male offspring. We link psychostimulant action specifically to serotonin signaling and implicate the sex-specific epigenetic reprogramming of serotonin-related gene regulatory networks upstream from the transcription factor Pet1/Fev as determinants of reduced insulin production.

Life-long epigenetic programming of cortical architecture by maternal 'Western' diet during pregnancy

The evolution of human diets led to preferences toward polyunsaturated fatty acid (PUFA) content with 'Western' diets enriched in ω-6 PUFAs. Mounting evidence points to ω-6 PUFA excess limiting metabolic and cognitive processes that define longevity in humans. When chosen during pregnancy, ω-6 PUFA-enriched 'Western' diets can reprogram maternal bodily metabolism with maternal nutrient supply precipitating the body-wide imprinting of molecular and cellular adaptations at the level of long-range intercellular signaling networks in the unborn fetus. Even though unfavorable neurological outcomes are amongst the most common complications of intrauterine ω-6 PUFA excess, cellular underpinnings of life-long modifications to brain architecture remain unknown. Here, we show that nutritional ω-6 PUFA-derived endocannabinoids desensitize CB₁ cannabinoid receptors, thus inducing epigenetic repression of transcriptional regulatory networks controlling neuronal differentiation. We found that cortical neurons lose their positional identity and axonal selectivity when mouse fetuses are exposed to excess ω-6 PUFAs in utero. Conversion of ω-6 PUFAs into endocannabinoids disrupted the temporal precision of signaling at neuronal CB₁ cannabinoid receptors, chiefly deregulating Stat3-dependent transcriptional cascades otherwise required to execute neuronal differentiation programs. Global proteomics identified the immunoglobulin family of cell adhesion molecules (IgCAMs) as direct substrates, with DNA methylation and chromatin accessibility profiling uncovering epigenetic reprogramming at >1400 sites in neurons after prolonged cannabinoid exposure. We found anxiety and depression-like behavioral traits to manifest in adult offspring, which is consistent with genetic models of reduced IgCAM expression, to suggest causality for cortical wiring defects. Overall, our data uncover a regulatory mechanism whose disruption by maternal food choices could limit an offspring's brain function for life.

Improving Stratification for Children With Late Bone Marrow B-Cell Acute Lymphoblastic Leukemia Relapses With Refined Response Classification and Integration of Genetics

PURPOSE

Minimal residual disease (MRD) helps to accurately assess when children with late bone marrow relapses of B-cell precursor (BCP) acute lymphoblastic leukemia (ALL) will benefit from allogeneic hematopoietic stem-cell transplantation (allo-HSCT). More detailed dissection of MRD response heterogeneity and the specific genetic aberrations could improve current practice.

PATIENTS AND METHODS

MRD was assessed after induction treatment and at different times during relapse treatment until allo-HSCT (indicated in poor responders to induction; MRD ≥ 10−3) for patients being treated for late BCP-ALL bone marrow relapses (n = 413; median follow-up, 9.4 years) in the ALL-REZ BFM 2002 trial/registry (ClinicalTrials.gov identifier: NCT00114348 ).

RESULTS

Patients with both good (MRD < 10−3) and poor responses to induction treatment reached excellent event-free survival (EFS; 72% v 65%) and overall survival (OS; 82% v 74%). Patients with MRD of 10−2 or greater after induction had reduced EFS (56%), and their MRD persisted until allo-HSCT more frequently than it did in patients with MRD of 10−3 or greater to less than 10−2 ( P = .037). Patients with 25% or more leukemic blasts after induction (early nonresponders) had the poorest prognosis (EFS, 22%). Interestingly, patients with MRD of 10−3 or greater before allo-HSCT (late nonresponders) still had an EFS of 50% and OS of 63%, which in principle justifies allo-HSCT in these patients. From a panel of selected candidate genes, TP53 alterations (frequency, 8%) were the only genetic alteration with independent prognostic value in any MRD-based response subgroup.

CONCLUSION

After induction treatment, MRD-based treatment stratification resulted in excellent survival in patients with late relapsed BCP-ALL. Prognosis could be further improved in very poor responders by intensifying treatment directly after induction. TP53 alterations can be defined as a novel genetic high-risk marker in all MRD response groups in late relapsed BCP-ALL. Here we identified early and late nonresponders to be considered as events in future trials.