Search

Extrachromosomal DNA (ecDNA) is a common mode of oncogene amplification but is challenging to analyze. Here, we adapt CRISPR-CATCH, in vitro CRISPR-Cas9 treatment and pulsed field gel electrophoresis of agarose-entrapped genomic DNA, previously developed for bacterial chromosome segments, to isolate megabase-sized human ecDNAs. We demonstrate strong enrichment of ecDNA molecules containing EGFR, FGFR2 and MYC from human cancer cells and NRAS ecDNA from human metastatic melanoma with acquired therapeutic resistance. Targeted enrichment of ecDNA versus chromosomal DNA enabled phasing of genetic variants, identified the presence of an EGFRvIII mutation exclusively on ecDNAs and supported an excision model of ecDNA genesis in a glioblastoma model. CRISPR-CATCH followed by nanopore sequencing enabled single-molecule ecDNA methylation profiling and revealed hypomethylation of the EGFR promoter on ecDNAs. We distinguished heterogeneous ecDNA species within the same sample by size and sequence with base-pair resolution and discovered functionally specialized ecDNAs that amplify select enhancers or oncogene-coding sequences.

Summary: Extrachromosomal DNA circles (ecDNA) are a common mechanism for oncogene amplification and are associated with worse clinical outcomes compared with other types of oncogene amplification. Several recent discoveries of ecDNA hubs—local congregations of ecDNAs in the nucleus—highlight unique features of ecDNA biology that may contribute to higher oncogene expression and rapid tumor evolution.

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).

T cells are the major arm of the immune system responsible for controlling and regressing cancers. To identify genes limiting T cell function, we conducted genome-wide CRISPR knockout screens in human chimeric antigen receptor (CAR) T cells. Top hits were MED12 and CCNC , components of the Mediator kinase module. Targeted MED12 deletion enhanced antitumor activity and sustained the effector phenotype in CAR- and T cell receptor–engineered T cells, and inhibition of CDK8/19 kinase activity increased expansion of nonengineered T cells. MED12 -deficient T cells manifested increased core Meditator chromatin occupancy at transcriptionally active enhancers—most notably for STAT and AP-1 transcription factors—and increased IL2RA expression and interleukin-2 sensitivity. These results implicate Mediator in T cell effector programming and identify the kinase module as a target for enhancing potency of antitumor T cell responses.

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.

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture 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 cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding–deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

Mitochondria are essential organelles in eukaryotic cells that provide critical support for energetic and metabolic homeostasis. Although the elimination of pathogenic mitochondrial DNA (mtDNA) mutations in somatic cells has been observed, the mechanisms to maintain proper functions despite their mtDNA mutation load are poorly understood. In this study, we analyzed somatic mtDNA mutations in more than 30,000 single human peripheral and bone marrow mononuclear cells. We observed a significant overrepresentation of homoplasmic mtDNA mutations in B, T, and natural killer (NK) lymphocytes. Intriguingly, their overall mutational burden was lower than that in hematopoietic progenitors and myeloid cells. This characteristic mtDNA mutational landscape indicates a genetic bottleneck during lymphoid development, as confirmed with single-cell datasets from multiple platforms and individuals. We further demonstrated that mtDNA replication lags behind cell proliferation in both pro-B and pre-B progenitor cells, thus likely causing the genetic bottleneck by diluting mtDNA copies per cell. Through computational simulations and approximate Bayesian computation (ABC), we recapitulated this lymphocyte-specific mutational landscape and estimated the minimal mtDNA copies as

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.

Advances in single-cell RNA sequencing (scRNA-seq) technologies in the past 10 years have had a transformative effect on biomedical research, enabling the profiling and analysis of the transcriptomes of single cells at unprecedented resolution and throughput. Specifically, scRNA-seq has facilitated the identification of novel or rare cell types, the analysis of single-cell trajectory construction and stem or progenitor cell differentiation, and the comparison of healthy and disease-related tissues at single-cell resolution. These applications have been critical in advances in cardiovascular research in the past decade as evidenced by the generation of cell atlases of mammalian heart and blood vessels and the elucidation of mechanisms involved in cardiovascular development and stem or progenitor cell differentiation. In this Review, we summarize the currently available scRNA-seq technologies and analytical tools and discuss the latest findings using scRNA-seq that have substantially improved our knowledge on the development of the cardiovascular system and the mechanisms underlying cardiovascular diseases. Furthermore, we examine emerging strategies that integrate multimodal single-cell platforms, focusing on future applications in cardiovascular precision medicine that use single-cell omics approaches to characterize cell-specific responses to drugs or environmental stimuli and to develop effective patient-specific therapeutics.

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.

BackgroundSystemic sclerosis (SSc) is characterised by autoimmune activation, tissue and vascular fibrosis in the skin and internal organs. Tissue fibrosis is driven by myofibroblasts, that are known to maintain their phenotype in vitro, which is associated with epigenetically driven trimethylation of lysine 27 of histone 3 (H3K27me3).MethodsFull-thickness skin biopsies were surgically obtained from the forearms of 12 adult patients with SSc of recent onset. Fibroblasts were isolated and cultured in monolayers and protein and RNA extracted. HOX transcript antisense RNA (HOTAIR) was expressed in healthy dermal fibroblasts by lentiviral induction employing a vector containing the specific sequence. Gamma secretase inhibitors were employed to block Notch signalling. Enhancer of zeste 2 (EZH2) was blocked with GSK126 inhibitor.ResultsSSc myofibroblasts in vitro and SSc skin biopsies in vivo display high levels of HOTAIR, a scaffold long non-coding RNA known to direct the histone methyltransferase EZH2 to induce H3K27me3 in specific target genes. Overexpression of HOTAIR in dermal fibroblasts induced EZH2-dependent increase in collagen and α-SMA expression in vitro, as well as repression of miRNA-34A expression and consequent NOTCH pathway activation. Consistent with these findings, we show that SSc dermal fibroblast display decreased levels of miRNA-34a in vitro. Further, EZH2 inhibition rescued miRNA-34a levels and mitigated the profibrotic phenotype of both SSc and HOTAIR overexpressing fibroblasts in vitro.ConclusionsOur data indicate that the EZH2-dependent epigenetic phenotype of myofibroblasts is driven by HOTAIR and is linked to miRNA-34a repression-dependent activation of NOTCH signalling.

Oncogene amplification on extrachromosomal DNA (ecDNA) is prevalent in human cancer and is associated with poor outcomes. Clonal, megabase-sized circular ecDNAs in cancer are distinct from nonclonal, small sub-kilobase-sized DNAs that may arise during normal tissue homeostasis. ecDNAs enable profound changes in gene regulation beyond copy-number gains. An emerging principle of ecDNA regulation is the formation of ecDNA hubs: micrometer-sized nuclear structures of numerous copies of ecDNAs tethered by proteins in spatial proximity. ecDNA hubs enable cooperative and intermolecular sharing of DNA regulatory elements for potent and combinatorial gene activation. The 3D context of ecDNA shapes its gene expression potential, selection for clonal heterogeneity among ecDNAs, distribution through cell division, and reintegration into chromosomes. Technologies for studying gene regulation and structure of ecDNA are starting to answer long-held questions on the distinct rules that govern cancer genes beyond chromosomes.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a common event, driving aggressive tumor growth, drug resistance and shorter survival. Currently, the impact of nonchromosomal oncogene inheritance—random identity by descent—is poorly understood. Also unclear is the impact of ecDNA on somatic variation and selection. Here integrating theoretical models of random segregation, unbiased image analysis, CRISPR-based ecDNA tagging with live-cell imaging and CRISPR-C, we demonstrate that random ecDNA inheritance results in extensive intratumoral ecDNA copy number heterogeneity and rapid adaptation to metabolic stress and targeted treatment. Observed ecDNAs benefit host cell survival or growth and can change within a single cell cycle. ecDNA inheritance can predict, a priori, some of the aggressive features of ecDNA-containing cancers. These properties are facilitated by the ability of ecDNA to rapidly adapt genomes in a way that is not possible through chromosomal oncogene amplification. These results show how the nonchromosomal random inheritance pattern of ecDNA contributes to poor outcomes for patients with cancer.

Chronic antigen exposure during viral infection or cancer promotes an exhausted T cell (Tex) state with reduced effector function. However, whether all antigen-specific T cell clones follow the same Tex differentiation trajectory remains unclear. Here, we generate a single-cell multiomic atlas of T cell exhaustion in murine chronic viral infection that redefines Tex phenotypic diversity, including two late-stage Tex subsets with either a terminal exhaustion (Tex

T cells mediate anti-tumor immune responses and are the key target of immune checkpoint therapy, but they can also promote immune tolerance. A clear understanding of the specific contributions and biology of different T cell subsets is required to fully harness the curative potential of immunotherapies. Experts discuss the state of the field and key challenges for moving forward.

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.

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.

Summary Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) structure probing techniques characterize the secondary structure of RNA molecules, which influence their functions and interactions. A variation of SHAPE, footprinting SHAPE (fSHAPE), probes RNA in the presence and absence of protein to identify RNA bases that hydrogen-bond with protein. SHAPE or fSHAPE coupled with enhanced crosslinking and immunoprecipitation (SHAPE-eCLIP or fSHAPE-eCLIP) pulls down RNAs bound by any protein of interest and returns their structure or protein interaction information, respectively. Here, we describe detailed protocols for SHAPE-eCLIP and fSHAPE-eCLIP and an analysis protocol for fSHAPE. For complete details on the use and execution of these protocols, please refer to Corley et al. (2020)., Graphical abstract, Highlights • Protocols for probing protein-RNA hydrogen bonds (fSHAPE) and RNA structure (SHAPE) • fSHAPE-eCLIP and SHAPE-eCLIP probe RNA regions bound by immunoprecipitated protein • Samples probed in parallel under different conditions is the basis of each protocol • Probed regions yield fSHAPE or SHAPE reactivities at nucleotide resolution, Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) structure probing techniques characterize the secondary structure of RNA molecules, which influence their functions and interactions. A variation of SHAPE, footprinting SHAPE (fSHAPE), probes RNA in the presence and absence of protein to identify RNA bases that hydrogen-bond with protein. SHAPE or fSHAPE coupled with enhanced crosslinking and immunoprecipitation (SHAPE-eCLIP or fSHAPE-eCLIP) pulls down RNAs bound by any protein of interest and returns their structure or protein interaction information, respectively. Here, we describe detailed protocols for SHAPE-eCLIP and fSHAPE-eCLIP and an analysis protocol for fSHAPE.

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.

Chimeric antigen receptor (CAR) T cells mediate anti-tumour effects in a small subset of patients with cancer1-3, but dysfunction due to T cell exhaustion is an important barrier to progress4-6. To investigate the biology of exhaustion in human T cells expressing CAR receptors, we used a model system with a tonically signaling CAR, which induces hallmark features of exhaustion6. Exhaustion was associated with a profound defect in the production of IL-2, along with increased chromatin accessibility of AP-1 transcription factor motifs and overexpression of the bZIP and IRF transcription factors that have been implicated in mediating dysfunction in exhausted T cells7-10. Here we show that CAR T cells engineered to overexpress the canonical AP-1 factor c-Jun have enhanced expansion potential, increased functional capacity, diminished terminal differentiation and improved anti-tumour potency in five different mouse tumour models in vivo. We conclude that a functional deficiency in c-Jun mediates dysfunction in exhausted human T cells, and that engineering CAR T cells to overexpress c-Jun renders them resistant to exhaustion, thereby addressing a major barrier to progress for this emerging class of therapeutic agents.

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.

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.

A systematic approach to understanding the noncoding genome in cancer promises to improve cancer diagnosis and therapy. Tracking the landscape of gene control in individual patients and single cells yields many insights for inherited and somatic mutations in cancer, extrachromosomal oncogene amplifications and cancer immunotherapy. New technologies and bold therapeutic approaches are paving the way to truly envisage personalized cancer medicine in the future. A systematic approach to understanding the noncoding genome in cancer promises to improve cancer diagnosis and therapy. New technologies and bold therapeutic approaches are paving the way to truly envisage personalized cancer medicine in the future.

Significance Simultaneous measurement of the DNA, RNA, and proteins of single cells can lead to important new insights but is experimentally challenging. This work introduces a deep learning algorithm that flexibly translates between chromatin, RNA, and protein profiles in single cells. This makes it possible to computationally synthesize matched multiomic measurements when only one modality is experimentally available. This algorithm complements experimental advances to efficiently achieve single-cell multiomic insight., Simultaneous profiling of multiomic modalities within a single cell is a grand challenge for single-cell biology. While there have been impressive technical innovations demonstrating feasibility—for example, generating paired measurements of single-cell transcriptome (single-cell RNA sequencing [scRNA-seq]) and chromatin accessibility (single-cell assay for transposase-accessible chromatin using sequencing [scATAC-seq])—widespread application of joint profiling is challenging due to its experimental complexity, noise, and cost. Here, we introduce BABEL, a deep learning method that translates between the transcriptome and chromatin profiles of a single cell. Leveraging an interoperable neural network model, BABEL can predict single-cell expression directly from a cell’s scATAC-seq and vice versa after training on relevant data. This makes it possible to computationally synthesize paired multiomic measurements when only one modality is experimentally available. Across several paired single-cell ATAC and gene expression datasets in human and mouse, we validate that BABEL accurately translates between these modalities for individual cells. BABEL also generalizes well to cell types within new biological contexts not seen during training. Starting from scATAC-seq of patient-derived basal cell carcinoma (BCC), BABEL generated single-cell expression that enabled fine-grained classification of complex cell states, despite having never seen BCC data. These predictions are comparable to analyses of experimental BCC scRNA-seq data for diverse cell types related to BABEL’s training data. We further show that BABEL can incorporate additional single-cell data modalities, such as protein epitope profiling, thus enabling translation across chromatin, RNA, and protein. BABEL offers a powerful approach for data exploration and hypothesis generation.

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.

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).

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 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

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.

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with 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 tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 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 change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

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.

We have known for decades that long noncoding RNAs (lncRNAs) can play essential functions across most forms of life. The maintenance of chromosome length requires an lncRNA (e.g., hTERC) and two lncRNAs in the ribosome that are required for protein synthesis. Thus, lncRNAs can represent powerful RNA machines. More recently, it has become clear that mammalian genomes encode thousands more lncRNAs. Thus, we raise the question: Which, if any, of these lncRNAs could also represent RNA-based machines? Here we synthesize studies that are beginning to address this question by investigating fundamental properties of lncRNA genes, revealing new insights into the RNA structure–function relationship, determining cis- and trans-acting lncRNAs in vivo, and generating new developments in high-throughput screening used to identify functional lncRNAs. Overall, these findings provide a context toward understanding the molecular grammar underlying lncRNA biology.

Summary Discovering the interaction mechanism and location of RNA-binding proteins (RBPs) on RNA is critical for understanding gene expression regulation. Here, we apply selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) on in vivo transcripts compared to protein-absent transcripts in four human cell lines to identify transcriptome-wide footprints (fSHAPE) on RNA. Structural analyses indicate that fSHAPE precisely detects nucleobases that hydrogen bond with protein. We demonstrate that fSHAPE patterns predict binding sites of known RBPs, such as iron response elements in both known loci and previously unknown loci in CDC34, SLC2A4RG, COASY, and H19. Furthermore, by integrating SHAPE and fSHAPE with crosslinking and immunoprecipitation (eCLIP) of desired RBPs, we interrogate specific RNA-protein complexes, such as histone stem-loop elements and their nucleotides that hydrogen bond with stem-loop-binding proteins. Together, these technologies greatly expand our ability to study and understand specific cellular RNA interactions in RNA-protein complexes.

/Summary SARS-CoV-2 genomic and subgenomic RNA (sgRNA) transcripts hijack the host cell's machinery. Subcellular localization of its viral RNA could thus play important roles in viral replication and host antiviral immune response. We perform computational modeling of SARS-CoV-2 viral RNA subcellular residency across eight subcellular neighborhoods. We compare hundreds of SARS-CoV-2 genomes to the human transcriptome and other coronaviruses. We predict the SARS-CoV-2 RNA genome and sgRNAs to be enriched towards the host mitochondrial matrix and nucleolus, and that the 5’ and 3’ viral untranslated regions contain the strongest, most distinct localization signals. We interpret the mitochondrial residency signal as an indicator of intracellular RNA trafficking with respect to double-membrane vesicles, a critical stage in the coronavirus life cycle. Our computational analysis serves as a hypothesis generation tool to suggest models for SARS-CoV-2 biology and inform experimental efforts to combat the virus. A record of this paper’s Transparent Peer Review process is included in the Supplemental Information., Graphical Abstract, Highlights ● Application of machine learning model of RNA subcellular localization to SARS-CoV-2 ● Viral RNAs show residency signal for host mitochondria and nucleolus ● Mitochondria prediction suggests viruses repurpose endogenous pathways ● Predictions may be linked to vesicle formation and viral-host protein interactions, Where the SARS-CoV-2 genome localizes inside human cells remains understudied but may regulate viral replication and host response. We use a machine learning model to predict subcellular residency of the SARS-CoV-2 genome and its encoded transcripts, as well as for other coronaviruses. Our predictions suggest new hypotheses for SARS-CoV-2 mechanisms.

Despite the fact that Otto H. Warburg discovered the Warburg effect almost one hundred years ago, why cancer cells waste most of the glucose carbon as lactate remains an enigma. Warburg proposed a connection between the Warburg effect and cell dedifferentiation. Hypoxia is a common tumor microenvironmental stress that induces the Warburg effect and blocks tumor cell differentiation. The underlying mechanism by which this occurs is poorly understood, and no effective therapeutic strategy has been developed to overcome this resistance to differentiation. Using a neuroblastoma differentiation model, we discovered that hypoxia repressed cell differentiation through reducing cellular acetyl-CoA levels, leading to reduction of global histone acetylation and chromatin accessibility. The metabolic switch triggering this global histone hypoacetylation was the induction of pyruvate dehydrogenase kinases (PDK1 and PDK3). Inhibition of PDKs using dichloroacetate (DCA) restored acetyl-CoA generation and histone acetylation under hypoxia. Knocking down PDK1 induced neuroblastoma cell differentiation, highlighting the critical role of PDK1 in cell fate control. Importantly, acetate or glycerol triacetate (GTA) supplementation restored differentiation markers expression and neuron differentiation under hypoxia. Moreover, ATAC-Seq analysis demonstrated that hypoxia treatment significantly reduced chromatin accessibility at RAR/RXR binding sites, which can be restored by acetate supplementation. In addition, hypoxia-induced histone hypermethylation by increasing 2-hydroxyglutarate (2HG) and reducing α-ketoglutarate (αKG). αKG supplementation reduced histone hypermethylation upon hypoxia, but did not restore histone acetylation or differentiation markers expression. Together, these findings suggest that diverting pyruvate flux away from acetyl-CoA generation to lactate production is the key mechanism that Warburg effect drives dedifferentiation and tumorigenesis. We propose that combining differentiation therapy with acetate/GTA supplementation might represent an effective therapy against neuroblastoma.

INTRODUCTION: The cerebral cortex is responsible for higher-order functions in the nervous system and has undergone substantial expansion in size in primates. The development of the forebrain, including the assembly of the expanded human cerebral cortex, is a lengthy process that involves the diversification and expansion of neural progenitors, the generation and positioning of layer-specific glutamatergic neurons, the cellular migration of γ-aminobutyric acid (GABA)-ergic neurons, and the formation and maturation of glial cells. Disruption of these cellular events by either genetic or environmental factors can lead to neurodevelopmental disease, including autism spectrum disorders and intellectual disability. RATIONALE: Human forebrain development is, to a large extent, inaccessible for cellular-level study, direct functional investigation, or manipulation. The lack of availability of primary brain tissue samples—in particular, at later stages—as well as the limitations of conventional in vitro cellular models have precluded a detailed mechanistic understanding of corticogenesis in healthy and disease states. Therefore, tracking epigenetic changes in specific forebrain cell lineages over long time periods, has the potential to unravel the molecular programs that underlie cell specification in the human cerebral cortex and, by temporally mapping disease risk onto these changes, to identify cell types and periods of increased disease susceptibility. RESULTS: We used three-dimensional (3D) directed differentiation of human pluripotent stem cells into dorsal and ventral forebrain domains and applied the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) in combination with RNA-sequencing (RNA-seq) to map the epigenetic and gene expression signatures of neuronal and glial cell lineages over 20 months in vitro. We show, through direct comparison with primary brain tissue from our study and several epigenetic datasets, that human stem cell-derived 3D forebrain organoids recapitulated in vivo chromatin accessibility patterns over time. We then integrated these data to discover putative enhancer-gene linkages and lineage-specific transcription factor regulators, including a diverse repertoire of factors that may control cortical specification. We validated protein expression of some of these transcription factors using immunofluorescence, confirming cellular and temporal dynamics in both primary tissue and forebrain organoids. Next, we used this resource to map genes and genetic variants associated with schizophrenia and autism spectrum disorders to distinct accessibility patterns to reveal cell types and periods of susceptibility. Last, we identified a wave of chromatin remodeling during cortical neurogenesis, during which a quarter of regulatory regions are active, and then highlighted transcription factors that may drive these developmental changes. CONCLUSION: Using long-term 3D neural differentiation of stem cells as well as primary brain tissue samples, we found that organoids intrinsically undergo chromatin state transitions in vitro that are closely related to human forebrain development in vivo. Leveraging this platform, we identified epigenetic alterations putatively driven by specific transcription factors and discovered a dynamic period of chromatin remodeling during human cortical neurogenesis. Finally, we nominated several key transcription factors that may coordinate over time to drive these changes and mapped cell type-specific disease-associated variation over time and in specific cell types.

CRISPR takes first steps in humans CRISPR-Cas9 is a revolutionary gene-editing technology that offers the potential to treat diseases such as cancer, but the effects of CRISPR in patients are currently unknown. Stadtmauer et al. report a phase 1 clinical trial to assess the safety and feasibility of CRISPR-Cas9 gene editing in three patients with advanced cancer (see the Perspective by Hamilton and Doudna). They removed immune cells called T lymphocytes from patients and used CRISPR-Cas9 to disrupt three genes ( TRAC, TRBC , and PDCD1 ) with the goal of improving antitumor immunity. A cancer-targeting transgene, NY-ESO-1, was also introduced to recognize tumors. The engineered cells were administered to patients and were well tolerated, with durable engraftment observed for the study duration. These encouraging observations pave the way for future trials to study CRISPR-engineered cancer immunotherapies. Science , this issue p. eaba7365 ; see also p. 976

Pediatric cancers often mimic fetal tissues and express proteins normally silenced postnatally that could serve as immune targets. We developed T cells expressing chimeric antigen receptors (CARs) targeting glypican-2 (GPC2), a fetal antigen expressed on neuroblastoma (NB) and several other solid tumors. CARs engineered using standard designs control NBs with transgenic GPC2 overexpression, but not those expressing clinically relevant GPC2 site density (∼5,000 molecules/cell, range 1-6 × 10

Human embryonic stem cell (hESC) differentiation promises advances in regenerative medicine1–3, yet conversion of hESCs into transplantable cells or tissues remains poorly understood. Using our keratinocyte differentiation system, we employ a multi-dimensional genomics approach to interrogate the contributions of inductive morphogens retinoic acid and bone morphogenetic protein 4 (BMP4) and the epidermal master regulator p63 (encoded by TP63)4,5 during surface ectoderm commitment. In contrast to other master regulators6–9, p63 effects major transcriptional changes only after morphogens alter chromatin accessibility, establishing an epigenetic landscape for p63 to modify. p63 distally closes chromatin accessibility and promotes accumulation of H3K27me3 (trimethylated histone H3 lysine 27). Cohesin HiChIP10 visualizations of chromosome conformation show that p63 and the morphogens contribute to dynamic long-range chromatin interactions, as illustrated by TFAP2C regulation11. Our study demonstrates the unexpected dependency of p63 on morphogenetic signaling and provides novel insights into how a master regulator can specify diverse transcriptional programs based on the chromatin landscape induced by exposure to specific morphogens.

Significance Biological samples are often heterogeneous mixtures of different types of cells. Suppose we have two single-cell datasets, each providing information on a different cellular feature and generated on a different sample from this mixture. Then, the clustering of cells in the two samples should be coupled as both clusterings are reflecting the underlying cell types in the same mixture. This “coupled clustering” problem is a new problem not covered by existing clustering methods. In this paper, we develop an approach for its solution based on the coupling of two nonnegative matrix factorizations. The method should be useful for integrative single-cell genomics analysis tasks such as the joint analysis of single-cell RNA-sequencing and single-cell ATAC-sequencing data., When different types of functional genomics data are generated on single cells from different samples of cells from the same heterogeneous population, the clustering of cells in the different samples should be coupled. We formulate this “coupled clustering” problem as an optimization problem and propose the method of coupled nonnegative matrix factorizations (coupled NMF) for its solution. The method is illustrated by the integrative analysis of single-cell RNA-sequencing (RNA-seq) and single-cell ATAC-sequencing (ATAC-seq) data.

Many craniofacial disorders are caused by heterozygous mutations in general regulators of housekeeping cellular functions such as transcription or ribosome biogenesis1,2. Although it is understood that many of these malformations are a consequence of defects in cranial neural crest cells, a cell type that gives rise to most of the facial structures during embryogenesis3,4, the mechanism underlying cell-type selectivity of these defects remains largely unknown. By exploring molecular functions of DDX21, a DEAD-box RNA helicase involved in control of both RNA polymerase (Pol) I- and II-dependent transcriptional arms of ribosome biogenesis5, we uncovered a previously unappreciated mechanism linking nucleolar dysfunction, ribosomal DNA (rDNA) damage, and craniofacial malformations. Here we demonstrate that genetic perturbations associated with Treacher Collins syndrome, a craniofacial disorder caused by heterozygous mutations in components of the Pol I transcriptional machinery or its cofactor TCOF1 (ref. 1), lead to relocalization of DDX21 from the nucleolus to the nucleoplasm, its loss from the chromatin targets, as well as inhibition of rRNA processing and downregulation of ribosomal protein gene transcription. These effects are cell-type-selective, cell-autonomous, and involve activation of p53 tumour-suppressor protein. We further show that cranial neural crest cells are sensitized to p53-mediated apoptosis, but blocking DDX21 loss from the nucleolus and chromatin rescues both the susceptibility to apoptosis and the craniofacial phenotypes associated with Treacher Collins syndrome. This mechanism is not restricted to cranial neural crest cells, as blood formation is also hypersensitive to loss of DDX21 functions. Accordingly, ribosomal gene perturbations associated with Diamond-Blackfan anaemia disrupt DDX21 localization. At the molecular level, we demonstrate that impaired rRNA synthesis elicits a DNA damage response, and that rDNA damage results in tissue-selective and dosage-dependent effects on craniofacial development. Taken together, our findings illustrate how disruption in general regulators that compromise nucleolar homeostasis can result in tissue-selective malformations.