Uncovering dark matter in cancer by identifying epigenetic drivers.

The complex relationship between chromatin accessibility, transcriptional regulation, and cancer transitions presents a daunting puzzle. Terekhanova et al. created a pan-cancer epigenetic and transcriptomic atlas at single-cell resolution, yielding important insights into the underlying chromatin architecture of cancer transitions and novel discoveries with the potential to advance precision medicine.

Definition of regulatory elements and transcription factors controlling porcine immune cell gene expression at single cell resolution using single nucleus ATAC-seq.

The transcriptome of porcine peripheral blood mononuclear cells (PBMC) at single cell (sc) resolution is well described, but little is understood about the cis-regulatory mechanism behind scPBMC gene expression. Here, we profiled the open chromatin landscape of porcine PBMC that define cis-regulatory elements and mechanism contributing to the transcription using single nucleus ATAC sequencing (snATAC-seq). Approximately 22 % of the identified peaks overlapped with annotated transcription start sites (TSS). Using clustering based on open chromatin pattern similarity, we demonstrate that cell type annotations using snATAC-seq are highly concordant to that reported for sc RNA sequencing (scRNA-seq). The differentially accessible peaks (DAPs) for each cell type were characterized and the pattern of accessibility of the DAPs near cell type markers across cell types was similar to that of the average gene expression level of corresponding marker genes. Additionally, we found that peaks identified in snATAC-seq have the potential power to predict the cell type specific transcription starting site (TSS). We identified both transcription factors (TFs) whose binding motif were enriched in cell type DAPs of multiple cell types and cell type specific TFs by conducting transcription factor binding motif (TFBM) analysis. Furthermore, we identified the putative enhancer or promoter regions bound by TFs for each differentially expressed gene (DEG) with a DAP that overlapped with its TSS by generating cis-co-accessibility networks (CCAN). To predict the regulators of such DEGs, TFBM analysis was performed for each CCAN. The regulator TF-target DEG pairs predicted in this way were largely consistent with the results reported in the ENCODE Transcription Factor Targets Dataset (TFTD). This snATAC-seq approach provides insights into the regulation of chromatin accessibility landscape of porcine PBMCs and enables discovery of TFs predicted to control DEG through binding regulatory elements whose chromatin accessibility correlates with the DEG promoter region.

Chromatin accessibility analysis and architectural profiling of human kidneys reveal key cell types and a regulator of diabetic kidney disease.

Diabetes is the leading cause of kidney disease that progresses to kidney failure. However, the key molecular and cellular pathways involved in diabetic kidney disease (DKD) pathogenesis are largely unknown. Here, we performed a comparative analysis of adult human kidneys by examining cell type-specific chromatin accessibility by single-nucleus ATAC-seq (snATAC-seq) and analyzing three-dimensional chromatin architecture via high-throughput chromosome conformation capture (Hi-C method) of paired samples. We mapped the cell type-specific and DKD-specific open chromatin landscape and found that genetic variants associated with kidney diseases were significantly enriched in the proximal tubule- (PT) and injured PT-specific open chromatin regions in samples from patients with DKD. BACH1 was identified as a core transcription factor of injured PT cells; its binding target genes were highly associated with fibrosis and inflammation, which were also key features of injured PT cells. Additionally, Hi-C analysis revealed global chromatin architectural changes in DKD, accompanied by changes in local open chromatin patterns. Combining the snATAC-seq and Hi-C data identified direct target genes of BACH1, and indicated that BACH1 binding regions showed increased chromatin contact frequency with promoters of their target genes in DKD. Thus, our multi-omics analysis revealed BACH1 target genes in injured PTs and highlighted the role of BACH1 as a novel regulator of tubular inflammation and fibrosis.

Single-nucleus chromatin accessibility profiling highlights distinct astrocyte signatures in progressive supranuclear palsy and corticobasal degeneration.

Tauopathies such as progressive supranuclear palsy (PSP) and corticobasal degeneration (CBD) exhibit characteristic neuronal and glial inclusions of hyperphosphorylated Tau (pTau). Although the astrocytic pTau phenotype upon neuropathological examination is the most guiding feature in distinguishing both diseases, regulatory mechanisms controlling their transitions into disease-specific states are poorly understood to date. Here, we provide accessible chromatin data of more than 45,000 single nuclei isolated from the frontal cortex of PSP, CBD, and control individuals. We found a strong association of disease-relevant molecular changes with astrocytes and demonstrate that tauopathy-relevant genetic risk variants are tightly linked to astrocytic chromatin accessibility profiles in the brains of PSP and CBD patients. Unlike the established pathogenesis in the secondary tauopathy Alzheimer disease, microglial alterations were relatively sparse. Transcription factor (TF) motif enrichments in pseudotime as well as modeling of the astrocytic TF interplay suggested a common pTau signature for CBD and PSP that is reminiscent of an inflammatory immediate-early response. Nonetheless, machine learning models also predicted discriminatory features, and we observed marked differences in molecular entities related to protein homeostasis between both diseases. Predicted TF involvement was supported by immunofluorescence analyses in postmortem brain tissue for their highly correlated target genes. Collectively, our data expand the current knowledge on risk gene involvement (e.g., MAPT, MAPK8, and NFE2L2) and molecular pathways leading to the phenotypic changes associated with CBD and PSP.

Transformation of brain myeloid cell populations by SIV in rhesus macaques revealed by multiomics.

The primary immune constituents in the brain, microglia and macrophages, are the target for HIV in people and simian immunodeficiency virus (SIV) in nonhuman primates. This infection can lead to neurological dysfunction, known as HIV-associated neurocognitive disorder (HAND). Given the gaps in our knowledge on how these cells respond in vivo to CNS infection, we performed single-cell multiomic sequencing, including gene expression and ATAC-seq, on myeloid cells from the brains of rhesus macaques with SIV-induced encephalitis (SIVE) as well as uninfected controls. We found that the myeloid cell populations were significantly changed by SIVE. In SIVE microglia-like cells express high levels of chemoattractants capable of recruiting highly activated CAM-like cells to the site of infection/inflammation. A unique population of microglia-like cells was found in which the chromatin accessibility of genes diverged from their RNA expression. Additionally, we observed a dramatic shift of upstream gene regulators and their targets in brain myeloid cells during SIVE. In summary, this study further uncovers the transcriptome, gene regulatory events and potential roles of different brain myeloid phenotypes in SIVE.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in ß cell function and the pathogenesis of type 2 diabetes.

Type 2 and type 1 diabetes (T2D, T1D) exhibit sex differences in insulin secretion, the mechanisms of which are unknown. We examined sex differences in human pancreatic islets from 52 donors with and without T2D combining single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), hormone secretion, and bioenergetics. In nondiabetic (ND) donors, sex differences in islet cells gene accessibility and expression predominantly involved sex chromosomes. Islets from T2D donors exhibited similar sex differences in sex chromosomes differentially expressed genes (DEGs), but also exhibited sex differences in autosomal genes. Comparing ß cells from T2D vs. ND donors, gene enrichment of female ß cells showed suppression in mitochondrial respiration, while male ß cells exhibited suppressed insulin secretion. Thus, although sex differences in gene accessibility and expression of ND ß cells predominantly affect sex chromosomes, the transition to T2D reveals sex differences in autosomes highlighting mitochondrial failure in females.

Single-nucleus multi-omics of Parkinson's disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks.

The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

Characterization and Optimization of Multiomic Single-Cell Epigenomic Profiling.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n = 18) and a solid tumor type, ovarian cancer (OC, n = 18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell type identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n = 6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Single nucleus multiome analysis of the prefrontal cortex from C9orf72 ALS/FTD patients illuminates pathways affected during disease progression.

Repeat expansions in the C9orf72 gene are the most common genetic cause of amyotrophic lateral sclerosis and familial frontotemporal dementia (ALS/FTD). To identify molecular defects that take place in the dorsolateral frontal cortex of patients with C9orf72 ALS/FTD, we compared healthy controls with C9orf72 ALS/FTD donor samples staged based on the levels of cortical phosphorylated TAR DNA binding protein (pTDP-43), a neuropathological hallmark of disease progression. We identified distinct molecular changes in different cell types that take place during disease progression. These alterations include downregulation of nuclear and mitochondrial ribosomal protein genes in early disease stages that become upregulated as the disease progresses. High ratios of premature oligodendrocytes expressing low levels of genes encoding major myelin protein components are characteristic of late disease stages and may represent a unique signature of C9orf72 ALS/FTD. Microglia with increased reactivity and astrocyte specific transcriptome changes in genes involved in glucose/glycogen metabolism are also associated with disease progression. Late stages of C9orf72 ALS/FTD correlate with sequential changes in the regulatory landscape of several genes in glial cells, namely MBP/MAG/MOG in oligodendrocytes, CD83/IRF8 in microglia, and GLUT1/GYS2/AGL in astrocytes. Only layer 2-3 cortical projection neurons with high expression of CUX2/LAMP5 are significantly reduced in C9orf72 ALS/FTD patients with respect to controls. Our findings reveal previously unknown progressive functional changes in cortical cells of C9orf72 ALS/FTD patients that shed light on the mechanisms underlying the pathology of this disease.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains.

BACKGROUND: The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD. METHODS: Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs). RESULTS: Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes. CONCLUSIONS: To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

Tumor heterogeneity and tumor-microglia interactions in primary and recurrent IDH1-mutant gliomas.

The isocitrate dehydrogenase (IDH) gene is recurrently mutated in adult diffuse gliomas. IDH-mutant gliomas are categorized into oligodendrogliomas and astrocytomas, each with unique pathological features. Here, we use single-nucleus RNA and ATAC sequencing to compare the molecular heterogeneity of these glioma subtypes. In addition to astrocyte-like, oligodendrocyte progenitor-like, and cycling tumor subpopulations, a tumor population enriched for ribosomal genes and translation elongation factors is primarily present in oligodendrogliomas. Longitudinal analysis of astrocytomas indicates that the proportion of tumor subpopulations remains stable in recurrent tumors. Analysis of tumor-associated microglia/macrophages (TAMs) reveals significant differences between oligodendrogliomas, with astrocytomas harboring inflammatory TAMs expressing phosphorylated STAT1, as confirmed by immunohistochemistry. Furthermore, inferred receptor-ligand interactions between tumor subpopulations and TAMs may contribute to TAM state diversity. Overall, our study sheds light on distinct tumor populations, TAM heterogeneity, TAM-tumor interactions in IDH-mutant glioma subtypes, and the relative stability of tumor subpopulations in recurrent astrocytomas.

Multi-omics analysis in developmental bone biology.

Single-cell omics and multi-omics have revolutionized our understanding of molecular and cellular biological processes at a single-cell level. In bone biology, the combination of single-cell RNA-sequencing analyses and in vivo lineage-tracing approaches has successfully identified multi-cellular diversity and dynamics of skeletal cells. This established a new concept that bone growth and regeneration are regulated by concerted actions of multiple types of skeletal stem cells, which reside in spatiotemporally distinct niches. One important subtype is endosteal stem cells that are particularly abundant in young bone marrow. The discovery of this new skeletal stem cell type has been facilitated by single-cell multi-omics, which simultaneously measures gene expression and chromatin accessibility. Using single-cell omics, it is now possible to computationally predict the immediate future state of individual cells and their differentiation potential. In vivo validation using histological approaches is the key to interpret the computational prediction. The emerging spatial omics, such as spatial transcriptomics and epigenomics, have major advantage in retaining the location of individual cells within highly complex tissue architecture. Spatial omics can be integrated with other omics to further obtain in-depth insights. Single-cell multi-omics are now becoming an essential tool to unravel intricate multicellular dynamics and intercellular interactions of skeletal cells.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise.

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

Deciphering the distinct transcriptomic and gene regulatory map in adult macaque basal ganglia cells.

BACKGROUND: The basal ganglia are a complex of interconnected subcortical structures located beneath the mammalian cerebral cortex. The degeneration of dopaminergic neurons in the basal ganglia is the primary pathological feature of Parkinson's disease. Due to a lack of integrated analysis of multiomics datasets across multiple basal ganglia brain regions, very little is known about the regulatory mechanisms of this area. FINDINGS: We utilized high-throughput transcriptomic and epigenomic analysis to profile over 270,000 single-nucleus cells to create a cellular atlas of the basal ganglia, characterizing the cellular composition of 4 regions of basal ganglia in adult macaque brain, including the striatum, substantia nigra (SN), globus pallidum, and amygdala. We found a distinct epigenetic regulation on gene expression of neuronal and nonneuronal cells across regions in basal ganglia. We identified a cluster of SN-specific astrocytes associated with neurodegenerative diseases and further explored the conserved and primate-specific transcriptomics in SN cell types across human, macaque, and mouse. Finally, we integrated our epigenetic landscape of basal ganglia cells with human disease heritability and identified a regulatory module consisting of candidate cis-regulatory elements that are specific to medium spiny neurons and associated with schizophrenia. CONCLUSIONS: In general, our macaque basal ganglia atlas provides valuable insights into the comprehensive transcriptome and epigenome of the most important and populous cell populations in the macaque basal ganglia. We have identified 49 cell types based on transcriptomic profiles and 47 cell types based on epigenomic profiles, some of which exhibit region specificity, and characterized the molecular relationships underlying these brain regions.

AMULET: a novel read count-based method for effective multiplet detection from single nucleus ATAC-seq data.

Detecting multiplets in single nucleus (sn)ATAC-seq data is challenging due to data sparsity and limited dynamic range. AMULET (ATAC-seq MULtiplet Estimation Tool) enumerates regions with greater than two uniquely aligned reads across the genome to effectively detect multiplets. We evaluate the method by generating snATAC-seq data in the human blood and pancreatic islet samples. AMULET has high precision, estimated via donor-based multiplexing, and high recall, estimated via simulated multiplets, compared to alternatives and identifies multiplets most effectively when a certain read depth of 25K median valid reads per nucleus is achieved.

Single-Nucleus Chromatin Accessibility Landscape Reveals Diversity in Regulatory Regions Across Distinct Adult Rat Cortex.

Rats have been widely used as an experimental organism in psychological, pharmacological, and behavioral studies by modeling human diseases such as neurological disorders. It is critical to identify and characterize cell fate determinants and their regulatory mechanisms in single-cell resolutions across rat brain regions. Here, we applied droplet-based single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq) to systematically profile the single-cell chromatin accessibility across four dissected brain areas in adult Sprague-Dawley (SD) rats with a total of 59,023 single nuclei and identified 16 distinct cell types. Interestingly, we found that different cortex regions exhibit diversity in both cellular compositions and gene regulatory regions. Several cell-type-specific transcription factors (TFs), including SPI1, KLF4, KLF6, and NEUROD2, have been shown to play important roles during the pathogenesis of various neurological diseases, such as Alzheimer's disease (AD), astrocytic gliomas, autism spectrum disorder (ASD), and intellectual disabilities. Therefore, our single-nucleus atlas of rat cortex could serve as an invaluable resource for dissecting the regulatory mechanisms underlying diverse cortex cell fates and further revealing the regulatory networks of neuropathogenesis.

Single-Cell Chromatin Analysis of Mammary Gland Development Reveals Cell-State Transcriptional Regulators and Lineage Relationships.

Technological improvements enable single-cell epigenetic analyses of organ development. We reasoned that high-resolution single-cell chromatin accessibility mapping would provide needed insight into the epigenetic reprogramming and transcriptional regulators involved in normal mammary gland development. Here, we provide a single-cell resource of chromatin accessibility for murine mammary development from the peak of fetal mammary stem cell (fMaSC) functional activity in late embryogenesis to the differentiation of adult basal and luminal cells. We find that the chromatin landscape within individual cells predicts both gene accessibility and transcription factor activity. The ability of single-cell chromatin profiling to separate E18 fetal mammary cells into clusters exhibiting basal-like and luminal-like chromatin features is noteworthy. Such distinctions were not evident in analyses of droplet-based single-cell transcriptomic data. We present a web application as a scientific resource for facilitating future analyses of the gene regulatory networks involved in mammary development.