Search

Human genes are arranged on 23 pairs of chromosomes, but in cancer, tumor-promoting genes and regulatory elements can free themselves from chromosomes and relocate to circular, extrachromosomal pieces of DNA (ecDNA). ecDNA, because of its nonchromosomal inheritance, drives high-copy-number oncogene amplification and enables tumors to evolve their genomes rapidly. Furthermore, the circular ecDNA architecture fundamentally alters gene regulation and transcription, and the higher-order organization of ecDNA contributes to tumor pathogenesis. Consequently, patients whose cancers harbor ecDNA have significantly shorter survival. Although ecDNA was first observed more than 50 years ago, its critical importance has only recently come to light. In this review, we discuss the current state of understanding of how ecDNAs form and function as well as how they contribute to drug resistance and accelerated cancer evolution. Expected final online publication date for the Annual Review of Pathology: Mechanisms of Disease, Volume 17 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The human transcriptome contains many types of noncoding RNAs, which rival the number of protein-coding species. From long noncoding RNAs (lncRNAs) that are over 200 nucleotides long to piwi-interacting RNAs (piRNAs) of only 20 nucleotides, noncoding RNAs play important roles in regulating transcription, epigenetic modifications, translation, and cell signaling. Roles for noncoding RNAs in disease mechanisms are also being uncovered, and several species have been identified as potential drug targets. On May 11-14, 2021, the Keystone eSymposium "Noncoding RNAs: Biology and Applications" brought together researchers working in RNA biology, structure, and technologies to accelerate both the understanding of RNA basic biology and the translation of those findings into clinical applications.

The chromodomain helicase DNA-binding protein CHD8 is among the most frequently found de-novo mutations in autism (1–3). Unlike most other autism-risk genes, CHD8 mutations appear to be fully penetrant (4). Despite its prominent disease involvement, little known about its molecular function. Based on sequence homology, CHD8 is believed to be a chromatin regulator, but mechanisms for its genomic targeting and its role on chromatin are unclear. Here, we developed a human cell model carrying conditional CHD8 loss-of-function alleles. Full knockout CHD8 was required for the viability of undifferentiated human embryonic stem (ES) cells, whereas postmitotic neurons survived following CHD8 depletion. However, chromatin accessibility maps and transcriptional profiling revealed that CHD8 is a potent general chromatin activator, enhancing transcription of its direct target genes, including a large group of autism genes. CHD8’s genomic binding sites in human neurons were significantly enriched for ELK1 (ETS) motifs. Moreover, positive CHD8-dependent chromatin remodeling was enhanced at ELK1 motif-containing CHD8 binding sites. ELK1 was the most prominent ETS factor expressed in human neurons and was necessary for CHD8 to target the sites that contained the ELK1 motif, demonstrating a cooperative interaction between ELK1 and CHD8 on chromatin. We also observed potential role of CHD8 in ELK1 localization on nuclear compartments in a transcription-stage-dependent manner. Finally, inhibition of ELK1 activity or ELK1 knockdown that enhances the neurogenesis from embryonic stem cells (ES) was dependent on the presence of CHD8. In summary, our results establish that CHD8 is a strong activator of chromatin accessibility and transcription in neurons and reveals a role in regulating many high-risk autism genes. Additionally, we show there is molecular and functional interdependence of CHD8 and ELK1 in chromatin binding of CHD8, nuclear interaction of ELK1, and neurogenesis enhancement. These data imply the involvement of the MAPK/ERK pathway effector ELK1 in pathogenesis of autism caused by CHD8 mutations (5).

Although alterations in chromatin structure are known to exist in tumors, how these alterations relate to molecular phenotypes in cancer remains to be demonstrated. Multi-omics profiling of human tumors can provide insight into how alterations in chromatin structure are propagated through the pathway of gene expression to result in malignant protein expression. We applied multi-omics profiling of chromatin accessibility, RNA abundance, and protein abundance to 36 human thyroid cancer primary tumors, metastases, and patient-match normal tissue. Through quantification of chromatin accessibility associated with active transcription units and global protein expression, we identify a local chromatin structure that is highly correlated with coordinated RNA and protein expression. In particular, we identify enhancers located within gene-bodies as predictive of correlated RNA and protein expression, that is independent of overall transcriptional activity. To demonstrate the generalizability of these findings we also identify similar results in an independent cohort of human breast cancers. Taken together, these analyses suggest that local enhancers, rather than distal enhancers, are likely most predictive of cancer gene expression phenotypes. This allows for identification of potential targets for cancer therapeutic approaches and reinforces the utility of multi-omics profiling as a methodology to understand human disease., Studies show the cancer transcriptome correlates poorly with the cancer proteome, questioning the role of chromatin regulation. Here the authors demonstrate proximal-gene-body chromatin elements and transcription predict abundances of differentially expressed proteins in thyroid and breast cancers.

Long noncoding RNAs (lncRNAs) are a diverse and pervasive class of genes. Recent studies in the mammalian brain have uncovered several novel mechanisms. LncRNA loci are often located in proximity to developmental transcriptional factors. The lncRNA product may act like a transcription factor to control distantly located genes, or in other instances, the lncRNA loci contain DNA regulatory elements that act locally on neighboring genes. Circular RNAs are covalently closed single-stranded RNAs that can control neuronal function by acting as microRNA sponges and additional mechanisms. LncRNAs can also engage in target-directed microRNA degradation to shape the pool of microRNAs and translation. Thus, diverse mechanisms allow lncRNAs to act in the nucleus and cytoplasm to control neuronal fate and function.

Genome-wide association studies (GWAS) of neurological diseases have identified thousands of variants associated with disease phenotypes. However, the majority of these variants do not alter coding sequences, making it difficult to assign their function. Here, we present a multi-omic epigenetic atlas of the adult human brain through profiling of single-cell chromatin accessibility landscapes and three-dimensional (3D) chromatin interactions of diverse adult brain regions across a cohort of cognitively healthy individuals. We developed a machine-learning classifier to integrate this multi-omic framework and predict dozens of functional single-nucleotide polymorphisms (SNPs) for Alzheimer’s disease (AD) and Parkinson’s disease (PD), nominating target genes and cell types for previously orphaned GWAS loci. Moreover, we dissected the complex inverted haplotype of the MAPT (encoding tau) PD risk locus, identifying putative ectopic regulatory interactions in neurons that may mediate this disease association. This work expands our understanding of inherited variation and provides a roadmap for the epigenomic dissection of causal regulatory variation in disease.

In celebration of the 20th anniversary of Nature Reviews Genetics, we asked 12 leading researchers to reflect on the key challenges and opportunities faced by the field of genetics and genomics. Keeping their particular research area in mind, they take stock of the current state of play and emphasize the work that remains to be done over the next few years so that, ultimately, the benefits of genetic and genomic research can be felt by everyone., To celebrate the first 20 years of Nature Reviews Genetics, we asked 12 leading scientists to reflect on the key challenges and opportunities faced by the field of genetics and genomics., The contributors Amy L. McGuire is the Leon Jaworski Professor of Biomedical Ethics and Director of the Center for Medical Ethics and Health Policy at Baylor College of Medicine. She has received numerous teaching awards at Baylor College of Medicine, was recognized by the Texas Executive Women as a Woman on the Move in 2016 and was invited to give a TedMed talk titled “There is No Genome for the Human Spirit” in 2014. In 2020, she was elected as a Hastings Center Fellow. Her research focuses on ethical and policy issues related to emerging technologies, with a particular focus on genomic research, personalized medicine and the clinical integration of novel neurotechnologies. Stacey Gabriel is the Senior Director of the Genomics Platform at the Broad Institute since 2012 and has led platform development, execution and operation since its founding. She is Chair of Institute Scientists and serves on the institute’s executive leadership team. She is widely recognized as a leader in genomic technology and project execution. She has led the Broad’s contributions to numerous flagship projects in human genetics, including the International HapMap Project, the 1000 Genomes Project, The Cancer Genome Atlas, the National Heart, Lung, and Blood Institute’s Exome Sequencing Project and the TOPMed programme. She is Principal Investigator of the Broad’s All of Us (AoU) Genomics Center and serves on the AoU Program Steering Committee. Sarah A. Tishkoff is the David and Lyn Silfen University Associate Professor in Genetics and Biology at the University of Pennsylvania, Philadelphia, USA, and holds appointments in the School of Medicine and the School of Arts and Sciences. She is a member of the US National Academy of Sciences and a recipient of an NIH Pioneer Award, a David and Lucile Packard Career Award, a Burroughs/Wellcome Fund Career Award and an American Society of Human Genetics Curt Stern Award. Her work focuses on genomic variation in Africa, human evolutionary history, the genetic basis of adaptation and phenotypic variation in Africa, and the genetic basis of susceptibility to infectious disease in Africa. Ambroise Wonkam is Professor of Medical Genetics, Director of GeneMAP (Genetic Medicine of African Populations Research Centre) and Deputy Dean Research in the Faculty of Health Sciences, University of Cape Town, South Africa. He has successfully led numerous NIH- and Wellcome Trust-funded projects over the past decade to investigate clinical variability in sickle cell disease, hearing impairment genetics and the return of individual findings in genetic research in Africa. He won the competitive Clinical Genetics Society International Award for 2014 from the British Society of Genetic Medicine. He is president of the African Society of Human Genetics. Aravinda Chakravarti is Director of the Center for Human Genetics and Genomics, the Muriel G. and George W. Singer Professor of Neuroscience and Physiology, and Professor of Medicine at New York University School of Medicine. He is an elected member of the US National Academy of Sciences, the US National Academy of Medicine and the Indian National Science Academy. He has been a key participant in the Human Genome Project, the International HapMap Project and the 1000 Genomes Project. His research attempts to understand the molecular basis of multifactorial disease. He was awarded the 2013 William Allan Award by the American Society of Human Genetics and the 2018 Chen Award by the Human Genome Organization. Eileen E. M. Furlong is Head of the Genome Biology Department at the European Molecular Biology Laboratory (EMBL) and a member of the EMBL Directorate. She is an elected member of the European Molecular Biology Organization (EMBO) and the Academia Europaea, and a European Research Council (ERC) advanced investigator. Her group dissects fundamental principles of how the genome is regulated and how it drives cell fate decisions during embryonic development, including how developmental enhancers are organized and function within the 3D nucleus. Her work combines genetics, (single-cell) genomics, imaging and computational approaches to understand these processes. Her research has advanced the development of genomic methods for use in complex multicellular organisms. Barbara Treutlein is Associate Professor of Quantitative Developmental Biology in the Department of Biosystems Science and Engineering of ETH Zurich in Basel, Switzerland. Her group uses and develops single-cell genomics approaches in combination with stem cell-based 2D and 3D culture systems to study how human organs develop and regenerate and how cell fate is regulated. For her work, Barbara has received multiple awards, including the Friedmund Neumann Prize of the Schering Foundation, the Dr. Susan Lim Award for Outstanding Young Investigator of the International Society of Stem Cell Research and the EMBO Young Investigator Award. Alexander Meissner is a scientific member of the Max Planck Society and currently Managing Director of the Max Planck Institute (MPI) for Molecular Genetics in Berlin, Germany. He heads the Department of Genome Regulation and is a visiting scientist in the Department of Stem Cell and Regenerative Biology at Harvard University. Before his move to the MPI, he was a tenured professor at Harvard University and a senior associate member of the Broad Institute, where he co-directed the epigenomics programme. In 2018, he was elected as an EMBO member. His laboratory uses genomic tools to study developmental and disease biology with a particular focus on epigenetic regulation. Howard Y. Chang is the Virginia and D. K. Ludwig Professor of Cancer Genomics at Stanford University and an investigator at the Howard Hughes Medical Institute. He is a physician–scientist who has focused on deciphering the hidden information in the non-coding genome. His laboratory is best known for studies of long non-coding RNAs in gene regulation and development of new epigenomic technologies. He is an elected member of the US National Academy of Sciences, the US National Academy of Medicine, and the American Academy of Arts and Sciences. Núria López-Bigas is ICREA research Professor at the Institute for Research in Biomedicine and Associate Professor at the University Pompeu Fabra. She obtained an ERC Consolidator Grant in 2015 and was elected as an EMBO member in 2016. Her work has been recognized with the prestigious Banc de Sabadell Award for Research in Biomedicine, the Catalan National Award for Young Research Talent and the Career Development Award from the Human Frontier Science Program. Her research focuses on the identification of cancer driver mutations, genes and pathways across tumour types and in understanding the mutational processes that lead to the accumulation of mutations in cancer cells. Eran Segal is Professor in the Department of Computer Science and Applied Mathematics at the Weizmann Institute of Science, heading a multidisciplinary laboratory with extensive experience in machine learning, computational biology and analysis of heterogeneous high-throughput genomic data. His research focuses on the microbiome, nutrition and genetics, and their effect on health and disease and aims to develop personalized medicine based on big data from human cohorts. He has published more than 150 publications and received several awards and honours for his work, including the Overton and the Michael Bruno awards. He was recently elected as an EMBO member and as a member of the Israel Young Academy. Jin-Soo Kim is Director of the Center for Genome Engineering in the Institute for Basic Science in Daejon, South Korea. He has received numerous awards, including the 2017 Asan Award in Medicine, the 2017 Yumin Award in Science and the 2019 Research Excellence Award (Federation of Asian and Oceanian Biochemists and Molecular Biologists). He was featured as one of ten Science Stars of East Asia in Nature (558, 502–510 (2018)) and has been recognized as a highly cited researcher by Clarivate Analytics since 2018. His work focuses on developing tools for genome editing in biomedical research.

Extrachromosomal DNA (ecDNA) amplification promotes intratumoral genetic heterogeneity and accelerated tumor evolution1–3, but its frequency and clinical impact are unclear. Here we show, using computational analysis of whole-genome sequencing data from 3,212 cancer patients, that ecDNA amplification frequently occurs in most cancer types, but not in blood or normal tissue. Oncogenes were highly enriched on amplified ecDNA and the most common recurrent oncogene amplifications arise on ecDNA. EcDNA amplifications resulted in higher levels of oncogene transcription compared to copy number matched linear DNA, coupled with enhanced chromatin accessibility and more frequently resulted in transcript fusions. Patients whose cancers carry ecDNAs have significantly shorter survival, even when controlled for tissue type, than do patients whose cancers are not driven by ecDNA-based oncogene amplification. The results presented here demonstrate that ecDNA-based oncogene amplification is common in cancer, is different from chromosomal amplification and drives poor outcome for patients across many cancer types.

The immune system is composed of a diverse array of cell types, each with a specialized role in orchestrating the immune response to pathogens or cancer. Even within a single cell 'type,' individual cells can access a wide spectrum of differentiation and activation states, which reflect the physiological response of each cell to the tissue environment and immune stimuli. Thus, the cellular diversity of the immune system is inherently quite complex and understanding this complexity has greatly benefited from technologies that measure immune responses at single-cell resolution, in addition to the systems-level response as a whole. In this Commentary, we focus on recent work at the interface of immunology and single-cell genomics and highlight advances in technologies and their application to immune cells. In particular, we highlight recent single-cell genomic profiling studies of T cells, since somatic rearrangements in the T cell receptor (TCR) loci enable the tracking of clonal T cell responses through space and time. Finally, we discuss opportunities for future use of these technologies in understanding vaccination and the basis for effective vaccine-induced immunity.

Genetic variations underlying susceptibility to complex autoimmune and allergic diseases are concentrated within noncoding regulatory elements termed enhancers1. The functions of a large majority of disease-associated enhancers are unknown, in part owing to their distance from the genes they regulate, a lack of understanding of the cell types in which they operate, and our inability to recapitulate the biology of immune diseases in vitro. Here, using shared synteny to guide loss-of-function analysis of homologues of human enhancers in mice, we show that the prominent autoimmune and allergic disease risk locus at chromosome 11q13.52–7 contains a distal enhancer that is functional in CD4+ regulatory T (Treg) cells and required for Treg-mediated suppression of colitis. The enhancer recruits the transcription factors STAT5 and NF-κB to mediate signal-driven expression of Lrrc32, which encodes the protein glycoprotein A repetitions predominant (GARP). Whereas disruption of the Lrrc32 gene results in early lethality, mice lacking the enhancer are viable but lack GARP expression in Foxp3+ Treg cells, which are unable to control colitis in a cell-transfer model of the disease. In human Treg cells, the enhancer forms conformational interactions with the promoter of LRRC32 and enhancer risk variants are associated with reduced histone acetylation and GARP expression. Finally, functional fine-mapping of 11q13.5 using CRISPR-activation (CRISPRa) identifies a CRISPRa-responsive element in the vicinity of risk variant rs11236797 capable of driving GARP expression. These findings provide a mechanistic basis for association of the 11q13.5 risk locus with immune-mediated diseases and identify GARP as a potential target in their therapy. Shared synteny guides loss-of-function analysis of human enhancer homologues in mice, identifying a distal enhancer at the autoimmune and allergic disease risk locus at chromosome 11q13.5 whose function in regulatory T cells provides a mechanistic basis for its role in disease.

The on-target pioneer factors Ascl1 and Myod1 are sequence-related but induce two developmentally unrelated lineages—that is, neuronal and muscle identities, respectively. It is unclear how these two basic helix–loop–helix (bHLH) factors mediate such fundamentally different outcomes. The chromatin binding of Ascl1 and Myod1 was surprisingly similar in fibroblasts, yet their transcriptional outputs were drastically different. We found that quantitative binding differences explained differential chromatin remodelling and gene activation. Although strong Ascl1 binding was exclusively associated with bHLH motifs, strong Myod1-binding sites were co-enriched with non-bHLH motifs, possibly explaining why Ascl1 is less context dependent. Finally, we observed that promiscuous binding of Myod1 to neuronal targets results in neuronal reprogramming when the muscle program is inhibited by Myt1l. Our findings suggest that chromatin access of on-target pioneer factors is primarily driven by the protein–DNA interaction, unlike ordinary context-dependent transcription factors, and that promiscuous transcription factor binding requires specific silencing mechanisms to ensure lineage fidelity. Lee, Mall et al. explore why the sequence-related transcription factors Myod1 and Ascl1 lead to different reprogramming outcomes when expressed in fibroblasts, despite binding similar targets.

To image the accessible genome at nanometer scale in situ, we developed three-dimensional assay for transposase-accessible chromatin-photoactivated localization microscopy (3D ATAC-PALM) that integrates an assay for transposase-accessible chromatin with visualization, PALM super-resolution imaging and lattice light-sheet microscopy. Multiplexed with oligopaint DNA-fluorescence in situ hybridization (FISH), RNA-FISH and protein fluorescence, 3D ATAC-PALM connected microscopy and genomic data, revealing spatially segregated accessible chromatin domains (ACDs) that enclose active chromatin and transcribed genes. Using these methods to analyze genetically perturbed cells, we demonstrated that genome architectural protein CTCF prevents excessive clustering of accessible chromatin and decompacts ACDs. These results highlight 3D ATAC-PALM as a useful tool to probe the structure and organizing mechanism of the genome.

SUMMARYThe p53 transcription factor, encoded by the most frequently mutated gene in human cancer, plays a critical role in tissue homeostasis in response to stress signals. The mechanisms through which p53 promotes downstream tumor suppressive gene expression programs remain, however, only superficially understood. Here, we used tandem affinity purification and mass spectrometry to reveal new components of the p53 response. This approach uncovered Mettl3, a component of the m6A RNA methyltransferase complex (MTC), as a p53-interacting protein. Analysis of Mettl3- deficient cells revealed that Mettl3 promotes p53 protein stabilization and target gene expression in response to DNA damage. Mettl3 acts in part by competing with the p53 negative regulator, Mdm2, for binding to the p53 transactivation domains to promote methyltransferase-independent stabilization of p53. In addition, Mettl3 relies on its catalytic activity to augment p53 responses, with p53 recruiting Mettl3 to p53 target genes to co-transcriptionally direct m6A modification of p53 pathway transcripts to enhance their expression. Mettl3 also promotes p53 activity downstream of oncogenic signals in vivo, in both allograft and autochthonous lung adenocarcinoma models, suggesting cooperative action of p53 and Mettl3 in tumor suppression. Accordingly, we found in diverse human cancers that mutations in MTC components perturb expression of p53 target genes and that MTC mutations are mutually exclusive with TP53 mutations, suggesting that the MTC enhances the p53 transcriptional program in human cancer. Together, these studies reveal a fundamental role for Mettl3 in amplifying p53 signaling through protein stabilization and epitranscriptome regulation.

Transcription factors bind DNA sequence motif vocabularies in cis-regulatory elements (CREs) to modulate chromatin state and gene expression during cell state transitions. A quantitative understanding of how motif lexicons influence dynamic regulatory activity has been elusive due to the combinatorial nature of the cis-regulatory code. To address this, we undertook multiomic data profiling of chromatin and expression dynamics across epidermal differentiation to identify 40,103 dynamic CREs associated with 3,609 dynamically expressed genes, then applied an interpretable deep-learning framework to model the cis-regulatory logic of chromatin accessibility. This analysis framework identified cooperative DNA sequence rules in dynamic CREs regulating synchronous gene modules with diverse roles in skin differentiation. Massively parallel reporter assay analysis validated temporal dynamics and cooperative cis-regulatory logic. Variants linked to human polygenic skin disease were enriched in these time-dependent combinatorial motif rules. This integrative approach shows the combinatorial cis-regulatory lexicon of epidermal differentiation and represents a general framework for deciphering the organizational principles of the cis-regulatory code of dynamic gene regulation.

Significance In the skin, tissue injury results in fibrosis in the form of a scar composed of dense extracellular matrix deposited by fibroblasts. Therapies that promote tissue regeneration rather than fibrosis remain elusive because principles of fibroblast programming and response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the study of cell populations in complex tissue, which has allowed us to characterize wound healing fibroblasts across both time and space. We identify functionally distinct fibroblast subpopulations and track cell fate during the response to wounding. We demonstrate that populations of fibroblasts migrate, proliferate, and differentiate in an adaptive response to disruption of their environment. These results illustrate fundamental principles underlying the cellular response to tissue injury., In the skin, tissue injury results in fibrosis in the form of scars composed of dense extracellular matrix deposited by fibroblasts. The therapeutic goal of regenerative wound healing has remained elusive, in part because principles of fibroblast programming and adaptive response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the comprehensive study of cell populations in complex tissue, which has allowed us to characterize the cells involved in wound healing across both time and space. We employ a stented wound model that recapitulates human tissue repair kinetics and multiple Rainbow transgenic lines to precisely track fibroblast fate during the physiologic response to skin injury. Through integrated analysis of single cell chromatin landscapes and gene expression states, coupled with spatial transcriptomic profiling, we are able to impute fibroblast epigenomes with temporospatial resolution. This has allowed us to reveal potential mechanisms controlling fibroblast fate during migration, proliferation, and differentiation following skin injury, and thereby reexamine the canonical phases of wound healing. These findings have broad implications for the study of tissue repair in complex organ systems.

Cells communicate with each other via receptor-ligand interactions. Here, we describe lentiviral-mediated cell entry by engineered receptor-ligand interaction (ENTER) to display ligand proteins, deliver payloads, and record receptor specificity. We optimize ENTER to decode interactions between T cell receptor (TCR)-MHC peptides, antibody-antigen, and other receptor-ligand pairs. A viral presentation strategy allows ENTER to capture interactions between B cell receptor and any antigen. We engineer ENTER to deliver genetic payloads to antigen-specific T or B cells to selectively modulate cellular behavior in mixed populations. Single-cell readout of ENTER by RNA sequencing (ENTER-seq) enables multiplexed enumeration of antigen specificities, TCR clonality, cell type, and states of individual T cells. ENTER-seq of CMV-seropositive patient blood samples reveals the viral epitopes that drive effector memory T cell differentiation and inter-clonal vs. intra-clonal phenotypic diversity targeting the same epitope. ENTER technology enables systematic discovery of receptor specificity, linkage to cell fates, and antigen-specific cargo delivery.

Mitochondria are essential organelles in eukaryotic cells that provide critical support for energetic and metabolic homeostasis. Mutations that accumulate in mitochondrial DNA (mtDNA) in somatic cells have been implicated in cancer, degenerative diseases, and the aging process. However, the mechanisms used by somatic cells to maintain proper functions despite their mtDNA mutation load are poorly understood. Here, we analyzed somatic mtDNA mutations in more than 30,000 human single peripheral and bone marrow mononuclear cells and observed a significant overrepresentation of homoplastic mtDNA mutations in B, T and NK lymphocytes despite their lower mutational burden than other hematopoietic cells. The characteristic mutational landscape of mtDNA in lymphocytes were validated with data from multiple platforms and individuals. Single-cell RNA-seq and computational modeling demonstrated a stringent mitochondrial bottleneck during lymphocyte development likely caused by lagging mtDNA replication relative to cell proliferation. These results illuminate a potential mechanism used by highly metabolically active immune cells for quality control of their mitochondrial genomes.

Most autosomal genes are thought to be expressed from both alleles, with some notable exceptions, including imprinted genes and genes showing random monoallelic expression (RME). The extent and nature of RME has been the subject of debate. Here we investigate the expression of several candidate RME genes in F1 hybrid mouse cells before and after differentiation, to define how they become persistently, monoallelically expressed. Clonal monoallelic expression is not present in embryonic stem cells, but we observe high frequencies of monoallelism in neuronal progenitor cells by assessing expression status in more than 200 clones. We uncover unforeseen modes of allelic expression that appear to be gene-specific and epigenetically regulated. This non-canonical allelic regulation has important implications for development and disease, including autosomal dominant disorders and opens up therapeutic perspectives., Some autosomal genes are expressed in a random monoallelic manner, but its extent and mechanisms have remained unclear. Here the authors show robust monoallelic expression in cell lines and mice, where the silent allele can be reexpressed using epidrugs. Further, they find these genes display various modalities of allelic expression with different degrees of allelic imbalance.

Pathologic skin scarring presents a vast economic and medical burden. Unfortunately, the molecular mechanisms underlying scar formation remain to be elucidated. We used a hypertrophic scarring (HTS) mouse model in which Jun is overexpressed globally or specifically in α-smooth muscle or collagen type I–expressing cells to cause excessive extracellular matrix deposition by skin fibroblasts in the skin after wounding. Jun overexpression triggered dermal fibrosis by modulating distinct fibroblast subpopulations within the wound, enhancing reticular fibroblast numbers, and decreasing lipofibroblasts. Analysis of human scars further revealed that JUN is highly expressed across the wide spectrum of scars, including HTS and keloids. CRISPR-Cas9–mediated JUN deletion in human HTS fibroblasts combined with epigenomic and transcriptomic analysis of both human and mouse HTS fibroblasts revealed that JUN initiates fibrosis by regulating CD36. Blocking CD36 with salvianolic acid B or CD36 knockout model counteracted JUN-mediated fibrosis efficacy in both human fibroblasts and mouse wounds. In summary, JUN is a critical regulator of pathological skin scarring, and targeting its downstream effector CD36 may represent a therapeutic strategy against scarring.

Identifying the causes of human diseases requires deconvolution of abnormal molecular phenotypes spanning DNA accessibility, gene expression and protein abundance1–3. We present a single-cell framework that integrates highly multiplexed protein quantification, transcriptome profiling and analysis of chromatin accessibility. Using this approach, we establish a normal epigenetic baseline for healthy blood development, which we then use to deconvolve aberrant molecular features within blood from patients with mixed-phenotype acute leukemia4,5. Despite widespread epigenetic heterogeneity within the patient cohort, we observe common malignant signatures across patients as well as patient-specific regulatory features that are shared across phenotypic compartments of individual patients. Integrative analysis of transcriptomic and chromatin-accessibility maps identified 91,601 putative peak-to-gene linkages and transcription factors that regulate leukemia-specific genes, such as RUNX1-linked regulatory elements proximal to the marker gene CD69. These results demonstrate how integrative, multiomic analysis of single cells within the framework of normal development can reveal both distinct and shared molecular mechanisms of disease from patient samples. Analyzing DNA accessibility, transcriptome and protein expression in single cells uncovers cancer-regulatory programs.

Oncogenes are commonly amplified on particles of extrachromosomal DNA (ecDNA) in cancer1,2, but our understanding of the structure of ecDNA and its effect on gene regulation is limited. Here, by integrating ultrastructural imaging, long-range optical mapping and computational analysis of whole-genome sequencing, we demonstrate the structure of circular ecDNA. Pan-cancer analyses reveal that oncogenes encoded on ecDNA are among the most highly expressed genes in the transcriptome of the tumours, linking increased copy number with high transcription levels. Quantitative assessment of the chromatin state reveals that although ecDNA is packaged into chromatin with intact domain structure, it lacks higher-order compaction that is typical of chromosomes and displays significantly enhanced chromatin accessibility. Furthermore, ecDNA is shown to have a significantly greater number of ultra-long-range interactions with active chromatin, which provides insight into how the structure of circular ecDNA affects oncogene function, and connects ecDNA biology with modern cancer genomics and epigenetics.

Induction of the transcription factor Irf8 in the common dendritic cell progenitor (CDP) is required for classical type 1 dendritic cell (cDC1) fate specification, but the mechanisms controlling this induction are unclear. In the present study Irf8 enhancers were identified via chromatin profiling of dendritic cells and CRISPR/Cas9 genome editing was used to assess their roles in Irf8 regulation. An enhancer 32 kilobases (kb) downstream of the Irf8 transcriptional start site (+32-kb Irf8) that was active in mature cDC1s was required for the development of this lineage, but not for its specification. Instead, a +41-kb Irf8 enhancer, previously thought to be active only in plasmacytoid dendritic cells, was found to also be transiently accessible in cDC1 progenitors, and deleting this enhancer prevented the induction of Irf8 in CDPs and abolished cDC1 specification. Thus, cryptic activation of the +41-kb Irf8 enhancer in dendritic cell progenitors is responsible for cDC1 fate specification. The transcription factor IRF8 is essential for classical type 1 dendritic cell (cDC1) development. Murphy and colleagues investigate in detail the molecular control of cDC1 fate specification by systematically unpicking the IRF8 enhancer regions.

Lamin A/C (LMNA) is one of the most frequently mutated genes associated with dilated cardiomyopathy (DCM). DCM related to mutations in LMNA is a common inherited cardiomyopathy that is associated with systolic dysfunction and cardiac arrhythmias. Here we modelled the LMNA-related DCM in vitro using patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Electrophysiological studies showed that the mutant iPSC-CMs displayed aberrant calcium homeostasis that led to arrhythmias at the single-cell level. Mechanistically, we show that the platelet-derived growth factor (PDGF) signalling pathway is activated in mutant iPSC-CMs compared to isogenic control iPSC-CMs. Conversely, pharmacological and molecular inhibition of the PDGF signalling pathway ameliorated the arrhythmic phenotypes of mutant iPSC-CMs in vitro. Taken together, our findings suggest that the activation of the PDGF pathway contributes to the pathogenesis of LMNA-related DCM and point to PDGF receptor-β (PDGFRB) as a potential therapeutic target.

The vast majority of polymorphisms for human dermatologic diseases fall in noncoding DNA regions, leading to difficulty interpreting their functional significance. Recent work using chromosome conformation capture technology in combination with chromatin immunoprecipitation (ChIP) has provided a systematic means of linking noncoding variants in active enhancer loci to putative gene targets. Here, we apply H3K27ac HiChIP high-resolution contact maps, generated from primary human T-cell subsets (CD4+ naive, T helper type 17, and regulatory T cells), to 21 dermatologic conditions associated with single nucleotide polymorphisms from 106 genome-wide association studies. This "enhancer connectome" identified 1,492 HiChIP gene targets from 542 noncoding SNPs (P ≤ 5.0 × 10–8). SNP-containing enhancers from inflammatory skin conditions were significantly enriched at the HLA locus and also targeted several key factors from the JAK-STAT signaling pathway, but nonimmune conditions were not. A focused profiling of systemic lupus erythematosus HiChIP genes identified enhancer interactions with factors important for effector CD4+ T-cell differentiation and function, including IRF8 and members of the Ikaros family of zinc-finger proteins. Our results show the ability of the enhancer connectome to nominate functionally relevant candidates from genome-wide association study–identified variants, representing a powerful tool to guide future studies into the genomic regulatory mechanisms underlying dermatologic diseases.

In the skin, tissue injury results in fibrosis in the form of scars composed of dense extracellular matrix deposited by fibroblasts. The therapeutic goal of regenerative wound healing has remained elusive in part because principles of fibroblast programming and adaptive response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the comprehensive study of cell populations in complex tissue, which has allowed us to characterize the cells involved in wound healing across both time and space. We employ a stented wound model that recapitulates human tissue repair kinetics and multiple Rainbow transgenic lines to precisely track fibroblast fate during the physiologic response to injury. Through integrated analysis of single cell chromatin landscapes and gene expression states, coupled with spatial transcriptomic profiling, we are able to impute fibroblast epigenomes with temporospatial resolution. This has allowed us to define the mechanisms controlling cell fate during migration, proliferation, and differentiation following tissue injury and thereby reexamine the canonical phases of wound healing. These findings have broad implications for the study of tissue repair in complex organ systems.

Metastasis is the leading cause of cancer-related deaths, enabling cancer cells to expand to secondary sites and compromise organ function1. Given that primary tumors and metastases often share the same constellation of driver mutations2–4, the mechanisms driving their distinct phenotypes are unclear. Here, we show that inactivation of the frequently mutated tumor suppressor gene, liver kinase B1 (LKB1), has evolving effects throughout lung cancer progression, leading to the differential epigenetic re-programming of early-stage primary tumors compared to late-stage metastases. By integrating genome-scale CRISPR/Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumors. Using an in vivo model of metastatic progression, we further reveal that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumors. SOX17 expression is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can appear to change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

SUMMARYThe long noncoding RNA (lncRNA) XIST establishes X chromosome inactivation (XCI) in female cells in early development and thereafter is thought to be largely dispensable. Here we show XIST is continually required in adult human B cells to silence a subset of X-linked immune genes such asTLR7. XIST-dependent genes lack promoter DNA methylation and require continual XIST-dependent histone deacetylation. XIST RNA-directed proteomics and CRISPRi screen reveal distinctive somatic cell-specific XIST complexes, and identify TRIM28 that mediates Pol II pausing at promoters of X-linked genes in B cells. XIST dysregylation, reflected by escape of XIST-dependent genes, occurs in CD11c+ atypical memory B cells across single-cell transcriptome data in patients with female-biased autoimmunity and COVID-19 infection. XIST inactivation with TLR7 agonism suffices to promote isotype-switched atypical B cells. These results suggest cell-type-specific diversification of lncRNA-protein complexes increase lncRNA functionalities, and expand roles for XIST in sex-differences in biology and medicine.HIGHLIGHTSXIST prevents escape of genes with DNA hypomethylated promoters in B cells.XIST maintains X-inactivation through continuous deacetylation of H3K27ac.XIST ChIRP-MS and allelic CRISPRi screen reveal a B cell-specific XIST cofactor TRIM28.XIST loss and TLR7 stimulation promotes CD11c+ atypical B cell formation.

RNA-binding proteins (RBPs) are essential for skeletal muscle regeneration and RBP dysfunction causes muscle degeneration and neuromuscular disease. How ubiquitously expressed RBPs orchestrate complex tissue regeneration and direct cell fate decisions in skeletal muscle remains poorly understood. Single cell RNA-sequencing of regenerating skeletal muscle reveals that RBP expression, including numerous neuromuscular disease-associated RBPs, is temporally regulated in skeletal muscle stem cells and correlates to stages of myogenic differentiation. By combining machine learning with RBP engagement scoring, we discover that the neuromuscular disease associated RBP Hnrnpa2b1 is a differentiation-specifying regulator of myogenesis controlling myogenic cell fate transitions during terminal differentiation. The timing of RBP expression specifies cell fate transitions by providing a layer of post-transcriptional regulation needed to coordinate stem cell fate decisions during complex tissue regeneration.

The advent of single-cell chromatin accessibility profiling has accelerated the ability to map gene regulatory landscapes but has outpaced the development of scalable software to rapidly extract biological meaning from these data. Here we present a software suite for single-cell analysis of regulatory chromatin in R (ArchR; https://www.archrproject.com/) that enables fast and comprehensive analysis of single-cell chromatin accessibility data. ArchR provides an intuitive, user-focused interface for complex single-cell analyses, including doublet removal, single-cell clustering and cell type identification, unified peak set generation, cellular trajectory identification, DNA element-to-gene linkage, transcription factor footprinting, mRNA expression level prediction from chromatin accessibility and multi-omic integration with single-cell RNA sequencing (scRNA-seq). Enabling the analysis of over 1.2 million single cells within 8 h on a standard Unix laptop, ArchR is a comprehensive software suite for end-to-end analysis of single-cell chromatin accessibility that will accelerate the understanding of gene regulation at the resolution of individual cells., ArchR is a software suite that enables efficient and end-to-end analysis of single-cell chromatin accessibility data (scATAC-seq).

Invariant natural killer T (iNKT) cells are a T-cell subset with potent immunomodulatory properties. Experimental evidence in mice and observational studies in humans indicate that iNKT cells have antitumor potential as well as the ability to suppress acute and chronic graft-versus-host-disease (GVHD). Murine iNKT cells differentiate during thymic development into iNKT1, iNKT2, and iNKT17 sublineages, which differ transcriptomically and epigenomically and have subset-specific developmental requirements. Whether distinct iNKT sublineages also differ in their antitumor effect and their ability to suppress GVHD is currently unknown. In this work, we generated highly purified murine iNKT sublineages, characterized their transcriptomic and epigenomic landscape, and assessed specific functions. We show that iNKT2 and iNKT17, but not iNKT1, cells efficiently suppress T-cell activation in vitro and mitigate murine acute GVHD in vivo. Conversely, we show that iNKT1 cells display the highest antitumor activity against murine B-cell lymphoma cells both in vitro and in vivo. Thus, we report for the first time that iNKT sublineages have distinct and different functions, with iNKT1 cells having the highest antitumor activity and iNKT2 and iNKT17 cells having immune-regulatory properties. These results have important implications for the translation of iNKT cell therapies to the clinic for cancer immunotherapy as well as for the prevention and treatment of GVHD.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Studying lncRNAs is challenging due to their low expression level, cell type-specific occurrence, poor sequence conservation between orthologs, and lack of information about RNA domains. LncRNAs direct the regulatory factors in the locations that are in cis to their transcription sites. We designed a model to predict if an lncRNA acts in cis based on its features and trained it using RNA-chromatin interaction data. The trained model is cell type-independent and does not require RNA-chromatin data. Combining RNA-chromatin and Hi-C data, we showed that lncRNA-chromatin binding sites are determined by chromosome conformation. For each lncRNA, the spatially proximal genes were identified as their potential targets by combining Hi-C and Cap Analysis Gene Expression (CAGE) data in 18 human cell types. RNA-protein and RNA-chromatin interaction data suggested that lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. We provide the data through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Genetic perturbations of cerebral cortical development can lead to neurodevelopmental disease, including autism spectrum disorder (ASD). To identify genomic regions crucial to corticogenesis, we mapped the activity of gene-regulatory elements generating a single-cell atlas of gene expression and chromatin accessibility both independently and jointly. This revealed waves of gene regulation by key transcription factors (TFs) across a nearly continuous differentiation trajectory into glutamatergic neurons, distinguished the expression programs of glial lineages, and identified lineage-determining TFs that exhibited strong correlation between linked gene-regulatory elements and expression levels. These highly connected genes adopted an active chromatin state in early differentiating cells, consistent with lineage commitment. Basepair-resolution neural network models identified strong cell-type specific enrichment of noncoding mutations predicted to be disruptive in a cohort of ASD subjects and identified frequently disrupted TF binding sites. This approach illustrates how cell-type specific mapping can provide insights into the programs governing human development and disease.

INTRODUCTION: The brains of extant animals have evolved over hundreds of millions of years from simple circuits. Cell types diversified, connections elaborated, and new brain regions emerged. Models for brain region evolution range from duplication of existing regions to splitting of previously multifunctional regions and de novo assembly from existing cell types. These models, however, have not been demonstrated in vertebrate brains at cell-type resolution. RATIONALE: We investigated brain region evolution using the cerebellar nuclei as a model system. The cerebellum is a major hindbrain structure in jawed vertebrates, comprising the cerebellar cortex and cerebellar nuclei. It is thought to act as a feedforward model for motor control and cognitive processes. The cerebellar cortex receives and processes inputs and sends outputs to the cerebellar nuclei, which route the results of cerebellar computations to the rest of the brain. Whereas the cerebellar cortex is well conserved across vertebrates, the cerebellar nuclei vary in number, with none in jawless vertebrates, one pair in cartilaginous fishes and amphibians, two pairs in reptiles and birds, and three pairs in mammals. This pattern suggests that extant cerebellar nuclei evolved from a single ancestral nucleus. Cerebellar nuclei thus provide a good model to interrogate brain region evolution. RESULTS: We characterized the cerebellar nuclei in mice, chickens, and humans using whole-brain and spinal cord projection mapping in cleared samples, single-nucleus RNA sequencing (snRNAseq), and spatially resolved transcript amplicon readout mapping (STARmap) analysis. We first compared the projection patterns of the three cerebellar nuclei of mice. Our data reveal broad projections of all nuclei, which in common target regions are shifted relative to each other. To understand the transcriptomic differences that underlie these shifting projections, we produced a cell-type atlas of the mouse cerebellar nuclei using snRNAseq. We discovered three region-invariant inhibitory cell classes and 15 region-specific excitatory cell types. Excitatory cell types fall into two classes with distinct gene expression and electrophysiological properties. Members of each class are present in every nucleus and are putative sister cell types. STARmap analysis in mice revealed that the organizational unit of the cerebellar nuclei is cytoarchitectonically distinguishable subnuclei, each of which contains the three inhibitory and two excitatory classes. To test whether this archetypal subnucleus is also the evolutionary unit of the cerebellar nuclei, we performed snRNAseq and STARmap on the chicken cerebellar nuclei. We identified four subnuclei, three of which had direct orthologs in mice. Each chicken subnucleus contained the same cell-type set of three inhibitory and two excitatory classes already identified in mice, confirming our hypothesis. Cerebellar nuclei vary in size across vertebrates. In particular, the human lateral nucleus is markedly expanded. To understand this expansion, we performed snRNAseq in humans. We found that the medial and interposed nuclei maintained the archetypal cerebellar nuclei composition. However, the lateral nucleus expanded one excitatory cell class at the expense of the other. Conditional tracing in the mouse lateral nucleus revealed that the cell class expanded in humans preferentially accesses lateral frontal cortices via specific intermediate thalamic nuclei. CONCLUSION: We identified a conserved cell-type set that forms an archetypal cerebellar nucleus as the unit of cerebellar nuclei organization and evolution. We propose that this archetypal nucleus was repeatedly duplicated during evolution, accompanied primarily by transcriptomic divergence of excitatory neurons and shifts in their projection patterns. Our data support a model of duplication-and-divergence of entire cell-type sets for brain region evolution.

ATAC-seq reveals unique neutrophil chromatin architecture changes in response to different stimuli before transcriptional activation, possibly regulating downstream gene expression., Sepsis, sequela of bloodstream infections and dysregulated host responses, is a leading cause of death globally. Neutrophils tightly regulate responses to pathogens to prevent organ damage. Profiling early host epigenetic responses in neutrophils may aid in disease recognition. We performed assay for transposase-accessible chromatin (ATAC)-seq of human neutrophils challenged with six toll-like receptor ligands and two organisms; and RNA-seq after Escherichia coli exposure for 1 and 4 h along with ATAC-seq. ATAC-seq of neutrophils facilitates detection of pathogen DNA. In addition, despite similarities in genomic distribution of differential chromatin changes across challenges, only a fraction overlaps between the challenges. Ligands depict shared signatures, but majority are unique in position, function, and challenge. Epigenomic changes are plastic, only ∼120 are shared by E. coli challenges over time, resulting in varied differential genes and associated processes. We identify three classes of gene regulation, chromatin access changes in the promoter; changes in the promoter and distal enhancers; and controlling expression through changes solely in distal enhancers. These and transcription factor footprinting reveal timely and challenge specific mechanisms of transcriptional regulation in neutrophils.

Extrachromosomal DNAs (ecDNAs) are prevalent in human cancers and mediate high oncogene expression through elevated copy number and altered gene regulation1. Gene expression typically involves distal enhancer DNA elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs, comprised of ~10-100 ecDNAs clustered in the nucleus of interphase cells, drive intermolecular enhancer input for amplified oncogene expression. Single-molecule sequencing, single-cell multiome, and 3D enhancer connectome reveal subspecies of MYC-PVT1 ecDNAs lacking enhancers that access intermolecular and ectopic enhancer-promoter interactions in ecDNA hubs. ecDNA hubs persist without transcription and are tethered by BET protein BRD4. BET inhibitor JQ1 disperses ecDNA hubs, preferentially inhibits ecDNA oncogene transcription, and kills ecDNA+ cancer cells. Two amplified oncogenes MYC and FGFR2 intermix in ecDNA hubs, engage in intermolecular enhancer-promoter interactions, and transcription is uniformly sensitive to JQ1. Thus, ecDNA hubs are nuclear bodies of many ecDNAs tethered by proteins and platforms for cooperative transcription, leveraging the power of oncogene diversification and combinatorial DNA interactions. We suggest ecDNA hubs, rather than individual ecDNAs, as units of oncogene function, cooperative evolution, and new targets for cancer therapy.

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs—clusters of around 10–100 ecDNAs within the nucleus—enable intermolecular enhancer–gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer–gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy. Extrachromosomal DNA (ecDNA) congregates in clusters called ecDNA hubs that promote intermolecular interactions between gene-regulatory regions and thereby amplify the expression of oncogenes such as MYC in cancer cell lines.

CD93 CAR T cells eliminate AML in preclinical models without targeting hematopoietic progenitor cells, and a NOT-gated CAR engineering strategy mitigates on-target, off-tumor toxicity to endothelial cells., Chimeric antigen receptor (CAR) T cells hold promise for the treatment of acute myeloid leukemia (AML), but optimal targets remain to be defined. We demonstrate that CD93 CAR T cells engineered from a novel humanized CD93-specific binder potently kill AML in vitro and in vivo but spare hematopoietic stem and progenitor cells (HSPC). No toxicity is seen in murine models, but CD93 is expressed on human endothelial cells, and CD93 CAR T cells recognize and kill endothelial cell lines. We identify other AML CAR T-cell targets with overlapping expression on endothelial cells, especially in the context of proinflammatory cytokines. To address the challenge of endothelial-specific cross-reactivity, we provide proof of concept for NOT-gated CD93 CAR T cells that circumvent endothelial cell toxicity in a relevant model system. We also identify candidates for combinatorial targeting by profiling the transcriptome of AML and endothelial cells at baseline and after exposure to proinflammatory cytokines. Significance: CD93 CAR T cells eliminate AML and spare HSPCs but exert on-target, off-tumor toxicity to endothelial cells. We show coexpression of other AML targets on endothelial cells, introduce a novel NOT-gated strategy to mitigate endothelial toxicity, and demonstrate use of high-dimensional transcriptomic profiling for rational design of combinatorial immunotherapies. See related commentary by Velasquez and Gottschalk, p. 559. This article is highlighted in the In This Issue feature, p. 549

Transcription factors (TFs) bind DNA sequence motif vocabularies in cis-regulatory elements (CREs) to modulate chromatin state and gene expression during cell state transitions. A quantitative understanding of how motif lexicons influence dynamic regulatory activity has been elusive due to the combinatorial nature of the cis-regulatory code. To address this, we undertook multi-omic data profiling of chromatin and expression dynamics across epidermal differentiation to identify 40,103 dynamic CREs associated with 3,609 dynamically expressed genes, then applied an interpretable deep learning framework to model the cis-regulatory logic of chromatin accessibility. This identified cooperative DNA sequence rules in dynamic CREs regulating synchronous gene modules with diverse roles in skin differentiation. Massively parallel reporter analysis validated temporal dynamics and cooperative cis-regulatory logic. Variants linked to human polygenic skin disease were enriched in these time-dependent combinatorial motif rules. This integrative approach reveals the combinatorial cis-regulatory lexicon of epidermal differentiation and represents a general framework for deciphering the organizational principles of the cis-regulatory code in dynamic gene regulation.HIGHLIGHTSAn integrative multi-omic resource profiling chromatin and expression dynamics across keratinocyte differentiationPredictive deep learning models of chromatin dynamics reveal a high-resolution cis-regulatory DNA motif lexicon of epidermal differentiationModel interpretation enables discovery of combinatorial cis-regulatory logic of homotypic and heterotypic motif combinationsMassively parallel reporter experiments validate temporal dynamics and cis-regulatory logic of the combinatorial motif lexicon

SARS-CoV-2 is the cause of a pandemic with growing global mortality. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with ChIRP-MS data from three other RNA viruses defined viral specificity of RNA-host protein interactions. Targeted CRISPR screens revealed that the majority of functional RNA-binding proteins protect the host from virus-induced cell death, and comparative CRISPR screens across seven RNA viruses revealed shared and SARS-specific antiviral factors. Finally, by combining the RNA-centric approach and functional CRISPR screens, we demonstrated a physical and functional connection between SARS-CoV-2 and mitochondria, highlighting this organelle as a general platform for antiviral activity. Altogether, these data provide a comprehensive catalog of functional SARS-CoV-2 RNA-host protein interactions, which may inform studies to understand the host-virus interface and nominate host pathways that could be targeted for therapeutic benefit., Graphical abstract, Interrogation of SARS-CoV-2 RNA-host protein interaction networks by ChIRP-MS and CRISPR screens, in comparison with other human viruses such as flaviviruses, picornavirus, and rhinovirus, identifies complexes specific to SARS-CoV-2 infection and highlights the role of mitochondria in mediating antiviral immunity.

B cells are critical for the production of antibodies and protective immunity to viruses. Here we show that patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) who develop coronavirus disease 2019 (COVID-19) display early recruitment of B cells expressing a limited subset of IGHV genes, progressing to a highly polyclonal response of B cells with broader IGHV gene usage and extensive class switching to IgG and IgA subclasses with limited somatic hypermutation in the initial weeks of infection. We identify convergence of antibody sequences across SARS-CoV-2-infected patients, highlighting stereotyped naive responses to this virus. Notably, sequence-based detection in COVID-19 patients of convergent B cell clonotypes previously reported in SARS-CoV infection predicts the presence of SARS-CoV/SARS-CoV-2 cross-reactive antibody titers specific for the receptor-binding domain. These findings offer molecular insights into shared features of human B cell responses to SARS-CoV-2 and SARS-CoV., Graphical Abstract, Highlights • Human IGH repertoire sequencing identifies SARS-CoV-2-specific B cell clones • Convergent virus-specific antibody sequences are shared between COVID-19 patients • SARS-CoV antibodies are detected by sequence in COVID-19 patients, B cells produce antibodies and provide protective immunity to viruses. In longitudinal data from COVID-19 patients, Nielsen et al. analyze the development of antibody gene repertoires responding to SARS-CoV-2. COVID-19 patients share subsets of similar antibody sequences that bind SARS-CoV-2 antigens, with rare antibodies that also recognize SARS-CoV.

SUMMARY Human dendritic cells (DCs) comprise subsets with distinct phenotypic and functional characteristics, but the transcriptional programs that dictate their identity remain elusive. Here, we analyze global chromatin accessibility profiles across resting and stimulated human DC subsets by means of the assay for transposase-accessible chromatin using sequencing (ATAC-seq). We uncover specific regions of chromatin accessibility for each subset and transcriptional regulators of DC function. By comparing plasmacytoid DC responses to IFN-I-producing and non-IFN-I-producing conditions, we identify genetic programs related to their function. Finally, by intersecting chromatin accessibility with genome-wide association studies, we recognize DC subset-specific enrichment of heritability in autoimmune diseases. Our results unravel the basis of human DC subset heterogeneity and provide a framework for their analysis in disease pathogenesis., In Brief Human dendritic cells (DCs) orchestrate immune responses by a division of labor between functionally specialized subsets; however, the transcriptional basis of this heterogeneity is poorly understood. Using ATAC-seq, Leylek et al. profile the chromatin landscape of human DC subsets, providing insight into the underlying regulatory mechanisms that modulate their function., Graphical Abstract

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

SUMMARYPolycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction1-5. This process is essential for embryonic stem cell (ESCs) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spread of H3K27me3 on the linear genome6. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remains unclear7-10. Here, we apply H3K27me3 HiChIP11-13, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture14 to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells7-10. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes, such as Hmx1, Wnt6 and Hoxa. Genetic deletion of H3K27me3 loop anchors revealed a coupling of Polycomb-associated genome architecture and H3K27me3 deposition evidenced by disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. Further, we find that global alterations in PRC2 occupancy resulting from an EZH2 mutant15 selectively deficient in RNA binding is accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 acts as a “genomic wormhole”, using RNA binding to enhance long range chromosome folding and H3K27me3 spreading. Additionally, developmental gene loci have novel roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

Physiologic laminar shear stress (LSS) induces an endothelial gene expression profile that is vasculo-protective. In this report, we delineate how LSS mediates changes in the epigenetic landscape to promote this beneficial response. We show that under LSS, KLF4 interacts with the SWI/SNF nucleosome remodeling complex to increase accessibility at enhancer sites that promote expression of homeostatic endothelial genes. By combining molecular and computational approaches we discovered enhancers that loop to promoters of known and novel KLF4- and LSS-responsive genes that stabilize endothelial cells and suppress inflammation, such as BMPR2 and DUSP5. By linking enhancers to genes that they regulate under physiologic LSS, our work establishes a foundation for interpreting how non-coding DNA variants in these regions might disrupt protective gene expression to influence vascular disease.

SUMMARYDuring virus infection B cells are critical for the production of antibodies and protective immunity. Here we show that the human B cell compartment in patients with diagnostically confirmed SARS-CoV-2 and clinical COVID-19 is rapidly altered with the early recruitment of B cells expressing a limited subset of IGHV genes, progressing to a highly polyclonal response of B cells with broader IGHV gene usage and extensive class switching to IgG and IgA subclasses with limited somatic hypermutation in the initial weeks of infection. We identify extensive convergence of antibody sequences across SARS-CoV-2 patients, highlighting stereotyped naïve responses to this virus. Notably, sequence-based detection in COVID-19 patients of convergent B cell clonotypes previously reported in SARS-CoV infection predicts the presence of SARS-CoV/SARS-CoV-2 cross-reactive antibody titers specific for the receptor-binding domain. These findings offer molecular insights into shared features of human B cell responses to SARS-CoV-2 and other zoonotic spillover coronaviruses.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

How have complex brains evolved from simple circuits? Here we investigated brain region evolution at cell type resolution in the cerebellar nuclei (CN), the output structures of the cerebellum. Using single-nucleus RNA sequencing in mice, chickens, and humans, as well as STARmap spatial transcriptomic analysis and whole-CNS projection tracing in mice, we identified a conserved cell type set containing two classes of region-specific excitatory neurons and three classes of region-invariant inhibitory neurons. This set constitutes an archetypal CN that was repeatedly duplicated to form new regions. Interestingly, the excitatory cell class that preferentially funnels information to lateral frontal cortices in mice becomes predominant in the massively expanded human Lateral CN. Our data provide the first characterization of CN transcriptomic cell types in three species and suggest a model of brain region evolution by duplication and divergence of entire cell type sets.