Targeted gene editing in hematopoietic stem cells (HSCs) is a promising treatment for several diseases. However, the limited efficiency of homology-directed repair (HDR) in HSCs and the unknown ...impact of the procedure on clonal composition and dynamics of transplantation have hampered clinical translation. Here, we apply a barcoding strategy to clonal tracking of edited cells (BAR-Seq) and show that editing activates p53, which substantially shrinks the HSC clonal repertoire in hematochimeric mice, although engrafted edited clones preserve multilineage and self-renewing capacity. Transient p53 inhibition restored polyclonal graft composition. We increased HDR efficiency by forcing cell-cycle progression and upregulating components of the HDR machinery through transient expression of the adenovirus 5 E4orf6/7 protein, which recruits the cell-cycle controller E2F on its target genes. Combined E4orf6/7 expression and p53 inhibition resulted in HDR editing efficiencies of up to 50% in the long-term human graft, without perturbing repopulation and self-renewal of edited HSCs. This enhanced protocol should broaden applicability of HSC gene editing and pave its way to clinical translation.
Single cell profiling has been proven to be a powerful tool in molecular biology to understand the complex behaviours of heterogeneous system. The definition of the properties of single cells is the ...primary endpoint of such analysis, cells are typically clustered to underpin the common determinants that can be used to describe functional properties of the cell mixture under investigation. Several approaches have been proposed to identify cell clusters; while this is matter of active research, one popular approach is based on community detection in neighbourhood graphs by optimisation of modularity. In this paper we propose an alternative and principled solution to this problem, based on Stochastic Block Models. We show that such approach not only is suitable for identification of cell groups, it also provides a solid framework to perform other relevant tasks in single cell analysis, such as label transfer. To encourage the use of Stochastic Block Models, we developed a python library, schist, that is compatible with the popular scanpy framework.
Transplantation of hematopoietic cells from a healthy individual (allogeneic hematopoietic cell transplantation (allo-HCT)) demonstrates that adoptive immunotherapy can cure blood cancers: still, ...post-transplantation relapses remain frequent. To explain their drivers, we analyzed the genomic and gene expression profiles of acute myeloid leukemia (AML) blasts purified from patients at serial time-points during their disease history. We identified a transcriptional signature specific for post-transplantation relapses and highly enriched in immune-related processes, including T cell costimulation and antigen presentation. In two independent patient cohorts we confirmed the deregulation of multiple costimulatory ligands on AML blasts at post-transplantation relapse (PD-L1, B7-H3, CD80, PVRL2), mirrored by concomitant changes in circulating donor T cells. Likewise, we documented the frequent loss of surface expression of HLA-DR, -DQ and -DP on leukemia cells, due to downregulation of the HLA class II regulator CIITA. We show that loss of HLA class II expression and upregulation of inhibitory checkpoint molecules represent alternative modalities to abolish AML recognition from donor-derived T cells, and can be counteracted by interferon-γ or checkpoint blockade, respectively. Our results demonstrate that the deregulation of pathways involved in T cell-mediated allorecognition is a distinctive feature and driver of AML relapses after allo-HCT, which can be rapidly translated into personalized therapies.
Follicular lymphoma is an incurable malignancy, with transformation to an aggressive subtype representing a critical event during disease progression. Here we performed whole-genome or whole-exome ...sequencing on 10 follicular lymphoma-transformed follicular lymphoma pairs followed by deep sequencing of 28 genes in an extension cohort, and we report the key events and evolutionary processes governing tumor initiation and transformation. Tumor evolution occurred through either a 'rich' or 'sparse' ancestral common progenitor clone (CPC). We identified recurrent mutations in linker histone, JAK-STAT signaling, NF-κB signaling and B cell developmental genes. Longitudinal analyses identified early driver mutations in chromatin regulator genes (CREBBP, EZH2 and KMT2D (MLL2)), whereas mutations in EBF1 and regulators of NF-κB signaling (MYD88 and TNFAIP3) were gained at transformation. Collectively, this study provides new insights into the genetic basis of follicular lymphoma and the clonal dynamics of transformation and suggests that personalizing therapies to target key genetic alterations in the CPC represents an attractive therapeutic strategy.
Gene editing by engineered nucleases has revolutionized the field of gene therapy by enabling targeted and precise modification of the genome. However, the limited availability of methods for clonal ...tracking of edited cells has resulted in a paucity of information on the diversity, abundance and behavior of engineered clones. Here we detail the wet laboratory and bioinformatic BAR-Seq pipeline, a strategy for clonal tracking of cells harboring homology-directed targeted integration of a barcoding cassette. We present the BAR-Seq web application, an online, freely available and easy-to-use software that allows performing clonal tracking analyses on raw sequencing data without any computational resources or advanced bioinformatic skills. BAR-Seq can be applied to most editing strategies, and we describe its use to investigate the clonal dynamics of human edited hematopoietic stem/progenitor cells in xenotransplanted hosts. Notably, BAR-Seq may be applied in both basic and translational research contexts to investigate the biology of edited cells and stringently compare editing protocols at a clonal level. Our BAR-Seq pipeline allows library preparation and validation in a few days and clonal analyses of edited cell populations in 1 week.
It is not clear how spontaneous DNA double-strand breaks (DSBs) form and are processed in normal cells, and whether they predispose to cancer-associated translocations. We show that DSBs in normal ...mammary cells form upon release of paused RNA polymerase II (Pol II) at promoters, 5' splice sites and active enhancers, and are processed by end-joining in the absence of a canonical DNA-damage response. Logistic and causal-association models showed that Pol II pausing at long genes is the main predictor and determinant of DSBs. Damaged introns with paused Pol II-pS5, TOP2B and XRCC4 are enriched in translocation breakpoints, and map at topologically associating domain boundary-flanking regions showing high interaction frequencies with distal loci. Thus, in unperturbed growth conditions, release of paused Pol II at specific loci and chromatin territories favors DSB formation, leading to chromosomal translocations.
Recent efforts have succeeded in surveying open chromatin at the single-cell level, but high-throughput, single-cell assessment of heterochromatin and its underlying genomic determinants remains ...challenging. We engineered a hybrid transposase including the chromodomain (CD) of the heterochromatin protein-1α (HP-1α), which is involved in heterochromatin assembly and maintenance through its binding to trimethylation of the lysine 9 on histone 3 (H3K9me3), and developed a single-cell method, single-cell genome and epigenome by transposases sequencing (scGET-seq), that, unlike single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), comprehensively probes both open and closed chromatin and concomitantly records the underlying genomic sequences. We tested scGET-seq in cancer-derived organoids and human-derived xenograft (PDX) models and identified genetic events and plasticity-driven mechanisms contributing to cancer drug resistance. Next, building upon the differential enrichment of closed and open chromatin, we devised a method, Chromatin Velocity, that identifies the trajectories of epigenetic modifications at the single-cell level. Chromatin Velocity uncovered paths of epigenetic reorganization during stem cell reprogramming and identified key transcription factors driving these developmental processes. scGET-seq reveals the dynamics of genomic and epigenetic landscapes underlying any cellular processes.
Background: Analysis of scATAC-seq data has been recently scaled to thousands of cells. While processing of other types of single cell data was boosted by the implementation of alignment-free ...techniques, pipelines available to process scATAC-seq data still require large computational resources. We propose here an approach based on pseudoalignment, which reduces the execution times and hardware needs at little cost for precision.
Methods: Public data for 10k PBMC were downloaded from 10x Genomics web site. Reads were aligned to various references derived from DNase I Hypersensitive Sites (DHS) using
kallisto and quantified with
bustools. We compared our results with the ones publicly available derived by
cellranger-atac. We subsequently tested our approach on scATAC-seq data for K562 cell line.
Results: We found that
kallisto does not introduce biases in quantification of known peaks; cells groups identified are consistent with the ones identified from standard method. We also found that cell identification is robust when analysis is performed using DHS-derived reference in place of
de novo identification of ATAC peaks. Lastly, we found that our approach is suitable for reliable quantification of gene activity based on scATAC-seq signal, thus allows for efficient labelling of cell groups based on marker genes.
Conclusions: Analysis of scATAC-seq data by means of
kallisto produces results in line with standard pipelines while being considerably faster; using a set of known DHS sites as reference does not affect the ability to characterize the cell populations.
The ribonuclease DIS3 is one of the most frequently mutated genes in the hematological cancer multiple myeloma, yet the basis of its tumor suppressor function in this disease remains unclear. Herein, ...exploiting the TCGA dataset, we found that DIS3 plays a prominent role in the DNA damage response. DIS3 inactivation causes genomic instability by increasing mutational load, and a pervasive accumulation of DNA:RNA hybrids that induces genomic DNA double‐strand breaks (DSBs). DNA:RNA hybrid accumulation also prevents binding of the homologous recombination (HR) machinery to double‐strand breaks, hampering DSB repair. DIS3‐inactivated cells become sensitive to PARP inhibitors, suggestive of a defect in homologous recombination repair. Accordingly, multiple myeloma patient cells mutated for DIS3 harbor an increased mutational burden and a pervasive overexpression of pro‐inflammatory interferon, correlating with the accumulation of DNA:RNA hybrids. We propose DIS3 loss in myeloma to be a driving force for tumorigenesis via DNA:RNA hybrid‐dependent enhanced genome instability and increased mutational rate. At the same time, DIS3 loss represents a liability that might be therapeutically exploited in patients whose cancer cells harbor DIS3 mutations.
Synopsis
The ribonuclease DIS3 is frequently mutated in the blood cancer multiple myeloma. Here, DIS3 inactivation is found to cause accumulation of DNA:RNA hybrids, as well as to increases interferon responses and reduce homologous recombination.
DIS3 loss triggers a genome‐wide increase in DNA:RNA hybrids, which in turn leads to DNA fragmentation and genomic instability.
Hybrids accumulation at the sites of DNA damage prevents BRCA1 binding to DNA, impairing homologous recombination‐based DNA repair.
DIS3 loss is associated with increased mutational rate both in vitro and in patient samples with DIS3 mutations.
Myeloma cells derived from patients presenting DIS3 mutations display an intense interferon response.
DIS3 mutation in hematological cancer causes reduced homologous recombination repair, increased mutational burden, and overactivation of inflammatory interferon responses.
We report the genome-wide mapping of ORC1 binding sites in mammals, by chromatin immunoprecipitation and parallel sequencing (ChIP-seq). ORC1 binding sites in HeLa cells were validated as active DNA ...replication origins (ORIs) using Repli-seq, a method that allows identification of ORI-containing regions by parallel sequencing of temporally ordered replicating DNA. ORC1 sites were universally associated with transcription start sites (TSSs) of coding or noncoding RNAs (ncRNAs). Transcription levels at the ORC1 sites directly correlated with replication timing, suggesting the existence of two classes of ORIs: those associated with moderate/high transcription levels (≥1 RNA copy/cell), firing in early S and mapping to the TSSs of coding RNAs; and those associated with low transcription levels (<1 RNA copy/cell), firing throughout the entire S and mapping to TSSs of ncRNAs. These findings are compatible with a scenario whereby TSS expression levels influence the efficiency of ORC1 recruitment at G(1) and the probability of firing during S.