Assay For Transposase-accessible Chromatin Sequencing Research Articles

Spatial multiomic landscape of the human placenta at molecular resolution.

Successful pregnancy relies directly on the placenta's complex, dynamic, gene-regulatory networks. Disruption of this vast collection of intercellular and intracellular programs leads to pregnancy complications and developmental defects. In the present study, we generated a comprehensive, spatially resolved, multimodal cell census elucidating the molecular architecture of the first trimester human placenta. We utilized paired single-nucleus (sn)ATAC (assay for transposase accessible chromatin) sequencing and RNA sequencing (RNA-seq), spatial snATAC-seq and RNA-seq, and in situ sequencing and hybridization mapping of transcriptomes at molecular resolution to spatially reconstruct the joint epigenomic and transcriptomic regulatory landscape. Paired analyses unraveled intricate tumor-like gene expression and transcription factor motif programs potentially sustaining the placenta in a hostile uterine environment; further investigation of gene-linked cis-regulatory elements revealed heightened regulatory complexity that may govern trophoblast differentiation and placental disease risk. Complementary spatial mapping techniques decoded these programs within the placental villous core and extravillous trophoblast cell column architecture while simultaneously revealing niche-establishing transcriptional elements and cell-cell communication. Finally, we computationally imputed genome-wide, multiomic single-cell profiles and spatially characterized the placental chromatin accessibility landscape. This spatially resolved, single-cell multiomic framework of the first trimester human placenta serves as a blueprint for future studies on early placental development and pregnancy.

Rewired chromatin structure and epigenetic gene dysregulation during HTLV-1 infection to leukemogenesis.

Human T-cell leukemia virus type 1 (HTLV-1) broadly impacts host genes, affecting the infected cell population and inducing the development of a disease with a poor prognosis, adult T-cell leukemia-lymphoma (ATL). This study aimed to provide a comprehensive epigenomic characterization of the infected cell population and evaluated the transcriptome and chromatin structures of peripheral blood cells in HTLV-1-infected individuals using RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin sequencing (ATAC-seq). The infected cells showed significant changes in gene expression patterns from the polyclonal stage and before ATL onset while demonstrating similarities to tumor-forming ATL cells. These similarities were a result of large-scale open chromatin changes, supporting the independent early formation of epigenomic aberrations as an underlying mechanism for later clonal propagation. This study also demonstrated that HTLV-1 Tax directly affects the host chromatin structure, thereby developing fundamental epigenomic characteristics. Several Tax target genes, including the RASGRP3-ERK pathway, were recognized, indicating an impact on signaling pathways. This genome-wide variability in chromatin structural property is a novel feature of HTLV-1 infection and may contribute to pathogenic mechanisms. In addition, it has crucial implications for better understanding the impact of HTLV-1 on the host genome and identifying novel therapeutic targets.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development.

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6weeks; n= 10) and mouse (post-natal day [P]0; n= 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

Single cell multi-omics analysis unveils a distinct mechanism of exercise-mediated cardioprotection in ischemic cardiomyopathy

Abstract Background Exercise is pivotal in heart disease prevention, with proven cardioprotective effects in clinical trials. However, its precise molecular mechanisms remain elusive. Recent advancements in single-cell omics have revealed cellular heterogeneity within heart tissues, offering a deeper understanding of the pathophysiology of cardiac diseases. Purpose To investigate the exercise-induced physiological changes in heart tissue at the single-cell level and elucidate the molecular mechanisms underlying its cardioprotective effect on ischemic cardiomyopathy. Methods Murine myocardial infarction (MI) model was established by ligating the left ascending coronary artery, and a six-week voluntary wheel running exercise regimen was implemented to investigate the effect of aerobic exercise on ischemic cardiomyopathy. Integrated multi-omics analyses of single-cell DNA methylation profiling, single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) and single-nucleus RNA sequencing (snRNA-seq) were conducted to clarify the molecular mechanism of exercise-mediated cardioprotection at single-cell resolution. Overexpression and short hairpin RNA-mediated knockdown of gene X were conducted for murine MI model using cardiomyocyte-specific Myo-AAV2a vector. snRNA-seq of heart tissues from heart failure patients was performed to validate the clinical significance of the expression of gene X in cardiomyocytes. Results Exercise attenuates pathological cardiac remodeling, inducing physiological hypertrophy in non-infarcted areas while suppressing pathological remodeling of cardiomyocytes in infarct border zones. Single-cell epigenetic analysis unveiled intricate molecular mechanisms underlying exercise-induced transcriptome changes, providing insights into the gene regulatory networks modulated by exercise in the heart. In addition, integrated with snRNA-seq, we identified exercise-specific cardiomyocyte clusters, highlighting the substantial role of gene X in the molecular mechanism of exercise. Functional studies confirmed the essential role of gene X in exercise-induced cardioprotection, as its overexpression suppressed the adverse remodeling after MI and knockdown abolished the beneficial effects of exercise. Furthermore, snRNA-seq of heart failure patients revealed that higher expression of gene X was associated with ventricular reverse remodeling and favorable clinical outcomes, underscoring its significance in the pathophysiology of patients with heart failure. Conclusion Our study elucidates exercise-specific gene expression regulation networks in ischemic cardiomyopathy by unraveling the cellular heterogeneity and molecular changes associated with exercise. We identified gene X as a central player in exercise-mediated cardioprotection and further confirmed its significance in clinical settings, providing a potential therapeutic target for heart failure.

Integrated temporal transcriptional and epigenetic single-cell analysis reveals the intrarenal immune characteristics in an early-stage model of IgA nephropathy during its acute injury.

Kidney inflammation plays a crucial role in the pathogenesis of IgA nephropathy (IgAN), yet the specific phenotypes of immune cells involved in disease progression remain incompletely understood. Utilizing joint profiling through longitudinal single-cell RNA-sequencing (scRNAseq) and single-cell assay for transposase-accessible chromatin sequencing (scATACseq) can provide a comprehensive framework for elucidating the development of cell subset diversity and how chromatin accessibility regulates transcription. We aimed to characterize the dynamic immune cellular landscape at a high resolution in an early IgAN mouse model with acute kidney injury (AKI). A murine model was utilized to mimic 3 immunological states -"immune stability (IS), immune activation (IA) and immune remission (IR)" in early human IgAN-associated glomerulopathy during AKI, achieved through lipopolysaccharide (LPS) injection. Urinary albumin to creatinine ratio (UACR) was measured to further validate the exacerbation and resolution of kidney inflammation during this course. Paired scRNAseq and scATACseq analysis was performed on CD45+ immune cells isolated from kidney tissues obtained from CTRL (healthy vehicle), IS, IA and IR (4 or 5 mice each). The analyses revealed 7 major cell types and 24 clusters based on 72304 single-cell transcriptomes, allowing for the identification and characterization of various immune cell types within each cluster. Our data offer an impartial depiction of the immunological characteristics, as the proportions of immune cell types fluctuated throughout different stages of the disease. Specifically, these analyses also revealed novel subpopulations, such as a macrophage subset (Nlrp1b Mac) with distinct epigenetic features and a unique transcription factor motif profile, potentially exerting immunoregulatory effects, as well as an early subset of Tex distinguished by their effector and cytolytic potential (CX3CR1-transTeff). Furthermore, in order to investigate the potential interaction between immune cells and renal resident cells, we conducted single-cell RNA sequencing on kidney cells obtained from a separate cohort of IS and IA mice without isolating immune cells. These findings underscored the diverse roles played by macrophages and CD8+ T cells in maintaining homeostasis of endothelial cells (ECs) under stress. This study presents a comprehensive analysis of the dynamic changes in immune cell profiles in a model of IgAN, identifying key cell types and their roles and interactions. These findings significantly contribute to the understanding of the pathogenesis of IgAN and may provide potential targets for therapeutic intervention.

Single-cell analysis of the multiple myeloma microenvironment after γ-secretase inhibition and CAR T-cell therapy

AbstractChimeric antigen receptor (CAR) T cells and bispecific antibodies targeting B-cell maturation antigen (BCMA) have significantly advanced the treatment of relapsed and refractory multiple myeloma. Resistance to BCMA-targeting therapies, nonetheless, remains a significant challenge. BCMA shedding by γ-secretase is a known resistance mechanism, and preclinical studies suggest that inhibition may improve anti-BCMA therapy. Leveraging a phase 1 clinical trial of the γ-secretase inhibitor (GSI), crenigacestat, with anti-BCMA CAR T cells (FCARH143), we used single-nuclei RNA sequencing and assay for transposase-accessible chromatin sequencing to characterize the effects of GSI on the tumor microenvironment. The most significant impacts of GSI involved effects on monocytes, which are known to promote tumor growth. In addition to observing a reduction in the frequency of nonclassical monocytes, we also detected significant changes in gene expression, chromatin accessibility, and inferred cell-cell interactions after exposure to GSI. Although many genes with altered expression are associated with γ-secretase–dependent signaling, such as Notch, other pathways were affected, indicating GSI has far-reaching effects. Finally, we detected monoallelic deletion of the BCMA locus in some patients with prior exposure to anti-BCMA therapy, which significantly correlated with reduced progression-free survival (PFS; median PFS, 57 vs 861 days). GSIs are being explored in combination with the full spectrum of BCMA-targeting agents, and our results reveal widespread effects of GSI on both tumor and immune cell populations, providing insight into mechanisms for enhancing BCMA-directed therapies.

Robust estimation of cancer and immune cell-type proportions from bulk tumor ATAC-Seq data.

Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq) is a widely used technique to explore gene regulatory mechanisms. For most ATAC-Seq data from healthy and diseased tissues such as tumors, chromatin accessibility measurement represents a mixed signal from multiple cell types. In this work, we derive reliable chromatin accessibility marker peaks and reference profiles for most non-malignant cell types frequently observed in the microenvironment of human tumors. We then integrate these data into the EPIC deconvolution framework (Racle et al., 2017) to quantify cell-type heterogeneity in bulk ATAC-Seq data. Our EPIC-ATAC tool accurately predicts non-malignant and malignant cell fractions in tumor samples. When applied to a human breast cancer cohort, EPIC-ATAC accurately infers the immune contexture of the main breast cancer subtypes.

Robust estimation of cancer and immune cell-type proportions from bulk tumor ATAC-Seq data

Assay for Transposase-Accessible Chromatin sequencing (ATAC-Seq) is a widely used technique to explore gene regulatory mechanisms. For most ATAC-Seq data from healthy and diseased tissues such as tumors, chromatin accessibility measurement represents a mixed signal from multiple cell types. In this work, we derive reliable chromatin accessibility marker peaks and reference profiles for most non-malignant cell types frequently observed in the microenvironment of human tumors. We then integrate these data into the EPIC deconvolution framework (Racle et al., 2017) to quantify cell-type heterogeneity in bulk ATAC-Seq data. Our EPIC-ATAC tool accurately predicts non-malignant and malignant cell fractions in tumor samples. When applied to a human breast cancer cohort, EPIC-ATAC accurately infers the immune contexture of the main breast cancer subtypes.

6524 Sex-Specific Single Cell Genetic Architecture Of Human Islets With And Without Type 2 Diabetes

Abstract Disclosure: M. Qadir: None. R.M. Elgamal: None. P. Kudtarkar: None. K.J. Gaulton: None. F. Mauvais-Jarvis: None. Type 2 diabetes (T2D) is a heterogeneous disease and biological sex affects its pathogenesis including the development of adiposity, insulin resistance, and β cell failure. For example, ketosis-prone diabetes is a phenotypically-defined form of T2D with acute β cell failure and severe insulin deficiency, predominantly observed in men. However, what mechanisms drive male predominance in β cell failure is unknown. To fully understand the sex-based pathogenesis of islet failure in diabetes, the study of sex differences in human islet biology and dysfunction is essential. To decipher sex differences in human islets genetic architecture and function, we performed an orthogonal series of experiments from 15 islets of non-diabetic (ND) donors, using single cell RNAseq (scRNAseq) and single nucleus assay for transposase-accessible chromatin sequencing (snATACseq). We also leveraged 37 ND and T2D donors for scRNAseq data from the human pancreas analysis program (HPAP) dataset. We studied over 200,000 single cells spanning the accessible genome and transcriptome. It is the largest study of sex-based single-cell genomics and transcriptomics of human islets to date. We observe that in cultured human ND islets, in the absence of the in vivo hormonal environment, islet cell transcriptome, and genome accessibility across coding regions are predominantly limited to sex chromosome genes. Of particular interest are sex differences in gene accessibility and expression of the X-linked KDM6A and Y-linked KDM5D demethylase genes in female (F) and male (M) islet cells respectively. Furthermore, analysis of T2D islets revealed a greater number of autosomal differentially expressed genes (DEGs) in F (721 upregulated and 1164 downregulated, FDR < 0.01) compared to M β cells (111 upregulated and 99 downregulated, FDR < 0.01). M β cell DEGs over-representation analysis revealed suppressed insulin secretion and processing pathways, while F β cell DEGs exhibited suppression of oxidative phosphorylation and electron transport chain pathways. Thus, although, sex differences in gene expression and accessibility of cultured ND human islets are restricted to sex chromosome genes, during the transition from ND to T2D, major sex differences in autosomal DEGs appear, with a greater number of autosomal DEGs in F compared to M T2D islets. Presentation: 6/3/2024

Dissecting the Single-Cell Diversity and Heterogeneity Underlying Cervical Precancerous Lesions and Cancer Tissues.

The underlying cellular diversity and heterogeneity from cervix precancerous lesions to cervical squamous cell carcinoma (CSCC) is investigated. Four single-cell datasets including normal tissues, normal adjacent tissues, precancerous lesions, and cervical tumors were integrated to perform disease stage analysis. Single-cell compositional data analysis (scCODA) was utilized to reveal the compositional changes of each cell type. Differentially expressed genes (DEGs) among cell types were annotated using BioCarta. An assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis was performed to correlate epigenetic alterations with gene expression profiles. Lastly, a logistic regression model was used to assess the similarity between the original and new cohort data (HRA001742). After global annotation, seven distinct cell types were categorized. Eight consensus-upregulated DEGs were identified in B cells among different disease statuses, which could be utilized to predict the overall survival of CSCC patients. Inferred copy number variation (CNV) analysis of epithelial cells guided disease progression classification. Trajectory and ATAC-seq integration analysis identified 95 key transcription factors (TF) and one immunohistochemistry (IHC) testified key-node TF (YY1) involved in epithelial cells from CSCC initiation to progression. The consistency of epithelial cell subpopulation markers was revealed with single-cell sequencing, bulk sequencing, and RT-qPCR detection. KRT8 and KRT15, markers of Epi6, showed progressively higher expression with disease progression as revealed by IHC detection. The logistic regression model testified the robustness of the resemblance of clusters among the various datasets utilized in this study. Valuable insights into CSCC cellular diversity and heterogeneity provide a foundation for future targeted therapy.

Excitatory neurons and oligodendrocyte precursor cells are vulnerable to focal cortical dysplasia type IIIa as suggested by single-nucleus multiomics.

Focal cortical dysplasia (FCD) is a heterogeneous group of cortical developmental malformations that constitute a common cause of medically intractable epilepsy. FCD type IIIa (FCD IIIa) refers to temporal neocortex alterations in architectural organisation or cytoarchitectural composition in the immediate vicinity of hippocampal sclerosis. Slight alterations in the temporal neocortex of FCD IIIa patients pose a challenge for the preoperative diagnosis and definition of the resection range. We have performed multimodal integration of single-nucleus RNA sequencing and single-nucleus assay for transposase-accessible chromatin sequencing in the epileptogenic cortex of four patients with FCD IIIa, and three relatively normal temporal neocortex were chosen as controls. Our study revealed that the most significant dysregulation occurred in excitatory neurons (ENs) and oligodendrocyte precursor cells (OPCs) in the epileptogenic cortex of FCD IIIa patients. In ENs, we constructed a transcription factor (TF)-hub gene regulatory network and found DAB1high ENs subpopulation mediates neuronal immunity characteristically in FCD IIIa. Western blotting and immunofluorescence were used to validate the changes in protein expression levels caused by some of the key genes. The OPCs were activated and exhibited aberrant phenotypes in FCD IIIa, and TFs regulating reconstructed pseudotime trajectory were identified. Finally, our results revealed aberrant intercellular communication between ENs and OPCs in FCD IIIa patients. Our study revealed significant and intricate alterations in the transcriptomes and epigenomes in ENs and OPCs of FCD IIIa patients, shedding light on their cell type-specific regulation and potential pathogenic involvement in this disorder. This work will help evaluate the pathogenesis of cortical dysplasia and epilepsy and explore potential therapeutic targets. Paired snRNA-seq and snATAC-seq data were intergrated and analysed to identify crucial subpopulations of ENs and OPCs in the epileptogenic cortex of FCD IIIa patients and explore their possible pathogenic role in the disease. A TF-hub gene regulatory network was constructed in ENs, and the DAB1high Ex-1 mediated neuronal immunity was characterstically in FCD IIIa patients. The OPCs were activated and exhibited aberrant phenotypes in FCD IIIa patients, and TFs regulating reconstructed pseudotime traectory were identified. Aberrant intercelluar communications between ENs and OPCs in FCD IIIa patients were identified.

Analysis of chromatin accessibility in peripheral blood mononuclear cells from patients with early-stage breast cancer.

Objective: This study was aimed at exploring a specific open region of chromatin in the peripheral blood mononuclear cells (PBMCs) of patients with breast cancer and evaluating its feasibility as a biomarker for diagnosing and predicting breast cancer prognosis. Methods: We obtained PBMCs from breast cancer patients and healthy people for the assay for transposase-accessible chromatin (ATAC) sequencing (n = 3) and obtained the GSE27562 chip sequencing data for secondary analyses. Through bioinformatics analysis, we mined the pattern changes for chromatin accessibility in the PBMCs of breast cancer patients. Results: A total of 1,906 differentially accessible regions (DARs) and 1,632 differentially expressed genes (DEGs) were identified via ATAC sequencing. The upregulated DEGs in the disease group were mainly distributed in the cells, organelles, and cell-intima-related structures and were mainly responsible for biological functions such as cell nitrogen complex metabolism, macromolecular metabolism, and cell communication, in addition to functions such as nucleic acid binding, enzyme binding, hydrolase reaction, and transferase activity. Combined with microarray data analysis, the following set of nine DEGs showed intersection between the ATAC and microarray data: JUN, MSL2, CDC42, TRIB1, SERTAD3, RAB14, RHOB, RAB40B, and PRKDC. HOMER predicted and identified five transcription factors that could potentially bind to these peak sites, namely NFY, Sp 2, GFY, NRF, and ELK 1. Conclusion: Chromatin accessibility analysis of the PBMCs from patients with early-stage breast cancer underscores its potential as a significant avenue for biomarker discovery in breast cancer diagnostics and treatment. By screening the transcription factors and DEGs related to breast cancer, this study provides a comprehensive theoretical foundation that is expected to guide future clinical applications and therapeutic developments.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function.

Kidney failure, the decrease of kidney function below a threshold necessary to support life, is a major cause of morbidity and mortality. We performed a genome-wide association study (GWAS) of 406,504 individuals in the UK Biobank, identifying 430 loci affecting kidney function in middle-aged adults. To investigate the cell types affected by these loci, we integrated the GWAS with human kidney candidate cis-regulatory elements (cCREs) identified using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq). Overall, 56% of kidney function heritability localized to kidney tubule epithelial cCREs and an additional 7% to kidney podocyte cCREs. Thus, most heritable differences in adult kidney function are a result of altered gene expression in these two cell types. Using enhancer assays, allele-specific scATAC-seq and machine learning, we found that many kidney function variants alter tubule epithelial cCRE chromatin accessibility and function. Using CRISPRi, we determined which genes some of these cCREs regulate, implicating NDRG1, CCNB1 and STC1 in human kidney function.

BACH2 regulates diversification of regulatory and proinflammatory chromatin states in TH17 cells.

Interleukin-17 (IL-17)-producing helper T (TH17) cells are heterogenous and consist of nonpathogenic TH17 (npTH17) cells that contribute to tissue homeostasis and pathogenic TH17 (pTH17) cells that mediate tissue inflammation. Here, we characterize regulatory pathways underlying TH17 heterogeneity and discover substantial differences in the chromatin landscape of npTH17 and pTH17 cells both in vitro and in vivo. Compared to other CD4+ T cell subsets, npTH17 cells share accessible chromatin configurations with regulatory T cells, whereas pTH17 cells exhibit features of both npTH17 cells and type 1 helper T (TH1) cells. Integrating single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq), we infer self-reinforcing and mutually exclusive regulatory networks controlling different cell states and predicted transcription factors regulating TH17 cell pathogenicity. We validate that BACH2 promotes immunomodulatory npTH17 programs and restrains proinflammatory TH1-like programs in TH17 cells in vitro and in vivo. Furthermore, human genetics implicate BACH2 in multiple sclerosis. Overall, our work identifies regulators of TH17 heterogeneity as potential targets to mitigate autoimmunity.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish

BackgroundCis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq).ResultsUsing the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions.ConclusionsThus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

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.

A single-cell chromatin accessibility dataset of human primed and naïve pluripotent stem cell-derived teratoma

Teratoma, due to its remarkable ability to differentiate into multiple cell lineages, is a valuable model for studying human embryonic development. The similarity of the gene expression and chromatin accessibility patterns in these cells to those observed in vivo further underscores its potential as a research tool. Notably, teratomas derived from human naïve (pre-implantation epiblast-like) pluripotent stem cells (PSCs) have larger embryonic cell diversity and contain extraembryonic lineages, making them more suitable to study developmental processes. However, the cell type-specific epigenetic profiles of naïve PSC teratomas have not been yet characterized. Using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we analyzed 66,384 cell profiles from five teratomas derived from human naïve PSCs and their post-implantation epiblast-like (primed) counterparts. We observed 17 distinct cell types from both embryonic and extraembryonic lineages, resembling the corresponding cell types in human fetal tissues. Additionally, we identified key transcription factors specific to different cell types. Our dataset provides a resource for investigating gene regulatory programs in a relevant model of human embryonic development.

Detecting Small Cell Transformation in Patients with Advanced EGFR Mutant Lung Adenocarcinoma through Epigenomic cfDNA Profiling.

Histologic transformation to small cell lung cancer (SCLC) is a mechanism of treatment resistance in patients with advanced oncogene-driven lung adenocarcinoma (LUAD) that currently requires histologic review for diagnosis. Herein, we sought to develop an epigenomic cell-free DNA (cfDNA)-based approach to noninvasively detect small cell transformation in patients with EGFR mutant (EGFRm) LUAD. To characterize the epigenomic landscape of transformed (t)SCLC relative to LUAD and de novo SCLC, we performed chromatin immunoprecipitation sequencing (ChIP-seq) to profile the histone modifications H3K27ac, H3K4me3, and H3K27me3; methylated DNA immunoprecipitation sequencing (MeDIP-seq); assay for transposase-accessible chromatin sequencing; and RNA sequencing on 26 lung cancer patient-derived xenograft (PDX) tumors. We then generated and analyzed H3K27ac ChIP-seq, MeDIP-seq, and whole genome sequencing cfDNA data from 1 mL aliquots of plasma from patients with EGFRm LUAD with or without tSCLC. Analysis of 126 epigenomic libraries from the lung cancer PDXs revealed widespread epigenomic reprogramming between LUAD and tSCLC, with a large number of differential H3K27ac (n = 24,424), DNA methylation (n = 3,298), and chromatin accessibility (n = 16,352) sites between the two histologies. Tumor-informed analysis of each of these three epigenomic features in cfDNA resulted in accurate noninvasive discrimination between patients with EGFRm LUAD versus tSCLC [area under the receiver operating characteristic curve (AUROC) = 0.82-0.87]. A multianalyte cfDNA-based classifier integrating these three epigenomic features discriminated between EGFRm LUAD versus tSCLC with an AUROC of 0.94. These data demonstrate the feasibility of detecting small cell transformation in patients with EGFRm LUAD through epigenomic cfDNA profiling of 1 mL of patient plasma.

Characterization of the DNA accessibility of chloroplast genomes in grasses

Although the chloroplast genome (cpDNA) of higher plants is known to exist as a large protein-DNA complex called ‘plastid nucleoid’, researches on its DNA state and regulatory elements are limited. In this study, we performed the assay for transposase-accessible chromatin sequencing (ATAC-seq) on five common tissues across five grasses, and found that the accessibility of different regions in cpDNA varied widely, with the transcribed regions being highly accessible and accessibility patterns around gene start and end sites varying depending on the level of gene expression. Further analysis identified a total of 3970 putative protein binding footprints on cpDNAs of five grasses. These footprints were enriched in intergenic regions and co-localized with known functional elements. Footprints and their flanking accessibility varied dynamically among tissues. Cross-species analysis showed that footprints in coding regions tended to overlap non-degenerate sites and contain a high proportion of highly conserved sites, indicating that they are subject to evolutionary constraints. Taken together, our results suggest that the accessibility of cpDNA has biological implications and provide new insights into the transcriptional regulation of chloroplasts.

Single-cell mitochondrial sequencing reveals low-frequency mitochondrial mutations in naturally aging mice.

Mitochondria play a crucial role in numerous biological processes; however, limited methods and research have focused on revealing mitochondrial heterogeneity at the single-cell level. In this study, we optimized the DNBelab C4 single-cell ATAC (assay for transposase-accessible chromatin) sequencing workflow for single-cell mitochondrial sequencing (C4_mtscATAC-seq). We validated the effectiveness of our C4_mtscATAC-seq protocol by sequencing the HEK-293T cell line with two biological replicates, successfully capturing both mitochondrial content (~68% of total sequencing data) and open chromatin status simultaneously. Subsequently, we applied C4_mtscATAC-seq to investigate two mouse tissues, spleen and bone marrow, obtained from two mice aged 2 months and two mice aged 23 months. Our findings revealed higher mitochondrial DNA (mtDNA) content in young tissues compared to more variable mitochondrial content in aged tissues, consistent with higher activity scores of nuclear genes associated with mitochondrial replication and transcription in young tissues. We detected a total of 22, 15, and 21 mtDNA mutations in the young spleen, aged spleen, and bone marrow, respectively, with most variant allele frequencies (VAF) below 1%. Moreover, we observed a higher number of mtDNA mutations with higher VAF in aged tissues compared to young tissues. Importantly, we identified three mtDNA variations (m.9821A>T, m.15219T>C, and m.15984C>T) with the highest VAF in both aged spleen and aged bone marrow. By comparing cells with and without these mtDNA variations, we analyzed differential open chromatin status to identify potential genes associated with these mtDNA variations, including transcription factors such as KLF15 and NRF1. Our study presents an alternative single-cell mitochondrial sequencing method and provides crude insights into age-related single-cell mitochondrial variations.