A human cell atlas of fetal chromatin accessibility

Advances in profiling chromatin architecture shed light on the regulatory dynamics underlying brain disorders

Understanding the exquisitely complex nature of the three-dimensional organization of the genome and how it affects gene regulation remains a central question in biology. Recent advances in sequencing- and imaging-based approaches in decoding the three-dimensional chromatin landscape have enabled a systematic characterization of gene regulatory architecture. In this review, we outline how chromatin architecture provides a reference atlas to predict the functional consequences of non-coding variants associated with human traits and disease. High-throughput perturbation assays such as massively parallel reporter assays (MPRA) and CRISPR-based genome engineering in combination with a reference atlas opened an avenue for going beyond observational studies to experimentally validating the regulatory principles of the genome. We conclude by providing a suggested path forward by calling attention to barriers that can be addressed for a more complete understanding of the regulatory landscape of the human brain.

Variant to function mapping at single-cell resolution through network propagation

Genome-wide association studies in combination with single-cell genomic atlases can provide insights into the mechanisms of disease-causal genetic variation. However, identification of disease-relevant or trait-relevant cell types, states and trajectories is often hampered by sparsity and noise, particularly in the analysis of single-cell epigenomic data. To overcome these challenges, we present SCAVENGE, a computational algorithm that uses network propagation to map causal variants to their relevant cellular context at single-cell resolution. We demonstrate how SCAVENGE can help identify key biological mechanisms underlying human genetic variation, applying the method to blood traits at distinct stages of human hematopoiesis, to monocyte subsets that increase the risk for severe Coronavirus Disease 2019 (COVID-19) and to intermediate lymphocyte developmental states that predispose to acute leukemia. Our approach not only provides a framework for enabling variant-to-function insights at single-cell resolution but also suggests a more general strategy for maximizing the inferences that can be made using single-cell genomic data.

Recapitulating human myogenesis ex vivo using human pluripotent stem cells

Human pluripotent stem cells (hPSCs) provide a human model for developmental myogenesis, disease modeling and development of therapeutics. Differentiation of hPSCs into muscle stem cells has the potential to provide a cell-based therapy for many skeletal muscle wasting diseases. This review describes the current state of hPSCs towards recapitulating human myogenesis ex vivo, considerations of stem cell and progenitor cell state as well as function for future use of hPSC-derived muscle cells in regenerative medicine.

Subtype-Independent ANP32E Reduction During Breast Cancer Progression in Accordance with Chromatin Relaxation

BACKGROUND

Chromatin state provides a clear decipherable blueprint for maintenance of transcriptional patterns, exemplifying a mitotically stable form of cellular programming in dividing cells. In this regard, genomic studies of chromatin states within cancerous tissues have the potential to uncover novel aspects of tumor biology and unique mechanisms associated with disease phenotypes and outcomes. The degree to which chromatin state differences occur in accordance with breast cancer features has not been established.

METHODS

We applied a series of unsupervised computational methods to identify chromatin and molecular differences associated with discrete physiologies across human breast cancer tumors.

RESULTS

Chromatin patterns alone are capable of stratifying tumors in association with cancer subtype and disease progression. Major differences occur at DNA motifs for the transcription factor FOXA1, in hormone receptor-positive tumors, and motifs for SOX9 in Basal-like tumors. We find that one potential driver of this effect, the histone chaperone ANP32E, is inversely correlated with tumor progression and relaxation of chromatin at FOXA1 binding sites. Tumors with high levels of ANP32E exhibit an immune response and proliferative gene expression signature, whereas tumors with low ANP32E levels appear programmed for differentiation.

CONCLUSIONS

Our results indicate that ANP32E may function through chromatin state regulation to control breast cancer differentiation and tumor plasticity. This study sets a precedent for future computational studies of chromatin changes in carcinogenesis.

Epigenetic loss of heterogeneity from low to high grade localized prostate tumours

Identifying precise molecular subtypes attributable to specific stages of localized prostate cancer has proven difficult due to high levels of heterogeneity. Bulk assays represent a population-average, which mask the heterogeneity that exists at the single-cell level. In this work, we sequence the accessible chromatin regions of 14,424 single-cells from 18 flash-frozen prostate tumours. We observe shared chromatin features among low-grade prostate cancer cells are lost in high-grade tumours. Despite this loss, high-grade tumours exhibit an enrichment for FOXA1, HOXB13 and CDX2 transcription factor binding sites, indicating a shared trans-regulatory programme. We identify two unique genes encoding neuronal adhesion molecules that are highly accessible in high-grade prostate tumours. We show NRXN1 and NLGN1 expression in epithelial, endothelial, immune and neuronal cells in prostate cancer using cyclic immunofluorescence. Our results provide a deeper understanding of the active gene regulatory networks in primary prostate tumours, critical for molecular stratification of the disease.

Intronic elements associated with insomnia and restless legs syndrome exhibit cell type-specific epigenetic features contributing to MEIS1 regulation

A highly evolutionarily conserved MEIS1 intronic region is strongly associated with restless legs syndrome (RLS) and insomnia. To understand its regulatory function, we dissected the region by analyzing chromatin accessibility, enhancer-promoter contacts, DNA methylation, and eQTLs in different human neural cell types and tissues. We observed specific activity with respect to cell type and developmental maturation, indicating a prominent role for distinct highly conserved intronic elements in forebrain inhibitory neuron differentiation. Two elements were hypomethylated in neural cells with higher MEIS1 expression, suggesting a role of enhancer demethylation in gene regulation. MEIS1 eQTLs showed a striking modular chromosomal distribution, with forebrain eQTLs clustering in intron 8/9. CRISPR interference targeting of individual elements in this region attenuated MEIS1 expression, revealing a complex regulatory interplay of distinct elements. In summary, we found that MEIS1 regulation is organized in a modular pattern. Disease-associated intronic regulatory elements control MEIS1 expression with cell type and maturation stage specificity, particularly in the inhibitory neuron lineage. The precise spatiotemporal activity of these elements likely contributes to the pathogenesis of insomnia and RLS.

Multi-omic insights into Parkinson's Disease: From genetic associations to functional mechanisms

Genome-Wide Association Studies (GWAS) have elucidated the genetic components of Parkinson's Disease (PD). However, because the vast majority of GWAS association signals fall within non-coding regions, translating these results into an interpretable, mechanistic understanding of the disease etiology remains a major challenge in the field. In this review, we provide an overview of the approaches to prioritize putative causal variants and genes as well as summarise the primary findings of previous studies. We then discuss recent efforts to integrate multi-omics data to identify likely pathogenic cell types and biological pathways implicated in PD pathogenesis. We have compiled full summary statistics of cell-type, tissue, and phentoype enrichment analyses from multiple studies of PD GWAS and provided them in a standardized format as a resource for the research community (https://github.com/RajLabMSSM/PD_omics_review). Finally, we discuss the experimental, computational, and conceptual advances that will be necessary to fully elucidate the effects of functional variants and genes on cellular dysregulation and disease risk.

Cell type-specific analysis by single-cell profiling identifies a stable mammalian tRNA-mRNA interface and increased translation efficiency in neurons

The correlation between codon and anticodon pools influences the efficiency of translation, but whether differences exist in these pools across individual cells is unknown. We determined that codon usage and amino acid demand are highly stable across different cell types using available mouse and human single-cell RNA-sequencing atlases. After showing the robustness of ATAC-sequencing measurements for the analysis of tRNA gene usage, we quantified anticodon usage and amino acid supply in both mouse and human single-cell ATAC-seq atlases. We found that tRNA gene usage is overall coordinated across cell types, except in neurons, which clustered separately from other cell types. Integration of these data sets revealed a strong and statistically significant correlation between amino acid supply and demand across almost all cell types. Neurons have an enhanced translation efficiency over other cell types, driven by an increased supply of tRNA^Ala (AGC) anticodons. This results in faster decoding of the Ala-GCC codon, as determined by cell type–specific ribosome profiling, suggesting that the reduction of tRNA^Ala (AGC) anticodon pools may be implicated in neurological pathologies. This study, the first such examination of codon usage, anticodon usage, and translation efficiency resolved at the cell-type level with single-cell information, identifies a conserved landscape of translation elongation across mammalian cellular diversity and evolution.

Genome-wide cancer-specific chromatin accessibility patterns derived from archival processed xenograft tumors

Chromatin accessibility states that influence gene expression and other nuclear processes can be altered in disease. The constellation of transcription factors and chromatin regulatory complexes in cells results in characteristic patterns of chromatin accessibility. The study of these patterns in tissues has been limited because existing chromatin accessibility assays are ineffective for archival formalin-fixed, paraffin-embedded (FFPE) tissues. We have developed a method to efficiently extract intact chromatin from archival tissue via enhanced cavitation with a nanodroplet reagent consisting of a lipid shell with a liquid perfluorocarbon core. Inclusion of nanodroplets during the extraction of chromatin from FFPE tissues enhances the recovery of intact accessible and nucleosome-bound chromatin. We show that the addition of nanodroplets to the chromatin accessibility assay formaldehyde-assisted isolation of regulatory elements (FAIRE), does not affect the accessible chromatin signal. Applying the technique to FFPE human tumor xenografts, we identified tumor-relevant regions of accessible chromatin shared with those identified in primary tumors. Further, we deconvoluted non-tumor signal to identify cellular components of the tumor microenvironment. Incorporation of this method of enhanced cavitation into FAIRE offers the potential for extending chromatin accessibility to clinical diagnosis and personalized medicine, while also enabling the exploration of gene regulatory mechanisms in archival samples.

High-content single-cell combinatorial indexing

Single-cell combinatorial indexing (sci) with transposase-based library construction increases the throughput of single-cell genomics assays but produces sparse coverage in terms of usable reads per cell. We develop symmetrical strand sci ('s3'), a uracil-based adapter switching approach that improves the rate of conversion of source DNA into viable sequencing library fragments following tagmentation. We apply this chemistry to assay chromatin accessibility (s3-assay for transposase-accessible chromatin, s3-ATAC) in human cortical and mouse whole-brain tissues, with mouse datasets demonstrating a six- to 13-fold improvement in usable reads per cell compared with other available methods. Application of s3 to single-cell whole-genome sequencing (s3-WGS) and to whole-genome plus chromatin conformation (s3-GCC) yields 148- and 14.8-fold improvements, respectively, in usable reads per cell compared with sci-DNA-sequencing and sci-HiC. We show that s3-WGS and s3-GCC resolve subclonal genomic alterations in patient-derived pancreatic cancer cell lines. We expect that the s3 platform will be compatible with other transposase-based techniques, including sci-MET or CUT&Tag.

Epigenetics and tissue immunity-Translating environmental cues into functional adaptations

There is an increasing appreciation that many innate and adaptive immune cell subsets permanently reside within non-lymphoid organs, playing a critical role in tissue homeostasis and defense. The best characterized are macrophages and tissue-resident T lymphocytes that work in concert with organ structural cells to generate appropriate immune responses and are functionally shaped by organ-specific environmental cues. The interaction of tissue epithelial, endothelial and stromal cells is also required to attract, differentiate, polarize and maintain organ immune cells in their tissue niche. All of these processes require dynamic regulation of cellular transcriptional programmes, with epigenetic mechanisms playing a critical role, including DNA methylation and post-translational histone modifications. A failure to appropriately regulate immune cell transcription inevitably results in inadequate or inappropriate immune responses and organ pathology. Here, with a focus on the mammalian kidney, an organ which generates differing regional environmental cues (including hypersalinity and hypoxia) due to its physiological functions, we will review the basic concepts of tissue immunity, discuss the technologies available to profile epigenetic modifications in tissue immune cells, including those that enable single-cell profiling, and consider how these mechanisms influence the development, phenotype, activation and function of different tissue immune cell subsets, as well as the immunological function of structural cells.

A single-cell atlas of chromatin accessibility in the human genome

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

Gene regulatory networks controlling temporal patterning, neurogenesis, and cell-fate specification in mammalian retina

Gene regulatory networks (GRNs), consisting of transcription factors and their target sites, control neurogenesis and cell-fate specification in the developing central nervous system. In this study, we use integrated single-cell RNA and single-cell ATAC sequencing (scATAC-seq) analysis in developing mouse and human retina to identify multiple interconnected, evolutionarily conserved GRNs composed of cell-type-specific transcription factors that both activate genes within their own network and inhibit genes in other networks. These GRNs control temporal patterning in primary progenitors, regulate transition from primary to neurogenic progenitors, and drive specification of each major retinal cell type. We confirm that NFI transcription factors selectively activate expression of genes promoting late-stage temporal identity in primary retinal progenitors and identify other transcription factors that regulate rod photoreceptor specification in postnatal retina. This study inventories cis- and trans-acting factors that control retinal development and can guide cell-based therapies aimed at replacing retinal neurons lost to disease.

New Approaches to Assess Mechanisms of Action of Selective Vitamin D Analogues

Recent studies of transcription have revealed an advanced set of overarching principles that govern vitamin D action on a genome-wide scale. These tenets of vitamin D transcription have emerged as a result of the application of now well-established techniques of chromatin immunoprecipitation coupled to next-generation DNA sequencing that have now been linked directly to CRISPR-Cas9 genomic editing in culture cells and in mouse tissues in vivo. Accordingly, these techniques have established that the vitamin D hormone modulates sets of cell-type specific genes via an initial action that involves rapid binding of the VDR-ligand complex to multiple enhancer elements at open chromatin sites that drive the expression of individual genes. Importantly, a sequential set of downstream events follows this initial binding that results in rapid histone acetylation at these sites, the recruitment of additional histone modifiers across the gene locus, and in many cases, the appearance of H3K36me3 and RNA polymerase II across gene bodies. The measured recruitment of these factors and/or activities and their presence at specific regions in the gene locus correlate with the emerging presence of cognate transcripts, thereby highlighting sequential molecular events that occur during activation of most genes both in vitro and in vivo. These features provide a novel approach to the study of vitamin D analogs and their actions in vivo and suggest that they can be used for synthetic compound evaluation and to select for novel tissue- and gene-specific features. This may be particularly useful for ligand activation of nuclear receptors given the targeting of these factors directly to genetic sites in the nucleus.

scEnhancer: a single-cell enhancer resource with annotation across hundreds of tissue/cell types in three species

Previous studies on enhancers and their target genes were largely based on bulk samples that represent ‘average’ regulatory activities from a large population of millions of cells, masking the heterogeneity and important effects from the sub-populations. In recent years, single-cell sequencing technology has enabled the profiling of open chromatin accessibility at the single-cell level (scATAC-seq), which can be used to annotate the enhancers and promoters in specific cell types. A comprehensive resource is highly desirable for exploring how the enhancers regulate the target genes at the single-cell level. Hence, we designed a single-cell database scEnhancer (http://enhanceratlas.net/scenhancer/), covering 14 527 776 enhancers and 63 658 600 enhancer-gene interactions from 1 196 906 single cells across 775 tissue/cell types in three species. An unsupervised learning method was employed to sort and combine tens or hundreds of single cells in each tissue/cell type to obtain the consensus enhancers. In addition, we utilized a cis-regulatory network algorithm to identify the enhancer-gene connections. Finally, we provided a user-friendly platform with seven useful modules to search, visualize, and browse the enhancers/genes. This database will facilitate the research community towards a functional analysis of enhancers at the single-cell level.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Chromatin-accessibility estimation from single-cell ATAC-seq data with scOpen

A major drawback of single-cell ATAC-seq (scATAC-seq) is its sparsity, i.e., open chromatin regions with no reads due to loss of DNA material during the scATAC-seq protocol. Here, we propose scOpen, a computational method based on regularized non-negative matrix factorization for imputing and quantifying the open chromatin status of regulatory regions from sparse scATAC-seq experiments. We show that scOpen improves crucial downstream analysis steps of scATAC-seq data as clustering, visualization, cis-regulatory DNA interactions, and delineation of regulatory features. We demonstrate the power of scOpen to dissect regulatory changes in the development of fibrosis in the kidney. This identifies a role of Runx1 and target genes by promoting fibroblast to myofibroblast differentiation driving kidney fibrosis.

Identification of LZTFL1 as a candidate effector gene at a COVID-19 risk locus

The severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) disease (COVID-19) pandemic has caused millions of deaths worldwide. Genome-wide association studies identified the 3p21.31 region as conferring a twofold increased risk of respiratory failure. Here, using a combined multiomics and machine learning approach, we identify the gain-of-function risk A allele of an SNP, rs17713054G>A, as a probable causative variant. We show with chromosome conformation capture and gene-expression analysis that the rs17713054-affected enhancer upregulates the interacting gene, leucine zipper transcription factor like 1 (LZTFL1). Selective spatial transcriptomic analysis of lung biopsies from patients with COVID-19 shows the presence of signals associated with epithelial-mesenchymal transition (EMT), a viral response pathway that is regulated by LZTFL1. We conclude that pulmonary epithelial cells undergoing EMT, rather than immune cells, are likely responsible for the 3p21.31-associated risk. Since the 3p21.31 effect is conferred by a gain-of-function, LZTFL1 may represent a therapeutic target.

Improving Chromatin-Interaction Prediction Using Single-Cell Open-Chromatin Profiles and Making Insight Into the Cis-Regulatory Landscape of the Human Brain

Single-cell open-chromatin profiles have the potential to reveal the pattern of chromatin-interaction in a cell type. However, currently available cis-regulatory network prediction methods using single-cell open-chromatin profiles focus more on local chromatin interactions despite the fact that long-range interactions among genomic sites play a significant role in gene regulation. Here, we propose a method that predicts both short and long-range interactions among genomic sites using single-cell open chromatin profiles. Our method, termed as single-cell epigenome based chromatin-interaction analysis (scEChIA) exploits signal imputation and refined L1 regularization. For a few single-cell open-chromatin profiles, scEChIA outperformed other tools even in terms of accuracy of prediction. Using scEChIA, we predicted almost 0.7 million interactions among genomic sites across seven cell types in the human brain. Further analysis revealed cell type for connection between genes and expression quantitative trait locus (eQTL) in the human brain and making insight about target genes of human-accelerated-elements and disease-associated mutations. Our analysis enabled by scEChIA also hints about the possible action of a few transcription factors (TFs), especially through long-range interaction in brain endothelial cells.

An integrated pipeline for mammalian genetic screening

With the recent advancements in genome editing, next-generation sequencing (NGS), and scalable cloning techniques, scientists can now conduct genetic screens at unprecedented levels of scale and precision. With such a multitude of technologies, there is a need for a simple yet comprehensive pipeline to enable systematic mammalian genetic screening. In this study, we develop unique algorithms for target identification and a toxin-less Gateway cloning tool, termed MegaGate, for library cloning which, when combined with existing genetic perturbation methods and NGS-coupled readouts, enable versatile engineering of relevant mammalian cell lines. Our integrated pipeline for sequencing-based target ascertainment and modular perturbation screening (STAMPScreen) can thus be utilized for a host of cell state engineering applications.

Rewiring of human neurodevelopmental gene regulatory programs by human accelerated regions

Human accelerated regions (HARs) are the fastest-evolving regions of the human genome, and many are hypothesized to function as regulatory elements that drive human-specific gene regulatory programs. We interrogate the in vitro enhancer activity and in vivo epigenetic landscape of more than 3,100 HARs during human neurodevelopment, demonstrating that many HARs appear to act as neurodevelopmental enhancers and that sequence divergence at HARs has largely augmented their neuronal enhancer activity. Furthermore, we demonstrate PPP1R17 to be a putative HAR-regulated gene that has undergone remarkable rewiring of its cell type and developmental expression patterns between non-primates and primates and between non-human primates and humans. Finally, we show that PPP1R17 slows neural progenitor cell cycle progression, paralleling the cell cycle length increase seen predominantly in primate and especially human neurodevelopment. Our findings establish HARs as key components in rewiring human-specific neurodevelopmental gene regulatory programs and provide an integrated resource to study enhancer activity of specific HARs.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

BACKGROUND

Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging.

METHODS

Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain.

RESULTS

Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism.

CONCLUSION

This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

Cell-type-specific Hox regulatory strategies orchestrate tissue identity

Hox proteins are homeodomain transcription factors that diversify serially homologous segments along the animal body axis, as revealed by the classic bithorax phenotype of Drosophila melanogaster, in which mutations in Ultrabithorax (Ubx) transform the third thoracic segment into the likeness of the second thoracic segment. To specify segment identity, we show that Ubx both increases and decreases chromatin accessibility, coinciding with its dual role as both an activator and repressor of transcription. However, the choice of transcriptional activity executed by Ubx is spatially regulated and depends on the availability of cofactors, with Ubx acting as a repressor in some populations and as an activator in others. Ubx-mediated changes to chromatin accessibility positively and negatively affect the binding of Scalloped (Sd), a transcription factor that is required for appendage development in both segments. These findings illustrate how a single Hox protein can modify complex gene regulatory networks to transform the identity of an entire tissue.

Novel gene discovery for hearing loss and other routes to increased diagnostic rates

Despite decades of research, there is much to be learned about the genetic landscape of sensorineural hearing loss. Novel genes for hearing loss remain to be identified while ‘secrets’ of the known genes need to be uncovered. These ‘secrets’ include regulatory mechanisms of gene activity and novel aspects of gene structure. To obtain a more complete picture of the genetics of hearing loss, the available experimental and bioinformatic tools need to be fully exploited. This is also true for data resources such as ENCODE. For the inner ear, however, such data resources and analytical tools need to be developed or extended. Collaborative studies provide opportunities to achieve this and to optimally use those tools and resources that are already available. This will accelerate the discoveries that are necessary for improving molecular genetic diagnostics and genetic counselling and for the development of therapeutic strategies.

Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics.

Single cell RNA sequencing approaches to cardiac development and congenital heart disease

The development of single cell RNA sequencing technologies has accelerated the ability of scientists to understand healthy and disease states of the cardiovascular system. Congenital heart defects occur in approximately 40,000 births each year and 1 out of 4 children are born with critical congenital heart disease requiring surgical interventions and a lifetime of monitoring. An understanding of how the normal heart develops and how each cell contributes to normal and pathological anatomy is an important goal in pediatric cardiovascular research. Single cell sequencing has provided the tools to increase the ability to discover rare cell types and novel genes involved in normal cardiac development. Knowledge of gene expression of single cells within cardiac tissue has contributed to the understanding of how each cell type contributes to the anatomic structures of the heart. In this review, we summarize how single cell RNA sequencing has been utilized to understand cardiac developmental processes and congenital heart disease. We discuss the advantages and disadvantages of whole cell versus single nuclei RNA sequencing and describe the approaches to analyze the interactomes, transcriptomes, and differentiation trajectory from single cell data. We summarize the currently available single cell RNA sequencing technologies and technical aspects of performing single cell analysis and how to overcome common obstacles. We also review data from the recently published human and mouse fetal heart atlases and advancements that have occurred within the field due to the application of these single cell tools. Finally we highlight the potential for single cell technologies to uncover novel mechanisms of disease pathogenesis by leveraging findings from genome wide association studies.

Tagmentation-based single-cell genomics

It has been just over 10 years since the initial description of transposase-based methods to prepare high-throughput sequencing libraries, or "tagmentation," in which a hyperactive transposase is used to simultaneously fragment target DNA and append universal adapter sequences. Tagmentation effectively replaced a series of processing steps in traditional workflows with one single reaction. It is the simplicity, coupled with the high efficiency of tagmentation, that has made it a favored means of sequencing library construction and fueled a diverse range of adaptations to assay a variety of molecular properties. In recent years, this has been centered in the single-cell space with a catalog of tagmentation-based assays that have been developed, covering a substantial swath of the regulatory landscape. To date, there have been a number of excellent reviews on single-cell technologies structured around the molecular properties that can be profiled. This review is instead framed around the central components and properties of tagmentation and how they have enabled the development of innovative molecular tools to probe the regulatory landscape of single cells. Furthermore, the granular specifics on cell throughput or richness of data provided by the extensive list of individual technologies are not discussed. Such metrics are rapidly changing and highly sample specific and are better left to studies that directly compare technologies for assays against one another in a rigorously controlled framework. The hope for this review is that, in laying out the diversity of molecular techniques at each stage of these assay platforms, new ideas may arise for others to pursue that will further advance the field of single-cell genomics.

Intricacies of single-cell multi-omics data integration

A wealth of single-cell protocols makes it possible to characterize different molecular layers at unprecedented resolution. Integrating the resulting multimodal single-cell data to find cell-to-cell correspondences remains a challenge. We argue that data integration needs to happen at a meaningful biological level of abstraction and that it is necessary to consider the inherent discrepancies between modalities to strike a balance between biological discovery and noise removal. A survey of current methods reveals that a distinction between technical and biological origins of presumed unwanted variation between datasets is not yet commonly considered. The increasing availability of paired multimodal data will aid the development of improved methods by providing a ground truth on cell-to-cell matches.

A roadmap for the Human Developmental Cell Atlas

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

Prioritization and functional analysis of GWAS risk loci for Barrett's esophagus and esophageal adenocarcinoma

Genome-wide association studies (GWAS) have identified ~ 20 genetic susceptibility loci for esophageal adenocarcinoma (EAC), and its precursor, Barrett's esophagus (BE). Despite such advances, functional/causal variants and gene targets at these loci remain undefined, hindering clinical translation. A key challenge is that most causal variants map to non-coding regulatory regions such as enhancers, and typically, numerous potential candidate variants at GWAS loci require testing. We developed a systematic informatics pipeline for prioritizing candidate functional variants via integrative functional potential scores consolidated from multi-omics annotations, and used this pipeline to identify two high-scoring variants for experimental interrogation: chr9q22.32/rs11789015 and chr19p13.11/rs10423674. Minimal candidate enhancer regions spanning these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two variants tested (rs10423674) exhibited allele-specific enhancer activity. CRISPR-mediated deletion of the putative enhancer region in EAC cell lines correlated with reduced expression of two genes-CREB-regulated transcription coactivator 1 (CRTC1) and Cartilage oligomeric matrix protein (COMP); expression of five other genes remained unchanged (CRLF1, KLHL26, TMEM59L, UBA52, RFXANK). Expression quantitative trait locus (eQTL) mapping indicated that rs10423674 genotype correlated with CRTC1 and COMP expression in normal esophagus. This study represents the first experimental effort to bridge GWAS associations to biology in BE/EAC, and supports the utility of functional potential scores to guide variant prioritization. Our findings reveal a functional variant and candidate risk enhancer at chr19p13.11, and implicate CRTC1 and COMP as putative gene targets, suggesting that altered expression of these genes may underlie the BE/EAC risk association.

EpiScanpy: integrated single-cell epigenomic analysis

EpiScanpy is a toolkit for the analysis of single-cell epigenomic data, namely single-cell DNA methylation and single-cell ATAC-seq data. To address the modality specific challenges from epigenomics data, epiScanpy quantifies the epigenome using multiple feature space constructions and builds a nearest neighbour graph using epigenomic distance between cells. EpiScanpy makes the many existing scRNA-seq workflows from scanpy available to large-scale single-cell data from other -omics modalities, including methods for common clustering, dimension reduction, cell type identification and trajectory learning techniques, as well as an atlas integration tool for scATAC-seq datasets. The toolkit also features numerous useful downstream functions, such as differential methylation and differential openness calling, mapping epigenomic features of interest to their nearest gene, or constructing gene activity matrices using chromatin openness. We successfully benchmark epiScanpy against other scATAC-seq analysis tools and show its outperformance at discriminating cell types.

Developmental and evolutionary dynamics of cis-regulatory elements in mouse cerebellar cells

Organ development is orchestrated by cell- and time-specific gene regulatory networks. In this study, we investigated the regulatory basis of mouse cerebellum development from early neurogenesis to adulthood. By acquiring snATAC-seq (single-nucleus assay for transposase accessible chromatin using sequencing) profiles for ~90,000 cells spanning 11 stages, we mapped cerebellar cell types and identified candidate cis-regulatory elements (CREs). We detected extensive spatiotemporal heterogeneity among progenitor cells and a gradual divergence in the regulatory programs of cerebellar neurons during differentiation. Comparisons to vertebrate genomes and snATAC-seq profiles for ∼20,000 cerebellar cells from the marsupial opossum revealed a shared decrease in CRE conservation during development and differentiation as well as differences in constraint between cell types. Our work delineates the developmental and evolutionary dynamics of gene regulation in cerebellar cells and provides insights into mammalian organ development.

Enhancers in disease: molecular basis and emerging treatment strategies

Enhancers are genomic sequences that play a key role in regulating tissue-specific gene expression levels. An increasing number of diseases are linked to impaired enhancer function through chromosomal rearrangement, genetic variation within enhancers, or epigenetic modulation. Here, we review how these enhancer disruptions have recently been implicated in congenital disorders, cancers, and common complex diseases and address the implications for diagnosis and treatment. Although further fundamental research into enhancer function, target genes, and context is required, enhancer-targeting drugs and gene editing approaches show great therapeutic promise for a range of diseases.

Mechanisms of Binding Specificity among bHLH Transcription Factors

The transcriptome of every cell is orchestrated by the complex network of interaction between transcription factors (TFs) and their binding sites on DNA. Disruption of this network can result in many forms of organism malfunction but also can be the substrate of positive natural selection. However, understanding the specific determinants of each of these individual TF-DNA interactions is a challenging task as it requires integrating the multiple possible mechanisms by which a given TF ends up interacting with a specific genomic region. These mechanisms include DNA motif preferences, which can be determined by nucleotide sequence but also by DNA’s shape; post-translational modifications of the TF, such as phosphorylation; and dimerization partners and co-factors, which can mediate multiple forms of direct or indirect cooperative binding. Binding can also be affected by epigenetic modifications of putative target regions, including DNA methylation and nucleosome occupancy. In this review, we describe how all these mechanisms have a role and crosstalk in one specific family of TFs, the basic helix-loop-helix (bHLH), with a very conserved DNA binding domain and a similar DNA preferred motif, the E-box. Here, we compile and discuss a rich catalog of strategies used by bHLH to acquire TF-specific genome-wide landscapes of binding sites.

Single-Cell Analysis for Whole-Organism Datasets

Cell atlases are essential companions to the genome as they elucidate how genes are used in a cell type-specific manner or how the usage of genes changes over the lifetime of an organism. This review explores recent advances in whole-organism single-cell atlases, which enable understanding of cell heterogeneity and tissue and cell fate, both in health and disease. Here we provide an overview of recent efforts to build cell atlases across species and discuss the challenges that the field is currently facing. Moreover, we propose the concept of having a knowledgebase that can scale with the number of experiments and computational approaches and a new feedback loop for development and benchmarking of computational methods that includes contributions from the users. These two aspects are key for community efforts in single-cell biology that will help produce a comprehensive annotated map of cell types and states with unparalleled resolution.

Molecular Profiling of the Cardiac Conduction System: the Dawn of a New Era

PURPOSE OF REVIEW

Recent technological advances have led to an increased ability to define the gene expression profile of the cardiac conduction system (CCS). Here, we review the most salient studies to emerge in recent years and discuss existing gaps in our knowledge as well as future areas of investigation.

RECENT FINDINGS

Molecular profiling of the CCS spans several decades. However, the advent of high-throughput sequencing strategies has allowed for the discovery of unique transcriptional programs of the many diverse CCS cell types. The CCS, a diverse structure with significant inter- and intra-component cellular heterogeneity, is essential to the normal function of the heart. Progress in transcriptomic profiling has improved the resolution and depth of characterization of these unique and clinically relevant CCS cell types. Future studies leveraging this big data will play a crucial role in improving our understanding of CCS development and function as well as translating these findings into tangible translational tools for the improved detection, prevention, and treatment of cardiac arrhythmias.

A Bayesian inference transcription factor activity model for the analysis of single-cell transcriptomes

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful experimental approach to study cellular heterogeneity. One of the challenges in scRNA-seq data analysis is integrating different types of biological data to consistently recognize discrete biological functions and regulatory mechanisms of cells, such as transcription factor activities and gene regulatory networks in distinct cell populations. We have developed an approach to infer transcription factor activities from scRNA-seq data that leverages existing biological data on transcription factor binding sites. The Bayesian inference transcription factor activity model (BITFAM) integrates ChIP-seq transcription factor binding information into scRNA-seq data analysis. We show that the inferred transcription factor activities for key cell types identify regulatory transcription factors that are known to mechanistically control cell function and cell fate. The BITFAM approach not only identifies biologically meaningful transcription factor activities, but also provides valuable insights into underlying transcription factor regulatory mechanisms.

A transcriptional switch governs fibroblast activation in heart disease

In diseased organs, stress-activated signalling cascades alter chromatin, thereby triggering maladaptive cell state transitions. Fibroblast activation is a common stress response in tissues that worsens lung, liver, kidney and heart disease, yet its mechanistic basis remains unclear1,2. Pharmacological inhibition of bromodomain and extra-terminal domain (BET) proteins alleviates cardiac dysfunction3-7, providing a tool to interrogate and modulate cardiac cell states as a potential therapeutic approach. Here we use single-cell epigenomic analyses of hearts dynamically exposed to BET inhibitors to reveal a reversible transcriptional switch that underlies the activation of fibroblasts. Resident cardiac fibroblasts demonstrated robust toggling between the quiescent and activated state in a manner directly correlating with BET inhibitor exposure and cardiac function. Single-cell chromatin accessibility revealed previously undescribed DNA elements, the accessibility of which dynamically correlated with cardiac performance. Among the most dynamic elements was an enhancer that regulated the transcription factor MEOX1, which was specifically expressed in activated fibroblasts, occupied putative regulatory elements of a broad fibrotic gene program and was required for TGFβ-induced fibroblast activation. Selective CRISPR inhibition of the single most dynamic cis-element within the enhancer blocked TGFβ-induced Meox1 activation. We identify MEOX1 as a central regulator of fibroblast activation associated with cardiac dysfunction and demonstrate its upregulation after activation of human lung, liver and kidney fibroblasts. The plasticity and specificity of BET-dependent regulation of MEOX1 in tissue fibroblasts provide previously unknown trans- and cis-targets for treating fibrotic disease.

Single-cell chromatin accessibility landscape of human umbilical cord blood in trisomy 18 syndrome

BACKGROUND

Trisomy 18 syndrome (Edwards syndrome, ES) is a type of aneuploidy caused by the presence of an extra chromosome 18. Aneuploidy is the leading cause of early pregnancy loss, intellectual disability, and multiple congenital anomalies. The research of trisomy 18 is progressing slowly, and the molecular characteristics of the disease mechanism and phenotype are still largely unclear.

RESULTS

In this study, we used the commercial Chromium platform (10× Genomics) to perform sc-ATAC-seq to measure chromatin accessibility in 11,611 single umbilical cord blood cells derived from one trisomy 18 syndrome patient and one healthy donor. We obtained 13 distinct major clusters of cells and identified them as 6 human umbilical cord blood mononuclear cell types using analysis tool. Compared with the NC group, the ES group had a lower ratio of T cells to NK cells, the ratio of monocytes/DC cell population did not change significantly, and the ratio of B cell nuclear progenitor and megakaryocyte erythroid cells was higher. The differential genes of ME-0 are enriched in Human T cell leukemia virus 1 infection pathway, and the differential peak genes of ME-1 are enriched in apopotosis pathway. We found that CCNB2 and MCM3 may be vital to the development of trisomy 18. CCNB2 and MCM3, which have been reported to be essential components of the cell cycle and chromatin.

CONCLUSIONS

We have identified 6 cell populations in cord blood. Disorder in megakaryocyte erythroid cells implicates trisomy 18 in perturbing fetal hematopoiesis. We identified a pathway in which the master differential regulatory pathway in the ME-0 cell population involves human T cell leukemia virus 1 infection, a pathway that is dysregulated in patients with trisomy 18 and which may increase the risk of leukemia in patients with trisomy 18. CCNB2 and MCM3 in progenitor may be vital to the development of trisomy 18. CCNB2 and MCM3, which have been reported to be essential components of the cell cycle and chromatin, may be related to chromosomal abnormalities in trisomy 18.

Shifting epigenetic contexts influence regulatory variation and disease risk

Epigenetic shifts are a hallmark of aging that impact transcriptional networks at regulatory level. These shifts may modify the effects of genetic regulatory variants during aging and contribute to disease pathomechanism. However, these shifts occur on the backdrop of epigenetic changes experienced throughout an individual's development into adulthood; thus, the phenotypic, and ultimately fitness, effects of regulatory variants subject to developmental- versus aging-related epigenetic shifts may differ considerably. Natural selection therefore may act differently on variants depending on their changing epigenetic context, which we propose as a novel lens through which to consider regulatory sequence evolution and phenotypic effects. Here, we define genomic regions subjected to altered chromatin accessibility as tissues transition from their fetal to adult forms, and subsequently from early to late adulthood. Based on these epigenomic datasets, we examine patterns of evolutionary constraint and potential functional impacts of sequence variation (e.g., genetic disease risk associations). We find that while the signals observed with developmental epigenetic changes are consistent with stronger fitness consequences (i.e., negative selection pressures), they tend to have weaker effects on genetic risk associations for aging-related diseases. Conversely, we see stronger effects of variants with increased local accessibility in adult tissues, strongest in young adult when compared to old. We propose a model for how epigenetic status of a region may influence the effects of evolutionary relevant sequence variation, and suggest that such a perspective on gene regulatory networks may elucidate our understanding of aging biology.

How to find genomic regions relevant for gene regulation

Genetic variants associated with human diseases are often located outside the protein coding regions of the genome. Identification and functional characterization of the regulatory elements in the non-coding genome is therefore of crucial importance for understanding the consequences of genetic variation and the mechanisms of disease. The past decade has seen rapid progress in high-throughput analysis and mapping of chromatin accessibility, looping, structure, and occupancy by transcription factors, as well as epigenetic modifications, all of which contribute to the proper execution of regulatory functions in the non-coding genome. Here, we review the current technologies for the definition and functional validation of non-coding regulatory regions in the genome.

North Carolina Macular Dystrophy: Phenotypic Variability and Computational Analysis of Disease-Associated Noncoding Variants

Purpose

North Carolina macular dystrophy (NCMD) is an autosomal dominant, congenital disorder affecting the central retina. Here, we report clinical and genetic findings in three families segregating NCMD and use epigenomic datasets from human tissues to gain insights into the effect of NCMD-implicated variants.

Methods

Clinical assessment and genetic testing were performed. Publicly available transcriptomic and epigenomic datasets were analyzed and the activity-by-contact method for scoring enhancer elements and linking them to target genes was used.

Results

A previously described, heterozygous, noncoding variant upstream of the PRDM13 gene was detected in all six affected study participants (chr6:100,040,987G>C [GRCh37/hg19]). Interfamilial and intrafamilial variability were observed; the visual acuity ranged from 0.0 to 1.6 LogMAR and fundoscopic findings ranged from visually insignificant, confluent, drusen-like macular deposits to coloboma-like macular lesions. Variable degrees of peripheral retinal spots (which were easily detected on widefield retinal imaging) were observed in all study subjects. Notably, a 6-year-old patient developed choroidal neovascularization and required treatment with intravitreal bevacizumab injections. Computational analysis of the five single nucleotide variants that have been implicated in NCMD revealed that these noncoding changes lie within two putative enhancer elements; these elements are predicted to interact with PRDM13 in the developing human retina. PRDM13 was found to be expressed in the fetal retina, with greatest expression in the amacrine precursor cell population.

Conclusions

We provide further evidence supporting the role of PRDM13 dysregulation in the pathogenesis of NCMD and highlight the usefulness of widefield retinal imaging in individuals suspected to have this condition.

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.

Defining the Emerging Blood System During Development at Single-Cell Resolution

Developmental hematopoiesis differs from adult and is far less described. In the developing embryo, waves of lineage-restricted blood precede the ultimate emergence of definitive hematopoietic stem cells (dHSCs) capable of maintaining hematopoiesis throughout life. During the last two decades, the advent of single-cell genomics has provided tools to circumvent previously impeding characteristics of embryonic hematopoiesis, such as cell heterogeneity and rare cell states, allowing for definition of lineage trajectories, cellular hierarchies, and cell-type specification. The field has rapidly advanced from microfluidic platforms and targeted gene expression analysis, to high throughput unbiased single-cell transcriptomic profiling, single-cell chromatin analysis, and cell tracing-offering a plethora of tools to resolve important questions within hematopoietic development. Here, we describe how these technologies have been implemented to address a wide range of aspects of embryonic hematopoiesis ranging from the gene regulatory network of dHSC formation via endothelial to hematopoietic transition (EHT) and how EHT can be recapitulated in vitro, to hematopoietic trajectories and cell fate decisions. Together, these studies have important relevance for regenerative medicine and for our understanding of genetic blood disorders and childhood leukemias.

Computational principles and challenges in single-cell data integration

The development of single-cell multimodal assays provides a powerful tool for investigating multiple dimensions of cellular heterogeneity, enabling new insights into development, tissue homeostasis and disease. A key challenge in the analysis of single-cell multimodal data is to devise appropriate strategies for tying together data across different modalities. The term 'data integration' has been used to describe this task, encompassing a broad collection of approaches ranging from batch correction of individual omics datasets to association of chromatin accessibility and genetic variation with transcription. Although existing integration strategies exploit similar mathematical ideas, they typically have distinct goals and rely on different principles and assumptions. Consequently, new definitions and concepts are needed to contextualize existing methods and to enable development of new methods.

Cardiac cell type-specific gene regulatory programs and disease risk association

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.