In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis.

Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis.

Gene therapy restores the transcriptional program of hematopoietic stem cells in Fanconi anemia.

Clinical trials have shown that lentiviral-mediated gene therapy can ameliorate bone marrow failure (BMF) in nonconditioned Fanconi anemia (FA) patients resulting from the proliferative advantage of corrected FA hematopoietic stem and progenitor cells (HSPC). However, it is not yet known if gene therapy can revert affected molecular pathways in diseased HSPC. Single-cell RNA sequencing was performed in chimeric populations of corrected and uncorrected HSPC co-existing in the BM of gene therapy-treated FA patients. Our study demonstrates that gene therapy reverts the transcriptional signature of FA HSPC, which then resemble the transcriptional program of healthy donor HSPC. This includes a down-regulated expression of TGF-ß and p21, typically up-regulated in FA HSPC, and upregulation of DNA damage response and telomere maintenance pathways. Our results show for the first time the potential of gene therapy to rescue defects in the HSPC transcriptional program from patients with inherited diseases; in this case, in FA characterized by BMF and cancer predisposition.

Uncovering perturbations in human hematopoiesis associated with healthy aging and myeloid malignancies at single-cell resolution.

Early hematopoiesis is a continuous process in which hematopoietic stem and progenitor cells (HSPCs) gradually differentiate toward specific lineages. Aging and myeloid malignant transformation are characterized by changes in the composition and regulation of HSPCs. In this study, we used single-cell RNA sequencing (scRNA-seq) to characterize an enriched population of human HSPCs obtained from young and elderly healthy individuals.Based on their transcriptional profile, we identified changes in the proportions of progenitor compartments during aging, and differences in their functionality, as evidenced by gene set enrichment analysis. Trajectory inference revealed that altered gene expression dynamics accompanied cell differentiation, which could explain aging-associated changes in hematopoiesis. Next, we focused on key regulators of transcription by constructing gene regulatory networks (GRNs) and detected regulons that were specifically active in elderly individuals. Using previous findings in healthy cells as a reference, we analyzed scRNA-seq data obtained from patients with myelodysplastic syndrome (MDS) and detected specific alterations of the expression dynamics of genes involved in erythroid differentiation in all patients with MDS such as TRIB2. In addition, the comparison between transcriptional programs and GRNs regulating normal HSPCs and MDS HSPCs allowed identification of regulons that were specifically active in MDS cases such as SMAD1, HOXA6, POU2F2, and RUNX1 suggesting a role of these transcription factors (TFs) in the pathogenesis of the disease.In summary, we demonstrate that the combination of single-cell technologies with computational analysis tools enable the study of a variety of cellular mechanisms involved in complex biological systems such as early hematopoiesis and can be used to dissect perturbed differentiation trajectories associated with perturbations such as aging and malignant transformation. Furthermore, the identification of abnormal regulatory mechanisms associated with myeloid malignancies could be exploited for personalized therapeutic approaches in individual patients.

Our blood contains many different types of cells; red blood cells carry oxygen through the body, platelets help to stop bleeding and a variety of white blood cells fight infections. All of these critical components come from a pool of immature cells in bone marrow, which can develop and specialise into any of these. However, as we get older, these immature cells can accumulate damage, including mutations in specific genes. This increases the risk of diseases such as myelodysplastic syndromes (MDS), a type of cancer in which the cells cannot develop and the patient does not have enough healthy mature blood cells. The changes in gene activity in the immature cells have previously been studied using samples from young and elderly people, as well as individuals with MDS. These studies examined large numbers of cells together, revealing differences between young and elderly people, and individuals with MDS. However, this does not describe how the different types alter their behaviour. To address this, Ainciburu, Ezponda et al. used a technique called single-cell RNA sequencing to study the gene activity in individual immature blood cells. This revealed changes associated with maturation that may account for the different combinations of cell populations in younger and older people. The results confirmed findings from previous studies and suggested new genes involved in ageing or MDS. Ainciburu, Ezponda et al. used these results to create an analytical system that highlights gene activity differences in individual MDS patients that are independent of age-related changes. These results provide new insights that could help further research into the development of MDS and the ageing process. In addition, scientists could study other diseases using this approach of analysing individual patients' gene activity. In future, this could help to personalise clinical decisions on diagnosis and treatment.

Integrated flow cytometry and sequencing to reconstruct evolutionary patterns from dysplasia to acute myeloid leukemia.

Clonal evolution in acute myeloid leukemia (AML) originates long before diagnosis and is a dynamic process that may affect survival. However, it remains uninvestigated during routine diagnostic workups. We hypothesized that the mutational status of bone marrow dysplastic cells and leukemic blasts, analyzed at the onset of AML using integrated multidimensional flow cytometry (MFC) immunophenotyping and fluorescence-activated cell sorting (FACS) with next-generation sequencing (NGS), could reconstruct leukemogenesis. Dysplastic cells were detected by MFC in 285 of 348 (82%) newly diagnosed patients with AML. Presence of dysplasia according to MFC and World Health Organization criteria had no prognostic value in older adults. NGS of dysplastic cells and blasts isolated at diagnosis identified 3 evolutionary patterns: stable (n = 12 of 21), branching (n = 4 of 21), and clonal evolution (n = 5 of 21). In patients achieving complete response (CR), integrated MFC and FACS with NGS showed persistent measurable residual disease (MRD) in phenotypically normal cell types, as well as the acquisition of genetic traits associated with treatment resistance. Furthermore, whole-exome sequencing of dysplastic and leukemic cells at diagnosis and of MRD uncovered different clonal involvement in dysplastic myelo-erythropoiesis, leukemic transformation, and chemoresistance. Altogether, we showed that it is possible to reconstruct leukemogenesis in ∼80% of patients with newly diagnosed AML, using techniques other than single-cell multiomics.

The transcription factor DDIT3 is a potential driver of dyserythropoiesis in myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in older individuals. Here we analyze the transcriptome of human HSCs purified from young and older healthy adults, as well as MDS patients, identifying transcriptional alterations following different patterns of expression. While aging-associated lesions seem to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detect the upregulation of the transcription factor DNA Damage Inducible Transcript 3 (DDIT3). Overexpression of DDIT3 in human healthy HSCs induces an MDS-like transcriptional state, and dyserythropoiesis, an effect associated with a failure in the activation of transcriptional programs required for normal erythroid differentiation. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia is able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythroid differentiation characterizing MDS patients.

Gemcitabine-mediated depletion of immunosuppressive dendritic cells enhances the efficacy of therapeutic vaccination.

Vaccination using optimized strategies may increase response rates to immune checkpoint inhibitors (ICI) in some tumors. To enhance vaccine potency and improve thus responses to ICI, we analyzed the gene expression profile of an immunosuppressive dendritic cell (DC) population induced during vaccination, with the goal of identifying druggable inhibitory mechanisms. RNAseq studies revealed targetable genes, but their inhibition did not result in improved vaccines. However, we proved that immunosuppressive DC had a monocytic origin. Thus, monocyte depletion by gemcitabine administration reduced the generation of these DC and increased vaccine-induced immunity, which rejected about 20% of LLC-OVA and B16-OVA tumors, which are non-responders to anti-PD-1. This improved efficacy was associated with higher tumor T-cell infiltration and overexpression of PD-1/PD-L1. Therefore, the combination of vaccine + gemcitabine with anti-PD-1 was superior to anti-PD-1 monotherapy in both models. B16-OVA tumors benefited from a synergistic effect, reaching 75% of tumor rejection, but higher levels of exhausted T-cells in LLC-OVA tumors co-expressing PD-1, LAG3 and TIM3 precluded similar levels of efficacy. Our results indicate that gemcitabine is a suitable combination therapy with vaccines aimed at enhancing PD-1 therapies by targeting vaccine-induced immunosuppressive DC.

CAR density influences antitumoral efficacy of BCMA CAR T cells and correlates with clinical outcome.

Identification of new markers associated with long-term efficacy in patients treated with CAR T cells is a current medical need, particularly in diseases such as multiple myeloma. In this study, we address the impact of CAR density on the functionality of BCMA CAR T cells. Functional and transcriptional studies demonstrate that CAR T cells with high expression of the CAR construct show an increased tonic signaling with up-regulation of exhaustion markers and increased in vitro cytotoxicity but a decrease in in vivo BM infiltration. Characterization of gene regulatory networks using scRNA-seq identified regulons associated to activation and exhaustion up-regulated in CARHigh T cells, providing mechanistic insights behind differential functionality of these cells. Last, we demonstrate that patients treated with CAR T cell products enriched in CARHigh T cells show a significantly worse clinical response in several hematological malignancies. In summary, our work demonstrates that CAR density plays an important role in CAR T activity with notable impact on clinical response.

Deconvolution of the hematopoietic stem cell microenvironment reveals a high degree of specialization and conservation.

Understanding the regulation of normal and malignant human hematopoiesis requires comprehensive cell atlas of the hematopoietic stem cell (HSC) regulatory microenvironment. Here, we develop a tailored bioinformatic pipeline to integrate public and proprietary single-cell RNA sequencing (scRNA-seq) datasets. As a result, we robustly identify for the first time 14 intermediate cell states and 11 stages of differentiation in the endothelial and mesenchymal BM compartments, respectively. Our data provide the most comprehensive description to date of the murine HSC-regulatory microenvironment and suggest a higher level of specialization of the cellular circuits than previously anticipated. Furthermore, this deep characterization allows inferring conserved features in human, suggesting that the layers of microenvironmental regulation of hematopoiesis may also be shared between species. Our resource and methodology is a stepping-stone toward a comprehensive cell atlas of the BM microenvironment.

Preneoplastic somatic mutations including MYD88L265P in lymphoplasmacytic lymphoma.

Normal cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B lymphocytes before lymphoma remains uncertain. We sequenced multiple stages of the B lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of the high prevalence of mutated MYD88. We observed similar accumulation of random mutations in B lineages from both cohorts and unexpectedly found MYD88L265P in normal precursor and mature B lymphocytes from patients with lymphoma. We uncovered genetic and transcriptional pathways driving malignant transformation and leveraged these to model lymphoplasmacytic lymphoma in mice, based on mutated MYD88 in B cell precursors and BCL2 overexpression. Thus, MYD88L265P is a preneoplastic event, which challenges the current understanding of lymphomagenesis and may have implications for early detection of B cell lymphomas.

TCR-induced FOXP3 expression by CD8+ T cells impairs their anti-tumor activity.

Adoptive cell transfer therapy using CD8+ T lymphocytes showed promising results eradicating metastatic malignancies. However, several regulatory mechanisms limit its efficacy. We studied the role of the expression of the transcription factor FOXP3 on CD8+ T cell function and anti-tumor immunity. Here we show that suboptimal T cell receptor stimulation of CD8+ T cells upregulates FOXP3 in vitro. Similarly, CD8 T cells transferred into tumor-bearing mice upregulate FOXP3 in vivo. Cell-intrinsic loss of FOXP3 by CD8+ T cells resulted in improved functionality after TCR stimulation and better antitumor responses in vivo. Inhibition of the FOXP3/NFAT interaction likewise improved CD8+ T cell functionality. Transcriptomic analysis of cells after TCR stimulation revealed an enrichment of genes implicated in the response to IFN-γ, IFN-α, inflammatory response, IL-6/JAK/STAT, G2M checkpoint and IL-2/STAT signaling in FOXP3-deficient CD8+ T cells with respect to FOXP3-wt CD8+ T cells. Our results suggest that transient expression of FOXP3 by CD8+ T cells in the tumor microenvironment restrains their anti-tumor activity, with clear implications for improving T cell responses during immunotherapy.

FlowCT for the analysis of large immunophenotypic data sets and biomarker discovery in cancer immunology.

Large-scale immune monitoring is becoming routinely used in clinical trials to identify determinants of treatment responsiveness, particularly to immunotherapies. Flow cytometry remains one of the most versatile and high throughput approaches for single-cell analysis; however, manual interpretation of multidimensional data poses a challenge when attempting to capture full cellular diversity and provide reproducible results. We present FlowCT, a semi-automated workspace empowered to analyze large data sets. It includes pre-processing, normalization, multiple dimensionality reduction techniques, automated clustering, and predictive modeling tools. As a proof of concept, we used FlowCT to compare the T-cell compartment in bone marrow (BM) with peripheral blood (PB) from patients with smoldering multiple myeloma (SMM), identify minimally invasive immune biomarkers of progression from smoldering to active MM, define prognostic T-cell subsets in the BM of patients with active MM after treatment intensification, and assess the longitudinal effect of maintenance therapy in BM T cells. A total of 354 samples were analyzed and immune signatures predictive of malignant transformation were identified in 150 patients with SMM (hazard ratio [HR], 1.7; P < .001). We also determined progression-free survival (HR, 4.09; P < .0001) and overall survival (HR, 3.12; P = .047) in 100 patients with active MM. New data also emerged about stem cell memory T cells, the concordance between immune profiles in BM and PB, and the immunomodulatory effect of maintenance therapy. FlowCT is a new open-source computational approach that can be readily implemented by research laboratories to perform quality control, analyze high-dimensional data, unveil cellular diversity, and objectively identify biomarkers in large immune monitoring studies. These trials were registered at www.clinicaltrials.gov as #NCT01916252 and #NCT02406144.

Complement C5a induces the formation of neutrophil extracellular traps by myeloid-derived suppressor cells to promote metastasis.

Myeloid-derived suppressor cells (MDSCs) play a major role in cancer progression. In this study, we investigated the mechanisms by which complement C5a increases the capacity of polymorphonuclear MDSCs (PMN-MDSCs) to promote tumor growth and metastatic spread. Stimulation of PMN-MDSCs with C5a favored the invasion of cancer cells via a process dependent on the formation of neutrophil extracellular traps (NETs). NETosis was dependent on the production of high mobility group box 1 (HMGB1) by cancer cells. Moreover, C5a induced the surface expression of the HMGB1 receptors TLR4 and RAGE in PMN-MDSCs. In a mouse lung metastasis model, inhibition of C5a, C5a receptor-1 (C5aR1) or NETosis reduced the number of circulating-tumor cells (CTCs) and the metastatic burden. In support of the translational relevance of these findings, C5a was able to stimulate migration and NETosis in PMN-MDSCs obtained from lung cancer patients. Furthermore, myeloperoxidase (MPO)-DNA complexes, as markers of NETosis, were elevated in lung cancer patients and significantly correlated with C5a levels. In conclusion, C5a induces the formation of NETs from PMN-MDSCs in the presence of cancer cells, which may facilitate cancer cell dissemination and metastasis.

Tumor cells in light-chain amyloidosis and myeloma show distinct transcriptional rewiring of normal plasma cell development.

Although light-chain amyloidosis (AL) and multiple myeloma (MM) are characterized by tumor plasma cell (PC) expansion in bone marrow (BM), their clinical presentation differs. Previous attempts to identify unique pathogenic mechanisms behind such differences were unsuccessful, and no studies have investigated the differentiation stage of tumor PCs in patients with AL and MM. We sought to define a transcriptional atlas of normal PC development in secondary lymphoid organs (SLOs), peripheral blood (PB), and BM for comparison with the transcriptional programs (TPs) of tumor PCs in AL, MM, and monoclonal gammopathy of undetermined significance (MGUS). Based on bulk and single-cell RNA sequencing, we observed 13 TPs during transition of normal PCs throughout SLOs, PB, and BM. We further noted the following: CD39 outperforms CD19 to discriminate newborn from long-lived BM-PCs; tumor PCs expressed the most advantageous TPs of normal PC differentiation; AL shares greater similarity to SLO-PCs whereas MM is transcriptionally closer to PB-PCs and newborn BM-PCs; patients with AL and MM enriched in immature TPs had inferior survival; and protein N-linked glycosylation-related TPs are upregulated in AL. Collectively, we provide a novel resource to understand normal PC development and the transcriptional reorganization of AL and other monoclonal gammopathies.

Immunological Biomarkers of Fatal COVID-19: A Study of 868 Patients.

Information on the immunopathobiology of coronavirus disease 2019 (COVID-19) is rapidly increasing; however, there remains a need to identify immune features predictive of fatal outcome. This large-scale study characterized immune responses to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection using multidimensional flow cytometry, with the aim of identifying high-risk immune biomarkers. Holistic and unbiased analyses of 17 immune cell-types were conducted on 1,075 peripheral blood samples obtained from 868 COVID-19 patients and on samples from 24 patients presenting with non-SARS-CoV-2 infections and 36 healthy donors. Immune profiles of COVID-19 patients were significantly different from those of age-matched healthy donors but generally similar to those of patients with non-SARS-CoV-2 infections. Unsupervised clustering analysis revealed three immunotypes during SARS-CoV-2 infection; immunotype 1 (14% of patients) was characterized by significantly lower percentages of all immune cell-types except neutrophils and circulating plasma cells, and was significantly associated with severe disease. Reduced B-cell percentage was most strongly associated with risk of death. On multivariate analysis incorporating age and comorbidities, B-cell and non-classical monocyte percentages were independent prognostic factors for survival in training (n=513) and validation (n=355) cohorts. Therefore, reduced percentages of B-cells and non-classical monocytes are high-risk immune biomarkers for risk-stratification of COVID-19 patients.

Characterization of complete lncRNAs transcriptome reveals the functional and clinical impact of lncRNAs in multiple myeloma.

Multiple myeloma (MM) is an incurable disease, whose clinical heterogeneity makes its management challenging, highlighting the need for biological features to guide improved therapies. Deregulation of specific long non-coding RNAs (lncRNAs) has been shown in MM, nevertheless, the complete lncRNA transcriptome has not yet been elucidated. In this work, we identified 40,511 novel lncRNAs in MM samples. lncRNAs accounted for 82% of the MM transcriptome and were more heterogeneously expressed than coding genes. A total of 10,351 overexpressed and 9,535 downregulated lncRNAs were identified in MM patients when compared with normal bone-marrow plasma cells. Transcriptional dynamics study of lncRNAs in the context of normal B-cell maturation revealed 989 lncRNAs with exclusive expression in MM, among which 89 showed de novo epigenomic activation. Knockdown studies on one of these lncRNAs, SMILO (specific myeloma intergenic long non-coding RNA), resulted in reduced proliferation and induction of apoptosis of MM cells, and activation of the interferon pathway. We also showed that the expression of lncRNAs, together with clinical and genetic risk alterations, stratified MM patients into several progression-free survival and overall survival groups. In summary, our global analysis of the lncRNAs transcriptome reveals the presence of specific lncRNAs associated with the biological and clinical behavior of the disease.

Deep MRD profiling defines outcome and unveils different modes of treatment resistance in standard- and high-risk myeloma.

Patients with multiple myeloma (MM) carrying standard- or high-risk cytogenetic abnormalities (CAs) achieve similar complete response (CR) rates, but the later have inferior progression-free survival (PFS). This questions the legitimacy of CR as a treatment endpoint and represents a biological conundrum regarding the nature of tumor reservoirs that persist after therapy in high-risk MM. We used next-generation flow (NGF) cytometry to evaluate measurable residual disease (MRD) in MM patients with standard- vs high-risk CAs (n = 300 and 90, respectively) enrolled in the PETHEMA/GEM2012MENOS65 trial, and to identify mechanisms that determine MRD resistance in both patient subgroups (n = 40). The 36-month PFS rates were higher than 90% in patients with standard- or high-risk CAs achieving undetectable MRD. Persistent MRD resulted in a median PFS of âˆ¼3 and 2 years in patients with standard- and high-risk CAs, respectively. Further use of NGF to isolate MRD, followed by whole-exome sequencing of paired diagnostic and MRD tumor cells, revealed greater clonal selection in patients with standard-risk CAs, higher genomic instability with acquisition of new mutations in high-risk MM, and no unifying genetic event driving MRD resistance. Conversely, RNA sequencing of diagnostic and MRD tumor cells uncovered the selection of MRD clones with singular transcriptional programs and reactive oxygen species-mediated MRD resistance in high-risk MM. Our study supports undetectable MRD as a treatment endpoint for patients with MM who have high-risk CAs and proposes characterizing MRD clones to understand and overcome MRD resistance. This trial is registered at www.clinicaltrials.gov as #NCT01916252.

Immunologic characterization of COVID-19 patients with hematological cancer.

Not available.

Id1 and PD-1 Combined Blockade Impairs Tumor Growth and Survival of KRAS-mutant Lung Cancer by Stimulating PD-L1 Expression and Tumor Infiltrating CD8+ T Cells.

The use of PD-1/PD-L1 checkpoint inhibitors in advanced NSCLC is associated with longer survival. However, many patients do not benefit from PD-1/PD-L1 blockade, largely because of immunosuppression. New immunotherapy-based combinations are under investigation in an attempt to improve outcomes. Id1 (inhibitor of differentiation 1) is involved in immunosuppression. In this study, we explored the potential synergistic effect of the combination of Id1 inhibition and pharmacological PD-L1 blockade in three different syngeneic murine KRAS-mutant lung adenocarcinoma models. TCGA analysis demonstrated a negative and statistically significant correlation between PD-L1 and Id1 expression levels. This observation was confirmed in vitro in human and murine KRAS-driven lung cancer cell lines. In vivo experiments in KRAS-mutant syngeneic and metastatic murine lung adenocarcinoma models showed that the combined blockade targeting Id1 and PD-1 was more effective than each treatment alone in terms of tumor growth impairment and overall survival improvement. Mechanistically, multiplex quantification of CD3+/CD4+/CD8+ T cells and flow cytometry analysis showed that combined therapy favors tumor infiltration by CD8+ T cells, whilst in vivo CD8+ T cell depletion led to tumor growth restoration. Co-culture assays using CD8+ cells and tumor cells showed that T cells present a higher antitumor effect when tumor cells lack Id1 expression. These findings highlight that Id1 blockade may contribute to a significant immune enhancement of antitumor efficacy of PD-1 inhibitors by increasing PD-L1 expression and harnessing tumor infiltration of CD8+ T lymphocytes.