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.

Chromatin activation as a unifying principle underlying pathogenic mechanisms in multiple myeloma.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

Cereblon enhancer methylation and IMiD resistance in multiple myeloma.

Cereblon is the direct binding target of the immunomodulatory drugs (IMiDs) that are commonly used to treat multiple myeloma (MM), the second most frequent hematologic malignancy. Patients respond well to initial treatment with IMiDs, but virtually all patients develop drug resistance over time, and the underlying mechanisms are poorly understood. We identified an as yet undescribed DNA hypermethylation in an active intronic CRBN enhancer. Differential hypermethylation in this region was found to be increased in healthy plasma cells, but was more pronounced in IMiD-refractory MM. Methylation significantly correlated with decreased CRBN expression levels. DNA methyltransferase inhibitor (DNTMi) in vitro experiments induced CRBN enhancer demethylation, and sensitizing effects on lenalidomide treatment were observed in 2 MM cell lines. Thus, we provide first evidence that aberrant CRBN DNA methylation is a novel mechanism of IMiD resistance in MM and may predict IMiD response prior to treatment.

A network-based approach to integrate nutrient microenvironment in the prediction of synthetic lethality in cancer metabolism.

Synthetic Lethality (SL) is currently defined as a type of genetic interaction in which the loss of function of either of two genes individually has limited effect in cell viability but inactivation of both genes simultaneously leads to cell death. Given the profound genomic aberrations acquired by tumor cells, which can be systematically identified with -omics data, SL is a promising concept in cancer research. In particular, SL has received much attention in the area of cancer metabolism, due to the fact that relevant functional alterations concentrate on key metabolic pathways that promote cellular proliferation. With the extensive prior knowledge about human metabolic networks, a number of computational methods have been developed to predict SL in cancer metabolism, including the genetic Minimal Cut Sets (gMCSs) approach. A major challenge in the application of SL approaches to cancer metabolism is to systematically integrate tumor microenvironment, given that genetic interactions and nutritional availability are interconnected to support proliferation. Here, we propose a more general definition of SL for cancer metabolism that combines genetic and environmental interactions, namely loss of gene functions and absence of nutrients in the environment. We extend our gMCSs approach to determine this new family of metabolic synthetic lethal interactions. A computational and experimental proof-of-concept is presented for predicting the lethality of dihydrofolate reductase (DHFR) inhibition in different environments. Finally, our approach is applied to identify extracellular nutrient dependences of tumor cells, elucidating cholesterol and myo-inositol depletion as potential vulnerabilities in different malignancies.

BOSO: A novel feature selection algorithm for linear regression with high-dimensional data.

With the frenetic growth of high-dimensional datasets in different biomedical domains, there is an urgent need to develop predictive methods able to deal with this complexity. Feature selection is a relevant strategy in machine learning to address this challenge. We introduce a novel feature selection algorithm for linear regression called BOSO (Bilevel Optimization Selector Operator). We conducted a benchmark of BOSO with key algorithms in the literature, finding a superior accuracy for feature selection in high-dimensional datasets. Proof-of-concept of BOSO for predicting drug sensitivity in cancer is presented. A detailed analysis is carried out for methotrexate, a well-studied drug targeting cancer metabolism.

Immunogenomic identification and characterization of granulocytic myeloid-derived suppressor cells in multiple myeloma.

Granulocytic myeloid-derived suppressor cells (G-MDSCs) promote tumor growth and immunosuppression in multiple myeloma (MM). However, their phenotype is not well established for accurate monitoring or clinical translation. We aimed to provide the phenotypic profile of G-MDSCs based on their prognostic significance in MM, immunosuppressive potential, and molecular program. The preestablished phenotype of G-MDSCs was evaluated in bone marrow samples from controls and MM patients using multidimensional flow cytometry; surprisingly, we found that CD11b+CD14-CD15+CD33+HLADR- cells overlapped with common eosinophils and neutrophils, which were not expanded in MM patients. Therefore, we relied on automated clustering to unbiasedly identify all granulocytic subsets in the tumor microenvironment: basophils, eosinophils, and immature, intermediate, and mature neutrophils. In a series of 267 newly diagnosed MM patients (GEM2012MENOS65 trial), only the frequency of mature neutrophils at diagnosis was significantly associated with patient outcome, and a high mature neutrophil/T-cell ratio resulted in inferior progression-free survival (P < .001). Upon fluorescence-activated cell sorting of each neutrophil subset, T-cell proliferation decreased in the presence of mature neutrophils (0.5-fold; P = .016), and the cytotoxic potential of T cells engaged by a BCMA×CD3-bispecific antibody increased notably with the depletion of mature neutrophils (fourfold; P = .0007). Most interestingly, RNA sequencing of the 3 subsets revealed that G-MDSC-related genes were specifically upregulated in mature neutrophils from MM patients vs controls because of differential chromatin accessibility. Taken together, our results establish a correlation between the clinical significance, immunosuppressive potential, and transcriptional network of well-defined neutrophil subsets, providing for the first time a set of optimal markers (CD11b/CD13/CD16) for accurate monitoring of G-MDSCs in MM.

Upregulated expression and function of the α4ß1 integrin in multiple myeloma cells resistant to bortezomib.

The interaction of multiple myeloma (MM) cells with the bone marrow (BM) microenvironment promotes MM cell retention, survival, and resistance to different anti-MM agents, including proteasome inhibitors (PIs) such as bortezomib (BTZ). The α4ß1 integrin is a main adhesion receptor mediating MM cell-stroma interactions and MM cell survival, and its expression and function are downregulated by BTZ, leading to inhibition of cell adhesion-mediated drug resistance (CAM-DR) and MM cell apoptosis. Whether decreased α4ß1 expression and activity are maintained or recovered upon development of resistance to BTZ represents an important question, as a potential rescue of α4ß1 function could boost MM cell survival and disease progression. Using BTZ-resistant MM cells, we found that they not only rescue their α4ß1 expression, but its levels were higher than in parental cells. Increased α4ß1 expression in resistant cells correlated with enhanced α4ß1-mediated cell lodging in the BM, and with disease progression. BTZ-resistant MM cells displayed enhanced NF-κB pathway activation relative to parental counterparts, which contributed to upregulated α4 expression and to α4ß1-dependent MM cell adhesion. These data emphasize the upregulation of α4ß1 expression and function as a key event during resistance to BTZ in MM, which might indirectly contribute to stabilize this resistance, as stronger MM cell attachment to BM stroma will regain CAM-DR and MM cell growth and survival. Finally, we found a strong correlation between high ITGB1 (integrin ß1) expression in MM and poor progression-free survival (PFS) and overall survival (OS) during treatment of MM patients with BTZ and IMIDs, and combination of high ITGB1 levels and presence of the high-risk genetic factor amp1q causes low PFS and OS. These results unravel a novel prognostic value for ITGB1 in myeloma. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Dual Targeting of Histone Methyltransferase G9a and DNA-Methyltransferase 1 for the Treatment of Experimental Hepatocellular Carcinoma.

Epigenetic modifications such as DNA and histone methylation functionally cooperate in fostering tumor growth, including that of hepatocellular carcinoma (HCC). Pharmacological targeting of these mechanisms may open new therapeutic avenues. We aimed to determine the therapeutic efficacy and potential mechanism of action of our dual G9a histone-methyltransferase and DNA-methyltransferase 1 (DNMT1) inhibitor in human HCC cells and their crosstalk with fibrogenic cells. The expression of G9a and DNMT1, along with that of their molecular adaptor ubiquitin-like with PHD and RING finger domains-1 (UHRF1), was measured in human HCCs (n = 268), peritumoral tissues (n = 154), and HCC cell lines (n = 32). We evaluated the effect of individual and combined inhibition of G9a and DNMT1 on HCC cell growth by pharmacological and genetic approaches. The activity of our lead compound, CM-272, was examined in HCC cells under normoxia and hypoxia, human hepatic stellate cells and LX2 cells, and xenograft tumors formed by HCC or combined HCC+LX2 cells. We found a significant and correlative overexpression of G9a, DNMT1, and UHRF1 in HCCs in association with poor prognosis. Independent G9a and DNMT1 pharmacological targeting synergistically inhibited HCC cell growth. CM-272 potently reduced HCC and LX2 cells proliferation and quelled tumor growth, particularly in HCC+LX2 xenografts. Mechanistically, CM-272 inhibited the metabolic adaptation of HCC cells to hypoxia and induced a differentiated phenotype in HCC and fibrogenic cells. The expression of the metabolic tumor suppressor gene fructose-1,6-bisphosphatase (FBP1), epigenetically repressed in HCC, was restored by CM-272. Conclusion: Combined targeting of G9a/DNMT1 with compounds such as CM-272 is a promising strategy for HCC treatment. Our findings also underscore the potential of differentiation therapy in HCC.

Extracellular vesicles in DLBCL provide abundant clues to aberrant transcriptional programming and genomic alterations.

The biological role of extracellular vesicles (EVs) in diffuse large B-cell lymphoma (DLBCL) initiation and progression remains largely unknown. We characterized EVs secreted by 5 DLBCL cell lines, a primary DLBCL tumor, and a normal control B-cell sample, optimized their purification, and analyzed their content. We found that DLBCLs secreted large quantities of CD63, Alix, TSG101, and CD81 EVs, which can be extracted using an ultracentrifugation-based method and traced by their cell of origin surface markers. We also showed that tumor-derived EVs can be exchanged between lymphoma cells, normal tonsillar cells, and HK stromal cells. We then examined the content of EVs, focusing on isolation of high-quality total RNA. We sequenced the total RNA and analyzed the nature of RNA species, including coding and noncoding RNAs. We compared whole-cell and EV-derived RNA composition in benign and malignant B cells and discovered that transcripts from EVs were involved in many critical cellular functions. Finally, we performed mutational analysis and found that mutations detected in EVs exquisitely represented mutations in the cell of origin. These results enhance our understanding and enable future studies of the role that EVs may play in the pathogenesis of DLBCL, particularly with regards to the exchange of genomic information. Current findings open a new strategy for liquid biopsy approaches in disease monitoring.

The hydroxymethylome of multiple myeloma identifies FAM72D as a 1q21 marker linked to proliferation.

Cell identity relies on the cross-talk between genetics and epigenetics and their impact on gene expression. Oxidation of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) is the first step of an active DNA demethylation process occurring mainly at enhancers and gene bodies and, as such, participates in processes governing cell identity in normal and pathological conditions. Although genetic alterations are well documented in multiple myeloma (MM), epigenetic alterations associated with this disease have not yet been thoroughly analyzed. To gain insight into the biology of MM, genome-wide 5hmC profiles were obtained and showed that regions enriched in this modified base overlap with MM enhancers and super enhancers and are close to highly expressed genes. Through the definition of a MM-specific 5hmC signature, we identified FAM72D as a poor prognostic gene located on 1q21, a region amplified in high risk myeloma. We further uncovered that FAM72D functions as part of the FOXM1 transcription factor network controlling cell proliferation and survival and we evidenced an increased sensitivity of cells expressing high levels of FOXM1 and FAM72 to epigenetic drugs targeting histone deacetylases and DNA methyltransferases.

Epigenomic profiling of myelofibrosis reveals widespread DNA methylation changes in enhancer elements and ZFP36L1 as a potential tumor suppressor gene that is epigenetically regulated.

In this study we interrogated the DNA methylome of myelofibrosis patients using high-density DNA methylation arrays. We detected 35,215 differentially methylated CpG, corresponding to 10,253 genes, between myelofibrosis patients and healthy controls. These changes were present both in primary and secondary myelofibrosis, which showed no differences between them. Remarkably, most differentially methylated CpG were located outside gene promoter regions and showed significant association with enhancer regions. This aberrant enhancer hypermethylation was negatively correlated with the expression of 27 genes in the myelofibrosis cohort. Of these, we focused on the ZFP36L1 gene and validated its decreased expression and enhancer DNA hypermethylation in an independent cohort of patients and myeloid cell-lines. In vitro reporter assay and 5'-azacitidine treatment confirmed the functional relevance of hyper-methylation of ZFP36L1 enhancer. Furthermore, in vitro rescue of ZFP36L1 expression had an impact on cell proliferation and induced apoptosis in SET-2 cell line indicating a possible role of ZFP36L1 as a tumor suppressor gene in myelofibrosis. Collectively, we describe the DNA methylation profile of myelofibrosis, identifying extensive changes in enhancer elements and revealing ZFP36L1 as a novel candidate tumor suppressor gene.

Whole-epigenome analysis in multiple myeloma reveals DNA hypermethylation of B cell-specific enhancers.

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

A DNA methylation fingerprint of 1628 human samples.

Most of the studies characterizing DNA methylation patterns have been restricted to particular genomic loci in a limited number of human samples and pathological conditions. Herein, we present a compromise between an extremely comprehensive study of a human sample population with an intermediate level of resolution of CpGs at the genomic level. We obtained a DNA methylation fingerprint of 1628 human samples in which we interrogated 1505 CpG sites. The DNA methylation patterns revealed show this epigenetic mark to be critical in tissue-type definition and stemness, particularly around transcription start sites that are not within a CpG island. For disease, the generated DNA methylation fingerprints show that, during tumorigenesis, human cancer cells underwent a progressive gain of promoter CpG-island hypermethylation and a loss of CpG methylation in non-CpG-island promoters. Although transformed cells are those in which DNA methylation disruption is more obvious, we observed that other common human diseases, such as neurological and autoimmune disorders, had their own distinct DNA methylation profiles. Most importantly, we provide proof of principle that the DNA methylation fingerprints obtained might be useful for translational purposes by showing that we are able to identify the tumor type origin of cancers of unknown primary origin (CUPs). Thus, the DNA methylation patterns identified across the largest spectrum of samples, tissues, and diseases reported to date constitute a baseline for developing higher-resolution DNA methylation maps and provide important clues concerning the contribution of CpG methylation to tissue identity and its changes in the most prevalent human diseases.

Unraveling a novel transcription factor code determining the human arterial-specific endothelial cell signature.

Endothelial cells (ECs) lining arteries and veins have distinct molecular/functional signatures. The underlying regulatory mechanisms are incompletely understood. Here, we established a specific fingerprint of freshly isolated arterial and venous ECs from human umbilical cord comprising 64 arterial and 12 venous genes, representing distinct functions/pathways. Among the arterial genes were 8 transcription factors (TFs), including Notch target HEY2, the current "gold standard" determinant for arterial EC (aEC) specification. Culture abrogated differential gene expression in part due to gradual loss of canonical Notch activity and HEY2 expression. Notably, restoring HEY2 expression or Delta-like4-induced Notch signaling in cultured ECs only partially reinstated the aEC gene signature, whereas combined overexpression of the 8 TFs restored this fingerprint more robustly. Whereas some TFs stimulated few genes, others boosted a large proportion of arterial genes. Although there was some overlap and cross-regulation, the TFs largely complemented each other in regulating the aEC gene profile. Finally, overexpression of the 8 TFs in human umbilical vein ECs conveyed an arterial-like behavior upon their implantation in a Matrigel plug in vivo. Thus, our study shows that Notch signaling determines only part of the aEC signature and identifies additional novel and complementary transcriptional players in the complex regulation of human arteriovenous EC identity.

Expression of MALT1 oncogene in hematopoietic stem/progenitor cells recapitulates the pathogenesis of human lymphoma in mice.

Chromosomal translocations involving the MALT1 gene are hallmarks of mucosa-associated lymphoid tissue (MALT) lymphoma. To date, targeting these translocations to mouse B cells has failed to reproduce human disease. Here, we induced MALT1 expression in mouse Sca1(+)Lin(-) hematopoietic stem/progenitor cells, which showed NF-κB activation and early lymphoid priming, being selectively skewed toward B-cell differentiation. These cells accumulated in extranodal tissues and gave rise to clonal tumors recapitulating the principal clinical, biological, and molecular genetic features of MALT lymphoma. Deletion of p53 gene accelerated tumor onset and induced transformation of MALT lymphoma to activated B-cell diffuse large-cell lymphoma (ABC-DLBCL). Treatment of MALT1-induced lymphomas with a specific inhibitor of MALT1 proteolytic activity decreased cell viability, indicating that endogenous Malt1 signaling was required for tumor cell survival. Our study shows that human-like lymphomas can be modeled in mice by targeting MALT1 expression to hematopoietic stem/progenitor cells, demonstrating the oncogenic role of MALT1 in lymphomagenesis. Furthermore, this work establishes a molecular link between MALT lymphoma and ABC-DLBCL, and provides mouse models to test MALT1 inhibitors. Finally, our results suggest that hematopoietic stem/progenitor cells may be involved in the pathogenesis of human mature B-cell lymphomas.

Down-regulated expression of hsa-miR-181c in Fanconi anemia patients: implications in TNFα regulation and proliferation of hematopoietic progenitor cells.

Fanconi anemia (FA) is an inherited genetic disorder associated with BM failure and cancer predisposition. In the present study, we sought to elucidate the role of microRNAs (miRNAs) in the hematopoietic defects observed in FA patients. Initial studies showed that 3 miRNAs, hsa-miR-133a, hsa-miR-135b, and hsa-miR-181c, were significantly down-regulated in lymphoblastoid cell lines and fresh peripheral blood cells from FA patients. In vitro studies with cells expressing the luciferase reporter fused to the TNFα 3'-untranslated region confirmed in silico predictions suggesting an interaction between hsa-miR-181c and TNFα mRNA. These observations were consistent with the down-regulated expression of TNFα mediated by hsa-miR-181c in cells from healthy donors and cells from FA patients. Because of the relevance of TNFα in the hematopoietic defects of FA patients, in the present study, we transfected BM cells from FA patients with hsa-miR-181c to evaluate the impact of this miRNA on their clonogenic potential. hsa-miR-181c markedly increased the number and size of the myeloid and erythroid colonies generated by BM cells from FA patients. Our results offer new clues toward understanding the biologic basis of BM failure in FA patients and open new possibilities for the treatment of the hematologic dysfunction in FA patients based on miRNA regulation.

Vav3 collaborates with p190-BCR-ABL in lymphoid progenitor leukemogenesis, proliferation, and survival.

Despite the introduction of tyrosine kinase inhibitor therapy, the prognosis for p190-BCR-ABL(+) acute lymphoblastic leukemia remains poor. In the present study, we present the cellular and molecular roles of the Rho GTPase guanine nucleotide exchange factor Vav in lymphoid leukemogenesis and explore the roles of Vav proteins in BCR-ABL-dependent signaling. We show that genetic deficiency of the guanine nucleotide exchange factor Vav3 delays leukemogenesis by p190-BCR-ABL and phenocopies the effect of Rac2 deficiency, a downstream effector of Vav3. Compensatory up-regulation of expression and activation of Vav3 in Vav1/Vav2-deficient B-cell progenitors increases the transformation ability of p190-BCR-ABL. Vav3 deficiency induces apoptosis of murine and human leukemic lymphoid progenitors, decreases the activation of Rho GTPase family members and p21-activated kinase, and is associated with increased Bad phosphorylation and up-regulation of Bax, Bak, and Bik. Finally, Vav3 activation only partly depends on ABL TK activity, and Vav3 deficiency collaborates with tyrosine kinase inhibitors to inhibit CrkL activation and impair leukemogenesis in vitro and in vivo. We conclude that Vav3 represents a novel specific molecular leukemic effector for multitarget therapy in p190-BCR-ABL-expressing acute lymphoblastic leukemia.

HDACi vorinostat protects muscle from degeneration after acute rotator cuff injury in mice.

Aims: Rotator cuff (RC) injuries are characterized by tendon rupture, muscle atrophy, retraction, and fatty infiltration, which increase injury severity and jeopardize adequate tendon repair. Epigenetic drugs, such as histone deacetylase inhibitors (HDACis), possess the capacity to redefine the molecular signature of cells, and they may have the potential to inhibit the transformation of the fibro-adipogenic progenitors (FAPs) within the skeletal muscle into adipocyte-like cells, concurrently enhancing the myogenic potential of the satellite cells. Methods: HDACis were added to FAPs and satellite cell cultures isolated from mice. The HDACi vorinostat was additionally administered into a RC injury animal model. Histological analysis was carried out on the isolated supra- and infraspinatus muscles to assess vorinostat anti-muscle degeneration potential. Results: Vorinostat, a HDACi compound, blocked the adipogenic transformation of muscle-associated FAPs in culture, promoting myogenic progression of the satellite cells. Furthermore, it protected muscle from degeneration after acute RC in mice in the earlier muscle degenerative stage after tenotomy. Conclusion: The HDACi vorinostat may be a candidate to prevent early muscular degeneration after RC injury.