T-bet suppresses proliferation of malignant B cells in chronic lymphocytic leukemia.

The T-box transcription factor T-bet is known as a master regulator of T-cell response but its role in malignant B cells is not sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with genetic knockout of TBX21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity induced by inflammatory signals provided by the microenvironment, triggered T-bet expression which impacted on promoter proximal and distal chromatin co-accessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling, and a negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of CLL patients. Our study uncovers a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling which has implications for stratification and therapy of CLL patients. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in inflammatory signaling pathways in CLL.

A multidimensional atlas of human glioblastoma-like organoids reveals highly coordinated molecular networks and effective drugs.

Recent advances in the genomics of glioblastoma (GBM) led to the introduction of molecular neuropathology but failed to translate into treatment improvement. This is largely attributed to the genetic and phenotypic heterogeneity of GBM, which are considered the major obstacle to GBM therapy. Here, we use advanced human GBM-like organoid (LEGO: Laboratory Engineered Glioblastoma-like Organoid) models and provide an unprecedented comprehensive characterization of LEGO models using single-cell transcriptome, DNA methylome, metabolome, lipidome, proteome, and phospho-proteome analysis. We discovered that genetic heterogeneity dictates functional heterogeneity across molecular layers and demonstrates that NF1 mutation drives mesenchymal signature. Most importantly, we found that glycerol lipid reprogramming is a hallmark of GBM, and several targets and drugs were discovered along this line. We also provide a genotype-based drug reference map using LEGO-based drug screen. This study provides new human GBM models and a research path toward effective GBM therapy.

Compartments in medulloblastoma with extensive nodularity are connected through differentiation along the granular precursor lineage.

Medulloblastomas with extensive nodularity are cerebellar tumors characterized by two distinct compartments and variable disease progression. The mechanisms governing the balance between proliferation and differentiation in MBEN remain poorly understood. Here, we employ a multi-modal single cell transcriptome analysis to dissect this process. In the internodular compartment, we identify proliferating cerebellar granular neuronal precursor-like malignant cells, along with stromal, vascular, and immune cells. In contrast, the nodular compartment comprises postmitotic, neuronally differentiated malignant cells. Both compartments are connected through an intermediate cell stage resembling actively migrating CGNPs. Notably, we also discover astrocytic-like malignant cells, found in proximity to migrating and differentiated cells at the transition zone between the two compartments. Our study sheds light on the spatial tissue organization and its link to the developmental trajectory, resulting in a more benign tumor phenotype. This integrative approach holds promise to explore intercompartmental interactions in other cancers with varying histology.

Nucleosome repositioning in chronic lymphocytic leukemia.

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.

Resolving the spatial architecture of myeloma and its microenvironment at the single-cell level.

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Resolving therapy resistance mechanisms in multiple myeloma by multiomics subclone analysis.

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them.

Transplanting rejuvenated blood stem cells extends lifespan of aged immunocompromised mice.

One goal of regenerative medicine is to rejuvenate tissues and extend lifespan by restoring the function of endogenous aged stem cells. However, evidence that somatic stem cells can be targeted in vivo to extend lifespan is still lacking. Here, we demonstrate that after a short systemic treatment with a specific inhibitor of the small RhoGTPase Cdc42 (CASIN), transplanting aged hematopoietic stem cells (HSCs) from treated mice is sufficient to extend the healthspan and lifespan of aged immunocompromised mice without additional treatment. In detail, we show that systemic CASIN treatment improves strength and endurance of aged mice by increasing the myogenic regenerative potential of aged skeletal muscle stem cells. Further, we show that CASIN modifies niche localization and H4K16ac polarity of HSCs in vivo. Single-cell profiling reveals changes in HSC transcriptome, which underlie enhanced lymphoid and regenerative capacity in serial transplantation assays. Overall, we provide proof-of-concept evidence that a short systemic treatment to decrease Cdc42 activity improves the regenerative capacity of different endogenous aged stem cells in vivo, and that rejuvenated HSCs exert a broad systemic effect sufficient to extend murine health- and lifespan.

Differentiation-related epigenomic changes define clinically distinct keratinocyte cancer subclasses.

Keratinocyte cancers (KC) are the most prevalent malignancies in fair-skinned populations, posing a significant medical and economic burden to health systems. KC originate in the epidermis and mainly comprise basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC). Here, we combined single-cell multi-omics, transcriptomics, and methylomics to investigate the epigenomic dynamics during epidermal differentiation. We identified ~3,800 differentially accessible regions between undifferentiated and differentiated keratinocytes, corresponding to regulatory regions associated with key transcription factors. DNA methylation at these regions defined AK/cSCC subtypes with epidermal stem cell- or keratinocyte-like features. Using cell-type deconvolution tools and integration of bulk and single-cell methylomes, we demonstrate that these subclasses are consistent with distinct cells-of-origin. Further characterization of the phenotypic traits of the subclasses and the study of additional unstratified KC entities uncovered distinct clinical features for the subclasses, linking invasive and metastatic KC cases with undifferentiated cells-of-origin. Our study provides a thorough characterization of the epigenomic dynamics underlying human keratinocyte differentiation and uncovers novel links between KC cells-of-origin and their prognosis.

Inflammatory exposure drives long-lived impairment of hematopoietic stem cell self-renewal activity and accelerated aging.

Hematopoietic stem cells (HSCs) mediate regeneration of the hematopoietic system following injury, such as following infection or inflammation. These challenges impair HSC function, but whether this functional impairment extends beyond the duration of inflammatory exposure is unknown. Unexpectedly, we observed an irreversible depletion of functional HSCs following challenge with inflammation or bacterial infection, with no evidence of any recovery up to 1 year afterward. HSCs from challenged mice demonstrated multiple cellular and molecular features of accelerated aging and developed clinically relevant blood and bone marrow phenotypes not normally observed in aged laboratory mice but commonly seen in elderly humans. In vivo HSC self-renewal divisions were absent or extremely rare during both challenge and recovery periods. The progressive, irreversible attrition of HSC function demonstrates that temporally discrete inflammatory events elicit a cumulative inhibitory effect on HSCs. This work positions early/mid-life inflammation as a mediator of lifelong defects in tissue maintenance and regeneration.

Subclone-specific microenvironmental impact and drug response in refractory multiple myeloma revealed by single-cell transcriptomics.

Virtually all patients with multiple myeloma become unresponsive to treatment over time. Relapsed/refractory multiple myeloma (RRMM) is accompanied by the clonal evolution of myeloma cells with heterogeneous genomic aberrations and profound changes of the bone marrow microenvironment (BME). However, the molecular mechanisms that drive drug resistance remain elusive. Here, we analyze the heterogeneous tumor cell population and its complex interaction network with the BME of 20 RRMM patients by single cell RNA-sequencing before/after treatment. Subclones with chromosome 1q-gain express a specific transcriptomic signature and frequently expand during treatment. Furthermore, RRMM cells shape an immune suppressive BME by upregulation of inflammatory cytokines and close interaction with the myeloid compartment. It is characterized by the accumulation of PD1+ γδ T-cells and tumor-associated macrophages as well as the depletion of hematopoietic progenitors. Thus, our study resolves transcriptional features of subclones in RRMM and mechanisms of microenvironmental reprogramming with implications for clinical decision-making.

The DNA methylation landscape of multiple myeloma shows extensive inter- and intrapatient heterogeneity that fuels transcriptomic variability.

BACKGROUND: Cancer evolution depends on epigenetic and genetic diversity. Historically, in multiple myeloma (MM), subclonal diversity and tumor evolution have been investigated mostly from a genetic perspective. METHODS: Here, we performed an analysis of 42 MM samples from 21 patients by using enhanced reduced representation bisulfite sequencing (eRRBS). We combined several metrics of epigenetic heterogeneity to analyze DNA methylation heterogeneity in MM patients. RESULTS: We show that MM is characterized by the continuous accumulation of stochastic methylation at the promoters of development-related genes. High combinatorial entropy change is associated with poor outcomes in our pilot study and depends predominantly on partially methylated domains (PMDs). These PMDs, which represent the major source of inter- and intrapatient DNA methylation heterogeneity in MM, are linked to other key epigenetic aberrations, such as CpG island (CGI)/transcription start site (TSS) hypermethylation and H3K27me3 redistribution as well as 3D organization alterations. In addition, transcriptome analysis revealed that intratumor methylation heterogeneity was associated with low-level expression and high variability. CONCLUSIONS: We propose that disrupted DNA methylation in MM is responsible for high epigenetic and transcriptomic instability allowing tumor cells to adapt to environmental changes by tapping into a pool of evolutionary trajectories.

Downregulation of TGR5 (GPBAR1) in biliary epithelial cells contributes to the pathogenesis of sclerosing cholangitis.

BACKGROUND & AIMS: Primary sclerosing cholangitis (PSC) is characterized by chronic inflammation and progressive fibrosis of the biliary tree. The bile acid receptor TGR5 (GPBAR1) is found on biliary epithelial cells (BECs), where it promotes secretion, proliferation and tight junction integrity. Thus, we speculated that changes in TGR5-expression in BECs may contribute to PSC pathogenesis. METHODS: TGR5-expression and -localization were analyzed in PSC livers and liver tissue, isolated bile ducts and BECs from Abcb4-/-, Abcb4-/-/Tgr5Tg and ursodeoxycholic acid (UDCA)- or 24-norursodeoxycholic acid (norUDCA)-fed Abcb4-/- mice. The effects of IL8/IL8 homologues on TGR5 mRNA and protein levels were studied. BEC gene expression was analyzed by single-cell transcriptomics (scRNA-seq) from distinct mouse models. RESULTS: TGR5 mRNA expression and immunofluorescence staining intensity were reduced in BECs of PSC and Abcb4-/- livers, in Abcb4-/- extrahepatic bile ducts, but not in intrahepatic macrophages. No changes in TGR5 BEC fluorescence intensity were detected in liver tissue of other liver diseases, including primary biliary cholangitis. Incubation of BECs with IL8/IL8 homologues, but not with other cytokines, reduced TGR5 mRNA and protein levels. BECs from Abcb4-/- mice had lower levels of phosphorylated Erk and higher expression levels of Icam1, Vcam1 and Tgfß2. Overexpression of Tgr5 abolished the activated inflammatory phenotype characteristic of Abcb4-/- BECs. NorUDCA-feeding restored TGR5-expression levels in BECs in Abcb4-/- livers. CONCLUSIONS: Reduced TGR5 levels in BECs from patients with PSC and Abcb4-/- mice promote development of a reactive BEC phenotype, aggravate biliary injury and thus contribute to the pathogenesis of sclerosing cholangitis. Restoration of biliary TGR5-expression levels represents a previously unknown mechanism of action of norUDCA. LAY SUMMARY: Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease-associated with progressive inflammation of the bile duct, leading to fibrosis and end-stage liver disease. Bile acid (BA) toxicity may contribute to the development and disease progression of PSC. TGR5 is a membrane-bound receptor for BAs, which is found on bile ducts and protects bile ducts from BA toxicity. In this study, we show that TGR5 levels were reduced in bile ducts from PSC livers and in bile ducts from a genetic mouse model of PSC. Our investigations indicate that lower levels of TGR5 in bile ducts may contribute to PSC development and progression. Furthermore, treatment with norUDCA, a drug currently being tested in a phase III trial for PSC, restored TGR5 levels in biliary epithelial cells.

Auto-aggressive CXCR6+ CD8 T cells cause liver immune pathology in NASH.

Nonalcoholic steatohepatitis (NASH) is a manifestation of systemic metabolic disease related to obesity, and causes liver disease and cancer1,2. The accumulation of metabolites leads to cell stress and inflammation in the liver3, but mechanistic understandings of liver damage in NASH are incomplete. Here, using a preclinical mouse model that displays key features of human NASH (hereafter, NASH mice), we found an indispensable role for T cells in liver immunopathology. We detected the hepatic accumulation of CD8 T cells with phenotypes that combined tissue residency (CXCR6) with effector (granzyme) and exhaustion (PD1) characteristics. Liver CXCR6+ CD8 T cells were characterized by low activity of the FOXO1 transcription factor, and were abundant in NASH mice and in patients with NASH. Mechanistically, IL-15 induced FOXO1 downregulation and CXCR6 upregulation, which together rendered liver-resident CXCR6+ CD8 T cells susceptible to metabolic stimuli (including acetate and extracellular ATP) and collectively triggered auto-aggression. CXCR6+ CD8 T cells from the livers of NASH mice or of patients with NASH had similar transcriptional signatures, and showed auto-aggressive killing of cells in an MHC-class-I-independent fashion after signalling through P2X7 purinergic receptors. This killing by auto-aggressive CD8 T cells fundamentally differed from that by antigen-specific cells, which mechanistically distinguishes auto-aggressive and protective T cell immunity.

Hematopoietic stem and progenitor cell-restricted Cdx2 expression induces transformation to myelodysplasia and acute leukemia.

The caudal-related homeobox transcription factor CDX2 is expressed in leukemic cells but not during normal blood formation. Retroviral overexpression of Cdx2 induces AML in mice, however the developmental stage at which CDX2 exerts its effect is unknown. We developed a conditionally inducible Cdx2 mouse model to determine the effects of in vivo, inducible Cdx2 expression in hematopoietic stem and progenitor cells (HSPCs). Cdx2-transgenic mice develop myelodysplastic syndrome with progression to acute leukemia associated with acquisition of additional driver mutations. Cdx2-expressing HSPCs demonstrate enrichment of hematopoietic-specific enhancers associated with pro-differentiation transcription factors. Furthermore, treatment of Cdx2 AML with azacitidine decreases leukemic burden. Extended scheduling of low-dose azacitidine shows greater efficacy in comparison to intermittent higher-dose azacitidine, linked to more specific epigenetic modulation. Conditional Cdx2 expression in HSPCs is an inducible model of de novo leukemic transformation and can be used to optimize treatment in high-risk AML.

Dissecting intratumour heterogeneity of nodal B-cell lymphomas at the transcriptional, genetic and drug-response levels.

Tumour heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is known to affect the efficacy of cancer therapy, most personalized treatment approaches do not account for intratumour heterogeneity. We addressed this issue by studying the heterogeneity of nodal B-cell lymphomas by single-cell RNA-sequencing and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subpopulations and compared their drug-response and genomic profiles. Malignant subpopulations from the same patient responded strikingly differently to anti-cancer drugs ex vivo, which recapitulated subpopulation-specific drug sensitivity during in vivo treatment. Infiltrating T cells represented the majority of non-malignant cells, whose gene-expression signatures were similar across all donors, whereas the frequencies of T-cell subsets varied significantly between the donors. Our data provide insights into the heterogeneity of nodal B-cell lymphomas and highlight the relevance of intratumour heterogeneity for personalized cancer therapy.

Methylome-based cell-of-origin modeling (Methyl-COOM) identifies aberrant expression of immune regulatory molecules in CLL.

BACKGROUND: In cancer, normal epigenetic patterns are disturbed and contribute to gene expression changes, disease onset, and progression. The cancer epigenome is composed of the epigenetic patterns present in the tumor-initiating cell at the time of transformation, and the tumor-specific epigenetic alterations that are acquired during tumor initiation and progression. The precise dissection of these two components of the tumor epigenome will facilitate a better understanding of the biological mechanisms underlying malignant transformation. Chronic lymphocytic leukemia (CLL) originates from differentiating B cells, which undergo extensive epigenetic programming. This poses the challenge to precisely determine the epigenomic ground state of the cell-of-origin in order to identify CLL-specific epigenetic aberrations. METHODS: We developed a linear regression model, methylome-based cell-of-origin modeling (Methyl-COOM), to map the cell-of-origin for individual CLL patients based on the continuum of epigenomic changes during normal B cell differentiation. RESULTS: Methyl-COOM accurately maps the cell-of-origin of CLL and identifies CLL-specific aberrant DNA methylation events that are not confounded by physiologic epigenetic B cell programming. Furthermore, Methyl-COOM unmasks abnormal action of transcription factors, altered super-enhancer activities, and aberrant transcript expression in CLL. Among the aberrantly regulated transcripts were many genes that have previously been implicated in T cell biology. Flow cytometry analysis of these markers confirmed their aberrant expression on malignant B cells at the protein level. CONCLUSIONS: Methyl-COOM analysis of CLL identified disease-specific aberrant gene regulation. The aberrantly expressed genes identified in this study might play a role in immune-evasion in CLL and might serve as novel targets for immunotherapy approaches. In summary, we propose a novel framework for in silico modeling of reference DNA methylomes and for the identification of cancer-specific epigenetic changes, a concept that can be broadly applied to other human malignancies.

Glioblastoma initiating cells are sensitive to histone demethylase inhibition due to epigenetic deregulation.

Tumor-initiating cells are a subpopulation of cells that have self-renewal capacity to regenerate a tumor. Here, we identify stem cell-like chromatin features in human glioblastoma initiating cells (GICs) and link them to a loss of the repressive histone H3 lysine 9 trimethylation (H3K9me3) mark. Increasing H3K9me3 levels by histone demethylase inhibition led to cell death in GICs but not in their differentiated counterparts. The induction of apoptosis was accompanied by a loss of the activating H3 lysine 9 acetylation (H3K9ac) modification and accumulation of DNA damage and downregulation of DNA damage response genes. Upon knockdown of histone demethylases, KDM4C and KDM7A both differentiation and DNA damage were induced. Thus, the H3K9me3-H3K9ac equilibrium is crucial for GIC viability and represents a chromatin feature that can be exploited to specifically target this tumor subpopulation.

Single-nucleus chromatin accessibility reveals intratumoral epigenetic heterogeneity in IDH1 mutant gliomas.

The presence of genome-wide DNA hypermethylation is a hallmark of lower grade gliomas (LGG) with isocitrate dehydrogenase (IDH) mutations. Further molecular classification of IDH mutant gliomas is defined by the presence (IDHmut-codel) or absence (IDHmut-noncodel) of hemizygous codeletion of chromosome arms 1p and 19q. Despite the DNA hypermethylation seen in bulk tumors, intra-tumoral heterogeneity at the epigenetic level has not been thoroughly analyzed. To address this question, we performed the first epigenetic profiling of single cells in a cohort of 5 gliomas with IDH1 mutation using single nucleus Assay for Transposase-Accessible Chromatin with high-throughput sequencing (snATAC-seq). Using the Fluidigm HT IFC microfluidics platform, we generated chromatin accessibility maps from 336 individual nuclei, and identified variable promoter accessibility of non-coding RNAs in LGGs. Interestingly, local chromatin structures of several non-coding RNAs are significant factors that contribute to heterogeneity, and show increased promoter accessibility in IDHmut-noncodel samples. As an example for clinical significance of this result, we identify CYTOR as a poor prognosis factor in gliomas with IDH mutation. Open chromatin assay points to differential accessibility of non-coding RNAs as an important source of epigenetic heterogeneity within individual tumors and between molecular subgroups. Rare populations of nuclei that resemble either IDH mutant molecular group co-exist within IDHmut-noncodel and IDHmut-codel groups, and along with non-coding RNAs may be an important issue to consider for future studies, as they may help guide predict treatment response and relapse.A web-based explorer for the data is available at shiny.turcanlab.org.

Haematopoietic stem cells in perisinusoidal niches are protected from ageing.

With ageing, intrinsic haematopoietic stem cell (HSC) activity decreases, resulting in impaired tissue homeostasis, reduced engraftment following transplantation and increased susceptibility to diseases. However, whether ageing also affects the HSC niche, and thereby impairs its capacity to support HSC function, is still widely debated. Here, by using in-vivo long-term label-retention assays we demonstrate that aged label-retaining HSCs, which are, in old mice, the most quiescent HSC subpopulation with the highest regenerative capacity and cellular polarity, reside predominantly in perisinusoidal niches. Furthermore, we demonstrate that sinusoidal niches are uniquely preserved in shape, morphology and number on ageing. Finally, we show that myeloablative chemotherapy can selectively disrupt aged sinusoidal niches in the long term, which is linked to the lack of recovery of endothelial Jag2 at sinusoids. Overall, our data characterize the functional alterations of the aged HSC niche and unveil that perisinusoidal niches are uniquely preserved and thereby protect HSCs from ageing.