Single-cell analysis of immune recognition in chronic myeloid leukemia patients following tyrosine kinase inhibitor discontinuation.

Immunological control of residual leukemia cells is thought to occur in patients with chronic myeloid leukemia (CML) that maintain treatment-free remission (TFR) following tyrosine kinase inhibitor (TKI) discontinuation. To study this, we analyzed 55 single-cell RNA and T cell receptor (TCR) sequenced samples (scRNA+TCRαß-seq) from patients with CML (n = 13, N = 25), other cancers (n = 28), and healthy (n = 7). The high number and active phenotype of natural killer (NK) cells in CML separated them from healthy and other cancers. Most NK cells in CML belonged to the active CD56dim cluster with high expression of GZMA/B, PRF1, CCL3/4, and IFNG, with interactions with leukemic cells via inhibitory LGALS9-TIM3 and PVR-TIGIT interactions. Accordingly, upregulation of LGALS9 was observed in CML target cells and TIM3 in NK cells when co-cultured together. Additionally, we created a classifier to identify TCRs targeting leukemia-associated antigen PR1 and quantified anti-PR1 T cells in 90 CML and 786 healthy TCRß-sequenced samples. Anti-PR1 T cells were more prevalent in CML, enriched in bone marrow samples, and enriched in the mature, cytotoxic CD8 + TEMRA cluster, especially in a patient maintaining TFR. Our results highlight the role of NK cells and anti-PR1 T cells in anti-leukemic immune responses in CML.

Tracking early mammalian organogenesis - prediction and validation of differentiation trajectories at whole organism scale.

Early organogenesis represents a key step in animal development, during which pluripotent cells diversify to initiate organ formation. Here, we sampled 300,000 single-cell transcriptomes from mouse embryos between E8.5 and E9.5 in 6-h intervals and combined this new dataset with our previous atlas (E6.5-E8.5) to produce a densely sampled timecourse of >400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction identified complex waves of blood and endothelial development, including a new programme for somite-derived endothelium. We also dissected the E7.5 primitive streak into four adjacent regions, performed scRNA-seq and predicted cell fates computationally. Finally, we defined developmental state/fate relationships by combining orthotopic grafting, microscopic analysis and scRNA-seq to transcriptionally determine cell fates of grafted primitive streak regions after 24 h of in vitro embryo culture. Experimentally determined fate outcomes were in good agreement with computationally predicted fates, demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in developmental biology and guide the in vitro production of cells for organ regeneration and repair.

Single-cell multiomics of human fetal hematopoiesis define a developmental-specific population and a fetal signature.

Knowledge of human fetal blood development and how it differs from adult blood is highly relevant to our understanding of congenital blood and immune disorders and childhood leukemia, of which the latter can originate in utero. Blood formation occurs in waves that overlap in time and space, adding to heterogeneity, which necessitates single-cell approaches. Here, a combined single-cell immunophenotypic and transcriptional map of first trimester primitive blood development is presented. Using CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing), the molecular profile of established immunophenotype-gated progenitors was analyzed in the fetal liver (FL). Classical markers for hematopoietic stem cells (HSCs), such as CD90 and CD49F, were largely preserved, whereas CD135 (FLT3) and CD123 (IL3R) had a ubiquitous expression pattern capturing heterogenous populations. Direct molecular comparison with an adult bone marrow data set revealed that the HSC state was less frequent in FL, whereas cells with a lymphomyeloid signature were more abundant. An erythromyeloid-primed multipotent progenitor cluster was identified, potentially representing a transient, fetal-specific population. Furthermore, differentially expressed genes between fetal and adult counterparts were specifically analyzed, and a fetal core signature was identified. The core gene set could separate subgroups of acute lymphoblastic leukemia by age, suggesting that a fetal program may be partially retained in specific subgroups of pediatric leukemia. Our detailed single-cell map presented herein emphasizes molecular and immunophenotypic differences between fetal and adult blood cells, which are of significance for future studies of pediatric leukemia and blood development in general.

Temporal multimodal single-cell profiling of native hematopoiesis illuminates altered differentiation trajectories with age.

Aging negatively affects hematopoiesis, with consequences for immunity and acquired blood cell disorders. Although impairments in hematopoietic stem cell (HSC) function contribute to this, the in vivo dynamics of such changes remain obscure. Here, we integrate extensive longitudinal functional assessments of HSC-specific lineage tracing with single-cell transcriptome and epitope profiling. In contrast to recent suggestions from single-cell RNA sequencing alone, our data favor a defined structure of HSC/progenitor differentiation that deviates substantially from HSC-derived hematopoiesis following transplantation. Native age-dependent attrition in HSC differentiation manifests as drastically reduced lymphoid output through an early lymphoid-primed progenitor (MPP Ly-I). While in vitro activation fails to rescue lymphoid differentiation from most aged HSCs, robust lymphopoiesis can be achieved by culturing elevated numbers of candidate HSCs. Therefore, our data position rare chronologically aged HSC clones, fully competent at producing lymphoid offspring, as a prime target for approaches aimed to improve lymphopoiesis in the elderly.

In vitro clonal multilineage differentiation of distinct murine hematopoietic progenitor populations.

Here we describe an in vitro co-culture system that can differentiate hematopoietic progenitor populations to all major hematopoietic lineages at clonal level. We present both a sensitive single-cell switch-culture system as well as a less laborious alternative barcoding protocol more convenient for larger cell numbers. Importantly, generation of all lineages from single long-term hematopoietic stem cells are described, following 21 days of culture. This protocol represents an efficient tool for validation experiments for single-cell genomics data. For complete details on the use and execution of this protocol, please refer to Safi et al. (2022).1.

Scarf enables a highly memory-efficient analysis of large-scale single-cell genomics data.

As the scale of single-cell genomics experiments grows into the millions, the computational requirements to process this data are beyond the reach of many. Herein we present Scarf, a modularly designed Python package that seamlessly interoperates with other single-cell toolkits and allows for memory-efficient single-cell analysis of millions of cells on a laptop or low-cost devices like single-board computers. We demonstrate Scarf's memory and compute-time efficiency by applying it to the largest existing single-cell RNA-Seq and ATAC-Seq datasets. Scarf wraps memory-efficient implementations of a graph-based t-stochastic neighbour embedding and hierarchical clustering algorithm. Moreover, Scarf performs accurate reference-anchored mapping of datasets while maintaining memory efficiency. By implementing a subsampling algorithm, Scarf additionally has the capacity to generate representative sampling of cells from a given dataset wherein rare cell populations and lineage differentiation trajectories are conserved. Together, Scarf provides a framework wherein any researcher can perform advanced processing, subsampling, reanalysis, and integration of atlas-scale datasets on standard laptop computers. Scarf is available on Github: https://github.com/parashardhapola/scarf .

Concurrent stem- and lineage-affiliated chromatin programs precede hematopoietic lineage restriction.

The emerging notion of hematopoietic stem and progenitor cells (HSPCs) as a low-primed cloud without sharply demarcated gene expression programs raises the question on how cellular-fate options emerge and at which stem-like stage lineage priming is initiated. Here, we investigate single-cell chromatin accessibility of Lineage-, cKit+, and Sca1+ (LSK) HSPCs spanning the early differentiation landscape. Application of a signal-processing algorithm to detect transition points corresponding to massive alterations in accessibility of 571 transcription factor motifs reveals a population of LSK FMS-like tyrosine kinase 3 (Flt3)intCD9high cells that concurrently display stem-like and lineage-affiliated chromatin signatures, pointing to a simultaneous gain of both lympho-myeloid and megakaryocyte-erythroid programs. Molecularly and functionally, these cells position between stem cells and committed progenitors and display multi-lineage capacity in vitro and in vivo but lack self-renewal activity. This integrative molecular analysis resolves the heterogeneity of cells along hematopoietic differentiation and permits investigation of chromatin-mediated transition between multipotency and lineage restriction.

Genome-wide association study on 13 167 individuals identifies regulators of blood CD34+cell levels.

Stem cell transplantation is a cornerstone in the treatment of blood malignancies. The most common method to harvest stem cells for transplantation is by leukapheresis, requiring mobilization of CD34+ hematopoietic stem and progenitor cells (HSPCs) from the bone marrow into the blood. Identifying the genetic factors that control blood CD34+ cell levels could reveal new drug targets for HSPC mobilization. Here we report the first large-scale, genome-wide association study on blood CD34+ cell levels. Across 13 167 individuals, we identify 9 significant and 2 suggestive associations, accounted for by 8 loci (PPM1H, CXCR4, ENO1-RERE, ITGA9, ARHGAP45, CEBPA, TERT, and MYC). Notably, 4 of the identified associations map to CXCR4, showing that bona fide regulators of blood CD34+ cell levels can be identified through genetic variation. Further, the most significant association maps to PPM1H, encoding a serine/threonine phosphatase never previously implicated in HSPC biology. PPM1H is expressed in HSPCs, and the allele that confers higher blood CD34+ cell levels downregulates PPM1H. Through functional fine-mapping, we find that this downregulation is caused by the variant rs772557-A, which abrogates an MYB transcription factor-binding site in PPM1H intron 1 that is active in specific HSPC subpopulations, including hematopoietic stem cells, and interacts with the promoter by chromatin looping. Furthermore, PPM1H knockdown increases the proportion of CD34+ and CD34+90+ cells in cord blood assays. Our results provide the first large-scale analysis of the genetic architecture of blood CD34+ cell levels and warrant further investigation of PPM1H as a potential inhibition target for stem cell mobilization.

CellexalVR: A virtual reality platform to visualize and analyze single-cell omics data.

Single-cell RNAseq is a routinely used method to explore heterogeneity within cell populations. Data from these experiments are often visualized using dimension reduction methods such as UMAP and tSNE, where each cell is projected in two or three dimensional space. Three-dimensional projections can be more informative for larger and complex datasets because they are less prone to merging and flattening similar cell-types/clusters together. However, visualizing and cross-comparing 3D projections using current software on conventional flat-screen displays is far from optimal as they are still essentially 2D, and lack meaningful interaction between the user and the data. Here we present CellexalVR (www.cellexalvr.med.lu.se), a feature-rich, fully interactive virtual reality environment for the visualization and analysis of single-cell experiments that allows researchers to intuitively and collaboratively gain an understanding of their data.

Loss of Cxcr5 alters neuroblast proliferation and migration in the aged brain.

Neurogenesis, the production of new neurons from neural stem cells, dramatically decreases during aging concomitantly with increased inflammation both systemically and in the brain. However, the precise role of inflammation and whether local or systemic factors drive the neurogenic decline during aging is poorly understood. Here, we identify CXCR5/5/CXCL13 signaling as a novel regulator of neurogenesis in the aged brain. The chemokine Cxcl13 was found to be upregulated in the brain during aging. Loss of its receptor, Cxcr5, led to increased proliferation and decreased numbers of neuroblasts in the aged subventricular zone (SVZ), together with accumulation of neuroblasts in the rostral migratory stream and olfactory bulb (OB), without increasing the amount of new mature neurons in the OB. The effect on proliferation and migration was specific to neuroblasts and likely mediated through increased levels of systemic IL-6 and local Cxcl12 expression in the SVZ. Our study raises the possibility of a new mechanism by which interplay between systemic and local alterations in inflammation regulates neurogenesis during aging.

BLM Potentiates c-Jun Degradation and Alters Its Function as an Oncogenic Transcription Factor.

Mutations in BLM helicase predispose Bloom syndrome (BS) patients to a wide spectrum of cancers. We demonstrate that MIB1-ubiquitylated BLM in G1 phase functions as an adaptor protein by enhancing the binding of transcription factor c-Jun and its E3 ligase, Fbw7α. BLM enhances the K48/K63-linked ubiquitylation on c-Jun, thereby enhancing the rate of its subsequent degradation. Functionally defective Fbw7α mutants prevalent in multiple human cancers are reactivated by BLM. However, BS patient-derived BLM mutants cannot potentiate Fbw7α-dependent c-Jun degradation. The decrease in the levels of c-Jun in cells expressing BLM prevents effective c-Jun binding to 2,584 gene promoters. This causes decreases in the transcript and protein levels of c-Jun targets in BLM-expressing cells, resulting in attenuated c-Jun-dependent effects during neoplastic transformation. Thus, BLM carries out its function as a tumor suppressor by enhancing c-Jun turnover and thereby preventing its activity as a proto-oncogene.

Novel internal regulators and candidate miRNAs within miR-379/miR-656 miRNA cluster can alter cellular phenotype of human glioblastoma.

Clustered miRNAs can affect functioning of downstream pathways due to possible coordinated function. We observed 78-88% of the miR-379/miR-656 cluster (C14MC) miRNAs were downregulated in three sub-types of diffuse gliomas, which was also corroborated with analysis from The Cancer Genome Atlas (TCGA) datasets. The miRNA expression levels decreased with increasing tumor grade, indicating this downregulation as an early event in gliomagenesis. Higher expression of the C14MC miRNAs significantly improved glioblastioma prognosis (Pearson's r = 0.62; p < 3.08e-22). ENCODE meta-data analysis, followed by reporter assays validated existence of two novel internal regulators within C14MC. CRISPR activation of the most efficient internal regulator specifically induced members of the downstream miRNA sub-cluster and apoptosis in glioblastoma cells. Luciferase assays validated novel targets for miR-134 and miR-485-5p, two miRNAs from C14MC with the most number of target genes relevant for glioma. Overexpression of miR-134 and miR-485-5p in human glioblastoma cells suppressed invasion and proliferation, respectively. Furthermore, apoptosis was induced by both miRs, individually and in combination. The results emphasize the tumor suppressive role of C14MC in diffuse gliomas, and identifies two specific miRNAs with potential therapeutic value and towards better disease management and therapy.

Transcription regulation of CDKN1A (p21/CIP1/WAF1) by TRF2 is epigenetically controlled through the REST repressor complex.

We observed extra-telomeric binding of the telomere repeat binding factor TRF2 within the promoter of the cyclin-dependent kinase CDKNIA (p21/CIP1/WAF1). This result in TRF2 induced transcription repression of p21. Interestingly, p21 repression was through engagement of the REST-coREST-LSD1-repressor complex and altered histone marks at the p21 promoter in a TRF2-dependent fashion. Furthermore, mutational analysis shows p21 repression requires interaction of TRF2 with a p21 promoter G-quadruplex. Physiologically, TRF2-mediated p21 repression attenuated drug-induced activation of cellular DNA damage response by evading G2/M arrest in cancer cells. Together these reveal for the first time role of TRF2 in REST- repressor complex mediated transcription repression.

ATF3 negatively regulates cellular antiviral signaling and autophagy in the absence of type I interferons.

Stringent regulation of antiviral signaling and cellular autophagy is critical for the host response to virus infection. However, little is known how these cellular processes are regulated in the absence of type I interferon signaling. Here, we show that ATF3 is induced following Japanese encephalitis virus (JEV) infection, and regulates cellular antiviral and autophagy pathways in the absence of type I interferons in mouse neuronal cells. We have identified new targets of ATF3 and show that it binds to the promoter regions of Stat1, Irf9, Isg15 and Atg5 thereby inhibiting cellular antiviral signaling and autophagy. Consistent with these observations, ATF3-depleted cells showed enhanced antiviral responses and induction of robust autophagy. Furthermore, we show that JEV replication was significantly reduced in ATF3-depleted cells. Our findings identify ATF3 as a negative regulator of antiviral signaling and cellular autophagy in mammalian cells, and demonstrate its important role in JEV life cycle.

Epigenetic suppression of human telomerase (hTERT) is mediated by the metastasis suppressor NME2 in a G-quadruplex-dependent fashion.

Transcriptional activation of the human telomerase reverse transcriptase (hTERT) gene, which remains repressed in adult somatic cells, is critical during tumorigenesis. Several transcription factors and the epigenetic state of the hTERT promoter are known to be important for tight control of hTERT in normal tissues, but the molecular mechanisms leading to hTERT reactivation in cancer are not well-understood. Surprisingly, here we found occupancy of the metastasis suppressor non-metastatic 2 (NME2) within the hTERT core promoter in HT1080 fibrosarcoma cells and HCT116 colon cancer cells and NME2-mediated transcriptional repression of hTERT in these cells. We also report that loss of NME2 results in up-regulated hTERT expression. Mechanistically, additional results indicated that the RE1-silencing transcription factor (REST)-lysine-specific histone demethylase 1 (LSD1) co-repressor complex associates with the hTERT promoter in an NME2-dependent way and that this assembly is required for maintaining repressive chromatin at the hTERT promoter. Interestingly, a G-quadruplex motif at the hTERT promoter was essential for occupancy of NME2 and the REST repressor complex on the hTERT promoter. In light of this mechanistic insight, we studied the effects of G-quadruplex-binding ligands on hTERT expression and observed that several of these ligands repressed hTERT expression. Together, our results support a mechanism of hTERT epigenetic control involving a G-quadruplex promoter motif, which potentially can be targeted by tailored small molecules.

QuadBase2: web server for multiplexed guanine quadruplex mining and visualization.

DNA guanine quadruplexes or G4s are non-canonical DNA secondary structures which affect genomic processes like replication, transcription and recombination. G4s are computationally identified by specific nucleotide motifs which are also called putative G4 (PG4) motifs. Despite the general relevance of these structures, there is currently no tool available that can allow batch queries and genome-wide analysis of these motifs in a user-friendly interface. QuadBase2 (quadbase.igib.res.in) presents a completely reinvented web server version of previously published QuadBase database. QuadBase2 enables users to mine PG4 motifs in up to 178 eukaryotes through the EuQuad module. This module interfaces with Ensembl Compara database, to allow users mine PG4 motifs in the orthologues of genes of interest across eukaryotes. PG4 motifs can be mined across genes and their promoter sequences in 1719 prokaryotes through ProQuad module. This module includes a feature that allows genome-wide mining of PG4 motifs and their visualization as circular histograms. TetraplexFinder, the module for mining PG4 motifs in user-provided sequences is now capable of handling up to 20 MB of data. QuadBase2 is a comprehensive PG4 motif mining tool that further expands the configurations and algorithms for mining PG4 motifs in a user-friendly way.

BreCAN-DB: a repository cum browser of personalized DNA breakpoint profiles of cancer genomes.

BreCAN-DB (http://brecandb.igib.res.in) is a repository cum browser of whole genome somatic DNA breakpoint profiles of cancer genomes, mapped at single nucleotide resolution using deep sequencing data. These breakpoints are associated with deletions, insertions, inversions, tandem duplications, translocations and a combination of these structural genomic alterations. The current release of BreCAN-DB features breakpoint profiles from 99 cancer-normal pairs, comprising five cancer types. We identified DNA breakpoints across genomes using high-coverage next-generation sequencing data obtained from TCGA and dbGaP. Further, in these cancer genomes, we methodically identified breakpoint hotspots which were significantly enriched with somatic structural alterations. To visualize the breakpoint profiles, a next-generation genome browser was integrated with BreCAN-DB. Moreover, we also included previously reported breakpoint profiles from 138 cancer-normal pairs, spanning 10 cancer types into the browser. Additionally, BreCAN-DB allows one to identify breakpoint hotspots in user uploaded data set. We have also included a functionality to query overlap of any breakpoint profile with regions of user's interest. Users can download breakpoint profiles from the database or may submit their data to be integrated in BreCAN-DB. We believe that BreCAN-DB will be useful resource for genomics scientific community and is a step towards personalized cancer genomics.