An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

KDM6A Regulates Immune Response Genes in Multiple Myeloma.

The histone H3K27 demethylase KDM6A is a tumor suppressor in multiple cancers, including multiple myeloma (MM). We created isogenic MM cells disrupted for KDM6A and tagged the endogenous protein to facilitate genome wide studies. KDM6A binds genes associated with immune recognition and cytokine signaling. Most importantly, KDM6A binds and activates NLRC5 and CIITA encoding regulators of Major Histocompatibility Complex (MHC) genes. Patient data indicate that NLRC5 and CIITA, are downregulated in MM with low KDM6A expression. Chromatin analysis shows that KDM6A binds poised and active enhancers and KDM6A loss led to decreased H3K27ac at enhancers, increased H3K27me3 levels in body of genes bound by KDM6A and decreased gene expression. Reestablishing histone acetylation with an HDAC3 inhibitor leads to upregulation of MHC expression, offering a strategy to restore immunogenicity of KDM6A deficient tumors. Loss of Kdm6a in murine RAS-transformed fibroblasts led to increased growth in vivo associated with decreased T cell infiltration. Statement of significance: We show that KDM6A participates in immune recognition of myeloma tumor cells by directly regulating the expression of the master regulators of MHC-I and II, NLRC5 and CIITA. The expression of these regulators can by rescued by the HDAC3 inhibitors in KDM6A-null cell lines.

Correction: Epigenetic, psychological, and EEG changes after a 1-week retreat based on mindfulness and compassion for stress reduction in healthy adults: Study protocol of a cross-over randomized controlled trial.

[This corrects the article DOI: 10.1371/journal.pone.0283169.].

Evaluation of clinical parameters and biomarkers in older, untreated mantle cell lymphoma patients receiving bendamustine-rituximab.

Mantle cell lymphoma (MCL) is clinically and biologically heterogeneous. While various prognostic features have been proposed, none currently impact therapy selection, particularly in older patients, for whom treatment is primarily dictated by age and comorbidities. Herein, we undertook a comprehensive comparison of clinicopathological features in a cohort of patients 60 years and older, uniformly treated with bendamustine and rituximab, with a median survival of >8 years. The strongest prognostic indicators in this cohort were a high-risk call by a simplified MCL international prognostic index (s-MIPI) (HR: 3.32, 95% CI: 1.65-6.68 compared to low risk), a high-risk call by MCL35 (HR: 10.34, 95% CI: 2.37-45.20 compared to low risk) and blastoid cytology (HR: 4.21, 95% CR: 1.92-9.22 compared to classic). Patients called high risk by both the s-MIPI and MCL35 had the most dismal prognosis (HR: 11.58, 95% CI: 4.10-32.72), while those with high risk by either had a moderate but clinically relevant prognosis (HR: 2.95, 95% CI: 1.49-5.82). A robust assay to assess proliferation, such as MCL35, along with stringent guidelines for cytological evaluation of MCL, in combination with MIPI, may be a strong path to risk-stratify older MCL patients in future clinical trials.

Epigenetic, psychological, and EEG changes after a 1-week retreat based on mindfulness and compassion for stress reduction in healthy adults: Study protocol of a cross-over randomized controlled trial.

INTRODUCTION: The main objective of the study will be to evaluate the effects of two widely used standardized mindfulness-based programs [Mindfulness-Based Stress Reduction (MBSR) and Compassion Cultivation Training (CCT)], on epigenetic, neurobiological, psychological, and physiological variables. METHODS: The programs will be offered in an intensive retreat format in a general population sample of healthy volunteer adults. During a 7-day retreat, participants will receive MBSR and CCT in a crossover design where participants complete both programs in random order. After finishing their first 3-day training with one of the two programs, participants will be assigned to the second 3-day training with the second program. The effects of the MBSR and CCT programs, and their combination, will be measured by epigenetic changes (i.e., DNA methylation biomarkers), neurobiological and psychophysiological measures (i.e., EEG resting state, EKG, respiration patterns, and diurnal cortisol slopes), self-report questionnaires belonging to different psychological domains (i.e., mindfulness, compassion, well-being, distress, and general functioning), and stress tasks (i.e., an Arithmetic Stress Test and the retrieval of negative autobiographical memories). These measures will be collected from both groups on the mornings of day 1 (pre-program), day 4 (after finishing the first program and before beginning the second program), and day 7 (post-second program). We will conduct a 3-month and a 12-month follow-up using only the set of self-report measures. DISCUSSION: This study aims to shed light on the neurobiological and psychological mechanisms linked to meditation and compassion in the general population. The protocol was registered at clinicaltrials.gov (Identifier: NCT05516355; August 23, 2022).

Genomic landscape of follicular lymphoma across a wide spectrum of clinical behaviors.

While some follicular lymphoma (FL) patients do not require treatment or experience prolonged responses, others relapse early, and little is known about genetic alterations specific to patients with a particular clinical behavior. We selected 56 grade 1-3A FL patients according to their need of treatment or timing of relapse: never treated (n = 7), non-relapsed (19), late relapse (14), early relapse or POD24 (11), and primary refractory (5). We analyzed 56 diagnostic and 12 paired relapse lymphoid tissue biopsies and performed copy number alteration (CNA) analysis and next generation sequencing (NGS). We identified six focal driver losses (1p36.32, 6p21.32, 6q14.1, 6q23.3, 9p21.3, 10q23.33) and 1p36.33 copy-neutral loss of heterozygosity (CN-LOH). By integrating CNA and NGS results, the most frequently altered genes/regions were KMT2D (79%), CREBBP (67%), TNFRSF14 (46%) and BCL2 (40%). Although we found that mutations in PIM1, FOXO1 and TMEM30A were associated with an adverse clinical behavior, definitive conclusions cannot be drawn, due to the small sample size. We identified common precursor cells harboring early oncogenic alterations of the KMT2D, CREBBP, TNFRSF14 and EP300 genes and 16p13.3-p13.2 CN-LOH. Finally, we established the functional consequences of mutations by means of protein modeling (CD79B, PLCG2, PIM1, MCL1 and IRF8). These data expand the knowledge on the genomics behind the heterogeneous FL population and, upon replication in larger cohorts, could contribute to risk stratification and the development of targeted therapies.

Genome-wide DNA methylation analysis in an antimigraine-treated preclinical model of cortical spreading depolarization.

BACKGROUND: Cortical spreading depolarization, the cause of migraine aura, is a short-lasting depolarization wave that moves across the brain cortex, transiently suppressing neuronal activity. Prophylactic treatments for migraine, such as topiramate or valproate, reduce the number of cortical spreading depression events in rodents. OBJECTIVE: To investigate whether cortical spreading depolarization with and without chronic treatment with topiramate or valproate affect the DNA methylation of the cortex. METHODS: Sprague-Dawley rats were intraperitoneally injected with saline, topiramate or valproate for four weeks when cortical spreading depolarization were induced and genome-wide DNA methylation was performed in the cortex of six rats per group. RESULTS: The DNA methylation profile of the cortex was significantly modified after cortical spreading depolarization, with and without topiramate or valproate. Interestingly, topiramate reduced by almost 50% the number of differentially methylated regions, whereas valproate increased them by 17%, when comparing to the non-treated group after cortical spreading depolarization induction. The majority of the differentially methylated regions lay within intragenic regions, and the analyses of functional group over-representation retrieved several enriched functions, including functions related to protein processing in the cortical spreading depolarization without treatment group; functions related to metabolic processes in the cortical spreading depolarization with topiramate group; and functions related to synapse and ErbB, MAPK or retrograde endocannabinoid signaling in the cortical spreading depolarization with valproate group. CONCLUSIONS: Our results may provide insights into the underlying physiological mechanisms of migraine with aura and emphasize the role of epigenetics in migraine susceptibility.

Tumorigenic role of Musashi-2 in aggressive mantle cell lymphoma.

SOX11 overexpression has been associated with aggressive behavior of mantle cell lymphomas (MCL). SOX11 is overexpressed in embryonic and cancer stem cells (CSC) of some tumors. Although CSC have been isolated from primary MCL, their relationship to SOX11 expression and contribution to MCL pathogenesis and clinical evolution remain unknown. Here, we observed enrichment in leukemic and hematopoietic stem cells gene signatures in SOX11+ compared to SOX11- MCL primary cases. Musashi-2 (MSI2) emerged as one of the most significant upregulated stem cell-related genes in SOX11+ MCLs. SOX11 is directly bound to the MSI2 promoter upregulating its expression in vitro. MSI2 intronic enhancers were strongly activated in SOX11+ MCL cell lines and primary cases. MSI2 upregulation was significantly associated with poor overall survival independently of other high-risk features of MCL. MSI2 knockdown decreased the expression of genes related to apoptosis and stem cell features and significantly reduced clonogenic growth, tumor cell survival and chemoresistance in MCL cells. MSI2-knockdown cells had reduced tumorigenic engraftment into mice bone marrow and spleen compared to control cells in xenotransplanted mouse models. Our results suggest that MSI2 might play a key role in sustaining stemness and tumor cell survival, representing a possible novel target for therapeutic interventions in MCL.

Integrative epigenomics in chronic lymphocytic leukaemia: Biological insights and clinical applications.

Chronic lymphocytic leukaemia (CLL) is not only characterised by driver genetic alterations but by extensive epigenetic changes. Over the last decade, epigenomic studies have described the DNA methylome, chromatin accessibility, histone modifications and the three-dimensional (3D) genome architecture of CLL. Beyond its regulatory role, the DNA methylome contains imprints of the cellular origin and proliferative history of CLL cells. These two aspects are strong independent prognostic factors. Integrative analyses of chromatin marks have uncovered novel regulatory elements and altered transcription factor networks as non-genetic means mediating gene deregulation in CLL. Additionally, CLL cells display a disease-specific pattern of 3D genome interactions. From the technological perspective, we are currently witnessing a transition from bulk omics to single-cell analyses. This review aims at summarising the major findings from the epigenomics field as well as providing a prospect of the present and future of single-cell analyses in CLL.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Remission of obesity and insulin resistance is not sufficient to restore mitochondrial homeostasis in visceral adipose tissue.

Metabolic plasticity is the ability of a biological system to adapt its metabolic phenotype to different environmental stressors. We used a whole-body and tissue-specific phenotypic, functional, proteomic, metabolomic and transcriptomic approach to systematically assess metabolic plasticity in diet-induced obese mice after a combined nutritional and exercise intervention. Although most obesity and overnutrition-related pathological features were successfully reverted, we observed a high degree of metabolic dysfunction in visceral white adipose tissue, characterized by abnormal mitochondrial morphology and functionality. Despite two sequential therapeutic interventions and an apparent global healthy phenotype, obesity triggered a cascade of events in visceral adipose tissue progressing from mitochondrial metabolic and proteostatic alterations to widespread cellular stress, which compromises its biosynthetic and recycling capacity. In humans, weight loss after bariatric surgery showed a transcriptional signature in visceral adipose tissue similar to our mouse model of obesity reversion. Overall, our data indicate that obesity prompts a lasting metabolic fingerprint that leads to a progressive breakdown of metabolic plasticity in visceral adipose tissue.

Comprehensive characterization of the epigenetic landscape in Multiple Myeloma.

Background: Human multiple myeloma (MM) cell lines (HMCLs) have been widely used to understand the molecular processes that drive MM biology. Epigenetic modifications are involved in MM development, progression, and drug resistance. A comprehensive characterization of the epigenetic landscape of MM would advance our understanding of MM pathophysiology and may attempt to identify new therapeutic targets. Methods: We performed chromatin immunoprecipitation sequencing to analyze histone mark changes (H3K4me1, H3K4me3, H3K9me3, H3K27ac, H3K27me3 and H3K36me3) on 16 HMCLs. Results: Differential analysis of histone modification profiles highlighted links between histone modifications and cytogenetic abnormalities or recurrent mutations. Using histone modifications associated to enhancer regions, we identified super-enhancers (SE) associated with genes involved in MM biology. We also identified promoters of genes enriched in H3K9me3 and H3K27me3 repressive marks associated to potential tumor suppressor functions. The prognostic value of genes associated with repressive domains and SE was used to build two distinct scores identifying high-risk MM patients in two independent cohorts (CoMMpass cohort; n = 674 and Montpellier cohort; n = 69). Finally, we explored H3K4me3 marks comparing drug-resistant and -sensitive HMCLs to identify regions involved in drug resistance. From these data, we developed epigenetic biomarkers based on the H3K4me3 modification predicting MM cell response to lenalidomide and histone deacetylase inhibitors (HDACi). Conclusions: The epigenetic landscape of MM cells represents a unique resource for future biological studies. Furthermore, risk-scores based on SE and repressive regions together with epigenetic biomarkers of drug response could represent new tools for precision medicine in MM.

PRC2 Inhibitors Overcome Glucocorticoid Resistance Driven by NSD2 Mutation in Pediatric Acute Lymphoblastic Leukemia.

Mutations in epigenetic regulators are common in relapsed pediatric acute lymphoblastic leukemia (ALL). Here, we uncovered the mechanism underlying the relapse of ALL driven by an activating mutation of the NSD2 histone methyltransferase (p.E1099K). Using high-throughput drug screening, we found that NSD2-mutant cells were specifically resistant to glucocorticoids. Correction of this mutation restored glucocorticoid sensitivity. The transcriptional response to glucocorticoids was blocked in NSD2-mutant cells due to depressed glucocorticoid receptor (GR) levels and the failure of glucocorticoids to autoactivate GR expression. Although H3K27me3 was globally decreased by NSD2 p.E1099K, H3K27me3 accumulated at the NR3C1 (GR) promoter. Pretreatment of NSD2 p.E1099K cell lines and patient-derived xenograft samples with PRC2 inhibitors reversed glucocorticoid resistance in vitro and in vivo. PRC2 inhibitors restored NR3C1 autoactivation by glucocorticoids, increasing GR levels and allowing GR binding and activation of proapoptotic genes. These findings suggest a new therapeutic approach to relapsed ALL associated with NSD2 mutation. SIGNIFICANCE: NSD2 histone methyltransferase mutations observed in relapsed pediatric ALL drove glucocorticoid resistance by repression of the GR and abrogation of GR gene autoactivation due to accumulation of K3K27me3 at its promoter. Pretreatment with PRC2 inhibitors reversed resistance, suggesting a new therapeutic approach to these patients with ALL.This article is highlighted in the In This Issue feature, p. 1.

Gene expression derived from alternative promoters improves prognostic stratification in multiple myeloma.

Clinical and genetic risk factors are currently used in multiple myeloma (MM) to stratify patients and to design specific therapies. However, these systems do not capture the heterogeneity of the disease supporting the development of new prognostic factors. In this study, we identified active promoters and alternative active promoters in 6 different B cell subpopulations, including bone-marrow plasma cells, and 32 MM patient samples, using RNA-seq data. We find that expression initiated at both regular and alternative promoters was specific of each B cell subpopulation or MM plasma cells, showing a remarkable level of consistency with chromatin-based promoter definition. Interestingly, using 595 MM patient samples from the CoMMpass dataset, we observed that the expression derived from some alternative promoters was associated with lower progression-free and overall survival in MM patients independently of genetic alterations. Altogether, our results define cancer-specific alternative active promoters as new transcriptomic features that can provide a new avenue for prognostic stratification possibilities in patients with MM.

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.

Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation.

To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.