The Concerted Action of E2-2 and HEB Is Critical for Early Lymphoid Specification.

The apparition of adaptive immunity in Gnathostomata correlates with the expansion of the E-protein family to encompass E2-2, HEB, and E2A. Within the family, E2-2 and HEB are more closely evolutionarily related but their concerted action in hematopoiesis remains to be explored. Here we show that the combined disruption of E2-2 and HEB results in failure to express the early lymphoid program in Common lymphoid precursors (CLPs) and a near complete block in B-cell development. In the thymus, Early T-cell progenitors (ETPs) were reduced and T-cell development perturbed, resulting in reduced CD4 T- and increased γδ T-cell numbers. In contrast, hematopoietic stem cells (HSCs), erythro-myeloid progenitors, and innate immune cells were unaffected showing that E2-2 and HEB are dispensable for the ancestral hematopoietic lineages. Taken together, this E-protein dependence suggests that the appearance of the full Gnathostomata E-protein repertoire was critical to reinforce the gene regulatory circuits that drove the emergence and expansion of the lineages constituting humoral immunity.

High-throughput ChIPmentation: freely scalable, single day ChIPseq data generation from very low cell-numbers.

BACKGROUND: Chromatin immunoprecipitation coupled to sequencing (ChIP-seq) is widely used to map histone modifications and transcription factor binding on a genome-wide level. RESULTS: We present high-throughput ChIPmentation (HT-ChIPmentation) that eliminates the need for DNA purification prior to library amplification and reduces reverse-crosslinking time from hours to minutes. CONCLUSIONS: The resulting workflow is easily established, extremely rapid, and compatible with requirements for very low numbers of FACS sorted cells, high-throughput applications and single day data generation.

Ibrutinib induces rapid down-regulation of inflammatory markers and altered transcription of chronic lymphocytic leukaemia-related genes in blood and lymph nodes.

In chronic lymphocytic leukaemia (CLL) patients, treatment with the Bruton tyrosine kinase inhibitor ibrutinib induces a rapid shift of tumour cells from lymph nodes (LN) to peripheral blood (PB). Here, we characterized in depth the dynamics of ibrutinib-induced inflammatory, transcriptional and cellular changes in different compartments immediately after treatment initiation in seven relapsed/refractory CLL patients. Serial PB and LN samples were taken before start and during the first 29 days of treatment. Changes in plasma inflammation-related biomarkers, CLL cell RNA expression, B-cell activation and migration markers expression, and PB mononuclear cell populations were assessed. A significant reduction of 10 plasma inflammation markers, the majority of which were chemokines and not CLL-derived, was observed within hours, and was paralleled by very early increase of CD19+ circulating cells. At the RNA level, significant and continuous changes in transcription factors and signalling molecules linked to B-cell receptor signalling and CLL biology was observed in both PB and LN CLL cells already after 2 days of treatment. In conclusion, ibrutinib seems to instantly shut off an ongoing inflammatory response and interfere with diverse sensitive pathways in the LN.

Active enhancer and chromatin accessibility landscapes chart the regulatory network of primary multiple myeloma.

Multiple myeloma (MM) is an aggressive cancer that originates from antibody-secreting plasma cells. Although genetically and transcriptionally well characterized, the aberrant gene regulatory networks that underpin this disease remain poorly understood. Here, we mapped regulatory elements, open chromatin, and transcription factor (TF) footprints in primary MM cells. In comparison with normal antibody-secreting cells, MM cells displayed consistent changes in enhancer activity that are connected to superenhancer (SE)-mediated deregulation of TF genes. MM cells also displayed widespread decompaction of heterochromatin that was associated with activation of regulatory elements and in a major subset of patients' deregulation of the cyclic adenosine monophosphate pathway. Finally, building SE-associated TF-based regulatory networks allowed identification of several novel TFs that are central to MM biology. Taken together, these findings significantly add to our understanding of the aberrant gene regulatory network that underpins MM.

Mice deficient of Myc super-enhancer region reveal differential control mechanism between normal and pathological growth.

The gene desert upstream of the MYC oncogene on chromosome 8q24 contains susceptibility loci for several major forms of human cancer. The region shows high conservation between human and mouse and contains multiple MYC enhancers that are activated in tumor cells. However, the role of this region in normal development has not been addressed. Here we show that a 538 kb deletion of the entire MYC upstream super-enhancer region in mice results in 50% to 80% decrease in Myc expression in multiple tissues. The mice are viable and show no overt phenotype. However, they are resistant to tumorigenesis, and most normal cells isolated from them grow slowly in culture. These results reveal that only cells whose MYC activity is increased by serum or oncogenic driver mutations depend on the 8q24 super-enhancer region, and indicate that targeting the activity of this element is a promising strategy of cancer chemoprevention and therapy.

Comprehensive mapping of the effects of azacitidine on DNA methylation, repressive/permissive histone marks and gene expression in primary cells from patients with MDS and MDS-related disease.

Azacitidine (Aza) is first-line treatment for patients with high-risk myelodysplastic syndromes (MDS), although its precise mechanism of action is unknown. We performed the first study to globally evaluate the epigenetic effects of Aza on MDS bone marrow progenitor cells assessing gene expression (RNA seq), DNA methylation (Illumina 450k) and the histone modifications H3K18ac and H3K9me3 (ChIP seq). Aza induced a general increase in gene expression with 924 significantly upregulated genes but this increase showed no correlation with changes in DNA methylation or H3K18ac, and only a weak association with changes in H3K9me3. Interestingly, we observed activation of transcripts containing 15 endogenous retroviruses (ERVs) confirming previous cell line studies. DNA methylation decreased moderately in 99% of all genes, with a median ß-value reduction of 0.018; the most pronounced effects seen in heterochromatin. Aza-induced hypomethylation correlated significantly with change in H3K9me3. The pattern of H3K18ac and H3K9me3 displayed large differences between patients and healthy controls without any consistent pattern induced by Aza. We conclude that the marked induction of gene expression only partly could be explained by epigenetic changes, and propose that activation of ERVs may contribute to the clinical effects of Aza in MDS.

miSTAR: miRNA target prediction through modeling quantitative and qualitative miRNA binding site information in a stacked model structure.

In microRNA (miRNA) target prediction, typically two levels of information need to be modeled: the number of potential miRNA binding sites present in a target mRNA and the genomic context of each individual site. Single model structures insufficiently cope with this complex training data structure, consisting of feature vectors of unequal length as a consequence of the varying number of miRNA binding sites in different mRNAs. To circumvent this problem, we developed a two-layered, stacked model, in which the influence of binding site context is separately modeled. Using logistic regression and random forests, we applied the stacked model approach to a unique data set of 7990 probed miRNA-mRNA interactions, hereby including the largest number of miRNAs in model training to date. Compared to lower-complexity models, a particular stacked model, named miSTAR (miRNA stacked model target prediction; www.mi-star.org), displays a higher general performance and precision on top scoring predictions. More importantly, our model outperforms published and widely used miRNA target prediction algorithms. Finally, we highlight flaws in cross-validation schemes for evaluation of miRNA target prediction models and adopt a more fair and stringent approach.

Ebf1 heterozygosity results in increased DNA damage in pro-B cells and their synergistic transformation by Pax5 haploinsufficiency.

Early B-cell factor 1 (Ebf1) is a transcription factor with documented dose-dependent functions in normal and malignant B-lymphocyte development. To understand more about the roles of Ebf1 in malignant transformation, we investigated the impact of reduced functional Ebf1 dosage on mouse B-cell progenitors. Gene expression analysis suggested that Ebf1 was involved in the regulation of genes important for DNA repair and cell survival. Investigation of the DNA damage in steady state, as well as after induction of DNA damage by UV light, confirmed that pro-B cells lacking 1 functional allele of Ebf1 display signs of increased DNA damage. This correlated to reduced expression of DNA repair genes including Rad51, and chromatin immunoprecipitation data suggested that Rad51 is a direct target for Ebf1. Although reduced dosage of Ebf1 did not significantly increase tumor formation in mice, a dramatic increase in the frequency of pro-B cell leukemia was observed in mice with combined heterozygous mutations in the Ebf1 and Pax5 genes, revealing a synergistic effect of combined dose reduction of these proteins. Our data suggest that Ebf1 controls DNA repair in a dose-dependent manner providing a possible explanation to the frequent involvement of EBF1 gene loss in human leukemia.

SNP-guided identification of monoallelic DNA-methylation events from enrichment-based sequencing data.

Monoallelic gene expression is typically initiated early in the development of an organism. Dysregulation of monoallelic gene expression has already been linked to several non-Mendelian inherited genetic disorders. In humans, DNA-methylation is deemed to be an important regulator of monoallelic gene expression, but only few examples are known. One important reason is that current, cost-affordable truly genome-wide methods to assess DNA-methylation are based on sequencing post-enrichment. Here, we present a new methodology based on classical population genetic theory, i.e. the Hardy-Weinberg theorem, that combines methylomic data from MethylCap-seq with associated SNP profiles to identify monoallelically methylated loci. Applied on 334 MethylCap-seq samples of very diverse origin, this resulted in the identification of 80 genomic regions featured by monoallelic DNA-methylation. Of these 80 loci, 49 are located in genic regions of which 25 have already been linked to imprinting. Further analysis revealed statistically significant enrichment of these loci in promoter regions, further establishing the relevance and usefulness of the method. Additional validation was done using both 14 whole-genome bisulfite sequencing data sets and 16 mRNA-seq data sets. Importantly, the developed approach can be easily applied to other enrichment-based sequencing technologies, like the ChIP-seq-based identification of monoallelic histone modifications.