Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells

Lactylation of histone by BRD4 regulates astrocyte polarization after experimental subarachnoid hemorrhage

Under subarachnoid hemorrhage (SAH) conditions, astrocytes undergo a marked intensification of glycolytic activity, resulting in the generation of substantial amounts of lactate to maintain the energy demand for neurons and other brain cells. Lactate has garnered increasing attention in recent years because of its emerging role in critical biological processes such as inflammation regulation and neuroprotection, particularly through its histone lactylation. Bromodomain-containing protein 4 (BRD4) plays a crucial role in maintaining neural development and promoting memory formation in the central nervous system. Nonetheless, the function and regulatory mechanism of BRD4 and histone lactylation in astrocytes following SAH remain elusive. Our findings indicate that BRD4, a crucial epigenetic regulator, plays a definitive role in histone lactylation. Both in vitro and in vivo, these results demonstrated that targeted silencing of BRD4 in astrocytes can significantly reduce H4K8la lactylation, thereby aggravating the A1 polarization of astrocytes and ultimately affecting the recovery of neural function and prognosis in mice after SAH. In summary, BRD4 plays a pivotal role in modulating astrocyte polarization following SAH via histone lactylation. Targeting this mechanism might offer an efficient therapeutic strategy for SAH.

Heterologous survey of 130 DNA transposons in human cells highlights their functional divergence and expands the genome engineering toolbox

Experimental studies on DNA transposable elements (TEs) have been limited in scale, leading to a lack of understanding of the factors influencing transposition activity, evolutionary dynamics, and application potential as genome engineering tools. We predicted 130 active DNA TEs from 102 metazoan genomes and evaluated their activity in human cells. We identified 40 active (integration-competent) TEs, surpassing the cumulative number (20) of TEs found previously. With this unified comparative data, we found that the Tc1/mariner superfamily exhibits elevated activity, potentially explaining their pervasive horizontal transfers. Further functional characterization of TEs revealed additional divergence in features such as insertion bias. Remarkably, in CAR-T therapy for hematological and solid tumors, Mariner2_AG (MAG), the most active DNA TE identified, largely outperformed two widely used vectors, the lentiviral vector and the TE-based vector SB100X. Overall, this study highlights the varied transposition features and evolutionary dynamics of DNA TEs and increases the TE toolbox diversity.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

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

Genomic Mosaicism of the Brain: Origin, Impact, and Utility

Genomic mosaicism describes the phenomenon where some but not all cells within a tissue harbor unique genetic mutations. Traditionally, research focused on the impact of genomic mosaicism on clinical phenotype-motivated by its involvement in cancers and overgrowth syndromes. More recently, we increasingly shifted towards the plethora of neutral mosaic variants that can act as recorders of cellular lineage and environmental exposures. Here, we summarize the current state of the field of genomic mosaicism research with a special emphasis on our current understanding of this phenomenon in brain development and homeostasis. Although the field of genomic mosaicism has a rich history, technological advances in the last decade have changed our approaches and greatly improved our knowledge. We will provide current definitions and an overview of contemporary detection approaches for genomic mosaicism. Finally, we will discuss the impact and utility of genomic mosaicism.

Effect of genomic and cellular environments on gene expression noise

BACKGROUND

Individual cells from isogenic populations often display large cell-to-cell differences in gene expression. This "noise" in expression derives from several sources, including the genomic and cellular environment in which a gene resides. Large-scale maps of genomic environments have revealed the effects of epigenetic modifications and transcription factor occupancy on mean expression levels, but leveraging such maps to explain expression noise will require new methods to assay how expression noise changes at locations across the genome.

RESULTS

To address this gap, we present Single-cell Analysis of Reporter Gene Expression Noise and Transcriptome (SARGENT), a method that simultaneously measures the noisiness of reporter genes integrated throughout the genome and the global mRNA profiles of individual reporter-gene-containing cells. Using SARGENT, we perform the first comprehensive genome-wide survey of how genomic locations impact gene expression noise. We find that the mean and noise of expression correlate with different histone modifications. We quantify the intrinsic and extrinsic components of reporter gene noise and, using the associated mRNA profiles, assign the extrinsic component to differences between the CD24+ "stem-like" substate and the more "differentiated" substate. SARGENT also reveals the effects of transgene integrations on endogenous gene expression, which will help guide the search for "safe-harbor" loci.

CONCLUSIONS

Taken together, we show that SARGENT is a powerful tool to measure both the mean and noise of gene expression at locations across the genome and that the data generatd by SARGENT reveals important insights into the regulation of gene expression noise genome-wide.

Human ASXL1-Mutant Hematopoiesis Is Driven by a Truncated Protein Associated with Aberrant Deubiquitination of H2AK119

Mutations in additional sex combs like 1 (ASXL1) confer poor prognosis both in myeloid malignancies and in premalignant clonal hematopoiesis (CH). However, the mechanisms by which these mutations contribute to disease initiation remain unresolved, and mutation-specific targeting has remained elusive. To address this, we developed a human disease model that recapitulates the disease trajectory from ASXL1-mutant CH to lethal myeloid malignancy. We demonstrate that mutations in ASXL1 lead to the expression of a functional, truncated protein and determine that truncated ASXL1 leads to global redistribution of the repressive chromatin mark H2AK119Ub, increased transposase-accessible chromatin, and activation of both myeloid and stem cell gene-expression programs. Finally, we demonstrate that H2AK119Ub levels are tied to truncated ASXL1 expression levels and leverage this observation to demonstrate that inhibition of the PRC1 complex might be an ASXL1-mutant-specific therapeutic vulnerability in both premalignant CH and myeloid malignancy.

SIGNIFICANCE

Mutant ASXL1 is a common driver of CH and myeloid malignancy. Using primary human HSPCs, we determine that truncated ASXL1 leads to redistribution of H2AK119Ub and may affect therapeutic vulnerability to PRC1 inhibition.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Pycallingcards: an integrated environment for visualizing, analyzing, and interpreting Calling Cards data

MOTIVATION

Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data.

RESULTS

We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data.

AVAILABILITY AND IMPLEMENTATION

Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards.

A Cre-dependent massively parallel reporter assay allows for cell-type specific assessment of the functional effects of non-coding elements in vivo

The function of regulatory elements is highly dependent on the cellular context, and thus for understanding the function of elements associated with psychiatric diseases these would ideally be studied in neurons in a living brain. Massively Parallel Reporter Assays (MPRAs) are molecular genetic tools that enable functional screening of hundreds of predefined sequences in a single experiment. These assays have not yet been adapted to query specific cell types in vivo in a complex tissue like the mouse brain. Here, using a test-case 3'UTR MPRA library with genomic elements containing variants from autism patients, we developed a method to achieve reproducible measurements of element effects in vivo in a cell type-specific manner, using excitatory cortical neurons and striatal medium spiny neurons as test cases. This targeted technique should enable robust, functional annotation of genetic elements in the cellular contexts most relevant to psychiatric disease.

Massively parallel in vivo Perturb-seq reveals cell type-specific transcriptional networks in cortical development

Systematic analysis of gene function across diverse cell types in vivo is hindered by two challenges: obtaining sufficient cells from live tissues and accurately identifying each cell's perturbation in high-throughput single-cell assays. Leveraging AAV's versatile cell type tropism and high labeling capacity, we expanded the resolution and scale of in vivo CRISPR screens: allowing phenotypic analysis at single-cell resolution across a multitude of cell types in the embryonic brain, adult brain, and peripheral nervous system. We undertook extensive tests of 86 AAV serotypes, combined with a transposon system, to substantially amplify labeling and accelerate in vivo gene delivery from weeks to days. Using this platform, we performed an in utero genetic screen as proof-of-principle and identified pleiotropic regulatory networks of Foxg1 in cortical development, including Layer 6 corticothalamic neurons where it tightly controls distinct networks essential for cell fate specification. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% (mediated by lentivirus), and achieve analysis of over 30,000 cells in one experiment, thus enabling massively parallel in vivo Perturb-seq. Compatible with various perturbation techniques (CRISPRa/i) and phenotypic measurements (single-cell or spatial multi-omics), our platform presents a flexible, modular approach to interrogate gene function across diverse cell types in vivo, connecting gene variants to their causal functions.

Calling Cards: A Customizable Platform to Longitudinally Record Protein-DNA Interactions Over Time in Cells and Tissues

Calling Cards is a platform technology to record a cumulative history of transient protein-DNA interactions in the genome of genetically targeted cell types. The record of these interactions is recovered by next-generation sequencing. Compared with other genomic assays, readouts of which provide a snapshot at the time of harvest, Calling Cards enables correlation of historical molecular states to eventual outcomes or phenotypes. To achieve this, Calling Cards uses the piggyBac transposase to insert self-reporting transposon "Calling Cards" into the genome, leaving permanent marks at interaction sites. Calling Cards can be deployed in a variety of in vitro and in vivo biological systems to study gene regulatory networks involved in development, aging, and disease. Out of the box, it assesses enhancer usage but can be adapted to profile-specific transcription factor (TF) binding with custom TF-piggyBac fusion proteins. The Calling Cards workflow has five main stages: delivery of Calling Cards reagents, sample preparation, library preparation, sequencing, and data analysis. Here, we first present a comprehensive guide for experimental design, reagent selection, and optional customization of the platform to study additional TFs. Then, we provide an updated protocol for the five steps, using reagents that improve throughput and decrease costs, including an overview of a newly deployed computational pipeline. This protocol is designed for users with basic molecular biology experience to process samples into sequencing libraries in 2 days. Familiarity with bioinformatic analysis and command line tools is required to set up the pipeline in a high-performance computing environment and to conduct downstream analyses. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Preparation and delivery of Calling Cards reagents Support Protocol 1: Next-generation sequencing quantification of barcode distribution within self-reporting transposon plasmid pool and adeno-associated virus genome Basic Protocol 2: Sample collection and RNA purification Support Protocol 2: Library density quantitative PCR Basic Protocol 3: Sequencing library preparation Basic Protocol 4: Library pooling and sequencing Basic Protocol 5: Data analysis.

The use of RNA-based treatments in the field of cancer immunotherapy

Over the past several decades, mRNA vaccines have evolved from a theoretical concept to a clinical reality. These vaccines offer several advantages over traditional vaccine techniques, including their high potency, rapid development, low-cost manufacturing, and safe administration. However, until recently, concerns over the instability and inefficient distribution of mRNA in vivo have limited their utility. Fortunately, recent technological advancements have mostly resolved these concerns, resulting in the development of numerous mRNA vaccination platforms for infectious diseases and various types of cancer. These platforms have shown promising outcomes in both animal models and humans. This study highlights the potential of mRNA vaccines as a promising alternative approach to conventional vaccine techniques and cancer treatment. This review article aims to provide a thorough and detailed examination of mRNA vaccines, including their mechanisms of action and potential applications in cancer immunotherapy. Additionally, the article will analyze the current state of mRNA vaccine technology and highlight future directions for the development and implementation of this promising vaccine platform as a mainstream therapeutic option. The review will also discuss potential challenges and limitations of mRNA vaccines, such as their stability and in vivo distribution, and suggest ways to overcome these issues. By providing a comprehensive overview and critical analysis of mRNA vaccines, this review aims to contribute to the advancement of this innovative approach to cancer treatment.

How chromosomal translocations arise to cause cancer: Gene proximity, trans-splicing, and DNA end joining

Chromosomal translocations (CTs) are a genetic hallmark of cancer. They could be identified as recurrent genetic aberrations in hemato-malignancies and solid tumors. More than 40% of all "cancer genes" were identified in recurrent CTs. Most of these CTs result in the production of oncofusion proteins of which many have been studied over the past decades. They influence signaling pathways and/or alter gene expression. However, a precise mechanism for how these CTs arise and occur in a nearly identical fashion in individuals remains to be elucidated. Here, we performed experiments that explain the onset of CTs: (1) proximity of genes able to produce prematurely terminated transcripts, which lead to the production of (2) trans-spliced fusion RNAs, and finally, the induction of (3) DNA double-strand breaks which are subsequently repaired via EJ repair pathways. Under these conditions, balanced chromosomal translocations could be specifically induced. The implications of these findings will be discussed.

Zinc cluster transcription factors frequently activate target genes using a non-canonical half-site binding mode

Gene expression changes are orchestrated by transcription factors (TFs), which bind to DNA to regulate gene expression. It remains surprisingly difficult to predict basic features of the transcriptional process, including in vivo TF occupancy. Existing thermodynamic models of TF function are often not concordant with experimental measurements, suggesting undiscovered biology. Here, we analyzed one of the most well-studied TFs, the yeast zinc cluster Gal4, constructed a Shea-Ackers thermodynamic model to describe its binding, and compared the results of this model to experimentally measured Gal4p binding in vivo. We found that at many promoters, the model predicted no Gal4p binding, yet substantial binding was observed. These outlier promoters lacked canonical binding motifs, and subsequent investigation revealed Gal4p binds unexpectedly to DNA sequences with high densities of its half site (CGG). We confirmed this novel mode of binding through multiple experimental and computational paradigms; we also found most other zinc cluster TFs we tested frequently utilize this binding mode, at 27% of their targets on average. Together, these results demonstrate a novel mode of binding where zinc clusters, the largest class of TFs in yeast, bind DNA sequences with high densities of half sites.

MYT1L is required for suppressing earlier neuronal development programs in the adult mouse brain

In vitro studies indicate the neurodevelopmental disorder gene myelin transcription factor 1-like (MYT1L) suppresses non-neuronal lineage genes during fibroblast-to-neuron direct differentiation. However, MYT1L's molecular and cellular functions in the adult mammalian brain have not been fully characterized. Here, we found that MYT1L loss leads to up-regulated deep layer (DL) gene expression, corresponding to an increased ratio of DL/UL neurons in the adult mouse cortex. To define potential mechanisms, we conducted Cleavage Under Targets & Release Using Nuclease (CUT&RUN) to map MYT1L binding targets and epigenetic changes following MYT1L loss in mouse developing cortex and adult prefrontal cortex (PFC). We found MYT1L mainly binds to open chromatin, but with different transcription factor co-occupancies between promoters and enhancers. Likewise, multiomic data set integration revealed that, at promoters, MYT1L loss does not change chromatin accessibility but increases H3K4me3 and H3K27ac, activating both a subset of earlier neuronal development genes as well as Bcl11b, a key regulator for DL neuron development. Meanwhile, we discovered that MYT1L normally represses the activity of neurogenic enhancers associated with neuronal migration and neuronal projection development by closing chromatin structures and promoting removal of active histone marks. Further, we showed that MYT1L interacts with HDAC2 and transcriptional repressor SIN3B in vivo, providing potential mechanisms underlying repressive effects on histone acetylation and gene expression. Overall, our findings provide a comprehensive map of MYT1L binding in vivo and mechanistic insights into how MYT1L loss leads to aberrant activation of earlier neuronal development programs in the adult mouse brain.

Tn5 DNA Transposase in Multi-Omics Research

Tn5 transposase use in biotechnology has substantially advanced the sequencing applications of genome-wide analysis of cells. This is mainly due to the ability of Tn5 transposase to efficiently transpose DNA essentially randomly into any target DNA without the aid of other factors. This concise review is focused on the advances in Tn5 applications in multi-omics technologies, genome-wide profiling, and Tn5 hybrid molecule creation. The possibilities of other transposase uses are also discussed.

A single-cell massively parallel reporter assay detects cell-type-specific gene regulation

Massively parallel reporter gene assays are key tools in regulatory genomics but cannot be used to identify cell-type-specific regulatory elements without performing assays serially across different cell types. To address this problem, we developed a single-cell massively parallel reporter assay (scMPRA) to measure the activity of libraries of cis-regulatory sequences (CRSs) across multiple cell types simultaneously. We assayed a library of core promoters in a mixture of HEK293 and K562 cells and showed that scMPRA is a reproducible, highly parallel, single-cell reporter gene assay that detects cell-type-specific cis-regulatory activity. We then measured a library of promoter variants across multiple cell types in live mouse retinas and showed that subtle genetic variants can produce cell-type-specific effects on cis-regulatory activity. We anticipate that scMPRA will be widely applicable for studying the role of CRSs across diverse cell types.

Transposase N-terminal phosphorylation and asymmetric transposon ends inhibit piggyBac transposition in mammalian cells

DNA transposon systems are widely used in mammalian cells for genetic modification experiments, but their regulation remains poorly understood. We used biochemical and cell-based assays together with AlphaFold modeling and rational protein redesign to evaluate aspects of piggyBac transposition including the previously unexplained role of the transposase N-terminus and the need for asymmetric transposon ends for cellular activity. We found that phosphorylation at predicted casein kinase II sites in the transposase N-terminus inhibits transposition, most likely by preventing transposase-DNA interactions. Deletion of the region containing these sites releases inhibition thereby enhancing activity. We also found that the N-terminal domain promotes transposase dimerization in the absence of transposon DNA. When the N-terminus is deleted, the transposase gains the ability to carry out transposition using symmetric transposon left ends. This novel activity is also conferred by appending a second C-terminal domain. When combined, these modifications together result in a transposase that is highly active when symmetric transposon ends are used. Our results demonstrate that transposase N-terminal phosphorylation and the requirement for asymmetric transposon ends both negatively regulate piggyBac transposition in mammalian cells. These novel insights into the mechanism and structure of the piggyBac transposase expand its potential use for genomic applications.

Transcriptomic Features in a Single Extracellular Vesicle via Single-Cell RNA Sequencing

Although many studies have investigated functional molecules in extracellular vesicles (EVs), the exact number of ribonucleic acid molecules in a single-EV is unknown. Therefore, it is critical to explore the transcriptomic features and heterogeneity at the level of a single-EV. Here, using the 10x Genomics platform, the RNA cargos are profiled in single EVs derived from human K562 and mesenchymal stem cells. The key steps are labeling intact EVs using calcein-AM, detecting the EV concentration via flow cytometry, and using the CB2 algorithm with adaptive thresholds to effectively distinguish real EVs from background. The gene number in a single-EV varied from 6 to 148, with a mean of 52. Ribosomal genes, mitochondrial genes, and eukaryotic translation elongation factor 1 alpha has a high EV percentage in all EV samples. Hemoglobin genes are uniquely highly expressed in K562-EVs, and cytoskeleton genes are enriched in MSC-EVs. Ten or more clusters with different marker genes in each single-EV dataset demonstrated EV heterogeneity. Moreover, integrating EVs and their parental cells reveal both EVs and cells in each cluster, indicating different cell origins of various EVs. To the best of the author's knowledge, this study provides the first high-throughput transcriptome at the single-EV level and improves the understanding of EVs.

Measuring transcription factor binding and gene expression using barcoded self-reporting transposon calling cards and transcriptomes

Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2 and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells.

Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Transcription factors (TFs) play a critical role in determining cell fate decisions by integrating developmental and environmental signals through binding to specific cis-regulatory modules and regulating spatio-temporal specificity of gene expression patterns. Precise identification of functional TF binding sites in time and space not only will revolutionize our understanding of regulatory networks governing cell fate decisions but is also instrumental to uncover how genetic variations cause morphological diversity or disease. In this review, we discuss recent advances in mapping TF binding sites and characterizing the various parameters underlying the complexity of binding site recognition by TFs.

MYT1L in the making: emerging insights on functions of a neurodevelopmental disorder gene

Large scale human genetic studies have shown that loss of function (LoF) mutations in MYT1L are implicated in neurodevelopmental disorders (NDDs). Here, we provide an overview of the growing number of published MYT1L patient cases, and summarize prior studies in cells, zebrafish, and mice, both to understand MYT1L's molecular and cellular role during brain development and consider how its dysfunction can lead to NDDs. We integrate the conclusions from these studies and highlight conflicting findings to reassess the current model of the role of MYT1L as a transcriptional activator and/or repressor based on the biological context. Finally, we highlight additional functional studies that are needed to understand the molecular mechanisms underlying pathophysiology and propose key questions to guide future preclinical studies.

scPCOR-seq enables co-profiling of chromatin occupancy and RNAs in single cells

Cell-to-cell variation in gene expression is a widespread phenomenon, which may play important roles in cellular differentiation, function, and disease development^1–9. Chromatin is implicated in contributing to the cellular heterogeneity in gene expression^10–16. Fully understanding the mechanisms of cellular heterogeneity requires simultaneous measurement of RNA and occupancy of histone modifications and transcription factors on chromatin due to their critical roles in transcriptional regulation^17,18. We generally term the occupancy of histone modifications and transcription factors as Chromatin occupancy. Here, we report a technique, termed scPCOR-seq (single-cell Profiling of Chromatin Occupancy and RNAs Sequencing), for simultaneously profiling genome-wide chromatin protein binding or histone modification marks and RNA expression in the same cell. We demonstrated that scPCOR-seq can profile either H3K4me3 or RNAPII and RNAs in a mixture of human H1, GM12878 and 293 T cells at a single-cell resolution and either H3K4me3, RNAPII, or RNA profile can correctly separate the cells. Application of scPCOR-seq to the in vitro differentiation of the erythrocyte precursor CD36 cells from human CD34 stem or progenitor cells revealed that H3K4me3 and RNA exhibit distinct properties in clustering cells during differentiation. Overall, our work provides a promising approach to understand the relationships among different omics layers.

scPCOR-seq is a single-cell sequencing technique that enables simultaneous profiling of genome-wide chromatin protein binding or histone modification marks and RNA expression in the same cell.

The immune checkpoint ICOSLG is a relapse-predicting biomarker and therapeutic target in infant t(4;11) acute lymphoblastic leukemia

The most frequent genetic aberration leading to infant ALL (iALL) is the chromosomal translocation t(4;11), generating the fusion oncogenes KMT2A:AFF1 and AFF1:KMT2A, respectively. KMT2A-r iALL displays a dismal prognosis through high relapse rates and relapse-associated mortality. Relapse occurs frequently despite ongoing chemotherapy and without the accumulation of secondary mutations. A rational explanation for the observed chemo-resistance and satisfactory treatment options remain to be elucidated. We found that elevated ICOSLG expression level at diagnosis was associated with inferior event free survival (EFS) in a cohort of 43 patients with t(4;-11) iALL and that a cohort of 18 patients with iALL at relapse displayed strongly increased ICOSLG expression. Furthermore, co-culturing t(4;11) ALL cells (ICOSLG^hi) with primary T-cells resulted in the development of T_regs. This was impaired through treatment with a neutralizing ICOSLG antibody. These findings imply ICOSLG (1) as a relapse-predicting biomarker, and (2) as a therapeutic target involved in a potential immune evasion relapse-mechanism of infant t(4;11) ALL.

•

Early growth response 3 (EGR3) is a direct transactivator of the immune checkpoint gene ICOSLG

•

high ICOSLG expression at diagnosis is predictive for ALL relapse

•

EGR3 and ICOSLG expressions are relapse-associated

•

expression of ICOSLG on t(4;11) ALL cells leads to the rapid expansion of Tregs

Metabolic and epigenetic dysfunctions underlie the arrest of in vitro fertilized human embryos in a senescent-like state

Around 60% of in vitro fertilized (IVF) human embryos irreversibly arrest before compaction between the 3- to 8-cell stage, posing a significant clinical problem. The mechanisms behind this arrest are unclear. Here, we show that the arrested embryos enter a senescent-like state, marked by cell cycle arrest, the down-regulation of ribosomes and histones and down-regulation of MYC and p53 activity. The arrested embryos can be divided into 3 types. Type I embryos fail to complete the maternal-zygotic transition, and Type II/III embryos have low levels of glycolysis and either high (Type II) or low (Type III) levels of oxidative phosphorylation. Treatment with the SIRT agonist resveratrol or nicotinamide riboside (NR) can partially rescue the arrested phenotype, which is accompanied by changes in metabolic activity. Overall, our data suggests metabolic and epigenetic dysfunctions underlie the arrest of human embryos.

Genome-wide protein-DNA interaction site mapping in bacteria using a double-stranded DNA-specific cytosine deaminase

DNA-protein interactions are central to fundamental cellular processes, yet widely implemented technologies for measuring these interactions on a genome scale in bacteria are laborious and capture only a snapshot of binding events. We devised a facile method for mapping DNA-protein interaction sites in vivo using the double-stranded DNA-specific cytosine deaminase toxin DddA. In 3D-seq (DddA-sequencing), strains containing DddA fused to a DNA-binding protein of interest accumulate characteristic mutations in DNA sequence adjacent to sites occupied by the DNA-bound fusion protein. High-depth sequencing enables detection of sites of increased mutation frequency in these strains, yielding genome-wide maps of DNA-protein interaction sites. We validated 3D-seq for four transcription regulators in two bacterial species, Pseudomonas aeruginosa and Escherichia coli. We show that 3D-seq offers ease of implementation, the ability to record binding event signatures over time and the capacity for single-cell resolution.

Single-Cell Multiomics Techniques: From Conception to Applications

Recent advances in methods for single-cell analyses and barcoding strategies have led to considerable progress in research. The development of multiplexed assays offers the possibility to conduct parallel analyses of multiple factors and processes for comprehensive characterization of cellular and molecular states in health and disease. These technologies have expanded extremely rapidly in the past years and constantly evolve and provide better specificity, precision and resolution. This review summarizes recent progress in single-cell multiomics approaches, and focuses, in particular, on the most innovative techniques that integrate genome, epigenome and transcriptome profiling. It describes the methodologies, discusses their advantages and limitations, and explains how they have been applied to studies on cell heterogeneity and differentiation, and epigenetic reprogramming.

Genomic context sensitivity of insulator function

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter–enhancer specificity.

A deep generative model for multi-view profiling of single-cell RNA-seq and ATAC-seq data

Here, we present a multi-modal deep generative model, the single-cell Multi-View Profiler (scMVP), which is designed for handling sequencing data that simultaneously measure gene expression and chromatin accessibility in the same cell, including SNARE-seq, sci-CAR, Paired-seq, SHARE-seq, and Multiome from 10X Genomics. scMVP generates common latent representations for dimensionality reduction, cell clustering, and developmental trajectory inference and generates separate imputations for differential analysis and cis-regulatory element identification. scMVP can help mitigate data sparsity issues with imputation and accurately identify cell groups for different joint profiling techniques with common latent embedding, and we demonstrate its advantages on several realistic datasets.

Genomic environments scale the activities of diverse core promoters

A classical model of gene regulation is that enhancers provide specificity whereas core promoters provide a modular site for the assembly of the basal transcriptional machinery. However, examples of core promoter specificity have led to an alternate hypothesis in which specificity is achieved by core promoters with different sequence motifs that respond differently to genomic environments containing different enhancers and chromatin landscapes. To distinguish between these models, we measured the activities of hundreds of diverse core promoters in four different genomic locations and, in a complementary experiment, six different core promoters at thousands of locations across the genome. Although genomic locations had large effects on expression, the intrinsic activities of different classes of promoters were preserved across genomic locations, suggesting that core promoters are modular regulatory elements whose activities are independently scaled up or down by different genomic locations. This scaling of promoter activities is nonlinear and depends on the genomic location and the strength of the core promoter. Our results support the classical model of regulation in which diverse core promoter motifs set the intrinsic strengths of core promoters, which are then amplified or dampened by the activities of their genomic environments.

In vivo targeted DamID identifies CHD8 genomic targets in fetal mouse brain

Genetic studies of autism have revealed causal roles for chromatin remodeling gene mutations. Chromodomain helicase DNA binding protein 8 (CHD8) encodes a chromatin remodeler with significant de novo mutation rates in sporadic autism. However, relationships between CHD8 genomic function and autism-relevant biology remain poorly elucidated. Published studies utilizing ChIP-seq to map CHD8 protein-DNA interactions have high variability, consistent with technical challenges and limitations associated with this method. Thus, complementary approaches are needed to establish CHD8 genomic targets and regulatory functions in developing brain. We used in utero CHD8 Targeted DamID followed by sequencing (TaDa-seq) to characterize CHD8 binding in embryonic mouse cortex. CHD8 TaDa-seq reproduced interaction patterns observed from ChIP-seq and further highlighted CHD8 distal interactions associated with neuronal loci. This study establishes TaDa-seq as a useful alternative for mapping protein-DNA interactions in vivo and provides insights into the regulatory targets of CHD8 and autism-relevant pathophysiology associated with CHD8 mutations.

Here and there: the double-side transgene localization

Random transgene integration is a powerful tool for developing new genome-wide screening approaches. These techniques have already been used for functional gene annotation by transposon-insertion sequencing, for identif ication of transcription factor binding sites and regulatory sequences, and for dissecting chromatin position effects. Precise localization of transgenes and accurate artifact f iltration are essential for this type of method. To date, many mapping assays have been developed, including Inverse-PCR, TLA, LAM-PCR, and splinkerette PCR. However, none of them is able to ensure localization of both transgene’s f lanking regions simultaneously, which would be necessary for some applications. Here we proposed a cheap and simple NGS-based approach that overcomes this limitation. The developed assay requires using intentionally designed vectors that lack recognition sites of one or a set of restriction enzymes used for DNA fragmentation. By looping and sequencing these DNA fragments, we obtain special data that allows us to link the two f lanking regions of the transposon. This can be useful for precise insertion mapping and for screening approaches in the f ield of chromosome engineering, where chromosomal recombination events between transgenes occur in a cell population. To demonstrate the method’s feasibility, we applied it for mapping SB transposon integration in the human HAP1 cell line. Our technique allowed us to eff iciently localize genomic transposon integrations, which was conf irmed via PCR analysis. For practical application of this approach, we proposed a set of recommendations and a normalization strategy. The developed method can be used for multiplex transgene localization and detection of rearrangements between them.

A single-cell guide to retinal development: Cell fate decisions of multipotent retinal progenitors in scRNA-seq

Recent advances in high throughput single-cell RNA sequencing (scRNA-seq) technology have enabled the simultaneous transcriptomic profiling of thousands of individual cells in a single experiment. To investigate the intrinsic process of retinal development, researchers have leveraged this technology to quantify gene expression in retinal cells across development, in multiple species, and from numerous important models of human disease. In this review, we summarize recent applications of scRNA-seq and discuss how these datasets have complemented and advanced our understanding of retinal progenitor cell competence, cell fate specification, and differentiation. Finally, we also highlight the outstanding questions in the field that advances in single-cell data generation and analysis will soon be able to answer.

Ultra-Fast Label-Free Serum Metabolic Diagnosis of Coronary Heart Disease via a Deep Stabilizer

Although mass spectrometry (MS) of metabolites has the potential to provide real-time monitoring of patient status for diagnostic purposes, the diagnostic application of MS is limited due to sample treatment and data quality/reproducibility. Here, the generation of a deep stabilizer for ultra-fast, label-free MS detection and the application of this method for serum metabolic diagnosis of coronary heart disease (CHD) are reported. Nanoparticle-assisted laser desorption/ionization-MS is used to achieve direct metabolic analysis of trace unprocessed serum in seconds. Furthermore, a deep stabilizer is constructed to map native MS results to high-quality results obtained by established methods. Finally, using the newly developed protocol and diagnosis variation characteristic surface to characterize sensitivity/specificity and variation, CHD is diagnosed with advanced accuracy in a high-throughput/speed manner. This work advances design of metabolic analysis tools for disease detection as it provides a direct label-free, ultra-fast, and stabilized platform for future protocol development in clinics.

The qBED track: a novel genome browser visualization for point processes

SUMMARY

Transposon calling cards is a genomic assay for identifying transcription factor binding sites in both bulk and single cell experiments. Here, we describe the qBED format, an open, text-based standard for encoding and analyzing calling card data. In parallel, we introduce the qBED track on the WashU Epigenome Browser, a novel visualization that enables researchers to inspect calling card data in their genomic context. Finally, through examples, we demonstrate that qBED files can be used to visualize non-calling card datasets, such as Combined Annotation-Dependent Depletion scores and GWAS/eQTL hits, and thus may have broad utility to the genomics community.

AVAILABILITY AND IMPLEMENTATION

The qBED track is available on the WashU Epigenome Browser (http://epigenomegateway.wustl.edu/browser), beginning with version 46. Source code for the WashU Epigenome Browser with qBED support is available on GitHub (http://github.com/arnavm/eg-react and http://github.com/lidaof/eg-react). A complete definition of the qBED format is available as part of the WashU Epigenome Browser documentation (https://eg.readthedocs.io/en/latest/tracks.html#qbed-track). We have also released a tutorial on how to upload qBED data to the browser (http://dx.doi.org/10.17504/protocols.io.bca8ishw).

Single-cell joint detection of chromatin occupancy and transcriptome enables higher-dimensional epigenomic reconstructions

Deciphering mechanisms in cell-fate decisions requires single-cell holistic reconstructions of multidimensional epigenomic states in transcriptional regulation. Here we develop CoTECH, a combinatorial barcoding method allowing high-throughput single-cell joint detection of chromatin occupancy and transcriptome. We used CoTECH to examine bivalent histone marks (H3K4me3 and H3K27me3) with transcription from naive to primed mouse embryonic stem cells. We also derived concurrent bivalent marks in pseudosingle cells using transcriptome as an anchor for resolving pseudotemporal bivalency trajectories and disentangling a context-specific interplay between H3K4me3/H3K27me3 and transcription level. Next, we revealed the regulatory basis of endothelial-to-hematopoietic transition in two waves of hematopoietic cells and distinctive enhancer-gene-linking schemes guiding hemogenic endothelial cell emergence, indicating a unique epigenetic control of transcriptional regulation for hematopoietic stem cell priming. CoTECH provides an efficient framework for single-cell coassay of chromatin occupancy and transcription, thus enabling higher-dimensional epigenomic reconstructions.

Single Cell Technologies: Beyond Microfluidics

Single-cell RNA-sequencing (scRNA-seq) has been widely adopted in recent years due to standardized protocols and automation, reliability, and standardized bioinformatic pipelines. The most widely adopted platform is the 10× Genomics solution. Although powerful, this system is limited by its high cost, moderate throughput, and the inability to customize due to fixed kit components. This study will cover new approaches that do not rely on microfluidics and thus have low entry costs, are highly customizable, and are within the reach of any laboratory possessing molecular biology expertise.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

Non-genetic mechanisms of therapeutic resistance in cancer

Therapeutic resistance continues to be an indominable foe in our ambition for curative cancer treatment. Recent insights into the molecular determinants of acquired treatment resistance in the clinical and experimental setting have challenged the widely held view of sequential genetic evolution as the primary cause of resistance and brought into sharp focus a range of non-genetic adaptive mechanisms. Notably, the genetic landscape of the tumour and the non-genetic mechanisms used to escape therapy are frequently linked. Remarkably, whereas some oncogenic mutations allow the cancer cells to rapidly adapt their transcriptional and/or metabolic programme to meet and survive the therapeutic pressure, other oncogenic drivers convey an inherent cellular plasticity to the cancer cell enabling lineage switching and/or the evasion of anticancer immunosurveillance. The prevalence and diverse array of non-genetic resistance mechanisms pose a new challenge to the field that requires innovative strategies to monitor and counteract these adaptive processes. In this Perspective we discuss the key principles of non-genetic therapy resistance in cancer. We provide a perspective on the emerging data from clinical studies and sophisticated cancer models that have studied various non-genetic resistance pathways and highlight promising therapeutic avenues that may be used to negate and/or counteract the non-genetic adaptive pathways.

Single-cell lineage tracing approaches in hematology research: technical considerations

The coordinated differentiation of hematopoietic stem and progenitor cells (HSPCs) into the various mature blood cell types is responsible for sustaining blood and immune system homeostasis. The cell fate decisions underlying this important biological process are made at the level of single cells. Methods to trace the fate of single cells are therefore essential for understanding hematopoietic system activity in health and disease and have had a major impact on how we understand and represent hematopoiesis. Here, we discuss the basic methodologies and technical considerations for three important clonal assays: single-cell transplantation, lentiviral barcoding, and Sleeping Beauty barcoding. This perspective is a synthesis of presentations and discussions from the 2019 International Society for Experimental Hematology (ISEH) Annual Meeting New Investigator Technology Session and the 2019 ISEH Winter Webinar.