Search

Palatine tonsils are secondary lymphoid organs representing the first line of immunological defense against inhaled or ingested pathogens. Here, we present a comprehensive census of cell types forming the human tonsil by applying single-cell transcriptome, epigenome, proteome and adaptive immune repertoire sequencing as well as spatial transcriptomics, resulting in an atlas of >357,000 cells. We provide a glossary of 121 annotated cell types and states, and disentangle gene regulatory mechanisms that drive cells through specialized lineage trajectories. Exemplarily, we stratify multiple tonsil-resident myeloid slancyte subtypes, establish a distant BCL6 superenhancer as locally active in both follicle-associated T and B cells, and describe SIX5 as a potentially novel transcriptional regulator of plasma cell maturation. Further, our atlas is a reference map to understand alterations observed in disease. Here, we discover immune-phenotype plasticity in tumoral cells and microenvironment shifts of mantle cell lymphomas (MCL). To facilitate such reference-based analysis, we develop HCATonsilData and SLOcatoR, a computational framework that provides programmatic and modular access to our dataset; and allows the straightforward annotation of future single-cell profiles from secondary lymphoid organs.

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF. This work was supported by the European Research Council under the 7th Framework Program FP7/2007-2013 (ERC grant agreement no. 609989 to M.A.M.-R.) and the Spanish Ministerio de Ciencia, Innovación y Universidades through nos. IJCI-2015-23352 to I.F. and BFU2017-85926-P and PID2020-115696RB-I00 to M.A.M.-R. CRG acknowledges support from ‘Centro de Excelencia Severo Ochoa 2013-2017’, SEV-2012-0208 and the CERCA Program/Generalitat de Catalunya, as well as support from the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III and the EMBL partnership, the Generalitat de Catalunya through Departament de Salut and Departament d’Empresa i Coneixement, and cofinancing with funds from the European Regional Development Fund by the Spanish Ministry of Science and Innovation corresponding to the Programa Opertaivo FEDER Plurirregional de España 2014–2020 and by the Secretaria d’Universitats i Recerca, Departament d’Empresa i Coneixement of the Generalitat de Catalunya corresponding to the program Operatiu FEDER Catalunya 2014–2020 and the NIH (to C.T. Wu no. R01HD091797 for supporting I.F.)

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., The dynamics of genome architecture during human cell differentiation and upon neoplastic transformation remain poorly characterized. Here, the authors integrate in situ Hi-C and nine additional omic layers to characterize the dynamic changes in 3D genome architecture during normal B cell differentiation and in neoplastic cells from chronic lymphocytic leukemia and mantle cell lymphoma patients.

Mètodes computacionals; Genòmica Métodos computacionales; Genómica Computational Methods; Genomics The rapid development of Chromosome Conformation Capture (3C-based techniques), as well as imaging together with bioinformatics analyses, has been fundamental for unveiling that chromosomes are organized into the so-called topologically associating domains or TADs. While TADs appear as nested patterns in the 3C-based interaction matrices, the vast majority of available TAD callers are based on the hypothesis that TADs are individual and unrelated chromatin structures. Here we introduce TADpole, a computational tool designed to identify and analyze the entire hierarchy of TADs in intra-chromosomal interaction matrices. TADpole combines principal component analysis and constrained hierarchical clustering to provide a set of significant hierarchical chromatin levels in a genomic region of interest. TADpole is robust to data resolution, normalization strategy and sequencing depth. Domain borders defined by TADpole are enriched in main architectural proteins (CTCF and cohesin complex subunits) and in the histone mark H3K4me3, while their domain bodies, depending on their activation-state, are enriched in either H3K36me3 or H3K27me3, highlighting that TADpole is able to distinguish functional TAD units. Additionally, we demonstrate that TADpole's hierarchical annotation, together with the new DiffT score, allows for detecting significant topological differences on Capture Hi-C maps between wild-type and genetically engineered mouse. European Research Council under the Seventh Framework Program FP7/2007-2013 [609989, in part]; European Union's Horizon 2020 Research and Innovation Programme [676556]; Spanish Ministry of Science and Innovation [BFU2013-47736-P, BFU2017-85926-P to M.A.M-R., IJCI-2015-23352 to I.F., BES-2014-070327 to P.S-V.]; ‘Centro de Excelencia Severo Ochoa 2013–2017’, SEV-2012-0208; CERCA Programme/Generalitat de Catalunya (to C.R.G.). Funding for open access charge: European Research Council under the Seventh Framework Program FP7/2007-2013 [609989]. We also acknowledge the support of the Spanish Ministry of Science and Innovation to the EMBL partnership, the ‘Centro de Excelencia Severo Ochoa 2013-2017’, SEV-2012-0208, the CERCA Programme/Generalitat de Catalunya, Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III, the Generalitat de Catalunya through Departament de Salut and Departament d’Empresa i Coneixement and the Co-financing by the Spanish Ministry of Science and Innovation with funds from the European Regional Development Fund (ERDF) corresponding to the 2014-2020 Smart Growth Operating Program to the CRG.

Despite recent advances, the dynamics of genome architecture and chromatin function during human cell differentiation and its potential reorganization upon neoplastic transformation remains poorly characterized. Here, we integrate in situ Hi-C and nine additional omic layers to define and biologically characterize the dynamic changes in three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia (CLL) and mantle cell lymphoma (MCL) patients. Beyond conventional active (A) and inactive (B) compartments, an integrative analysis of Hi-C data reveals the presence of a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, we detect that 28% of the compartments change at defined maturation stages and mostly involve the intermediate compartment. The transition from naive to germinal center B cells is associated with widespread chromatin activation, which mostly reverts into the naive state upon further maturation of germinal center cells into memory B cells. The analysis of CLL and MCL neoplastic cells points both to entity and subtype-specific alterations in chromosome organization. Remarkably, we observe that large chromatin blocks containing key disease-specific genes alter their 3D genome organization. These include the inactivation of a 2Mb region containing the EBF1 gene in CLL and the activation of a 6.1Mb region containing the SOX11 gene in clinically aggressive MCL. 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.

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. This research was funded by the European Union's Seventh Framework Programme through the Blueprint Consortium (grant agreement 282510), Fundació La Marató de TV3, Instituto de Salud Carlos III (ISCIII) PI14/01867, PI16/02024, and PI17/00701, TRASCAN (EPICA), MINECO Explora (RTHALMY), Departamento de Salud del Gobierno de Navarra 40/2016, Gilead Fellowship Program (GLD16/00142), Multiple Myeloma Research Foundation Networks of excellence, the International Myeloma Foundation (Brian van Novis) and the Qatar National Research Fund award 7-916-3-237. Furthermore, the authors acknowledge the support of the Generalitat de Catalunya Suport Grups de Recerca AGAUR 2017-SGR-736 and 2017-SGR-1142, TRASCAN-iMMunocelland European Research Council starting grant (MYELOMANEXT), CIBERONC (CB16/12/00489, CB16/12/00369 and CB16/12/00225), co-financed with FEDER funds, the Accelerator award CRUK/AIRC/AECC joint funder-partnership, as well as NCI R01 CA180475 and a MMRF collaborative grant. R.O. was supported by a FPU Fellowship of the Spanish Government, M.K. by an AOI grant of the Spanish Association Against Cancer and N.R. by the Acció instrumental d'incorporació de scientífics i tecnòlegs PERIS 2016 from the Generalitat de Catalunya. This work was partially developed at the Centro Esther Koplowitz (CEK, Barcelona, Spain). We particularly acknowledge the patients for their participation and the Biobank of the University of Navarra for its collaboration.

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method–integrative modeling of genomic regions (IMGR)–to increase the genomic resolution of our traces to 10 kb., Author summary Questions regarding the impact of chromosome structure on genome function are focusing increasingly on the manner in which chromosomes are organized within the nucleus. In fact, studies of processes as diverse as gene activation and repression as well as genome repair and stability are all querying how the 3D organization of chromosomal DNA may be a major player. Here, we apply our strategy for tracing chromosomes at super-resolution, traversing over 8 megabases of human chromosome 19 while visualizing genomic features ranging in size from kilobases to megabases. This technology has enabled exploration of the physical nature of a genomic feature called the compartment; compartments are widely hypothesized to reflect the partitioning of genomes into relatively more and less active regions. Excitingly, we find that compartments are, indeed, physical structures, that they are sometimes distinct and other times entangled. We also find that the maternally-derived and the paternally-derived homologous regions can differ more than would be expected by chance. Finally, by integrating image data with information regarding the frequency with which genomic segments contact each other, we produce a 3D model of how the 8 megabases we have imaged may be organized within a single nucleus, achieving 10 kb genomic resolution.