Search

The presence of HIV in sequestered reservoirs is a central impediment to a functional cure, allowing HIV to persist despite life-long antiretroviral therapy (ART), and driving a variety of comorbid conditions. Our understanding of the latent HIV reservoir in the central nervous system is incomplete, because of difficulties in accessing human central nervous system tissues. Microglia contribute to HIV reservoirs, but the molecular phenotype of HIV-infected microglia is poorly understood. We leveraged the unique "Last Gift" rapid autopsy program, in which people with HIV are closely followed until days or even hours before death. Microglial populations were heterogeneous regarding their gene expression profiles but showed similar chromatin accessibility landscapes. Despite ART, we detected occasional microglia containing cell-associated HIV RNA and HIV DNA integrated into open regions of the host's genome (∼0.005%). Microglia with detectable HIV RNA showed an inflammatory phenotype. These results demonstrate a distinct myeloid cell reservoir in the brains of people with HIV despite suppressive ART. Strategies for curing HIV and neurocognitive impairment will need to consider the myeloid compartment to be successful.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Human cell line models, including the neuronal precursor line LUHMES, are important for investigating developmental transcriptional dynamics within imprinted regions, particularly the 15q11-q13 Angelman (AS) and Prader-Willi (PWS) syndrome locus. AS results from loss of maternal UBE3A in neurons, where the paternal allele is silenced by a convergent antisense transcript UBE3A-ATS, a lncRNA that terminates at PWAR1 in non-neurons. qRT-PCR analysis confirmed the exclusive and progressive increase in UBE3A-ATS in differentiating LUHMES neurons, validating their use for studying UBE3A silencing. Genome-wide transcriptome analyses revealed changes to 11 834 genes during neuronal differentiation, including the upregulation of most genes within the 15q11-q13 locus. To identify dynamic changes in chromatin loops linked to transcriptional activity, we performed a HiChIP validated by 4C, which identified two neuron-specific CTCF loops between MAGEL2-SNRPN and PWAR1-UBE3A. To determine if allele-specific differentially methylated regions (DMR) may be associated with CTCF loop anchors, whole genome long-read nanopore sequencing was performed. We identified a paternally hypomethylated DMR near the SNRPN upstream loop anchor exclusive to neurons and a paternally hypermethylated DMR near the PWAR1 CTCF anchor exclusive to undifferentiated cells, consistent with increases in neuronal transcription. Additionally, DMRs near CTCF loop anchors were observed in both cell types, indicative of allele-specific differences in chromatin loops regulating imprinted transcription. These results provide an integrated view of the 15q11-q13 epigenetic landscape during LUHMES neuronal differentiation, underscoring the complex interplay of transcription, chromatin looping, and DNA methylation. They also provide insights for future therapeutic approaches for AS and PWS.

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 h period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

In many animal species, the oocyte meiotic spindle, which is required for chromosome segregation, forms without centrosomes. In some systems, Ran-GEF on chromatin initiates spindle assembly. We found that in Caenorhabditis elegans oocytes, endogenously-tagged Ran-GEF dissociates from chromatin during spindle assembly but re-associates during meiotic anaphase. Meiotic spindle assembly occurred after auxin-induced degradation of Ran-GEF, but anaphase I was faster than controls and extrusion of the first polar body frequently failed. In search of a possible alternative pathway for spindle assembly, we found that soluble tubulin concentrates in the nuclear volume during germinal vesicle breakdown. We found that the concentration of soluble tubulin in the metaphase spindle region is enclosed by ER sheets which exclude cytoplasmic organelles including mitochondria and yolk granules. Measurement of the volume occupied by yolk granules and mitochondria indicated that volume exclusion would be sufficient to explain the concentration of tubulin in the spindle volume. We suggest that this concentration of soluble tubulin may be a redundant mechanism promoting spindle assembly near chromosomes.

Long noncoding RNAs (LncRNAs) are key regulators of gene expression and can mediate their effects in both the nucleus and cytoplasm. Some of the best-characterized lncRNAs are localized within the nucleus, where they modulate the nuclear architecture and influence gene expression. In this review, we discuss the role of lncRNAs in nuclear architecture in the context of their gene regulatory functions in innate immunity. Here, we discuss various approaches to functionally characterize nuclear-localized lncRNAs and the challenges faced in the field.

EWS fusion oncoproteins underlie several human malignancies including Desmoplastic Small Round Cell Tumor (DSRCT), an aggressive cancer driven by EWS-WT1 fusion proteins. Here we combine chromatin occupancy and 3D profiles to identify EWS-WT1-dependent gene regulation networks and target genes. We show that EWS-WT1 is a powerful chromatin activator controlling an oncogenic gene expression program that characterizes primary tumors. Similar to wild type WT1, EWS-WT1 has two isoforms that differ in their DNA binding domain and we find that they have distinct DNA binding profiles and are both required to generate viable tumors that resemble primary DSRCT. Finally, we identify candidate EWS-WT1 target genes with potential therapeutic implications, including CCND1, whose inhibition by the clinically-approved drug Palbociclib leads to marked tumor burden decrease in DSRCT PDXs in vivo. Taken together, our studies identify gene regulation programs and therapeutic vulnerabilities in DSRCT and provide a mechanistic understanding of the complex oncogenic activity of EWS-WT1.

Although microglia are macrophages of the central nervous system, their involvement is not limited to immune functions. The roles of microglia during development in humans remain poorly understood due to limited access to fetal tissue. To understand how microglia can impact human neurodevelopment, the methyl-CpG binding protein 2 (MECP2) gene was knocked out in human microglia-like cells (MGLs). Disruption of the MECP2 in MGLs led to transcriptional and functional perturbations, including impaired phagocytosis. The co-culture of healthy MGLs with MECP2-knockout (KO) neurons rescued synaptogenesis defects, suggesting a microglial role in synapse formation. A targeted drug screening identified ADH-503, a CD11b agonist, restored phagocytosis and synapse formation in spheroid-MGL co-cultures, significantly improved disease progression, and increased survival in MeCP2-null mice. These results unveil a MECP2-specific regulation of human microglial phagocytosis and identify a novel therapeutic treatment for MECP2-related conditions.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

BACKGROUND: A number of deep learning models have been developed to predict epigenetic features such as chromatin accessibility from DNA sequence. Model evaluations commonly report performance genome-wide; however, cis regulatory elements (CREs), which play critical roles in gene regulation, make up only a small fraction of the genome. Furthermore, cell type-specific CREs contain a large proportion of complex disease heritability. RESULTS: We evaluate genomic deep learning models in chromatin accessibility regions with varying degrees of cell type specificity. We assess two modeling directions in the field: general purpose models trained across thousands of outputs (cell types and epigenetic marks) and models tailored to specific tissues and tasks. We find that the accuracy of genomic deep learning models, including two state-of-the-art general purpose models-Enformer and Sei-varies across the genome and is reduced in cell type-specific accessible regions. Using accessibility models trained on cell types from specific tissues, we find that increasing model capacity to learn cell type-specific regulatory syntax-through single-task learning or high capacity multi-task models-can improve performance in cell type-specific accessible regions. We also observe that improving reference sequence predictions does not consistently improve variant effect predictions, indicating that novel strategies are needed to improve performance on variants. CONCLUSIONS: Our results provide a new perspective on the performance of genomic deep learning models, showing that performance varies across the genome and is particularly reduced in cell type-specific accessible regions. We also identify strategies to maximize performance in cell type-specific accessible regions.

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we have compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improve the models predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores reveal the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models predictions. Furthermore, we identify non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several new RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT1, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

The chromatin-associated protein WDR5 has been widely studied due to its role in histone modification and its potential as a pharmacological target for the treatment of cancer. In humans, the protein with highest sequence homology to WDR5 is encoded by the retrogene WDR5B, which remains unexplored. Here, we used CRISPR-Cas9 genome editing to generate WDR5B knockout and WDR5B-FLAG knock-in cell lines for further characterization. In contrast to WDR5, WDR5B exhibits low expression in pluripotent cells and is upregulated upon neural differentiation. Loss or shRNA depletion of WDR5B impairs cell growth and increases the fraction of non-viable cells in proliferating retinal pigment epithelial (RPE) cultures. CUT&RUN chromatin profiling in RPE and neural progenitors indicates minimal WDR5B enrichment at established WDR5 binding sites. These results suggest that WDR5 and WDR5B exhibit several divergent biological properties despite sharing a high degree of sequence homology.

During postnatal cardiac hypertrophy, cardiomyocytes undergo mitotic exit, relying on DNA replication-independent mechanisms of histone turnover to maintain chromatin organization and gene transcription. In other tissues, circadian oscillations in nucleosome occupancy influence clock-controlled gene expression, suggesting a role for the circadian clock in temporal control of histone turnover and coordinated cardiomyocyte gene expression. We sought to elucidate roles for the master circadian transcription factor, Bmal1, in histone turnover, chromatin organization, and myocyte-specific gene expression and cell growth in the neonatal period. Bmal1 knockdown in neonatal rat ventricular myocytes decreased myocyte size, total cellular protein synthesis, and transcription of the fetal hypertrophic gene Nppb after treatment with serum or the α-adrenergic agonist phenylephrine. Depletion of Bmal1 decreased the expression of clock-controlled genes Per2 and Tcap, as well as Sik1, a Bmal1 target upregulated in adult versus embryonic hearts. Bmal1 knockdown impaired Per2 and Sik1 promoter accessibility as measured by micrococcal nuclease-quantitative PCR and impaired histone turnover as measured by metabolic labeling of acid-soluble chromatin fractions. Sik1 knockdown in turn decreased myocyte size, while simultaneously inhibiting natriuretic peptide B transcription and activating Per2 transcription. Linking these changes to chromatin remodeling, depletion of the replication-independent histone variant H3.3a inhibited myocyte hypertrophy and prevented phenylephrine-induced changes in clock-controlled gene transcription. Bmal1 is required for neonatal myocyte growth, replication-independent histone turnover, and chromatin organization at the Sik1 promoter. Sik1 represents a novel clock-controlled gene that coordinates myocyte growth with hypertrophic and clock-controlled gene transcription. Replication-independent histone turnover is required for transcriptional remodeling of clock-controlled genes in cardiac myocytes in response to growth stimuli.

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

The copper reductase activity of histone H3 suggests undiscovered characteristics within the protein. Here, we investigated the function of leucine 126 (H3L126), which occupies an axial position relative to the copper binding. Typically found as methionine or leucine in copper-binding proteins, the axial ligand influences the reduction potential of the bound ion, modulating its tendency to accept or yield electrons. We found that mutation of H3L126 to methionine (H3L126M) enhanced the enzymatic activity of native yeast nucleosomes in vitro and increased intracellular levels of Cu1+, leading to improved copper-dependent activities including mitochondrial respiration and growth in oxidative media with low copper. Conversely, H3L126 to histidine (H3L126H) mutation decreased nucleosome's enzymatic activity and adversely affected copper-dependent activities in vivo. Our findings demonstrate that H3L126 fine-tunes the copper reductase activity of nucleosomes and highlights the utility of nucleosome enzymatic activity as a novel paradigm to uncover previously unnoticed features of histones.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and diseases. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 open chromatin regions, including thousands of sequences with cell type-specific accessibility and variants associated with brain gene regulation. In primary cells, we identified 46,802 active enhancer sequences and 164 variants that alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimers disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Pluss utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Amyotrophic Lateral Sclerosis (ALS), like many other neurodegenerative diseases, is highly heritable, but with only a small fraction of cases explained by monogenic disease alleles. To better understand sporadic ALS, we report epigenomic profiles, as measured by ATAC-seq, of motor neuron cultures derived from a diverse group of 380 ALS patients and 80 healthy controls. We find that chromatin accessibility is heavily influenced by sex, the iPSC cell type of origin, ancestry, and the inherent variance arising from sequencing. Once these covariates are corrected for, we are able to identify ALS-specific signals in the data. Additionally, we find that the ATAC-seq data is able to predict ALS disease progression rates with similar accuracy to methods based on biomarkers and clinical status. These results suggest that iPSC-derived motor neurons recapitulate important disease-relevant epigenomic changes.

The mechanisms by which brain insults lead to subsequent epilepsy remain unclear. Insults including trauma, stroke, infections, and long seizures (status epilepticus, SE) increase the nuclear expression and chromatin binding of the neuron-restrictive silencing factor/RE-1 silencing transcription factor (NRSF/REST). REST/NRSF orchestrates major disruption of the expression of key neuronal genes, including ion channels and neurotransmitter receptors, potentially contributing to epileptogenesis. Accordingly, transient interference with REST/NRSF chromatin binding after an epilepsy-provoking SE suppressed spontaneous seizures for the 12 d duration of a prior study. However, whether the onset of epileptogenesis was suppressed or only delayed has remained unresolved. The current experiments determined if transient interference with REST/NRSF chromatin binding prevented epileptogenesis enduringly or, alternatively, slowed epilepsy onset. Epileptogenesis was elicited in adult male rats via systemic kainic acid-induced SE (KA-SE). We then determined if decoy, NRSF-binding-motif oligodeoxynucleotides (NRSE-ODNs), given twice following KA-SE (1) prevented REST/NRSF binding to chromatin, using chromatin immunoprecipitation, or (2) prevented the onset of spontaneous seizures, measured with chronic digital video-electroencephalogram. Blocking NRSF function transiently after KA-SE significantly lengthened the latent period to a first spontaneous seizure. Whereas this intervention did not influence the duration and severity of spontaneous seizures, total seizure number and seizure burden were lower in the NRSE-ODN compared with scrambled-ODN cohorts. Transient interference with REST/NRSF function after KA-SE delays and moderately attenuates insult-related hippocampal epilepsy, but does not abolish it. Thus, the anticonvulsant and antiepileptogenic actions of NRSF are but one of the multifactorial mechanisms generating epilepsy in the adult brain.

Plasticity among cell lineages is a fundamental, but poorly understood, property of regenerative tissues. In the gut tube, the small intestine absorbs nutrients, whereas the colon absorbs electrolytes. In a striking display of inherent plasticity, adult colonic mucosa lacking the chromatin factor SATB2 is converted to small intestine. Using proteomics and CRISPR-Cas9 screening, we identify MTA2 as a crucial component of the molecular machinery that, together with SATB2, restrains colonic plasticity. MTA2 loss in the adult mouse colon activated lipid absorptive genes and functional lipid uptake. Mechanistically, MTA2 co-occupies DNA with HNF4A, an activating pan-intestinal transcription factor (TF), on colonic chromatin. MTA2 loss leads to HNF4A release from colonic chromatin, and accumulation on small intestinal chromatin. SATB2 similarly restrains colonic plasticity through an HNF4A-dependent mechanism. Our study provides a generalizable model of lineage plasticity in which broadly-expressed TFs are retained on tissue-specific enhancers to maintain cell identity and prevent activation of alternative lineages, and their release unleashes plasticity.

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Meiosis is a key step during germ cell differentiation, accompanied by the activation of thousands of genes through germline-specific chromatin reorganization. The chromatin remodeling mechanisms underpinning early meiotic stages remain poorly understood. Here we focus on the function of one of the major autism genes, CHD8, in spermatogenesis, based on the epidemiological association between autism and low fertility rates. Specific ablation of Chd8 in germ cells results in gradual depletion of undifferentiated spermatogonia and the failure of meiotic double-strand break (DSB) formation, leading to meiotic prophase I arrest and cell death. Transcriptional analyses demonstrate that CHD8 is required for extensive activation of spermatogenic genes in spermatogonia, necessary for spermatogonial proliferation and meiosis. CHD8 directly binds and regulates genes crucial for meiosis, including H3K4me3 histone methyltransferase genes, meiotic cohesin genes, HORMA domain-containing genes, synaptonemal complex genes, and DNA damage response genes. We infer that CHD8 contributes to meiotic DSB formation and subsequent meiotic progression through combined regulation of these meiosis-related genes. Our study uncovers an essential role of CHD8 in the proliferation of undifferentiated spermatogonia and the successful progression of meiotic prophase I.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as chromatinopathies. We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an older transcriptomic signature and Alzheimers disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Vocal production learning (vocal learning) is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P

The meiotic chromosome axis coordinates chromosome organization and interhomolog recombination in meiotic prophase and is essential for fertility. In S. cerevisiae, the HORMAD protein Hop1 mediates the enrichment of axis proteins at nucleosome-rich islands through a central chromatin-binding region (CBR). Here, we use cryoelectron microscopy to show that the Hop1 CBR directly recognizes bent nucleosomal DNA through a composite interface in its PHD and winged helix-turn-helix domains. Targeted disruption of the Hop1 CBR-nucleosome interface causes a localized reduction of axis protein binding and meiotic DNA double-strand breaks (DSBs) in axis islands and leads to defects in chromosome synapsis. Synthetic effects with mutants of the Hop1 regulator Pch2 suggest that nucleosome binding delays a conformational switch in Hop1 from a DSB-promoting, Pch2-inaccessible state to a DSB-inactive, Pch2-accessible state to regulate the extent of meiotic DSB formation. Phylogenetic analyses of meiotic HORMADs reveal an ancient origin of the CBR, suggesting that the mechanisms we uncover are broadly conserved.

Tau is a microtubule-associated protein often found in neurofibrillary tangles (NFTs) in the brains of patients with Alzheimers disease. Beyond this context, mounting evidence suggests that tau localizes into the nucleus, where it may play a role in DNA protection and heterochromatin regulation. The molecular mechanisms behind these observations are currently unclear. Using in vitro biophysical experiments, here we demonstrate that tau can undergo liquid-liquid phase separation (LLPS) with DNA, mononucleosomes, and reconstituted nucleosome arrays under low salt conditions. Low concentrations of tau promote chromatin compaction and protect DNA from digestion. While the material state of samples at physiological salt is dominated by chromatin oligomerization, tau can still associate strongly and reversibly with nucleosome arrays. These properties are driven by taus strong interactions with linker and nucleosomal DNA. In addition, tau co-localizes into droplets formed by nucleosome arrays and phosphorylated HP1α, a key heterochromatin constituent thought to function through an LLPS mechanism. Importantly, LLPS and chromatin interactions are disrupted by aberrant tau hyperphosphorylation. These biophysical properties suggest that tau may directly impact DNA and chromatin accessibility and that loss of these interactions could contribute to the aberrant nuclear effects seen in tau pathology.

XIST noncoding RNA promotes the initiation of X chromosome silencing by recruiting the protein SPEN to one X chromosome in female mammals. The SPEN protein is also called SHARP (SMRT and HDAC-associated repressor protein) and MINT (Msx-2 interacting nuclear target) in humans. SPEN recruits N-CoR2 and HDAC3 to initiate histone deacetylation on the X chromosome, leading to the formation of repressive chromatin marks and silencing gene expression. We dissected the contributions of different RNA and protein regions to the formation of a human XIST-SPEN complex in vitro and identified novel sequence and structure determinants that may contribute to X chromosome silencing initiation. Binding of SPEN to XIST RNA requires RRM 4 of the protein, in contrast to the requirement of RRM 3 and RRM 4 for specific binding to SRA RNA. Measurements of SPEN binding to full-length, dimeric, trimeric, or other truncated versions of the A-repeat region revealed that high-affinity binding of XIST to SPEN in vitro requires a minimum of four A-repeat segments. SPEN binding to XIST A-repeat RNA changes the accessibility of the RNA at specific nucleotide sequences, as indicated by changes in RNA reactivity through chemical structure probing. Based on computational modeling, we found that inter-repeat duplexes formed by multiple A-repeats can present an unpaired adenosine in the context of a double-stranded region of RNA. The presence of this specific combination of sequence and structural motifs correlates with high-affinity SPEN binding in vitro. These data provide new information on the molecular basis of the XIST and SPEN interaction.

Relationships between form and function have been the subject of many morphological and evolutionary studies (e.g., Dean et al., 2007; Kolmann et al., 2014). For example, elongated rostra have appeared [...]

We investigate how matrix stiffness regulates chromatin reorganization and cell reprogramming and find that matrix stiffness acts as a biphasic regulator of epigenetic state and fibroblast-to-neuron conversion efficiency, maximized at an intermediate stiffness of 20 kPa. ATAC sequencing analysis shows the same trend of chromatin accessibility to neuronal genes at these stiffness levels. Concurrently, we observe peak levels of histone acetylation and histone acetyltransferase (HAT) activity in the nucleus on 20 kPa matrices, and inhibiting HAT activity abolishes matrix stiffness effects. G-actin and cofilin, the cotransporters shuttling HAT into the nucleus, rises with decreasing matrix stiffness; however, reduced importin-9 on soft matrices limits nuclear transport. These two factors result in a biphasic regulation of HAT transport into nucleus, which is directly demonstrated on matrices with dynamically tunable stiffness. Our findings unravel a mechanism of the mechano-epigenetic regulation that is valuable for cell engineering in disease modeling and regenerative medicine applications.

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidneys active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

SNCAIP duplication may promote Group 4 medulloblastoma via induction of PRDM6, a poorly characterized member of the PRDF1 and RIZ1 homology domain-containing (PRDM) family of transcription factors. Here, we investigated the function of PRDM6 in human hindbrain neuroepithelial stem cells and tested PRDM6 as a driver of Group 4 medulloblastoma. We report that human PRDM6 localizes predominantly to the nucleus, where it causes widespread repression of chromatin accessibility and complex alterations of gene expression patterns. Genome-wide mapping of PRDM6 binding reveals that PRDM6 binds to chromatin regions marked by histone H3 lysine 27 trimethylation that are located within, or proximal to, genes. Moreover, we show that PRDM6 expression in neuroepithelial stem cells promotes medulloblastoma. Surprisingly, medulloblastomas derived from PRDM6-expressing neuroepithelial stem cells match human Group 3, but not Group 4, medulloblastoma. We conclude that PRDM6 expression has oncogenic potential but is insufficient to drive Group 4 medulloblastoma from neuroepithelial stem cells. We propose that both PRDM6 and additional factors, such as specific cell-of-origin features, are required for Group 4 medulloblastoma. Given the lack of PRDM6 expression in normal tissues and its oncogenic potential shown here, we suggest that PRDM6 inhibition may have therapeutic value in PRDM6-expressing medulloblastomas.

The master DNA damage repair histone protein, H2AX, is essential for orchestrating the recruitment of downstream mediator and effector proteins at damaged chromatin. The phosphorylation of H2AX at S139, γH2AX, is well-studied for its DNA repair function. However, the extended C-terminal tail is not characterized. Here, we define the minimal motif on H2AX for the canonical function in activating the MDC1-RNF8-RNF168 phosphorylation-ubiquitination pathway that is important for recruiting repair proteins, such as 53BP1 and BRCA1. Interestingly, H2AX recruits 53BP1 independently from the MDC1-RNF8-RNF168 pathway through its evolved C-terminal linker region with S139 phosphorylation. Mechanistically, 53BP1 recruitment to damaged chromatin is mediated by the interaction between the H2AX C-terminal tail and the 53BP1 Oligomerization-Tudor domains. Moreover, γH2AX-linker mediated 53BP1 recruitment leads to camptothecin resistance in H2AX knockout cells. Overall, our study uncovers an evolved mechanism within the H2AX C-terminal tail for regulating DNA repair proteins at damaged chromatin.

The conformational dynamics of nucleosome arrays generate a diverse spectrum of microscopic states, posing challenges to their structural determination. Leveraging cryogenic electron tomography (cryo-ET), we determine the three-dimensional (3D) structures of individual mononucleosomes and arrays comprising di-, tri-, and tetranucleosomes. By slowing the rate of condensation through a reduction in ionic strength, we probe the intra-array structural transitions that precede inter-array interactions and liquid droplet formation. Under these conditions, the arrays exhibite irregular zig-zag conformations with loose packing. Increasing the ionic strength promoted intra-array compaction, yet we do not observe the previously reported regular 30-nanometer fibers. Interestingly, the presence of H1 do not induce array compaction; instead, one-third of the arrays display nucleosomes invaded by foreign DNA, suggesting an alternative role for H1 in chromatin network construction. We also find that the crucial parameter determining the structure adopted by chromatin arrays is the angle between the entry and exit of the DNA and the corresponding tangents to the nucleosomal disc. Our results provide insights into the initial stages of intra-array compaction, a critical precursor to condensation in the regulation of chromatin organization.

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Background: The role of histone methyltransferase SETDB1 in renal ischemia–reperfusion (I/R) injury has not been explored yet. This study aims to investigate the potential mechanism of SETDB1 in regulating renal I/R injury and its impact on mitochondrial damage and oxidative stress. Methods: The in vivo model of renal I/R in mice and the in vitro model of hypoxia/reoxygenation (H/R) in human renal tubular epithelial cells (HK-2) were constructed to detect the expression of SETDB1. Next, the specific inhibitor (R,R)-59 and knockdown viruses were used to inhibit SETDB1 and verify its effects on mitochondrial damage and oxidative stress. Chromatin immunoprecipitation (ChIP) and coimmunoprecipitation (CoIP) were implemented to explore the in-depth mechanism of SETDB1 regulating renal I/R injury. Results: The study found that SETDB1 had a regulatory role in mitochondrial damage and oxidative stress during renal I/R injury. Notably, SESN2 was identified as a target of SETDB1, and its expression was under the influence of SETDB1. Besides, SESN2 mediated the regulation of SETDB1 on renal I/R injury. Through deeper mechanistic studies, we uncovered that SETDB1 collaborates with heterochromatin HP1β, facilitating the labeling of H3K9me3 on the SESN2 promoter and impeding SESN2 expression. Conclusions: The SETDB1/HP1β-SESN2 axis emerges as a potential therapeutic strategy for mitigating renal I/R injury. [ABSTRACT FROM AUTHOR]

The organization of DNA into nucleosomes is a ubiquitous and ancestral feature that was once thought to be exclusive to the eukaryotic domain of life. Intriguingly, several representatives of the Nucleocytoplasmic Large DNA Viruses (NCLDV) encode histone-like proteins that in Melbournevirus were shown to form nucleosome-like particles. Medusavirus medusae (MM), a distantly related giant virus, encodes all four core histone proteins and, unique amongst most giant viruses, a putative acidic protein with two domains resembling eukaryotic linker histone H1. Here, we report the structure of nucleosomes assembled with MM histones and highlight similarities and differences with eukaryotic and Melbournevirus nucleosomes. Our structure provides insight into how variations in histone tail and loop lengths are accommodated within the context of the nucleosome. We show that MM-histones assemble into tri-nucleosome arrays, and that the putative linker histone H1 does not function in chromatin compaction. These findings expand our limited understanding of chromatin organization by virus-encoded histones. Some large DNA viruses encode histone-like proteins similar to those in eukaryotes, which organize DNA into nucleosome-like particles. Here the authors find that Medusavirus histones form nucleosomes and arrays, but its unique linker histone H1 does not interact with chromatin. This reveals insights into how viruses organize their genetic material. [ABSTRACT FROM AUTHOR]