Principles of nucleosome organization revealed by single-cell micrococcal nuclease sequencing

NucMap 2.0: An Updated Database of Genome-wide Nucleosome Positioning Maps Across Species

Nucleosome dynamics plays important roles in many biological processes, such as DNA replication and gene expression. NucMap (https://ngdc.cncb.ac.cn/nucmap) is the first database of genome-wide nucleosome positioning maps across species. Here, we present an updated version, NucMap 2.0, by incorporating more species and MNase-seq samples. In addition, we integrate other related omics data for each MNase-seq sample to provide a comprehensive view of nucleosome positioning, such as gene expression, transcription factor binding sites, histone modifications and DNA methylation. In particular, NucMap 2.0 integrates and pre-analyzes RNA-seq data and ChIP-seq data of human-related samples, which facilitates the interpretation of nucleosome positioning in humans. All processed data are integrated into an in-built genome browser, and users can make comprehensive side-by-side analyses. In addition, more online analytical functions are developed, which allows researchers to identify differential nucleosome regions and explore potential gene regulatory regions. All resources are open access with a user-friendly web interface.

Epigenetic regulation of innate lymphoid cells

Innate lymphoid cells (ILCs) lack antigen-specific receptors and are considered the innate arm of the immune system, phenotypically and functionally mirroring CD4+ helper T cells. ILCs are categorized into groups 1, 2, and 3 based on transcription factors and cytokine expression. ILCs predominantly reside in mucosal tissues and play important roles in regional immune responses. The development and function of ILC subsets are controlled by both transcriptional and epigenetic mechanisms, which have been extensively studied in recent years. Epigenetic regulation refers to inheritable changes in gene expression that occur without affecting DNA sequences. This mainly includes chromatin status, histone modifications, and DNA methylation. In this review, we summarize recent discoveries on epigenetic mechanisms regulating ILC development and function, and how these regulations affect disease progression under pathological conditions. Although the ablation of specific epigenetic regulators can cause global changes in corresponding epigenetic modifications to the chromatin, only partial genes with altered epigenetic modifications change their mRNA expression, resulting in specific outcomes in cell differentiation and function. Therefore, elucidating epigenetic mechanisms underlying the regulation of ILCs will provide potential targets for the diagnosis and treatment of inflammatory diseases.

Cell-type-specific consequences of mosaic structural variants in hematopoietic stem and progenitor cells

The functional impact and cellular context of mosaic structural variants (mSVs) in normal tissues is understudied. Utilizing Strand-seq, we sequenced 1,133 single-cell genomes from 19 human donors of increasing age, and discovered the heterogeneous mSV landscapes of hematopoietic stem and progenitor cells. While mSVs are continuously acquired throughout life, expanded subclones in our cohort are confined to individuals >60. Cells already harboring mSVs are more likely to acquire additional somatic structural variants, including megabase-scale segmental aneuploidies. Capitalizing on comprehensive single-cell micrococcal nuclease digestion with sequencing reference data, we conducted high-resolution cell-typing for eight hematopoietic stem and progenitor cells. Clonally expanded mSVs disrupt normal cellular function by dysregulating diverse cellular pathways, and enriching for myeloid progenitors. Our findings underscore the contribution of mSVs to the cellular and molecular phenotypes associated with the aging hematopoietic system, and establish a foundation for deciphering the molecular links between mSVs, aging and disease susceptibility in normal tissues.

Acute depletion of BRG1 reveals its primary function as an activator of transcription

The mammalian SWI/SNF-like BAF complexes play critical roles during animal development and pathological conditions. Previous gene deletion studies and characterization of human gene mutations implicate that the complexes both repress and activate a large number of genes. However, the direct function of the complexes in cells remains largely unclear due to the relatively long-term nature of gene deletion or natural mutation. Here we generate a mouse line by knocking in the auxin-inducible degron tag (AID) to the Smarca4 gene, which encodes BRG1, the essential ATPase subunit of the BAF complexes. We show that the tagged BRG1 can be efficiently depleted by osTIR1 expression and auxin treatment for 6 to 10 h in CD4 + T cells, hepatocytes, and fibroblasts isolated from the knock-in mice. The acute depletion of BRG1 leads to decreases in nascent RNAs and RNA polymerase II binding at a large number of genes, which are positively correlated with the loss of BRG1. Further, these changes are correlated with diminished accessibility at DNase I Hypersensitive Sites (DHSs) and p300 binding. The acute BRG1 depletion results in three major patterns of nucleosome shifts leading to narrower nucleosome spacing surrounding transcription factor motifs and at enhancers and transcription start sites (TSSs), which are correlated with loss of BRG1, decreased chromatin accessibility and decreased nascent RNAs. Acute depletion of BRG1 severely compromises the Trichostatin A (TSA) -induced histone acetylation, suggesting a substantial interplay between the chromatin remodeling activity of BRG1 and histone acetylation. Our data suggest BRG1 mainly plays a direct positive role in chromatin accessibility, RNAPII binding, and nascent RNA production by regulating nucleosome positioning and facilitating transcription factor binding to their target sites.

Context-specific functions of chromatin remodellers in development and disease

Chromatin remodellers were once thought to be highly redundant and nonspecific in their actions. However, recent human genetic studies demonstrate remarkable biological specificity and dosage sensitivity of the thirty-two adenosine triphosphate (ATP)-dependent chromatin remodellers encoded in the human genome. Mutations in remodellers produce many human developmental disorders and cancers, motivating efforts to investigate their distinct functions in biologically relevant settings. Exquisitely specific biological functions seem to be an emergent property in mammals, and in many cases are based on the combinatorial assembly of subunits and the generation of stable, composite surfaces. Critical interactions between remodelling complex subunits, the nucleosome and other transcriptional regulators are now being defined from structural and biochemical studies. In addition, in vivo analyses of remodellers at relevant genetic loci have provided minute-by-minute insights into their dynamics. These studies are proposing new models for the determinants of remodeller localization and function on chromatin.

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Hypobaric hypoxia modulated structural characteristics of circulating cell-free DNA in high-altitude pulmonary edema

The utility of cell-free (cf) DNA has extended as a surrogate or clinical biomarker for various diseases. However, a more profound and expanded understanding of the diverse cfDNA population and its correlation with physiological phenotypes and environmental factors is imperative for using its full potential. The high-altitude (HA; altitude > 2,500 m above sea level) environment characterized by hypobaric hypoxia offers an observational case-control design to study the differential cfDNA profile in patients with high-altitude pulmonary edema (HAPE) (number of subjects, n = 112) and healthy HA sojourners (n = 111). The present study investigated cfDNA characteristics such as concentration, fragment length size, degree of integrity, and subfractions reflecting mitochondrial-cfDNA copies in the two groups. The total cfDNA level was significantly higher in patients with HAPE, and the level increased with increasing HAPE severity (P = 0.0036). A lower degree of cfDNA integrity of 0.346 in patients with HAPE (P = 0.001) indicated the prevalence of shorter cfDNA fragments in circulation in patients compared with the healthy HA sojourners. A significant correlation of cfDNA characteristics with the peripheral oxygen saturation levels in the patient group demonstrated the translational relevance of cfDNA molecules. The correlation was further supported by multivariate logistic regression and receiver operating characteristic curve. To our knowledge, our study is the first to highlight the association of higher cfDNA concentration, a lower degree of cfDNA integrity, and increased mitochondrial-derived cfDNA population with HAPE disease severity. Further deep profiling of cfDNA fragments, which preserves cell-type specific genetic and epigenetic features, can provide dynamic physiological responses to hypoxia.NEW & NOTEWORTHY This study observed altered cell-free (cf) DNA fragment patterns in patients with high-altitude pulmonary edema and the significant correlation of these patterns with peripheral oxygen saturation levels. This suggests deep profiling of cfDNA fragments in the future may identify genetic and epigenetic mechanisms underlying physiological and pathophysiological responses to hypoxia.

Nucleosome Dynamics Derived at the Single-Molecule Level Bridges Its Structures and Functions

Nucleosome, the building block of chromatin, plays pivotal roles in all DNA-related processes. While cryogenic-electron microscopy (cryo-EM) has significantly advanced our understanding of nucleosome structures, the emerging field of single-molecule force spectroscopy is illuminating their dynamic properties. This technique is crucial for revealing how nucleosome behavior is influenced by chaperones, remodelers, histone variants, and post-translational modifications, particularly in their folding and unfolding mechanisms under tension. Such insights are vital for deciphering the complex interplay in nucleosome assembly and structural regulation, highlighting the nucleosome's versatility in response to DNA activities. In this Perspective, we aim to consolidate the latest advancements in nucleosome dynamics, with a special focus on the revelations brought forth by single-molecule manipulation. Our objective is to highlight the insights gained from studying nucleosome dynamics through this innovative approach, emphasizing the transformative impact of single-molecule manipulation techniques in the field of chromatin research.

Clonally heritable gene expression imparts a layer of diversity within cell types

Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information.

Histone globular domain epigenetic modifications: The regulators of chromatin dynamics in malaria parasite

Plasmodium species adapt a complex lifecycle with multiple phenotypes to survive inside various cell types of humans and mosquitoes. Stage-specific gene expression in the developmental stages of parasites is tightly controlled in Plasmodium species; however, the underlying mechanisms have yet to be explored. Genome organization and gene expression for each stage of the malaria parasite need to be better characterized. Recent studies indicated that epigenetic modifications of histone proteins play a vital role in chromatin plasticity. Like other eukaryotes, Plasmodium species N-terminal tail modifications form a distinct "histone code," which creates the docking sites for histone reader proteins, including gene activator/repressor complexes, to regulate gene expression. The emerging research findings shed light on various unconventional epigenetic changes in histone proteins' core/globular domain regions, which might contribute to the chromatin organization in different developmental stages of the malaria parasite. The malaria parasite lost many transcription factors during evolution, and it is proposed that the nature of local chromatin structure essentially regulates the stage-specific gene expression. This review highlights recent discoveries of unconventional histone globular domain epigenetic modifications and their functions in regulating chromatin structure dynamics in various developmental stages of malaria parasites.

Chromatin organization of muscle stem cell

The proper functioning of skeletal muscles is essential throughout life. A crucial crosstalk between the environment and several cellular mechanisms allows striated muscles to perform successfully. Notably, the skeletal muscle tissue reacts to an injury producing a completely functioning tissue. The muscle's robust regenerative capacity relies on the fine coordination between muscle stem cells (MuSCs or "satellite cells") and their specific microenvironment that dictates stem cells' activation, differentiation, and self-renewal. Critical for the muscle stem cell pool is a fine regulation of chromatin organization and gene expression. Acquiring a lineage-specific 3D genome architecture constitutes a crucial modulator of muscle stem cell function during development, in the adult stage, in physiological and pathological conditions. The context-dependent relationship between genome structure, such as accessibility and chromatin compartmentalization, and their functional effects will be analysed considering the improved 3D epigenome knowledge, underlining the intimate liaison between environmental encounters and epigenetics.

Progress in single-cell multimodal sequencing and multi-omics data integration

With the rapid advance of single-cell sequencing technology, cell heterogeneity in various biological processes was dissected at different omics levels. However, single-cell mono-omics results in fragmentation of information and could not provide complete cell states. In the past several years, a variety of single-cell multimodal omics technologies have been developed to jointly profile multiple molecular modalities, including genome, transcriptome, epigenome, and proteome, from the same single cell. With the availability of single-cell multimodal omics data, we can simultaneously investigate the effects of genomic mutation or epigenetic modification on transcription and translation, and reveal the potential mechanisms underlying disease pathogenesis. Driven by the massive single-cell omics data, the integration method of single-cell multi-omics data has rapidly developed. Integration of the massive multi-omics single-cell data in public databases in the future will make it possible to construct a cell atlas of multi-omics, enabling us to comprehensively understand cell state and gene regulation at single-cell resolution. In this review, we summarized the experimental methods for single-cell multimodal omics data and computational methods for multi-omics data integration. We also discussed the future development of this field.

Structure of the ISW1a complex bound to the dinucleosome

Nucleosomes are basic repeating units of chromatin and form regularly spaced arrays in cells. Chromatin remodelers alter the positions of nucleosomes and are vital in regulating chromatin organization and gene expression. Here we report the cryo-EM structure of chromatin remodeler ISW1a complex from Saccharomyces cerevisiae bound to the dinucleosome. Each subunit of the complex recognizes a different nucleosome. The motor subunit binds to the mobile nucleosome and recognizes the acidic patch through two arginine residues, while the DNA-binding module interacts with the entry DNA at the nucleosome edge. This nucleosome-binding mode provides the structural basis for linker DNA sensing of the motor. Notably, the Ioc3 subunit recognizes the disk face of the adjacent nucleosome through interacting with the H4 tail, the acidic patch and the nucleosomal DNA, which plays a role in the spacing activity in vitro and in nucleosome organization and cell fitness in vivo. Together, these findings support the nucleosome spacing activity of ISW1a and add a new mode of nucleosome remodeling in the context of a chromatin environment.

A single fiber view of the nucleosome organization in eukaryotic chromatin

Eukaryotic cells are thought to arrange nucleosomes into extended arrays with evenly spaced nucleosomes phased at genomic landmarks. Here we tested to what extent this stereotypic organization describes the nucleosome organization in Saccharomyces cerevisiae using Fiber-Seq, a long-read sequencing technique that maps entire nucleosome arrays on individual chromatin fibers in a high throughput manner. With each fiber coming from a different cell, Fiber-Seq uncovers cell-to-cell heterogeneity. The long reads reveal the nucleosome architecture even over repetitive DNA such as the ribosomal DNA repeats. The absolute nucleosome occupancy, a parameter that is difficult to obtain with conventional sequencing approaches, is a direct readout of Fiber-Seq. We document substantial deviations from the stereotypical nucleosome organization with unexpectedly long linker DNAs between nucleosomes, gene bodies missing entire nucleosomes, cell-to-cell heterogeneity in nucleosome occupancy, heterogeneous phasing of arrays and irregular nucleosome spacing. Nucleosome array structures are indistinguishable throughout the gene body and with respect to the direction of transcription arguing against transcription promoting array formation. Acute nucleosome depletion destroyed most of the array organization indicating that nucleosome remodelers cannot efficiently pack nucleosomes under those conditions. Given that nucleosomes are cis-regulatory elements, the cell-to-cell heterogeneity uncovered by Fiber-Seq provides much needed information to understand chromatin structure and function.

A genome-wide comprehensive analysis of nucleosome positioning in yeast

In eukaryotic cells, the one-dimensional DNA molecules need to be tightly packaged into the spatially constraining nucleus. Folding is achieved on its lowest level by wrapping the DNA around nucleosomes. Their arrangement regulates other nuclear processes, such as transcription and DNA repair. Despite strong efforts to study nucleosome positioning using Next Generation Sequencing (NGS) data, the mechanism of their collective arrangement along the gene body remains poorly understood. Here, we classify nucleosome distributions of protein-coding genes in Saccharomyces cerevisiae according to their profile similarity and analyse their differences using functional Principal Component Analysis. By decomposing the NGS signals into their main descriptive functions, we compared wild type and chromatin remodeler-deficient strains, keeping position-specific details preserved whilst considering the nucleosome arrangement as a whole. A correlation analysis with other genomic properties, such as gene size and length of the upstream Nucleosome Depleted Region (NDR), identified key factors that influence the nucleosome distribution. We reveal that the RSC chromatin remodeler-which is responsible for NDR maintenance-is indispensable for decoupling nucleosome arrangement within the gene from positioning outside, which interfere in rsc8-depleted conditions. Moreover, nucleosome profiles in chd1Δ strains displayed a clear correlation with RNA polymerase II presence, whereas wild type cells did not indicate a noticeable interdependence. We propose that RSC is pivotal for global nucleosome organisation, whilst Chd1 plays a key role for maintaining local arrangement.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

Oleic acid availability impacts thymocyte preprogramming and subsequent peripheral Treg cell differentiation

The nature of activation signals is essential in determining T cell subset differentiation; however, the features that determine T cell subset preference acquired during intrathymic development remain elusive. Here we show that naive CD4+ T cells generated in the mouse thymic microenvironment lacking Scd1, encoding the enzyme catalyzing oleic acid (OA) production, exhibit enhanced regulatory T (Treg) cell differentiation and attenuated development of experimental autoimmune encephalomyelitis. Scd1 deletion in K14+ thymic epithelia recapitulated the enhanced Treg cell differentiation phenotype of Scd1-deficient mice. The dearth of OA permitted DOT1L to increase H3K79me2 levels at the Atp2a2 locus of thymocytes at the DN2-DN3 transition stage. Such epigenetic modification persisted in naive CD4+ T cells and facilitated Atp2a2 expression. Upon T cell receptor activation, ATP2A2 enhanced the activity of the calcium-NFAT1-Foxp3 axis to promote naive CD4+ T cells to differentiate into Treg cells. Therefore, OA availability is critical for preprogramming thymocytes with Treg cell differentiation propensities in the periphery.

Single-cell sequencing technology applied to epigenetics for the study of tumor heterogeneity

BACKGROUND

Previous studies have traditionally attributed the initiation of cancer cells to genetic mutations, considering them as the fundamental drivers of carcinogenesis. However, recent research has shed light on the crucial role of epigenomic alterations in various cell types present within the tumor microenvironment, suggesting their potential contribution to tumor formation and progression. Despite these significant findings, the progress in understanding the epigenetic mechanisms regulating tumor heterogeneity has been impeded over the past few years due to the lack of appropriate technical tools and methodologies.

RESULTS

The emergence of single-cell sequencing has enhanced our understanding of the epigenetic mechanisms governing tumor heterogeneity by revealing the distinct epigenetic layers of individual cells (chromatin accessibility, DNA/RNA methylation, histone modifications, nucleosome localization) and the diverse omics (transcriptomics, genomics, multi-omics) at the single-cell level. These technologies provide us with new insights into the molecular basis of intratumoral heterogeneity and help uncover key molecular events and driving mechanisms in tumor development.

CONCLUSION

This paper provides a comprehensive review of the emerging analytical and experimental approaches of single-cell sequencing in various omics, focusing specifically on epigenomics. These approaches have the potential to capture and integrate multiple dimensions of individual cancer cells, thereby revealing tumor heterogeneity and epigenetic features. Additionally, this paper outlines the future trends of these technologies and their current technical limitations.

Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation

High-throughput sequencing technology is central to our current understanding of the human methylome. The vast majority of studies use chemical conversion to analyse bulk-level patterns of DNA methylation across the genome from a population of cells. While this technology has been used to probe single-molecule methylation patterns, such analyses are limited to short reads of a few hundred basepairs. DNA methylation can also be directly detected using Nanopore sequencing which can generate reads measuring megabases in length. However, thus far these analyses have largely focused on bulk-level assessment of DNA methylation. Here, we analyse DNA methylation in single Nanopore reads from human lymphoblastoid cells, to show that bulk-level metrics underestimate large-scale heterogeneity in the methylome. We use the correlation in methylation state between neighbouring sites to quantify single-molecule heterogeneity and find that heterogeneity varies significantly across the human genome, with some regions having heterogeneous methylation patterns at the single-molecule level and others possessing more homogeneous methylation patterns. By comparing the genomic distribution of the correlation to epigenomic annotations, we find that the greatest heterogeneity in single-molecule patterns is observed within heterochromatic partially methylated domains (PMDs). In contrast, reads originating from euchromatic regions and gene bodies have more ordered DNA methylation patterns. By analysing the patterns of single molecules in more detail, we show the existence of a nucleosome-scale periodicity in DNA methylation that accounts for some of the heterogeneity we uncover in long single-molecule DNA methylation patterns. We find that this periodic structure is partially masked in bulk data and correlates with DNA accessibility as measured by nanoNOMe-seq, suggesting that it could be generated by nucleosomes. Our findings demonstrate the power of single-molecule analysis of long-read data to understand the structure of the human methylome.

Nucleosome density shapes kilobase-scale regulation by a mammalian chromatin remodeler

Nearly all essential nuclear processes act on DNA packaged into arrays of nucleosomes. However, our understanding of how these processes (for example, DNA replication, RNA transcription, chromatin extrusion and nucleosome remodeling) occur on individual chromatin arrays remains unresolved. Here, to address this deficit, we present SAMOSA-ChAAT: a massively multiplex single-molecule footprinting approach to map the primary structure of individual, reconstituted chromatin templates subject to virtually any chromatin-associated reaction. We apply this method to distinguish between competing models for chromatin remodeling by the essential imitation switch (ISWI) ATPase SNF2h: nucleosome-density-dependent spacing versus fixed-linker-length nucleosome clamping. First, we perform in vivo single-molecule nucleosome footprinting in murine embryonic stem cells, to discover that ISWI-catalyzed nucleosome spacing correlates with the underlying nucleosome density of specific epigenomic domains. To establish causality, we apply SAMOSA-ChAAT to quantify the activities of ISWI ATPase SNF2h and its parent complex ACF on reconstituted nucleosomal arrays of varying nucleosome density, at single-molecule resolution. We demonstrate that ISWI remodelers operate as density-dependent, length-sensing nucleosome sliders, whose ability to program DNA accessibility is dictated by single-molecule nucleosome density. We propose that the long-observed, context-specific regulatory effects of ISWI complexes can be explained in part by the sensing of nucleosome density within epigenomic domains. More generally, our approach promises molecule-precise views of the essential processes that shape nuclear physiology.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

New genetic and epigenetic insights into the chemokine system: the latest discoveries aiding progression toward precision medicine

Over the past thirty years, the importance of chemokines and their seven-transmembrane G protein-coupled receptors (GPCRs) has been increasingly recognized. Chemokine interactions with receptors trigger signaling pathway activity to form a network fundamental to diverse immune processes, including host homeostasis and responses to disease. Genetic and nongenetic regulation of both the expression and structure of chemokines and receptors conveys chemokine functional heterogeneity. Imbalances and defects in the system contribute to the pathogenesis of a variety of diseases, including cancer, immune and inflammatory diseases, and metabolic and neurological disorders, which render the system a focus of studies aiming to discover therapies and important biomarkers. The integrated view of chemokine biology underpinning divergence and plasticity has provided insights into immune dysfunction in disease states, including, among others, coronavirus disease 2019 (COVID-19). In this review, by reporting the latest advances in chemokine biology and results from analyses of a plethora of sequencing-based datasets, we outline recent advances in the understanding of the genetic variations and nongenetic heterogeneity of chemokines and receptors and provide an updated view of their contribution to the pathophysiological network, focusing on chemokine-mediated inflammation and cancer. Clarification of the molecular basis of dynamic chemokine-receptor interactions will help advance the understanding of chemokine biology to achieve precision medicine application in the clinic.

Restraint of IFN-γ expression through a distal silencer CNS-28 for tissue homeostasis

Interferon-γ (IFN-γ) is a key cytokine in response to viral or intracellular bacterial infection in mammals. While a number of enhancers are described to promote IFN-γ responses, to the best of our knowledge, no silencers for the Ifng gene have been identified. By examining H3K4me1 histone modification in naive CD4+ T cells within Ifng locus, we identified a silencer (CNS-28) that restrains Ifng expression. Mechanistically, CNS-28 maintains Ifng silence by diminishing enhancer-promoter interactions within Ifng locus in a GATA3-dependent but T-bet-independent manner. Functionally, CNS-28 restrains Ifng transcription in NK cells, CD4+ cells, and CD8+ T cells during both innate and adaptive immune responses. Moreover, CNS-28 deficiency resulted in repressed type 2 responses due to elevated IFN-γ expression, shifting Th1 and Th2 paradigm. Thus, CNS-28 activity ensures immune cell quiescence by cooperating with other regulatory cis elements within the Ifng gene locus to minimize autoimmunity.

Recent revelations and future directions using single-cell technologies in chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is a clinically and biologically heterogeneous disease with varying outcomes. In the last decade, the application of next-generation sequencing technologies has allowed extensive mapping of disease-specific genomic, epigenomic, immunogenetic, and transcriptomic signatures linked to CLL pathogenesis. These technologies have improved our understanding of the impact of tumor heterogeneity and evolution on disease outcome, although they have mostly been performed on bulk preparations of nucleic acids. As a further development, new technologies have emerged in recent years that allow high-resolution mapping at the single-cell level. These include single-cell RNA sequencing for assessment of the transcriptome, both of leukemic and non-malignant cells in the tumor microenvironment; immunogenetic profiling of B and T cell receptor rearrangements; single-cell sequencing methods for investigation of methylation and chromatin accessibility across the genome; and targeted single-cell DNA sequencing for analysis of copy-number alterations and single nucleotide variants. In addition, concomitant profiling of cellular subpopulations, based on protein expression, can also be obtained by various antibody-based approaches. In this review, we discuss different single-cell sequencing technologies and how they have been applied so far to study CLL onset and progression, also in response to treatment. This latter aspect is particularly relevant considering that we are moving away from chemoimmunotherapy to targeted therapies, with a potentially distinct impact on clonal dynamics. We also discuss new possibilities, such as integrative multi-omics analysis, as well as inherent limitations of the different single-cell technologies, from sample preparation to data interpretation using available bioinformatic pipelines. Finally, we discuss future directions in this rapidly evolving field.

Rig1 receptor plays a critical role in cardiac reprogramming via YY1 signaling

We discovered that innate immunity plays an important role in the reprogramming of fibroblasts into cardiomyocytes. In this report, we define the role of a novel retinoic acid-inducible gene 1 Yin Yang 1 (Rig1:YY1) pathway. We found that fibroblast to cardiomyocyte reprogramming efficacy was enhanced by specific Rig1 activators. To understand the mechanism of action, we performed various transcriptomic, nucleosome occupancy, and epigenomic approaches. Analysis of the datasets indicated that Rig1 agonists had no effect on reprogramming-induced changes in nucleosome occupancy or loss of inhibitory epigenetic motifs. Instead, Rig1 agonists were found to modulate cardiac reprogramming by promoting the binding of YY1 specifically to cardiac genes. To conclude, these results show that the Rig1:YY1 pathway plays a critical role in fibroblast to cardiomyocyte reprogramming.

Live imaging reveals chromatin compaction transitions and dynamic transcriptional bursting during stem cell differentiation in vivo

Stem cell differentiation requires dramatic changes in gene expression and global remodeling of chromatin architecture. How and when chromatin remodels relative to the transcriptional, behavioral, and morphological changes during differentiation remain unclear, particularly in an intact tissue context. Here, we develop a quantitative pipeline which leverages fluorescently-tagged histones and longitudinal imaging to track large-scale chromatin compaction changes within individual cells in a live mouse. Applying this pipeline to epidermal stem cells, we reveal that cell-to-cell chromatin compaction heterogeneity within the stem cell compartment emerges independent of cell cycle status, and instead is reflective of differentiation status. Chromatin compaction state gradually transitions over days as differentiating cells exit the stem cell compartment. Moreover, establishing live imaging of Keratin-10 (K10) nascent RNA, which marks the onset of stem cell differentiation, we find that Keratin-10 transcription is highly dynamic and largely precedes the global chromatin compaction changes associated with differentiation. Together, these analyses reveal that stem cell differentiation involves dynamic transcriptional states and gradual chromatin rearrangement.

Rare diseases of epigenetic origin: Challenges and opportunities

Rare diseases (RDs), more than 80% of which have a genetic origin, collectively affect approximately 350 million people worldwide. Progress in next-generation sequencing technology has both greatly accelerated the pace of discovery of novel RDs and provided more accurate means for their diagnosis. RDs that are driven by altered epigenetic regulation with an underlying genetic basis are referred to as rare diseases of epigenetic origin (RDEOs). These diseases pose unique challenges in research, as they often show complex genetic and clinical heterogeneity arising from unknown gene-disease mechanisms. Furthermore, multiple other factors, including cell type and developmental time point, can confound attempts to deconvolute the pathophysiology of these disorders. These challenges are further exacerbated by factors that contribute to epigenetic variability and the difficulty of collecting sufficient participant numbers in human studies. However, new molecular and bioinformatics techniques will provide insight into how these disorders manifest over time. This review highlights recent studies addressing these challenges with innovative solutions. Further research will elucidate the mechanisms of action underlying unique RDEOs and facilitate the discovery of treatments and diagnostic biomarkers for screening, thereby improving health trajectories and clinical outcomes of affected patients.

Factors and Methods for the Detection of Gene Expression Regulation

Gene-expression regulation involves multiple processes and a range of regulatory factors. In this review, we describe the key factors that regulate gene expression, including transcription factors (TFs), chromatin accessibility, histone modifications, DNA methylation, and RNA modifications. In addition, we also describe methods that can be used to detect these regulatory factors.

Matching queried single-cell open-chromatin profiles to large pools of single-cell transcriptomes and epigenomes for reference supported analysis

The true benefits of large single-cell transcriptome and epigenome data sets can be realized only with the development of new approaches and search tools for annotating individual cells. Matching a single-cell epigenome profile to a large pool of reference cells remains a major challenge. Here, we present scEpiSearch, which enables searching, comparison, and independent classification of single-cell open-chromatin profiles against a large reference of single-cell expression and open-chromatin data sets. Across performance benchmarks, scEpiSearch outperformed multiple methods in accuracy of search and low-dimensional coembedding of single-cell profiles, irrespective of platforms and species. Here we also demonstrate the unconventional utilities of scEpiSearch by applying it on single-cell epigenome profiles of K562 cells and samples from patients with acute leukaemia to reveal different aspects of their heterogeneity, multipotent behavior, and dedifferentiated states. Applying scEpiSearch on our single-cell open-chromatin profiles from embryonic stem cells (ESCs), we identified ESC subpopulations with more activity and poising for endoplasmic reticulum stress and unfolded protein response. Thus, scEpiSearch solves the nontrivial problem of amalgamating information from a large pool of single cells to identify and study the regulatory states of cells using their single-cell epigenomes.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100 bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Single-Cell Joint Profiling of Open Chromatin and Transcriptome by Paired-Seq

Simultaneous detection of chromatin accessibility and transcription from the same cells promises to greatly facilitate the dissection of cell-type-specific gene regulatory programs in complex tissues. Paired-seq enables joint analysis of open chromatin and nuclear transcriptome from up to a million cells in parallel. It achieves ultra-high-throughput single-cell multiomics with the use of a combinatorial barcoding strategy involving sequential ligation of multiplexed DNA barcodes to chromatin DNA fragments and reverse transcription products, followed by high-throughput DNA sequencing of the resulting DNA libraries and deconvolution of single-cell multiomic maps based on cell-specific barcodes.

Comprehensive mapping of alternative polyadenylation site usage and its dynamics at single-cell resolution

Alternative polyadenylation (APA) plays an important role in posttranscriptional gene regulation such as transcript stability and translation efficiency. However, our knowledge about APA dynamics at the single-cell level is largely unexplored. Here, we developed single-cell polyadenylation sequencing, a strand-specific approach for sequencing the 3' end of transcripts, to investigate the landscape of APA at the single-cell level. By analyzing several cell lines, we found many genes using multiple polyA sites in bulk data are prone to use only one polyA site in each single cell. Interestingly, cell cycle genes were significantly enriched in genes with high variation in polyA site usages. Furthermore, the 414 genes showing a polyA site usage switch after cell synchronization enriched cell cycle genes, while the differentially expressed genes after cell synchronization did not enrich cell cycle genes. We further identified 812 genes showing polyA site usage changes between neighboring cell cycles, which were grouped into six clusters, with cell phase-specific functional categories enriched in each cluster. Deletion of one polyA site in MSL1 and SCCPDH results in slower and faster cell cycle progression, respectively, supporting polyA site usage switch played an important role in cell cycle. These results indicate that APA is an important layer for cell cycle regulation.

A framework for clinical cancer subtyping from nucleosome profiling of cell-free DNA

Cell-free DNA (cfDNA) has the potential to inform tumor subtype classification and help guide clinical precision oncology. Here we develop Griffin, a framework for profiling nucleosome protection and accessibility from cfDNA to study the phenotype of tumors using as low as 0.1x coverage whole genome sequencing data. Griffin employs a GC correction procedure tailored to variable cfDNA fragment sizes, which generates a better representation of chromatin accessibility and improves the accuracy of cancer detection and tumor subtype classification. We demonstrate estrogen receptor subtyping from cfDNA in metastatic breast cancer. We predict estrogen receptor subtype in 139 patients with at least 5% detectable circulating tumor DNA with an area under the receive operator characteristic curve (AUC) of 0.89 and validate performance in independent cohorts (AUC = 0.96). In summary, Griffin is a framework for accurate tumor subtyping and can be generalizable to other cancer types for precision oncology applications.

Bidirectional Functional Effects of Staphylococcus on Carcinogenesis

As a Gram-positive cocci existing in nature, Staphylococcus has a variety of species, such as Staphylococcus aureus and Staphylococcus epidermidis, etc. Growing evidence reveals that Staphylococcus is closely related to the occurrence and development of various cancers. On the one hand, cancer patients are more likely to suffer from bacterial infection and antibiotic-resistant strain infection compared to healthy controls. On the other hand, there exists an association between staphylococcal infection and carcinogenesis. Staphylococcus often plays a pathogenic role and evades the host immune system through surface adhesion molecules, α-hemolysin, PVL (Panton-Valentine leukocidin), SEs (staphylococcal enterotoxins), SpA (staphylococcal protein A), TSST-1 (Toxic shock syndrom toxin-1) and other factors. Staphylococcal nucleases (SNases) are extracellular nucleases that serve as genomic markers for Staphylococcus aureus. Interestingly, a human homologue of SNases, SND1 (staphylococcal nuclease and Tudor domain-containing 1), has been recognized as an oncoprotein. This review is the first to summarize the reported basic and clinical evidence on staphylococci and neoplasms. Investigations on the correlation between Staphylococcus and the occurrence, development, diagnosis and treatment of breast, skin, oral, colon and other cancers, are made from the perspectives of various virulence factors and SND1.

Functional analysis of structural variants in single cells using Strand-seq

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Lineage tracing reveals B cell antibody class switching is stochastic, cell-autonomous, and tuneable

To optimize immunity to pathogens, B lymphocytes generate plasma cells with functionally diverse antibody isotypes. By lineage tracing single cells within differentiating B cell clones, we identified the heritability of discrete fate controlling mechanisms to inform a general mathematical model of B cell fate regulation. Founder cells highly influenced clonal plasma-cell fate, whereas class switch recombination (CSR) was variegated within clones. In turn, these CSR patterns resulted from independent all-or-none expression of both activation-induced cytidine deaminase (AID) and IgH germline transcription (GLT), with the latter being randomly re-expressed after each cell division. A stochastic model premised on these molecular transition rules accurately predicted antibody switching outcomes under varied conditions in vitro and during an immune response in vivo. Thus, the generation of functionally diverse antibody types follows rules of autonomous cellular programming that can be adapted and modeled for the rational control of antibody classes for potential therapeutic benefit.

Transcription factors TCF-1 and GATA3 are key factors for the epigenetic priming of early innate lymphoid progenitors toward distinct cell fates

The differentiation of innate lymphoid cells (ILCs) from hematopoietic stem cells needs to go through several multipotent progenitor stages. However, it remains unclear whether the fates of multipotent progenitors are predefined by epigenetic states. Here, we report the identification of distinct accessible chromatin regions in all lymphoid progenitors (ALPs), EILPs, and ILC precursors (ILCPs). Single-cell MNase-seq analyses revealed that EILPs contained distinct subpopulations epigenetically primed toward either dendritic cell lineages or ILC lineages. We found that TCF-1 and GATA3 co-bound to the lineage-defining sites for ILCs (LDS-Is), whereas PU.1 binding was enriched in the LDSs for alternative dendritic cells (LDS-As). TCF-1 and GATA3 were indispensable for the epigenetic priming of LDSs at the EILP stage. Our results suggest that the multipotency of progenitor cells is defined by the existence of a heterogeneous population of cells epigenetically primed for distinct downstream lineages, which are regulated by key transcription factors.

NucleoMap: A computational tool for identifying nucleosomes in ultra-high resolution contact maps

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.

Single cell cancer epigenetics

Bulk sequencing methodologies have allowed us to make great progress in cancer research. Unfortunately, these techniques lack the resolution to fully unravel the epigenetic mechanisms that govern tumor heterogeneity. Consequently, many novel single cell-sequencing methodologies have been developed over the past decade, allowing us to explore the epigenetic components that regulate different aspects of cancer heterogeneity, namely: clonal heterogeneity, tumor microenvironment (TME), spatial organization, intratumoral differentiation programs, metastasis, and resistance mechanisms. In this review, we explore the different sequencing techniques that enable researchers to study different aspects of epigenetics (DNA methylation, chromatin accessibility, histone modifications, DNA-protein interactions, and chromatin 3D architecture) at the single cell level, their potential applications in cancer, and their current technical limitations.

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.

Single-cell sequencing: expansion, integration and translation

With the rapid advancement in sequencing technologies, the concept of omics has revolutionized our understanding of cellular behaviors. Conventional omics investigation approaches measure the averaged behaviors of multiple cells, which may easily hide signals represented by a small-cell cohort, urging for the development of techniques with enhanced resolution. Single-cell RNA sequencing, investigating cell transcriptomics at the resolution of a single cell, has been rapidly expanded to investigate other omics such as genomics, proteomics and metabolomics since its invention. The requirement for comprehensive understanding of complex cellular behavior has led to the integration of multi-omics and single-cell sequencing data with other layers of information such as spatial data and the CRISPR screening technique towards gained knowledge or innovative functionalities. The development of single-cell sequencing in both dimensions has rendered it a unique field that offers us a versatile toolbox to delineate complex diseases, including cancers.

Nucleosome-Omics: A Perspective on the Epigenetic Code and 3D Genome Landscape

Genetic information is loaded on chromatin, which involves DNA sequence arrangement and the epigenetic landscape. The epigenetic information including DNA methylation, nucleosome positioning, histone modification, 3D chromatin conformation, and so on, has a crucial impact on gene transcriptional regulation. Out of them, nucleosomes, as basal chromatin structural units, play an important central role in epigenetic code. With the discovery of nucleosomes, various nucleosome-level technologies have been developed and applied, pushing epigenetics to a new climax. As the underlying methodology, next-generation sequencing technology has emerged and allowed scientists to understand the epigenetic landscape at a genome-wide level. Combining with NGS, nucleosome-omics (or nucleosomics) provides a fresh perspective on the epigenetic code and 3D genome landscape. Here, we summarized and discussed research progress in technology development and application of nucleosome-omics. We foresee the future directions of epigenetic development at the nucleosome level.

Histone H1 binding to nucleosome arrays depends on linker DNA length and trajectory

Throughout the genome, nucleosomes often form regular arrays that differ in nucleosome repeat length (NRL), occupancy of linker histone H1 and transcriptional activity. Here, we report cryo-EM structures of human H1-containing tetranucleosome arrays with four physiologically relevant NRLs. The structures show a zig-zag arrangement of nucleosomes, with nucleosomes 1 and 3 forming a stack. H1 binding to stacked nucleosomes depends on the NRL, whereas H1 always binds to the non-stacked nucleosomes 2 and 4. Short NRLs lead to altered trajectories of linker DNA, and these altered trajectories sterically impair H1 binding to the stacked nucleosomes in our structures. As the NRL increases, linker DNA trajectories relax, enabling H1 contacts and binding. Our results provide an explanation for why arrays with short NRLs are depleted of H1 and suited for transcription, whereas arrays with long NRLs show full H1 occupancy and can form transcriptionally silent heterochromatin regions.

Cryo-EM structures of human H1-containing tetranucleosome arrays with distinct, physiological nucleosome repeat lengths reveal that nucleosomes assume a zig-zag arrangement and H1 binds to stacked nucleosomes with longer linker DNA.

Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools—software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

Differences in nanoscale organization of regulatory active and inactive human chromatin

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.

Experimental development of the epigenomic library construction method to elucidate the epigenetic diversity and causal relationship between epigenome and transcriptome at a single-cell level

The method of single-cell RNA sequencing has been rapidly developed, and numerous experiments have been conducted over the past decade. Their results allow us to recognize various subpopulations and rare cell states in tissues, tumors, and immune systems that are previously unidentified, and guide us to understand fundamental biological processes that determine cell identity based on single-cell gene expression profiles. However, it is still challenging to understand the principle of comprehensive gene regulation that determines the cell fate only with transcriptome, a consequential output of the gene expression program. To elucidate the mechanisms related to the origin and maintenance of comprehensive single-cell transcriptome, we require a corresponding single-cell epigenome, which is a differentiated information of each cell with an identical genome. This review deals with the current development of single-cell epigenomic library construction methods, including multi-omics tools with crucial factors and additional requirements in the future focusing on DNA methylation, chromatin accessibility, and histone post-translational modifications. The study of cellular differentiation and the disease occurrence at a single-cell level has taken the first step with single-cell transcriptome and is now taking the next step with single-cell epigenome.

Computational Methods for Single-Cell Imaging and Omics Data Integration

Integrating single cell omics and single cell imaging allows for a more effective characterisation of the underlying mechanisms that drive a phenotype at the tissue level, creating a comprehensive profile at the cellular level. Although the use of imaging data is well established in biomedical research, its primary application has been to observe phenotypes at the tissue or organ level, often using medical imaging techniques such as MRI, CT, and PET. These imaging technologies complement omics-based data in biomedical research because they are helpful for identifying associations between genotype and phenotype, along with functional changes occurring at the tissue level. Single cell imaging can act as an intermediary between these levels. Meanwhile new technologies continue to arrive that can be used to interrogate the genome of single cells and its related omics datasets. As these two areas, single cell imaging and single cell omics, each advance independently with the development of novel techniques, the opportunity to integrate these data types becomes more and more attractive. This review outlines some of the technologies and methods currently available for generating, processing, and analysing single-cell omics- and imaging data, and how they could be integrated to further our understanding of complex biological phenomena like ageing. We include an emphasis on machine learning algorithms because of their ability to identify complex patterns in large multidimensional data.