Quiescence Modulates Stem Cell Maintenance and Regenerative Capacity in the Aging Brain.

The function of somatic stem cells declines with age. Understanding the molecular underpinnings of this decline is key to counteract age-related disease. Here, we report a dramatic drop in the neural stem cells (NSCs) number in the aging murine brain. We find that this smaller stem cell reservoir is protected from full depletion by an increase in quiescence that makes old NSCs more resistant to regenerate the injured brain. Once activated, however, young and old NSCs show similar proliferation and differentiation capacity. Single-cell transcriptomics of NSCs indicate that aging changes NSCs minimally. In the aging brain, niche-derived inflammatory signals and the Wnt antagonist sFRP5 induce quiescence. Indeed, intervention to neutralize them increases activation of old NSCs during homeostasis and following injury. Our study identifies quiescence as a key feature of old NSCs imposed by the niche and uncovers ways to activate NSCs to repair the aging brain.

Genome-wide DNA-methylation landscape defines specialization of regulatory T cells in tissues.

Regulatory T cells (Treg cells) perform two distinct functions: they maintain self-tolerance, and they support organ homeostasis by differentiating into specialized tissue Treg cells. We found that epigenetic modifications defined the molecular characteristics of tissue Treg cells. Tagmentation-based whole-genome bisulfite sequencing revealed more than 11,000 regions that were methylated differentially in pairwise comparisons of tissue Treg cell populations and lymphoid T cells. Similarities in the epigenetic landscape led to the identification of a common tissue Treg cell population that was present in many organs and was characterized by gain and loss of DNA methylation that included many gene sites associated with the TH2 subset of helper T cells, such as the gene encoding cytokine IL-33 receptor ST2, as well as the production of tissue-regenerative factors. Furthermore, the ST2-expressing population was dependent on the transcriptional regulator BATF and could be expanded by IL-33. Thus, tissue Treg cells integrate multiple waves of epigenetic reprogramming that define their tissue-restricted specialization.

Corrigendum: Genome-wide DNA-methylation landscape defines specialization of regulatory T cells in tissues.

This corrects the article DOI: 10.1038/ni.3799.

Precursors for Nonlymphoid-Tissue Treg Cells Reside in Secondary Lymphoid Organs and Are Programmed by the Transcription Factor BATF.

Specialized regulatory T (Treg) cells accumulate and perform homeostatic and regenerative functions in nonlymphoid tissues. Whether common precursors for nonlymphoid-tissue Treg cells exist and how they differentiate remain elusive. Using transcription factor nuclear factor, interleukin 3 regulated (Nfil3) reporter mice and single-cell RNA-sequencing (scRNA-seq), we identified two precursor stages of interleukin 33 (IL-33) receptor ST2-expressing nonlymphoid tissue Treg cells, which resided in the spleen and lymph nodes. Global chromatin profiling of nonlymphoid tissue Treg cells and the two precursor stages revealed a stepwise acquisition of chromatin accessibility and reprogramming toward the nonlymphoid-tissue Treg cell phenotype. Mechanistically, we identified and validated the transcription factor Batf as the driver of the molecular tissue program in the precursors. Understanding this tissue development program will help to harness regenerative properties of tissue Treg cells for therapy.

Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating 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.

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.

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.

HMGB2 Loss upon Senescence Entry Disrupts Genomic Organization and Induces CTCF Clustering across Cell Types.

Processes like cellular senescence are characterized by complex events giving rise to heterogeneous cell populations. However, the early molecular events driving this cascade remain elusive. We hypothesized that senescence entry is triggered by an early disruption of the cells' three-dimensional (3D) genome organization. To test this, we combined Hi-C, single-cell and population transcriptomics, imaging, and in silico modeling of three distinct cells types entering senescence. Genes involved in DNA conformation maintenance are suppressed upon senescence entry across all cell types. We show that nuclear depletion of the abundant HMGB2 protein occurs early on the path to senescence and coincides with the dramatic spatial clustering of CTCF. Knocking down HMGB2 suffices for senescence-induced CTCF clustering and for loop reshuffling, while ectopically expressing HMGB2 rescues these effects. Our data suggest that HMGB2-mediated genomic reorganization constitutes a primer for the ensuing senescent program.

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.

scSNPdemux: a sensitive demultiplexing pipeline using single nucleotide polymorphisms for improved pooled single-cell RNA sequencing analysis.

BACKGROUND: Here we present scSNPdemux, a sample demultiplexing pipeline for single-cell RNA sequencing data using natural genetic variations in humans. The pipeline requires alignment files from Cell Ranger (10× Genomics), a population SNP database and genotyped single nucleotide polymorphisms (SNPs) per sample. The tool works on sparse genotyping data in VCF format for sample identification. RESULTS: The pipeline was tested on both single-cell and single-nuclei based RNA sequencing datasets and showed superior demultiplexing performance over the lipid-based CellPlex and Multi-seq sample multiplexing technique which incurs additional single cell library preparation steps. Specifically, our pipeline demonstrated superior sensitivity and specificity in cell-identity assignment over CellPlex, especially on immune cell types with low RNA content. CONCLUSIONS: We designed a streamlined pipeline for single-cell sample demultiplexing, aiming to overcome common problems in multiplexing samples using single cell libraries which might affect data quality and can be costly.

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.

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.

Suv4-20h2 mediates chromatin compaction and is important for cohesin recruitment to heterochromatin.

Cohesin plays an important role in chromatid cohesion and has additional functions in higher-order chromatin organization and in transcriptional regulation. The binding of cohesin to euchromatic regions is largely mediated by CTCF or the mediator complex. However, it is currently unknown how cohesin is recruited to pericentric heterochromatin in mammalian cells. Here we define the histone methyltransferase Suv4-20h2 as a major structural constituent of heterochromatin that mediates chromatin compaction and cohesin recruitment. Suv4-20h2 stably associates with pericentric heterochromatin through synergistic interactions with multiple heterochromatin protein 1 (HP1) molecules, resulting in compaction of heterochromatic regions. Suv4-20h mutant cells display an overall reduced chromatin compaction and an altered chromocenter organization in interphase referred to as "chromocenter scattering." We found that Suv4-20h-deficient cells display chromosome segregation defects during mitosis that coincide with reduced sister chromatid cohesion. Notably, cohesin subunits interact with Suv4-20h2 both in vitro and in vivo. This interaction is necessary for cohesin binding to heterochromatin, as Suv4-20h mutant cells display substantially reduced cohesin levels at pericentric heterochromatin. This defect is most prominent in G0-phase cells, where cohesin is virtually lost from heterochromatin, suggesting that Suv4-20h2 is involved in the initial loading or maintenance of cohesion subunits. In summary, our data provide the first compelling evidence that Suv4-20h2 plays essential roles in regulating nuclear architecture and ensuring proper chromosome segregation.

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.

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks.

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

Retrieving Chromatin Patterns from Deep Sequencing Data Using Correlation Functions.

Epigenetic modifications and other chromatin features partition the genome on multiple length scales. They define chromatin domains with distinct biological functions that come in sizes ranging from single modified DNA bases to several megabases in the case of heterochromatic histone modifications. Due to chromatin folding, domains that are well separated along the linear nucleosome chain can form long-range interactions in three-dimensional space. It has now become a routine task to map epigenetic marks and chromatin structure by deep sequencing methods. However, assessing and comparing the properties of chromatin domains and their positional relationships across data sets without a priori assumptions remains challenging. Here, we introduce multiscale correlation evaluation (MCORE), which uses the fluctuation spectrum of mapped sequencing reads to quantify and compare chromatin patterns over a broad range of length scales in a model-independent manner. We applied MCORE to map the chromatin landscape in mouse embryonic stem cells and differentiated neural cells. We integrated sequencing data from chromatin immunoprecipitation, RNA expression, DNA methylation, and chromosome conformation capture experiments into network models that reflect the positional relationships among these features on different genomic scales. Furthermore, we used MCORE to compare our experimental data to models for heterochromatin reorganization during differentiation. The application of correlation functions to deep sequencing data complements current evaluation schemes and will support the development of quantitative descriptions of chromatin networks.

Nucleosome repositioning links DNA (de)methylation and differential CTCF binding during stem cell development.

During differentiation of embryonic stem cells, chromatin reorganizes to establish cell type-specific expression programs. Here, we have dissected the linkages between DNA methylation (5mC), hydroxymethylation (5hmC), nucleosome repositioning, and binding of the transcription factor CTCF during this process. By integrating MNase-seq and ChIP-seq experiments in mouse embryonic stem cells (ESC) and their differentiated counterparts with biophysical modeling, we found that the interplay between these factors depends on their genomic context. The mostly unmethylated CpG islands have reduced nucleosome occupancy and are enriched in cell type-independent binding sites for CTCF. The few remaining methylated CpG dinucleotides are preferentially associated with nucleosomes. In contrast, outside of CpG islands most CpGs are methylated, and the average methylation density oscillates so that it is highest in the linker region between nucleosomes. Outside CpG islands, binding of TET1, an enzyme that converts 5mC to 5hmC, is associated with labile, MNase-sensitive nucleosomes. Such nucleosomes are poised for eviction in ESCs and become stably bound in differentiated cells where the TET1 and 5hmC levels go down. This process regulates a class of CTCF binding sites outside CpG islands that are occupied by CTCF in ESCs but lose the protein during differentiation. We rationalize this cell type-dependent targeting of CTCF with a quantitative biophysical model of competitive binding with the histone octamer, depending on the TET1, 5hmC, and 5mC state.

Environment-induced epigenetic reprogramming in genomic regulatory elements in smoking mothers and their children.

Epigenetic mechanisms have emerged as links between prenatal environmental exposure and disease risk later in life. Here, we studied epigenetic changes associated with maternal smoking at base pair resolution by mapping DNA methylation, histone modifications, and transcription in expectant mothers and their newborn children. We found extensive global differential methylation and carefully evaluated these changes to separate environment associated from genotype-related DNA methylation changes. Differential methylation is enriched in enhancer elements and targets in particular "commuting" enhancers having multiple, regulatory interactions with distal genes. Longitudinal whole-genome bisulfite sequencing revealed that DNA methylation changes associated with maternal smoking persist over years of life. Particularly in children prenatal environmental exposure leads to chromatin transitions into a hyperactive state. Combined DNA methylation, histone modification, and gene expression analyses indicate that differential methylation in enhancer regions is more often functionally translated than methylation changes in promoters or non-regulatory elements. Finally, we show that epigenetic deregulation of a commuting enhancer targeting c-Jun N-terminal kinase 2 (JNK2) is linked to impaired lung function in early childhood.