Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

Non-small cell lung carcinoma (NSCLC) is a complex malignancy that owing to its heterogeneity and poor prognosis poses many challenges to diagnosis, prognosis and patient treatment. DNA methylation is an important mechanism of epigenetic regulation involved in normal development and cancer. It is a very stable and specific modification and therefore in principle a very suitable marker for epigenetic phenotyping of tumors. Here we present a genome-wide DNA methylation analysis of NSCLC samples and paired lung tissues, where we combine MethylCap and next generation sequencing (MethylCap-seq) to provide comprehensive DNA methylation maps of the tumor and paired lung samples. The MethylCap-seq data were validated by bisulfite sequencing and methyl-specific polymerase chain reaction of selected regions.

RiceGeneThresher is a public online resource for mining genes underlying genome regions of interest or quantitative trait loci (QTL) in rice genome. It is a compendium of rice genomic resources consisting of genetic markers, genome annotation, expressed sequence tags (ESTs), protein domains, gene ontology, plant stress-responsive genes, metabolic pathways and prediction of protein-protein interactions. RiceGeneThresher system integrates these diverse data sources and provides powerful web-based applications, and flexible tools for delivering customized set of biological data on rice. Its system supports whole-genome gene mining for QTL by querying using DNA marker intervals or genomic loci. RiceGeneThresher provides biologically supported evidences that are essential for targeting groups or networks of genes involved in controlling traits underlying QTL. Users can use it to discover and to assign the most promising candidate genes in preparation for the further gene function validation analysis. The web-based application is freely available at http://rice.kps.ku.ac.th.

Alternative splicing is an important source of heterogeneity underlying gene expression between individual cells but remains an understudied area due to the paucity of computational tools to analyze splicing dynamics at single-cell resolution. Here, we present MARVEL, a comprehensive R package for single-cell splicing analysis applicable to RNA sequencing generated from the plate- and droplet-based methods. We performed extensive benchmarking of MARVEL against available tools and demonstrated its utility by analyzing multiple publicly available datasets in diverse cell types, including in disease. MARVEL enables systematic and integrated splicing and gene expression analysis of single cells to characterize the splicing landscape and reveal biological insights.

Myelofibrosis is a severe myeloproliferative neoplasm characterized by increased numbers of abnormal bone marrow megakaryocytes that induce fibrosis, destroying the hematopoietic microenvironment. To determine the cellular and molecular basis for aberrant megakaryopoiesis in myelofibrosis, we performed single-cell transcriptome profiling of 135,929 CD34+ lineage- hematopoietic stem and progenitor cells (HSPCs), single-cell proteomics, genomics, and functional assays. We identified a bias toward megakaryocyte differentiation apparent from early multipotent stem cells in myelofibrosis and associated aberrant molecular signatures. A sub-fraction of myelofibrosis megakaryocyte progenitors (MkPs) are transcriptionally similar to healthy-donor MkPs, but the majority are disease specific, with distinct populations expressing fibrosis- and proliferation-associated genes. Mutant-clone HSPCs have increased expression of megakaryocyte-associated genes compared to wild-type HSPCs, and we provide early validation of G6B as a potential immunotherapy target. Our study paves the way for selective targeting of the myelofibrosis clone and illustrates the power of single-cell multi-omics to discover tumor-specific therapeutic targets and mediators of tissue fibrosis.

Acute erythroid leukemia (AEL) commonly involves both myeloid and erythroid lineage transformation. However, the mutations that cause AEL and the cell(s) that sustain the bilineage leukemia phenotype remain unknown. We here show that combined biallelic Cebpa and Gata2 zinc finger-1 (ZnF1) mutations cooperatively induce bilineage AEL, and that the major leukemia-initiating cell (LIC) population has a neutrophil-monocyte progenitor (NMP) phenotype. In pre-leukemic NMPs Cebpa and Gata2 mutations synergize by increasing erythroid transcription factor (TF) expression and erythroid TF chromatin access, respectively, thereby installing ectopic erythroid potential. This erythroid-permissive chromatin conformation is retained in bilineage LICs. These results demonstrate that synergistic transcriptional and epigenetic reprogramming by leukemia-initiating mutations can generate neomorphic pre-leukemic progenitors, defining the lineage identity of the resulting leukemia.

CD19-directed immunotherapies have revolutionized the treatment of advanced B-cell acute lymphoblastic leukemia (B-ALL). Despite initial impressive rates of complete remission (CR) many patients ultimately relapse. Patients with B-ALL successfully treated with CD19-directed T cells eventually relapse, which, coupled with the early onset of CD22 expression during B-cell development, suggests that preexisting CD34+CD22+CD19- (pre)-leukemic cells represent an "early progenitor origin-related" mechanism underlying phenotypic escape to CD19-directed immunotherapies. We demonstrate that CD22 expression precedes CD19 expression during B-cell development. CD34+CD19-CD22+ cells are found in diagnostic and relapsed bone marrow samples of ∼70% of patients with B-ALL, and their frequency increases twofold in patients with B-ALL in CR after CD19 CAR T-cell therapy. The median of CD34+CD19-CD22+ cells before treatment was threefold higher in patients in whom B-ALL relapsed after CD19-directed immunotherapy (median follow-up, 24 months). Fluorescence in situ hybridization analysis in flow-sorted cell populations and xenograft modeling revealed that CD34+CD19-CD22+ cells harbor the genetic abnormalities present at diagnosis and initiate leukemogenesis in vivo. Our data suggest that preleukemic CD34+CD19-CD22+ progenitors underlie phenotypic escape after CD19-directed immunotherapies and reinforce ongoing clinical studies aimed at CD19/CD22 dual targeting as a strategy for reducing CD19- relapses. The implementation of CD34/CD19/CD22 immunophenotyping in clinical laboratories for initial diagnosis and subsequent monitoring of patients with B-ALL during CD19-targeted therapy is encouraged.

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr+ perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis.

Current therapies for myeloproliferative neoplasms (MPNs) improve symptoms but have limited effect on tumor size. In preclinical studies, tamoxifen restored normal apoptosis in mutated hematopoietic stem/progenitor cells (HSPCs). TAMARIN Phase-II, multicenter, single-arm clinical trial assessed tamoxifen's safety and activity in patients with stable MPNs, no prior thrombotic events and mutated JAK2V617F, CALRins5 or CALRdel52 peripheral blood allele burden ≥20% (EudraCT 2015-005497-38). 38 patients were recruited over 112w and 32 completed 24w-treatment. The study's A'herns success criteria were met as the primary outcome ( ≥ 50% reduction in mutant allele burden at 24w) was observed in 3/38 patients. Secondary outcomes included ≥25% reduction at 24w (5/38), ≥50% reduction at 12w (0/38), thrombotic events (2/38), toxicities, hematological response, proportion of patients in each IWG-MRT response category and ELN response criteria. As exploratory outcomes, baseline analysis of HSPC transcriptome segregates responders and non-responders, suggesting a predictive signature. In responder HSPCs, longitudinal analysis shows high baseline expression of JAK-STAT signaling and oxidative phosphorylation genes, which are downregulated by tamoxifen. We further demonstrate in preclinical studies that in JAK2V617F+ cells, 4-hydroxytamoxifen inhibits mitochondrial complex-I, activates integrated stress response and decreases pathogenic JAK2-signaling. These results warrant further investigation of tamoxifen in MPN, with careful consideration of thrombotic risk.

Regulation of immunoglobulin (Ig) V(D)J gene rearrangement is dependent on higher-order chromatin organization. Here, we studied the in vivo function of the DNA-binding zinc-finger protein CTCF, which regulates interactions between enhancers and promoters. By conditional deletion of the Ctcf gene in the B cell lineage, we demonstrate that loss of CTCF allowed Ig heavy chain recombination, but pre-B cell proliferation and differentiation was severely impaired. In the absence of CTCF, the Igκ light chain locus showed increased proximal and reduced distal Vκ usage. This was associated with enhanced proximal Vκ and reduced Jκ germline transcription. Chromosome conformation capture experiments demonstrated that CTCF limits interactions of the Igκ enhancers with the proximal V(κ) gene region and prevents inappropriate interactions between these strong enhancers and elements outside the Igκ locus. Thus, although Ig gene recombination can occur in the absence of CTCF, it is a critical factor determining Vκ segment choice for recombination.

Thymic T cell development is initiated from bone-marrow-derived multi potent thymus-seeding progenitors. During the early stages of thymocyte differentiation, progenitors become T cell restricted. However, the cellular environments supporting these critical initial stages of T cell development within the thymic cortex are not known. Here we use the dependence of early, c-Kit-expressing thymic progenitors on Kit ligand (KitL) to show that CD4(-)CD8(-)c-Kit(+)CD25(-) DN1-stage progenitors associate with, and depend on, the membrane-bound form of KitL (mKitL) provided by a cortex-specific KitL-expressing vascular endothelial cell (VEC) population. In contrast, the subsequent CD4(-)CD8(-)c-Kit(+)CD25(+) DN2-stage progenitors associate selectively with cortical thymic epithelial cells (cTECs) and depend on cTEC-presented mKitL. These results show that the dynamic process of early thymic progenitor differentiation is paralleled by migration-dependent change to the supporting niche, and identify VECs as a thymic niche cell, with mKitL as a critical ligand.

The coupling of chromosome conformation capture (3C) with next-generation sequencing technologies enables the high-throughput detection of long-range genomic interactions, via the generation of ligation products between DNA sequences, which are closely juxtaposed in vivo. These interactions involve promoter regions, enhancers and other regulatory and structural elements of chromosomes and can reveal key details of the regulation of gene expression. 3C-seq is a variant of the method for the detection of interactions between one chosen genomic element (viewpoint) and the rest of the genome. We present r3Cseq, an R/Bioconductor package designed to perform 3C-seq data analysis in a number of different experimental designs. The package reads a common aligned read input format, provides data normalization, allows the visualization of candidate interaction regions and detects statistically significant chromatin interactions, thus greatly facilitating hypothesis generation and the interpretation of experimental results. We further demonstrate its use on a series of real-world applications.

The neural crest (NC) is a transient embryonic stem cell-like population characterized by its multipotency and broad developmental potential. Here, we perform NC-specific transcriptional and epigenomic profiling of foxd3-mutant cells in vivo to define the gene regulatory circuits controlling NC specification. Together with global binding analysis obtained by foxd3 biotin-ChIP and single cell profiles of foxd3-expressing premigratory NC, our analysis shows that, during early steps of NC formation, foxd3 acts globally as a pioneer factor to prime the onset of genes regulating NC specification and migration by re-arranging the chromatin landscape, opening cis-regulatory elements and reshuffling nucleosomes. Strikingly, foxd3 then gradually switches from an activator to its well-described role as a transcriptional repressor and potentially uses differential partners for each role. Taken together, these results demonstrate that foxd3 acts bimodally in the neural crest as a switch from "permissive" to "repressive" nucleosome and chromatin organization to maintain multipotency and define cell fates.

Somatic mutations in acute myeloid leukemia are acquired sequentially and hierarchically. First, pre-leukemic mutations, such as t(8;21) that encodes AML1-ETO, are acquired within the hematopoietic stem cell (HSC) compartment, while signaling pathway mutations, including KRAS activating mutations, are late events acquired during transformation of leukemic progenitor cells and are rarely detectable in HSC. This raises the possibility that signaling pathway mutations are detrimental to clonal expansion of pre-leukemic HSC. To address this hypothesis, we used conditional genetics to introduce Aml1-ETO and K-RasG12D into murine HSC, either individually or in combination. In the absence of activated Ras, Aml1-ETO-expressing HSC conferred a competitive advantage. However, activated K-Ras had a marked detrimental effect on Aml1-ETO-expressing HSC, leading to loss of both phenotypic and functional HSC. Cell cycle analysis revealed a loss of quiescence in HSC co-expressing Aml1-ETO and K-RasG12D, accompanied by an enrichment in E2F and Myc target gene expression and depletion of HSC self-renewal-associated gene expression. These findings provide a mechanistic basis for the observed absence of KRAS signaling mutations in the pre-malignant HSC compartment.

Although previous studies suggested that the expression of FMS-like tyrosine kinase 3 (Flt3) initiates downstream of mouse hematopoietic stem cells (HSCs), FLT3 internal tandem duplications (FLT3 ITDs) have recently been suggested to intrinsically suppress HSCs. Herein, single-cell interrogation found Flt3 mRNA expression to be absent in the large majority of phenotypic HSCs, with a strong negative correlation between Flt3 and HSC-associated gene expression. Flt3-ITD knock-in mice showed reduced numbers of phenotypic HSCs, with an even more severe loss of long-term repopulating HSCs, likely reflecting the presence of non-HSCs within the phenotypic HSC compartment. Competitive transplantation experiments established that Flt3-ITD compromises HSCs through an extrinsically mediated mechanism of disrupting HSC-supporting bone marrow stromal cells, with reduced numbers of endothelial and mesenchymal stromal cells showing increased inflammation-associated gene expression. Tumor necrosis factor (TNF), a cell-extrinsic potent negative regulator of HSCs, was overexpressed in bone marrow niche cells from FLT3-ITD mice, and anti-TNF treatment partially rescued the HSC phenotype. These findings, which establish that Flt3-ITD-driven myeloproliferation results in cell-extrinsic suppression of the normal HSC reservoir, are of relevance for several aspects of acute myeloid leukemia biology.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

Understanding the genetic and nongenetic determinants of tumor protein 53 (TP53)-mutation-driven clonal evolution and subsequent transformation is a crucial step toward the design of rational therapeutic strategies. Here we carry out allelic resolution single-cell multi-omic analysis of hematopoietic stem/progenitor cells (HSPCs) from patients with a myeloproliferative neoplasm who transform to TP53-mutant secondary acute myeloid leukemia (sAML). All patients showed dominant TP53 'multihit' HSPC clones at transformation, with a leukemia stem cell transcriptional signature strongly predictive of adverse outcomes in independent cohorts, across both TP53-mutant and wild-type (WT) AML. Through analysis of serial samples, antecedent TP53-heterozygous clones and in vivo perturbations, we demonstrate a hitherto unrecognized effect of chronic inflammation, which suppressed TP53 WT HSPCs while enhancing the fitness advantage of TP53-mutant cells and promoted genetic evolution. Our findings will facilitate the development of risk-stratification, early detection and treatment strategies for TP53-mutant leukemia, and are of broad relevance to other cancer types.

How transcription factors (TFs) cooperate within large protein complexes to allow rapid modulation of gene expression during development is still largely unknown. Here we show that the key haematopoietic LIM-domain-binding protein-1 (LDB1) TF complex contains several activator and repressor components that together maintain an erythroid-specific gene expression programme primed for rapid activation until differentiation is induced. A combination of proteomics, functional genomics and in vivo studies presented here identifies known and novel co-repressors, most notably the ETO2 and IRF2BP2 proteins, involved in maintaining this primed state. The ETO2-IRF2BP2 axis, interacting with the NCOR1/SMRT co-repressor complex, suppresses the expression of the vast majority of archetypical erythroid genes and pathways until its decommissioning at the onset of terminal erythroid differentiation. Our experiments demonstrate that multimeric regulatory complexes feature a dynamic interplay between activating and repressing components that determines lineage-specific gene expression and cellular differentiation.

The key haematopoietic regulator Myb is essential for coordinating proliferation and differentiation. ChIP-Sequencing and Chromosome Conformation Capture (3C)-Sequencing were used to characterize the structural and protein-binding dynamics of the Myb locus during erythroid differentiation. In proliferating cells expressing Myb, enhancers within the Myb-Hbs1l intergenic region were shown to form an active chromatin hub (ACH) containing the Myb promoter and first intron. This first intron was found to harbour the transition site from transcription initiation to elongation, which takes place around a conserved CTCF site. Upon erythroid differentiation, Myb expression is downregulated and the ACH destabilized. We propose a model for Myb activation by distal enhancers dynamically bound by KLF1 and the GATA1/TAL1/LDB1 complex, which primarily function as a transcription elongation element through chromatin looping.

The Generation Challenge programme (GCP) is a global crop research consortium directed toward crop improvement through the application of comparative biology and genetic resources characterization to plant breeding. A key consortium research activity is the development of a GCP crop bioinformatics platform to support GCP research. This platform includes the following: (i) shared, public platform-independent domain models, ontology, and data formats to enable interoperability of data and analysis flows within the platform; (ii) web service and registry technologies to identify, share, and integrate information across diverse, globally dispersed data sources, as well as to access high-performance computational (HPC) facilities for computationally intensive, high-throughput analyses of project data; (iii) platform-specific middleware reference implementations of the domain model integrating a suite of public (largely open-access/-source) databases and software tools into a workbench to facilitate biodiversity analysis, comparative analysis of crop genomic data, and plant breeding decision making.