UBR7 acts as a histone chaperone for post-nucleosomal histone H3.

Histone chaperones modulate the stability of histones beginning from histone synthesis, through incorporation into DNA, and during recycling during transcription and replication. Following histone removal from DNA, chaperones regulate histone storage and degradation. Here, we demonstrate that UBR7 is a histone H3.1 chaperone that modulates the supply of pre-existing post-nucleosomal histone complexes. We demonstrate that UBR7 binds to post-nucleosomal H3K4me3 and H3K9me3 histones via its UBR box and PHD. UBR7 binds to the non-nucleosomal histone chaperone NASP. In the absence of UBR7, the pool of NASP-bound post-nucleosomal histones accumulate and chromatin is depleted of H3K4me3-modified histones. We propose that the interaction of UBR7 with NASP and histones opposes the histone storage functions of NASP and that UBR7 promotes reincorporation of post-nucleosomal H3 complexes.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease.

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSß0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with pain-associated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the ß-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

Chromatin Hyperacetylation Impacts Chromosome Folding by Forming a Nuclear Subcompartment.

Delineating how chromosomes fold at length scales beyond one megabase remains obscure relative to smaller-scale folding into TADs, loops, and nucleosomes. We find that rather than simply unfolding chromatin, histone hyperacetylation results in interactions between distant genomic loci separated by tens to hundreds of megabases, even in the absence of transcription. These hyperacetylated "megadomains" are formed by the BRD4-NUT fusion oncoprotein, interact both within and between chromosomes, and form a specific nuclear subcompartment that has elevated gene activity with respect to other subcompartments. Pharmacological degradation of BRD4-NUT results in collapse of megadomains and attenuation of the interactions between them. In contrast, these interactions persist and contacts between newly acetylated regions are formed after inhibiting RNA polymerase II initiation. Our structure-function approach thus reveals that broad chromatin domains of identical biochemical composition, independent of transcription, form nuclear subcompartments, and also indicates the potential of altering chromosome structure for treating human disease.

Latent cellular analysis robustly reveals subtle diversity in large-scale single-cell RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for characterizing the cell-to-cell variation and cellular dynamics in populations which appear homogeneous otherwise in basic and translational biological research. However, significant challenges arise in the analysis of scRNA-seq data, including the low signal-to-noise ratio with high data sparsity, potential batch effects, scalability problems when hundreds of thousands of cells are to be analyzed among others. The inherent complexities of scRNA-seq data and dynamic nature of cellular processes lead to suboptimal performance of many currently available algorithms, even for basic tasks such as identifying biologically meaningful heterogeneous subpopulations. In this study, we developed the Latent Cellular Analysis (LCA), a machine learning-based analytical pipeline that combines cosine-similarity measurement by latent cellular states with a graph-based clustering algorithm. LCA provides heuristic solutions for population number inference, dimension reduction, feature selection, and control of technical variations without explicit gene filtering. We show that LCA is robust, accurate, and powerful by comparison with multiple state-of-the-art computational methods when applied to large-scale real and simulated scRNA-seq data. Importantly, the ability of LCA to learn from representative subsets of the data provides scalability, thereby addressing a significant challenge posed by growing sample sizes in scRNA-seq data analysis.

Resolving medulloblastoma cellular architecture by single-cell genomics.

Medulloblastoma is a malignant childhood cerebellar tumour type that comprises distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. Here we used single-cell transcriptomics to investigate intra- and intertumoral heterogeneity in 25 medulloblastomas spanning all molecular subgroups. WNT, SHH and Group 3 tumours comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumours consisted exclusively of differentiated neuronal-like neoplastic cells. SHH tumours closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumours exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, the relative proportions of which distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

Correction to: H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

The original article can be found online.

H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

Histone H3 K27M mutation is the defining molecular feature of the devastating pediatric brain tumor, diffuse intrinsic pontine glioma (DIPG). The prevalence of histone H3 K27M mutations indicates a critical role in DIPGs, but the contribution of the mutation to disease pathogenesis remains unclear. We show that knockdown of this mutation in DIPG xenografts restores K27M-dependent loss of H3K27me3 and delays tumor growth. Comparisons of matched DIPG xenografts with and without K27M knockdown allowed identification of mutation-specific effects on the transcriptome and epigenome. The resulting transcriptional changes recapitulate expression signatures from K27M primary DIPG tumors and are strongly enriched for genes associated with nervous system development. Integrated analysis of ChIP-seq and expression data showed that genes upregulated by the mutation are overrepresented in apparently bivalent promoters. Many of these targets are associated with more immature differentiation states. Expression profiles indicate K27M knockdown decreases proliferation and increases differentiation within lineages represented in DIPG. These data suggest that K27M-mediated loss of H3K27me3 directly regulates a subset of genes by releasing poised promoters, and contributes to tumor phenotype and growth by limiting differentiation. The delayed tumor growth associated with knockdown of H3 K27M provides evidence that this highly recurrent mutation is a relevant therapeutic target.

Metabolic heterogeneity underlies reciprocal fates of TH17 cell stemness and plasticity.

A defining feature of adaptive immunity is the development of long-lived memory T cells to curtail infection. Recent studies have identified a unique stem-like T-cell subset amongst exhausted CD8-positive T cells in chronic infection1-3, but it remains unclear whether CD4-positive T-cell subsets with similar features exist in chronic inflammatory conditions. Amongst helper T cells, TH17 cells have prominent roles in autoimmunity and tissue inflammation and are characterized by inherent plasticity4-7, although how such plasticity is regulated is poorly understood. Here we demonstrate that TH17 cells in a mouse model of autoimmune disease are functionally and metabolically heterogeneous; they contain a subset with stemness-associated features but lower anabolic metabolism, and a reciprocal subset with higher metabolic activity that supports transdifferentiation into TH1-like cells. These two TH17-cell subsets are defined by selective expression of the transcription factors TCF-1 and T-bet, and by discrete levels of CD27 expression. We also identify signalling via the kinase complex mTORC1 as a central regulator of TH17-cell fate decisions by coordinating metabolic and transcriptional programmes. TH17 cells with disrupted mTORC1 signalling or anabolic metabolism fail to induce autoimmune neuroinflammation or to develop into TH1-like cells, but instead upregulate TCF-1 expression and acquire stemness-associated features. Single-cell RNA sequencing and experimental validation reveal heterogeneity in fate-mapped TH17 cells, and a developmental arrest in the TH1 transdifferentiation trajectory upon loss of mTORC1 activity or metabolic perturbation. Our results establish that the dichotomy of stemness and effector function underlies the heterogeneous TH17 responses and autoimmune pathogenesis, and point to previously unappreciated metabolic control of plasticity in helper T cells.

A Single-Cell Transcriptional Atlas of the Developing Murine Cerebellum.

The cerebellum develops from a restricted number of cell types that precisely organize to form the circuitry that controls sensory-motor coordination and some higher-order cognitive processes. To acquire an enhanced understanding of the molecular processes that mediate cerebellar development, we performed single-cell RNA-sequencing of 39,245 murine cerebellar cells at twelve critical developmental time points. Using recognized lineage markers, we confirmed that the single-cell data accurately recapitulate cerebellar development. We then followed distinct populations from emergence through migration and differentiation, and determined the associated transcriptional cascades. After identifying key lineage commitment decisions, focused analyses uncovered waves of transcription factor expression at those branching points. Finally, we created Cell Seek, a flexible online interface that facilitates exploration of the dataset. Our study provides a transcriptional summarization of cerebellar development at single-cell resolution that will serve as a valuable resource for future investigations of cerebellar development, neurobiology, and disease.

High-resolution transcriptional dissection of in vivo Atoh1-mediated hair cell conversion in mature cochleae identifies Isl1 as a co-reprogramming factor.

In vivo direct conversion of differentiated cells holds promise for regenerative medicine; however, improving the conversion efficiency and producing functional target cells remain challenging. Ectopic Atoh1 expression in non-sensory supporting cells (SCs) in mouse cochleae induces their partial conversion to hair cells (HCs) at low efficiency. Here, we performed single-cell RNA sequencing of whole mouse sensory epithelia harvested at multiple time points after conditional overexpression of Atoh1. Pseudotemporal ordering revealed that converted HCs (cHCs) are present along a conversion continuum that correlates with both endogenous and exogenous Atoh1 expression. Bulk sequencing of isolated cell populations and single-cell qPCR confirmed 51 transcription factors, including Isl1, are differentially expressed among cHCs, SCs and HCs. In transgenic mice, co-overexpression of Atoh1 and Isl1 enhanced the HC conversion efficiency. Together, our study shows how high-resolution transcriptional profiling of direct cell conversion can identify co-reprogramming factors required for efficient conversion.

Metabolic signaling directs the reciprocal lineage decisions of αß and γδ T cells.

The interaction between extrinsic factors and intrinsic signal strength governs thymocyte development, but the mechanisms linking them remain elusive. We report that mechanistic target of rapamycin complex 1 (mTORC1) couples microenvironmental cues with metabolic programs to orchestrate the reciprocal development of two fundamentally distinct T cell lineages, the αß and γδ T cells. Developing thymocytes dynamically engage metabolic programs including glycolysis and oxidative phosphorylation, as well as mTORC1 signaling. Loss of RAPTOR-mediated mTORC1 activity impairs the development of αß T cells but promotes γδ T cell generation, associated with disrupted metabolic remodeling of oxidative and glycolytic metabolism. Mechanistically, we identify mTORC1-dependent control of reactive oxygen species production as a key metabolic signal in mediating αß and γδ T cell development, and perturbation of redox homeostasis impinges upon thymocyte fate decisions and mTORC1-associated phenotypes. Furthermore, single-cell RNA sequencing and genetic dissection reveal that mTORC1 links developmental signals from T cell receptors and NOTCH to coordinate metabolic activity and signal strength. Our results establish mTORC1-driven metabolic signaling as a decisive factor for reciprocal αß and γδ T cell development and provide insight into metabolic control of cell signaling and fate decisions.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.