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, but there are no studies investigating the differentiation stage of tumor PCs in patients with AL and MM. We sought to define a transcriptional atlas of normal PC development (n=11) in secondary lymphoid organs (SLO), peripheral blood (PB) and BM for comparison with the transcriptional programs (TPs) of tumor PCs in AL (n=37), MM (n=46) and MGUS (n=6). Based on bulk and single-cell RNAseq, we observed thirteen TPs during transition of normal PCs throughout SLO, PB and BM; that CD39 outperforms CD19 to discriminate new-born from long-lived BM-PCs; that tumor PCs expressed the most advantageous TPs of normal PC differentiation; that AL shares greater similarity to SLO-PCs whereas MM is transcriptionally closer to PB-PCs and new-born BM-PCs; that AL and MM patients enriched in immature TPs had inferior survival; and that TPs related with protein N-linked glycosylation 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.

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.

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.

BACKGROUND: Cardiac fibroblasts (CFs) have a central role in the ventricular remodeling process associated with different types of fibrosis. Recent studies have shown that fibroblasts do not respond homogeneously to heart injury. Because of the limited set of bona fide fibroblast markers, a proper characterization of fibroblast population heterogeneity in response to cardiac damage is lacking. The purpose of this study was to define CF heterogeneity during ventricular remodeling and the underlying mechanisms that regulate CF function. METHODS: Collagen1 alpha 1-GFP (green fluorescent protein)-positive CFs were characterized after myocardial infarction (MI) by single-cell and bulk RNA sequencing, assay for transposase-accessible chromatin sequencing, and functional assays. Swine and patient samples were studied using bulk RNA sequencing. RESULTS: We identified and characterized a unique CF subpopulation that emerges after MI in mice. These activated fibroblasts exhibit a clear profibrotic signature, express high levels of Cthrc1 (collagen triple helix repeat containing 1), and localize into the scar. Noncanonical transforming growth factor-beta signaling and different transcription factors including SOX9 are important regulators mediating their response to cardiac injury. Absence of CTHRC1 results in pronounced lethality attributable to ventricular rupture. A population of CFs with a similar transcriptome was identified in a swine model of MI and in heart tissue from patients with MI and dilated cardiomyopathy. CONCLUSIONS: We report CF heterogeneity and their dynamics during the course of MI and redefine the CFs that respond to cardiac injury and participate in myocardial remodeling. Our study identifies CTHRC1 as a novel regulator of the healing scar process and a target for future translational studies.

BACKGROUND: Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide. FINDINGS: We performed H3K27me3 chromatin immunoprecipitation followed by high-throughput sequencing and transcriptomic analysis by 3'-end RNA sequencing from B. rapa leaves and inflorescences. To analyze these data we developed a Reproducible Epigenomic Analysis pipeline using Galaxy and Jupyter, packaged into Docker images to facilitate transparency and reuse. We found that H3K27me3 covers roughly one-third of all B. rapa protein-coding genes and its presence correlates with low transcript levels. The comparative analysis between leaves and inflorescences suggested that the expression of various floral regulatory genes during development depends on H3K27me3. To demonstrate the importance of H3K27me3 for B. rapa development, we characterized a mutant line deficient in the H3K27 methyltransferase activity. We found that braA.clf mutant plants presented pleiotropic alterations, e.g., curly leaves due to increased expression and reduced H3K27me3 levels at AGAMOUS-like loci. CONCLUSIONS: We characterized the epigenetic mark H3K27me3 at genome-wide levels and provide genetic evidence for its relevance in B. rapa development. Our work reveals the epigenomic landscape of H3K27me3 in B. rapa and provides novel genomics datasets and bioinformatics analytical resources. We anticipate that this work will lead the way to further epigenomic studies in the complex genome of Brassica crops.