Benchmarking of 4C-seq pipelines based on real and simulated data

"jats:sec" "jats:title"Motivation"/jats:title" "jats:p"With its capacity for high-resolution data output in one region of interest, chromosome conformation capture combined with high-throughput sequencing (4C-seq) is a state-of-the-art next-generation sequencing technique that provides epigenetic insights, and regularly advances current medical research. However, 4C-seq data are complex and prone to biases, and while specialized programs exist, an unbiased, extensive benchmarking is still lacking. Furthermore, neither substantial datasets with fully characterized ground truth, nor simulation programs for realistic 4C-seq data have been published."/jats:p" "/jats:sec" "jats:sec" "jats:title"Results"/jats:title" "jats:p"We conducted a benchmarking study on 66 4C-seq samples from 20 datasets, and developed a novel 4C-seq simulation software, Basic4CSim, to allow for detailed comparisons of 4C-seq algorithms on 50 simulated datasets with 10–120 samples each. Simulations and benchmarking were adapted to address different characteristics of 4C-seq data. Simulated data were compared with published samples to validate simulation settings. We identified differences between 4C-seq algorithms in terms of precision, recall, interaction structure, and run time, and observed general trends. Novel differential pipeline versions of single-sample based 4C-seq algorithms were included in the benchmarking. While no single tool was optimally suited for both near-cis and far-cis, and both single-sample and differential analyses, choosing a high-performing algorithm variant did improve results considerably. For near-cis scenarios, r3Cseq, peakC and FourCSeq offered high precision, while fourSig demonstrated high overall F1 scores in far-cis analyses. Finally, 4C-seq simulations may aid in the development of improved analysis algorithms."/jats:p" "/jats:sec" "jats:sec" "jats:title"Availability and implementation"/jats:title" "jats:p"Basic4CSim is available at https://github.com/walter–ca/Basic4CSim."/jats:p" "/jats:sec" "jats:sec" "jats:title"Supplementary information"/jats:title" "jats:p"Supplementary data are available at Bioinformatics online."/jats:p" "/jats:sec Document type: Articl

Hemojuvelin-Mediated Hepcidin Induction Requires Both Bone Morphogenetic Protein Type I Receptors ALK2 and ALK3

Hemojuvelin (HJV) is a glycosylphosphatidylinositol-anchored protein of the repulsive guidance molecule family acting as a bone morphogenetic protein (BMP) coreceptor to induce the hepatic iron regulatory protein hepcidin. Hepcidin causes ubiquitination and degradation of the sole known iron exporter ferroportin, thereby limiting iron availability. The detailed signaling mechanism of HJV in vivo has yet to be investigated. In the current manuscript, we used an established model of adeno-associated virus (AAV)-mediated liver-specific overexpression of HJV in murine models of hepatocyte-specific deficiency of the BMP type I receptors Alk2 or Alk3. In control mice, HJV overexpression increased hepatic Hamp messenger RNA (mRNA) levels, soluble HJV (sHJV), splenic iron content (SIC), as well as phosphorylated small mothers against decapentaplegic protein (pSMAD1/5/8) levels. In contrast, in Alk2fl/fl;Alb-Cre and Alk3fl/fl;Alb-Cre mice, which present with moderate and severe iron overload, respectively, the administration of AAV-HJV induced HJV and sHJV. However, it did not rescue the iron overload phenotypes of those mice. Serum iron levels were induced in Alk2fl/fl;Alb-Cre mice after HJV overexpression. In phosphate-buffered saline-injected Alk3fl/fl;Alb-Cre mice, serum iron levels and the expression of duodenal ferroportin remained high, whereas Hamp mRNA levels were decreased to 1% to 5% of the levels detected in controls. This was reduced even further by AAV-HJV overexpression. SIC remained low in mice with hepatocyte-specific Alk2 or Alk3 deficiency, reflecting disturbed iron homeostasis with high serum iron levels and transferrin saturation and an inability to induce hepcidin by HJV overexpression. The data indicate that ALK2 and ALK3 are both required in vivo for the HJV-mediated induction of hepcidin

A new source of suprathermal O+ions near the dayside polar cap boundary

A new dayside source of O+ ions for the polar magnetosphere is described, and a statistical survey presented of upward flows of O+ ions using 2 years of data from the retarding ion mass spectrometer (RIMS) experiment on board DE 1, at geocentric distances below 3 RE and invariant latitudes above 40°. The flows are classified according to their spin angle distributions. It is believed that the spacecraft potential near perigee is generally less than +2 V, in which case the entire O+ population at energies below about 60 eV is sampled. Examples are given of field-aligned flow and of transversely accelerated “core” O+ ions; in the latter events a large fraction of the total O+ ion population has been transversely accelerated, and in some extreme cases all the observed ions (of all ion species) have been accelerated, and no residual cold population is observed (“toroidal” distributions). However, by far the most common type of O+ upflow seen by DE RIMS lies near the dayside polar cap boundary (particularly in the prenoon sector) and displays an asymmetric spin angle distribution. In such events the ions carry an upward heat flux, and strong upflow of all species is present (H+, He+, O+, O++, and N+ have all been observed with energies up to about 30 eV, but with the majority of ions below about 2 eV); hence, these have been termed upwelling ion events. The upwelling ions are embedded in larger regions of classical light ion polar wind and are persistently found under the following conditions: at geocentric distances greater than 1.4 RE; at all Kp in summer, but only at high Kp in winter. Low-energy conical ions (<30 eV) are only found near the equatorial edge of the events, the latitude of which moves equatorward with increasing Kp and is highly correlated with the location of field-aligned currents. The RIMS data are fully consistent with a “mass spectrometer effect,” whereby light ions and the more energetic O+ ions flow into the lobes and mantle and hence the far-tail plasma sheet, but lower-energy O+ is swept across the polar cap by the convection electric field, potentially acting as a source for the nightside auroral acceleration regions. The occurrence probability of upwelling ion events, as compared to those of low-altitude transversely accelerated core ions and of field-aligned flow, suggests this could be the dominant mechanism for supplying the nightside auroral acceleration region, and subsequently the ring current and near-earth plasma sheet, with ionospheric O+ ions. It is shown that the total rate of O+ outflow in upwelling ion events (greater than 10^25 s^{−1}) is sufficient for the region near the dayside polar cap boundary to be an important ionospheric heavy ion source

Cross talk between Wnt/β-catenin and Irf8 in leukemia progression and drug resistance

Progression and disease relapse of chronic myeloid leukemia (CML) depends on leukemia-initiating cells (LIC) that resist treatment. Using mouse genetics and a BCR-ABL model of CML, we observed cross talk between Wnt/β-catenin signaling and the interferon-regulatory factor 8 (Irf8). In normal hematopoiesis, activation of β-catenin results in up-regulation of Irf8, which in turn limits oncogenic β-catenin functions. Self-renewal and myeloproliferation become dependent on β-catenin in Irf8-deficient animals that develop a CML-like disease. Combined Irf8 deletion and constitutive β-catenin activation result in progression of CML into fatal blast crisis, elevated leukemic potential of BCR-ABL–induced LICs, and Imatinib resistance. Interestingly, activated β-catenin enhances a preexisting Irf8-deficient gene signature, identifying β-catenin as an amplifier of progression-specific gene regulation in the shift of CML to blast crisis. Collectively, our data uncover Irf8 as a roadblock for β-catenin–driven leukemia and imply both factors as targets in combinatorial therapy

Phosphorylation of Serine 248 of C/EBPα Is Dispensable for Myelopoiesis but Its Disruption Leads to a Low Penetrant Myeloid Disorder with Long Latency

BACKGROUND: Transcription factors play a key role in lineage commitment and differentiation of stem cells into distinct mature cells. In hematopoiesis, they regulate lineage-specific gene expression in a stage-specific manner through various physical and functional interactions with regulatory proteins that are simultanously recruited and activated to ensure timely gene expression. The transcription factor CCAAT/enhancer binding protein α (C/EBPα) is such a factor and is essential for the development of granulocytic/monocytic cells. The activity of C/EBPα is regulated on several levels including gene expression, alternative translation, protein interactions and posttranslational modifications, such as phosphorylation. In particular, the phosphorylation of serine 248 of the transactivation domain has been shown to be of crucial importance for granulocytic differentiation of 32Dcl3 cells in vitro. METHODOLOGY/PRINCIPAL FINDINGS: Here, we use mouse genetics to investigate the significance of C/EBPα serine 248 in vivo through the construction and analysis of Cebpa(S248A/S248A) knock-in mice. Surprisingly, 8-week old Cebpa(S248A/S248A) mice display normal steady-state hematopoiesis including unaltered development of mature myeloid cells. However, over time some of the animals develop a hematopoietic disorder with accumulation of multipotent, megakaryocytic and erythroid progenitor cells and a mild impairment of differentiation along the granulocytic-monocytic lineage. Furthermore, BM cells from Cebpa(S248A/S248A) animals display a competitive advantage compared to wild type cells in a transplantation assay. CONCLUSIONS/SIGNIFICANCE: Taken together, our data shows that the substitution of C/EBPα serine 248 to alanine favors the selection of the megakaryocytic/erythroid lineage over the monocytic/granulocytic compartment in old mice and suggests that S248 phosphorylation may be required to maintain proper hematopoietic homeostasis in response to changes in the wiring of cellular signalling networks. More broadly, the marked differences between the phenotype of the S248A variant in vivo and in vitro highlight the need to exert caution when extending in vitro phenotypes to the more appropriate in vivo context

Uncovering a new cellular origin for acute myeloid leukemia with lineage plasticity

Although acute myeloid leukaemia (AML) is assumed to be driven by transformed haematopoietic stem and progenitor cells, we have described recently a new pathway leading to AML from T-cell progenitors. Furthermore, we could identify a subgroup of human AML with gene expression profile suggesting T-lymphoid origin and potentially novel treatment