Assay For Transposase-accessible Chromatin Research Articles

CloudATAC: a cloud-based framework for ATAC-Seq data analysis.

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) generates genome-wide chromatin accessibility profiles, providing valuable insights into epigenetic gene regulation at both pooled-cell and single-cell population levels. Comprehensive analysis of ATAC-seq data involves the use of various interdependent programs. Learning the correct sequence of steps needed to process the data can represent a major hurdle. Selecting appropriate parameters at each stage, including pre-analysis, core analysis, and advanced downstream analysis, is important to ensure accurate analysis and interpretation of ATAC-seq data. Additionally, obtaining and working within a limited computational environment presents a significant challenge to non-bioinformatic researchers. Therefore, we present Cloud ATAC, an open-source, cloud-based interactive framework with a scalable, flexible, and streamlined analysis framework based on the best practices approach for pooled-cell and single-cell ATAC-seq data. These frameworks use on-demand computational power and memory, scalability, and a secure and compliant environment provided by the Google Cloud. Additionally, we leverage Jupyter Notebook's interactive computing platform that combines live code, tutorials, narrative text, flashcards, quizzes, and custom visualizations to enhance learning and analysis. Further, leveraging GPU instances has significantly improved the run-time of the single-cell framework. The source codes and data are publicly available through NIH Cloud lab https://github.com/NIGMS/ATAC-Seq-and-Single-Cell-ATAC-Seq-Analysis. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

TGF-β1 mediates hypoxia-preconditioned olfactory mucosa mesenchymal stem cells improved neural functional recovery in Parkinson’s disease models and patients

BackgroundParkinson’s disease (PD) is a neurodegenerative disorder characterized by the degeneration of dopaminergic neurons in the substantia nigra (SN). Activation of the neuroinflammatory response has a pivotal role in PD. Mesenchymal stem cells (MSCs) have emerged as a promising therapeutic approach for various nerve injuries, but there are limited reports on their use in PD and the underlying mechanisms remain unclear.MethodsWe investigated the effects of clinical-grade hypoxia-preconditioned olfactory mucosa (hOM)-MSCs on neural functional recovery in both PD models and patients, as well as the preventive effects on mouse models of PD. To assess improvement in neuroinflammatory response and neural functional recovery induced by hOM-MSCs exposure, we employed single-cell RNA sequencing (scRNA-seq), assay for transposase accessible chromatin with high-throughput sequencing (ATAC-seq) combined with full-length transcriptome isoform-sequencing (ISO-seq), and functional assay. Furthermore, we present the findings from an initial cohort of patients enrolled in a phase I first-in-human clinical trial evaluating the safety and efficacy of intraspinal transplantation of hOM-MSC transplantation into severe PD patients.ResultsA functional assay identified that transforming growth factor-β1 (TGF-β1), secreted from hOM-MSCs, played a critical role in modulating mitochondrial function recovery in dopaminergic neurons. This effect was achieved through improving microglia immune regulation and autophagy homeostasis in the SN, which are closely associated with neuroinflammatory responses. Mechanistically, exposure to hOM-MSCs led to an improvement in neuroinflammation and neural function recovery partially mediated by TGF-β1 via activation of the anaplastic lymphoma kinase/phosphatidylinositol-3-kinase/protein kinase B (ALK/PI3K/Akt) signaling pathway in microglia located in the SN of PD patients. Furthermore, intraspinal transplantation of hOM-MSCs improved the recovery of neurologic function and regulated the neuroinflammatory response without any adverse reactions observed in patients with PD.ConclusionsThese findings provide compelling evidence for the involvement of TGF-β1 in mediating the beneficial effects of hOM-MSCs on neural functional recovery in PD. Treatment and prevention of hOM-MSCs could be a promising and effective neuroprotective strategy for PD. Additionally, TGF-β1 may be used alone or combined with hOM-MSCs therapy for treating PD.

Characterization and integrated analysis of extrachromosomal DNA amplification in hematological malignancies

The study of extrachromosomal DNA (ecDNA), an element existing beyond classical chromosomes, contributes to creating a more comprehensive map of the cancer genome. In hematological malignancies, research on ecDNA has lacked comprehensive investigation into its frequency, structure, function, and mechanisms of formation. We re-analyzed WGS data from 208 hematological cancer samples across 11 types, focusing on ecDNA characteristics. Amplification of ecDNA was observed in 7 of these cancer types, with no instances found in normal blood cells. Patients with leukemia carrying ecDNA showed a low induction therapy remission rate (<30 %), a high relapse rate (75 %) among those who achieved complete remission, and a significantly lower survival rate compared to the general leukemia population, even those with complex chromosomal karyotypes. Among the 55 identified ecDNA amplicons, 268 genes were detected, of which 38 are known cancer-related genes exhibiting significantly increased copy numbers. By integrating RNA-Seq data, we discovered that the increased copy number, resulting in a higher amount of available DNA templates, indeed leads to the elevated expression of genes encoded on ecDNA. Additionally, through the integration of H3K4me3/H3K27ac chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin with sequencing, and high-throughput chromosome conformation capture data, we identified that ecDNA amplifications can also facilitate efficient, copy number-independent amplification of oncogenes. This process is linked to active histone modifications, improved chromatin accessibility, and enhancer hijacking, all of which are effects of ecDNA amplification. Mechanistically, chromothripsis and dysfunction of the DNA repair pathway can, to some extent, explain the origin of ecDNA.

Epigenomic, cistromic, and transcriptomic profiling of primary kidney tubular cells

Spatiotemporal regulation of gene expression is essential for maintaining cellular homeostasis throughout kidney development and disease progression. Transcription factors (TFs) and epigenetic modifications play pivotal roles in controlling gene expression. Profiling chromatin modifications across the genome, along with the distribution and target regulation by TFs in specific kidney cell types, is crucial for understanding the dynamic changes in gene expression. Here, we presented a comprehensive workflow for epigenomic, cistromic, and transcriptomic analyses of primary kidney tubular cells. Specifically, our methodologies included the isolation of primary kidney tubular epithelial cells, RNA extraction, assay for transposase-accessible chromatin using sequencing, ultra-low-input micrococcal nuclease-based native chromatin immunoprecipitation, cleavage under targets and release using nuclease, and subsequent bioinformatic analysis. This protocol provides a methodological framework for investigating the roles of TFs and epigenetic modifications in kidney development and diseases.

Light-responsive transcription factors VvHYH and VvGATA24 mediate wax terpenoid biosynthesis in Vitis vinifera.

The cuticular wax that covers the surfaces of plants is the first barrier against environmental stresses and increasingly accumulates with light exposure. However, the molecular basis of light-responsive wax biosynthesis remains elusive. In grape (Vitis vinifera), light exposure resulted in higher wax terpenoid content and lower decay and abscission rates than controls kept in darkness. Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq data were integrated to draw the chromatin accessibility and cis-elements regulatory map to identify the potential action sites. Terpenoid synthase 12 (VvTPS12) and 3-Hydroxy-3-methylglutaryl-CoA reductase 2 (VvHMGR2) were identified as grape wax biosynthesis targets, while VvHYH and VvGATA24 were identified as terpenoid biosynthesis activators, as more abundant wax crystals and higher wax terpenoid content were observed in transiently overexpressed grape berries and Nicotiana benthamiana leaves. The interaction between VvHYH and the open chromatin of VvTPS12 was confirmed qualitatively using a dual luciferase assay and quantitatively using surface plasma resonance, with an equilibrium dissociation constant of 2.81 nM identified via the latter approach. Molecular docking simulation implied the structural nature of this interaction, indicating that 24 amino acid residues of VvHYH, including Arg106A, could bind to the VvTPS12 G-box cis-element. VvGATA24 directly bound to the open chromatin of VvHMGR2, with an equilibrium dissociation constant of 8.59 nM. 12 amino acid residues of VvGATA24, including Pro218B, interacted with the VvHMGR2 GATA-box cis-element. Our work characterizes the mechanism underlying light-mediated wax terpenoid biosynthesis and provides gene targets for future molecular breeding.

Chromatin accessibility and transcriptional landscape in PK-15 cells during early exposure to Aflatoxin B1

Aflatoxin B1 (AFB1) not only causes significant losses in livestock production but also poses a serious threat to human health. It is the most carcinogenic among known chemicals. Pigs are more susceptible to AFB1 and experience a higher incidence. However, the molecular mechanism of the toxic effect of AFB1 remains unclear. In this study, we used assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA-seq to uncover chromatin accessibility and gene expression dynamics in PK-15 cells during early exposure to AFB1. We observed that the toxic effects of AFB1 involve signaling pathways such as p53, PI3K-AKT, Hippo, MAPK, TLRs, apoptosis, autophagy, and cancer pathways. Basic leucine zipper (bZIP) transcription factors (TFs), including AP-1, Fos, JunB, and Fra2, play a crucial role in regulating the biological processes involved in AFB1 challenge. Several new TFs, such as BORIS, HNF1b, Atf1, and KNRNPH2, represent potential targets for the toxic mechanism of AFB1. In addition, it is crucial to focus on the concentration of intracellular zinc ions. These findings will contribute to a better understanding of the mechanisms underlying AFB1-induced nephrotoxicity and offer new molecular targets.

Establishing an optimized ATAC-seq protocol for the maize.

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Inhibition of HDAC2 sensitises antitumour therapy by promoting NLRP3/GSDMD-mediated pyroptosis in colorectal cancer.

Although numerous studies have indicated that activated pyroptosis can enhance the efficacy of antitumour therapy in several tumours, the precise mechanism of pyroptosis in colorectal cancer (CRC) remains unclear. Pyroptosis in CRC cells treated with antitumour agents was assessed using various techniques, including Western blotting, lactate dehydrogenase release assay and microscopy analysis. To uncover the epigenetic mechanisms that regulate NLRP3, chromatin changes and NLRP3 promoter histone modifications were assessed using Assay for Transposase-Accessible Chromatin using sequencing and RNA sequencing. Chromatin immunoprecipitation‒quantitative polymerase chain reaction was used to investigate the NLRP3 transcriptional regulatory mechanism. Additionally, xenograft and patient-derived xenograft models were constructed to validate the effects of the drug combinations. As the core molecule of the inflammasome, NLRP3 expression was silenced in CRC, thereby limiting gasdermin D (GSDMD)-mediated pyroptosis. Supplementation with NLRP3 can rescue pyroptosis induced by antitumour therapy. Overexpression of HDAC2 in CRC silences NLRP3 via epigenetic regulation. Mechanistically, HDAC2 suppressed chromatin accessibility by eliminating H3K27 acetylation. HDAC2 knockout promotes H3K27ac-mediated recruitment of the BRD4-p-P65 complex to enhance NLRP3 transcription. Inhibiting HDAC2 by Santacruzamate A in combination with classic antitumour agents (5-fluorouracil or regorafenib) in CRC xenograft-bearing animals markedly activated pyroptosis and achieved a significant therapeutic effect. Clinically, HDAC2 is inversely correlated with H3K27ac/p-P65/NLRP3 and is a prognostic factor for CRC patients. Collectively, our data revealed a crucial role for HDAC2 in inhibiting NLRP3/GSDMD-mediated pyroptosis in CRC cells and highlighted HDAC2 as a potential therapeutic target for antitumour therapy. Silencing of NLRP3 limits the GSDMD-dependent pyroptosis in colorectal cancer. HDAC2-mediated histone deacetylation leads to epigenetic silencing of NLRP3. HDAC2 suppresses the NLRP3 transcription by inhibiting the formation of H3K27ac/BRD4/p-P65 complex. Targeting HDAC2 activates pyroptosis and enhances therapeutic effect.

Clock-dependent chromatin accessibility rhythms regulate circadian transcription.

Chromatin organization plays a crucial role in gene regulation by controlling the accessibility of DNA to transcription machinery. While significant progress has been made in understanding the regulatory role of clock proteins in circadian rhythms, how chromatin organization affects circadian rhythms remains poorly understood. Here, we employed ATAC-seq (Assay for Transposase-Accessible Chromatin with Sequencing) on FAC-sorted Drosophila clock neurons to assess genome-wide chromatin accessibility at dawn and dusk over the circadian cycle. We observed significant oscillations in chromatin accessibility at promoter and enhancer regions of hundreds of genes, with enhanced accessibility either at dusk or dawn, which correlated with their peak transcriptional activity. Notably, genes with enhanced accessibility at dusk were enriched with E-box motifs, while those more accessible at dawn were enriched with VRI/PDP1-box motifs, indicating that they are regulated by the core circadian feedback loops, PER/CLK and VRI/PDP1, respectively. Further, we observed a complete loss of chromatin accessibility rhythms in per01 null mutants, with chromatin consistently accessible at both dawn and dusk, underscoring the critical role of Period protein in driving chromatin compaction during the repression phase at dawn. Together, this study demonstrates the significant role of chromatin organization in circadian regulation, revealing how the interplay between clock proteins and chromatin structure orchestrates the precise timing of biological processes throughout the day. This work further implies that variations in chromatin accessibility might play a central role in the generation of diverse circadian gene expression patterns in clock neurons.

Dependency-aware deep generative models for multitasking analysis of spatial omics data.

Spatially resolved transcriptomics (SRT) technologies have significantly advanced biomedical research, but their data analysis remains challenging due to the discrete nature of the data and the high levels of noise, compounded by complex spatial dependencies. Here, we propose spaVAE, a dependency-aware, deep generative spatial variational autoencoder model that probabilistically characterizes count data while capturing spatial correlations. spaVAE introduces a hybrid embedding combining a Gaussian process prior with a Gaussian prior to explicitly capture spatial correlations among spots. It then optimizes the parameters of deep neural networks to approximate the distributions underlying the SRT data. With the approximated distributions, spaVAE can contribute to several analytical tasks that are essential for SRT data analysis, including dimensionality reduction, visualization, clustering, batch integration, denoising, differential expression, spatial interpolation, resolution enhancement and identification of spatially variable genes. Moreover, we have extended spaVAE to spaPeakVAE and spaMultiVAE to characterize spatial ATAC-seq (assay for transposase-accessible chromatin using sequencing) data and spatial multi-omics data, respectively.

P12 Transcriptional and epigenetic changes drive increased cellular plasticity in keloid fibroblasts

Abstract Introduction and aims Keloids are fibroproliferative scars that progressively grow beyond the original wound site. They occur as a result of genetic, epigenetic and environmental abnormalities. Fibroblasts are best characterized for their role in depositing and remodelling the extracellular matrix (ECM) and are believed to have a crucial role in keloid formation. However, the detailed molecular signatures that promote this pathogenic phenotype are still unknown. Methods To investigate the epigenetic and genetic changes in keloid fibroblasts, we performed assay for transposase-accessible chromatin (ATAC)-Seq and RNA-Seq on a patient-matched set of keloid and normal fibroblasts. ATAC-Seq is a high-throughput sequencing method that defines chromatin accessible regions in which the binding of sites of transcription factors that regulate gene expression can be mapped. Identified candidates are validated in keloid fibroblasts via quantitative polymerase chain reaction and their functional role is analysed in various in vitro assays. Results ATAC-Seq analysis revealed more accessibility in coding regions of keloid fibroblasts compared with normal fibroblasts. Pathways enriched in these regions include transcriptional regulation, embryonic development, cell differentiation, calcium-binding and collagen remodelling. Intriguingly, several developmental epigenetic modulators were identified that may alter the transcription and chromatin landscape of keloid fibroblasts to a more plastic state and have the capacity to influence cell differentiation choices, potentially promoting chondrogenic and adipogenic cell features. Currently, we are investigating the ability of keloid fibroblasts to be more readily differentiated into other cell types including induced pluripotent stem cells, and the contribution of those candidates on fibrotic features to the cell are also being explored. Conclusions Our analysis reveals that the chromatin of keloid fibroblasts is more accessible, especially in genomic regions associated with embryonic development and keloid fibroblasts show an increased plasticity in vitro. These findings may support the development of fibroblast-targeted therapies for keloid scars and possibly other connective tissue disorders.

O16 Understanding the pleiotropic effects of the Rho-kinase inhibitor Y-27632 on human keratinocytes

Abstract Introduction and aims Successful epidermal transplant generation depends on sufficiently high numbers of self-renewing stem cells that give rise to clonal colonies termed ‘holoclones’. Treatment with Y-27632, a nonspecific Rho-kinase 1/2 inhibitor, conditionally immortalizes epidermal keratinocyte cultures. However, the pleiotropic signalling mechanisms downstream of Y-27632 remain largely unexplored. We have therefore characterized the effects of continuous long-term Y-27632 treatment on chromatin regulation, and investigated whether Y-27632 reprogrammes a holoclone-like state. Methods Neonatal foreskin-derived human epidermal keratinocytes were cultured from passage p5 to p30, with or without Y treatment. Samples were collected for RNA sequencing analysis from progenitor-enriched middle passage cells (p5, MPK), progenitor-depleted late-passage cells (p8, LPK) and Y-treated cells at p8 (LPK-Y). Samples were also collected for assay for transposase-accessible chromatin (ATAC)-Seq analysis to assess open and closed chromatin regions. Results Untreated cultures became depleted of proliferating cells by p11, whereas Y maintained a highly proliferative culture up to p30 and beyond. Y withdrawal at p11 resulted in rapid loss of proliferative capacity within a single passage. ATAC-Seq identified 196 differentially regulated chromatin sites in MPK vs. LPK. Y treatment maintained accessibility in approximately 40% of these MPK-enriched chromatin regions. Approximately 3000 differentially expressed genes (DEGs) were identified in MPK vs. LPK, while Y treatment induced approximately 11 000 DEGs. Activated genes in LPK-Y were enriched for a holoclone signature, with no evidence for activation of a precancerous transcriptome. Y treatment altered the expression of various genes associated with amino acid metabolism, enabling cells to overcome serine/glycine autotrophy. Conclusions Y acts by preserving accessibility of progenitor-enriched chromatin sites and further reprogrammes the keratinocyte transcriptome towards a more holoclone-like state. However, keratinocytes retain a ‘memory’ of their doubling number, which is conditionally suppressed by Y. The observed global transcriptomic changes are indicative of the highly pleiotropic nature of the Y-induced signalling changes, warranting further deconvolution.

The single-cell chromatin landscape in gonadal cell lineage specification

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

The enriched expression of aquaporin-2 in the renal inner medulla is mediated by increased chromatin accessibility

Introduction: The osmotic gradient between renal cortex (CTX) and renal inner medulla (IM) is important for the generation of a concentrated urine. In addition, the unique hypertonic environment in the IM induces a specific gene expression pattern. For example the expression of aquaporin-2 (Aqp2) is enriched in the IM compared to CTX. We postulated that these changes must be associated with epigenetic mechanisms that could regulate the spatial expression pattern in CTX and IM. We performed assay for transposase-accessible chromatin (ATAC) sequencing to identify genomic regions with differences in accessibility in CTX and IM. Methods: We used kidneys from control mice and collecting duct principal cell specific nuclear factor of activated t cells (Nfat5)-KO mice and dissected CTX and IM. These samples were subjected for ATAC-seq analysis. Accessible genomic regions were identified with MACS2. We next compared differentially accessible regions between CTX and IM in control mice and between control and Nfat5-KO mice. The results were correlated with corresponding gene expression profiles that have been analyzed in a former study. Results: There are massive differences in the accessibility of genomic regions between CTX and IM. In CTX, more than 100,000 regions are differentially accessible compared to IM, with 48,698 regions having higher accessibility and 55,453 regions having lower accessibility. Genes with higher accessibility in CTX show also higher or even specific gene expression in CTX and vice versa. Since loss of NFAT5 is associated with a massive change in gene expression in kidneys from Nfat5-KO mice, we analyzed if this is affected by differences in chromatin accessibility. Our data shows that not only spatial localization in the kidney affects chromatin accessibility, but also the function of NFAT5. Loss of NFAT5 results in 10,456 regions with higher accessibility and 3,858 regions with closed chromatin in the CTX and 16,452 regions losing accessibility and 9,494 regions with increased accessibility in the IM. The accessibility of the Aqp2 gene is higher in IM compared with CTX. The open chromatin of the Aqp2 promoter region correlates with higher expression of Aqp2 in the IM. This suggests that hypertonicity leads to an improvement in the accessibility of the promoter region of Aqp2. In Nfat5-KO IM the accessibility of this region was decreased, which could explain the reduced Aqp2 expression in Nfat5-KO mice. Conclusions: To date, studying hypertonicity-dependent regulation of chromatin accessibility and associated gene expression in CTX and IM in the kidney were missing. Here we show that the differences in gene expression between CTX and IM might be due to changes in accessibility of genomic regions. We also show that the renal expression of Aqp2 might be a result of an epigenetic mechanism. In this context, NFAT5 was found to be an important regulator of chromatin accessibility and gene expression in the kidney. The authors have nothing to disclose. This project was funded by the Deutsche Forschungsgemeinschaft(ED181/9-3). This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Abstract 3064: Distinct Regulatory Programs Govern Smooth Muscle Cell Fate In Atherosclerosis

Introduction: During atherosclerosis, some vascular smooth muscle cells (SMCs) undergo phenotypic modulation, de-differentiating and adopting a fibroblast-like (FMC) or a chondrocyte/osteoblast-like phenotype (CMC). Although FMCs and CMCs appear to occur in succession in single-cell RNAseq data, they occupy distinct locations within the lesion and are associated with protective and deleterious plaque features, respectively. Hypothesis: Smooth muscle cell transitions to FMC and CMC occur via discrete transcriptional pathways during atherosclerosis. Methods: We performed single-cell RNA-sequencing (scRNAseq) and assay for transposase-accessible chromatin with sequencing (ATAC-seq) on lineage-traced vascular SMCs isolated from the aortic root of ApoE mice (n=3). Our analysis included pseudotime and RNA velocity analyses of the scRNAseq data, which independently identified a subgroup of Vcam+Cd74+ modulated SMCs as a potential cell state from which FMCs and CMCs appear to originate, each following distinct trajectories. We explored lineage-specific gene expression using the TradeSeq R package and conducted Ingenuity Pathway Analysis on differentially expressed genes to identify relevant upstream transcription factors, which were assayed for lineage-specific motif accessibility. Results: Our analysis showed profound changes in gene expression between the two trajectories as well as differential activation and chromatin accessibility of specific transcription factors at the point where FMCs diverge from CMCs. Within the CMC lineage, we noted heightened activation and chromatin accessibility of transcription factors associated with cellular stress and inflammatory responses, including Tp53, Foxo4, and Irf1. In stark contrast, these factors exhibited diminished activity within the FMC lineage during the same pseudotime interval. Furthermore, coronary artery disease GWAS transcription factors Twist1 and Tcf21 displayed differential motif accessibility over the FMC and CMC trajectories. Conclusions: Our study provides valuable insights into the regulatory mechanisms governing SMC phenotypic transitions in atherosclerosis, and further underscores the importance of GWAS genes in influencing these transitions.

Histone β-hydroxybutyrylation-mediated Chromatin Remodeling Upregulates β-oxidation of Fatty Acids and Attenuates Metabolic Syndrome

Center for Hypertension and Precision Medicine, Department of Physiology and Pharmacology, University of Toledo, Toledo, OH Introduction: Metabolic syndrome (MetS) is a cluster of metabolic dysregulations including high blood pressure (BP), dyslipidemia, central obesity, hyperglycemia, and insulin resistance. Regular exercise is one of the lifestyle modifications known to prevent MetS by the generation of ketone bodies. Precisely how ketone bodies operate to lower MetS is largely unknown. Previously, we demonstrated that the ketone body, β-hydroxybutyrate, is a potent anti-hypertensive metabolite. β-hydroxybutyrate has been recently recognized to epigenetically modify histones by β-hydroxybutyrylation. Using assay for transposase-accessible chromatin with sequencing and chromatin immunoprecipitation assays, we found that histone β-hydroxybutyrylation causes chromatin remodeling and exposes regions for transcription of genes for β-oxidation of fatty acids. Therefore, we hypothesized that β-hydroxybutyrate alleviates MetS via histone β-hydroxybutyrylation-mediated transcriptional upregulation of β-oxidation of fatty acids. To test this hypothesis, we examined the epigenetic effect of increasing systemic β-hydroxybutyrate in a genetic rat model of MetS inbred in our laboratory, the low-capacity runner (LCR/Tol) rats. Methods: Randomly selected two groups of six months old female LCR/Tol rats (n=5/group) were surgically implanted with radiotelemetry transmitters for BP measurement. The control group was given water, whereas the experimental group was given water containing 20% v/v of 1,3-butanediol (1,3-BD), which is a precursor of β-hydroxybutyrate for 10 weeks. Following euthanasia, sera and tissues were harvested. Histones were isolated and subjected to western blotting for quantitation of β-hydroxybutyrylation. Results: Compared to the control group, rats given 1,3-BD had a significantly increased circulating β-hydroxybutyrate (mean 1.98 mM vs 0.71 mM, p&lt;0.05), lower body weight (mean 231.6 gm vs 248.4 gm, p&lt;0.05), lower fasting glucose (mean 79.60 mg/dl vs 100.8 mg/dl, p&lt;0.05) and lower 24-hour average systolic BP (116±2 mmHg vs. 120±2 mmHg, p&lt;0.05). Interestingly, histones isolated from the experimental group with higher circulating β-hydroxybutyrate also displayed a marked increase in β-hydroxybutyrylation (mean 0.82 vs 0.22, p&lt;0.01). Further, the target transcriptional genes of histone β-hydroxybutyrylation, 3-hydroxy-3-methylglutaryl-CoA synthase 2 (Hmgcs2) and cytochrome P450 family 2 subfamily d polypeptide 4 (Cyp2d4) were both significantly upregulated in the experimental group (p&lt;0.01). However, markers of liver injury, alanine transaminase, alkaline phosphatase, aspartate transaminase, and total bile acid levels, were comparable between the two groups, indicating no liver damage. Conclusion: This is the first study to demonstrate that enhanced systemic β-hydroxybutyrate attenuates MetS by histone β-hydroxybutyrylation-mediated chromatin remodeling to upregulate the transcription of Hmgcs2 and Cyp2d4 to promote β-oxidation of fatty acids. Further, our results demonstrate that for individuals with MetS who may be unable to exercise, utilizing 1,3-BD could serve as an alternate source to combat the MetS ameliorating benefits of the ketone body β-hydroxybutyrate. Funding Sources: Authors acknowledge the National Institutes of Health for funding support to Bina Joe (R01HL1430820). This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

MMGAT: a graph attention network framework for ATAC-seq motifs finding.

Motif finding in Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data is essential to reveal the intricacies of transcription factor binding sites (TFBSs) and their pivotal roles in gene regulation. Deep learning technologies including convolutional neural networks (CNNs) and graph neural networks (GNNs), have achieved success in finding ATAC-seq motifs. However, CNN-based methods are limited by the fixed width of the convolutional kernel, which makes it difficult to find multiple transcription factor binding sites with different lengths. GNN-based methods has the limitation of using the edge weight information directly, makes it difficult to aggregate the neighboring nodes' information more efficiently when representing node embedding. To address this challenge, we developed a novel graph attention network framework named MMGAT, which employs an attention mechanism to adjust the attention coefficients among different nodes. And then MMGAT finds multiple ATAC-seq motifs based on the attention coefficients of sequence nodes and k-mer nodes as well as the coexisting probability of k-mers. Our approach achieved better performance on the human ATAC-seq datasets compared to existing tools, as evidenced the highest scores on the precision, recall, F1_score, ACC, AUC, and PRC metrics, as well as finding 389 higher quality motifs. To validate the performance of MMGAT in predicting TFBSs and finding motifs on more datasets, we enlarged the number of the human ATAC-seq datasets to 180 and newly integrated 80 mouse ATAC-seq datasets for multi-species experimental validation. Specifically on the mouse ATAC-seq dataset, MMGAT also achieved the highest scores on six metrics and found 356 higher-quality motifs. To facilitate researchers in utilizing MMGAT, we have also developed a user-friendly web server named MMGAT-S that hosts the MMGAT method and ATAC-seq motif finding results. The advanced methodology MMGAT provides a robust tool for finding ATAC-seq motifs, and the comprehensive server MMGAT-S makes a significant contribution to genomics research. The open-source code of MMGAT can be found at https://github.com/xiaotianr/MMGAT , and MMGAT-S is freely available at https://www.mmgraphws.com/MMGAT-S/ .

Integrated analysis of ATAC-seq and transcriptomic reveals the ScDof3-ScproC molecular module regulating the cold acclimation capacity of potato

Low temperature severely affects the geographical distribution and production of potato, which may incur cold damage in early spring or winter. Cultivated potatoes, mainly derived from Solanum tuberosum, are sensitive to freezing stress, but wild species of potato such as S. commersonii exhibit both constitutive freezing tolerance and/or cold acclimation tolerance. Hence, such wild species could assist in cold hardiness breeding. Yet the key transcription factors and their downstream functional genes that confer freezing tolerance are far from clear, hindering the breeding process. Here, we used ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) alongside RNA-seq to investigate the variation in chromatin accessibility and patterns of gene expression in freezing-tolerant CMM5 (S. commersonii), before and after its cold treatment. Our results suggest that after exposure to cold, transcription factors including Dof3, ABF2, PIF4, and MYB4 were predicted to further control the genes active in the synthetic/metabolic pathways of plant hormones, namely abscisic acid, polyamine, and reductive glutathione (among others). This suggests these transcription factors could regulate freezing tolerance of CMM5 leaves. In particular, ScDof3 was proven to regulate the expression of ScproC (pyrroline-5-carboxylate reductase, P5CR) according to dual-LUC assays. Overexpressing ScDof3 in Nicotiana benthamiana leaves led to an increase in both the proline content and expression level of NbproC (homolog of ScproC). These results demonstrate the ScDof3-ScproC module regulates the proline content and thus promotes freezing tolerance in potato. Our research provides valuable genetic resources to further study the molecular mechanisms underpinning cold tolerance in potato.

Single-cell RNA sequencing and ATAC sequencing identify novel biomarkers for bicuspid aortic valve-associated thoracic aortic aneurysm.

Bicuspid aortic valve (BAV) is the most prevalent congenital cardiovascular defect and known to cause thoracic aortic aneurysms (TAAs). To improve our understanding of BAV pathogenesis, we characterized the cellular composition of BAV tissues and identified molecular changes in each cell population. Tissue samples from two patients with BAV and two heart transplant donors were analyzed using single-cell RNA sequencing, assay for transposase-accessible chromatin using sequencing, and weighted gene coexpression network analysis for differential gene analysis. TAA-related changes were evaluated by comparing the proportion of each cell type and gene expression profiles between TAA and control tissues. Further, by combining our single-cell RNA sequencing data with publicly available data from genome-wide association studies, we determined critical genes for BAV. We found 20 cell subpopulations in TAA tissues, including multiple subtypes of smooth muscle cells, fibroblasts, macrophages, and T lymphocytes. This result suggested that these cells play multiple functional roles in BAV development. Several differentially expressed genes, including CD9, FHL1y, HSP90AA1, GAS6, PALLD, and ACTA2, were identified. We believe that this comprehensive assessment of the cellular composition of TAA tissues and the insights into altered gene expression patterns can facilitate identification of novel diagnostic biomarkers and therapeutic targets for BAV-associated TAA.

Timeless noncoding DNA contains cell-type preferential enhancers important for proper Drosophila circadian regulation

To address the contribution of transcriptional regulation to Drosophila clock gene expression and to behavior, we generated a series of CRISPR-mediated deletions within two regions of the circadian gene timeless (tim), an intronic E-box region and an upstream E-box region that are both recognized by the key transcription factor Clock (Clk) and its heterodimeric partner Cycle. The upstream deletions but not an intronic deletion dramatically impact tim expression in fly heads; the biggest upstream deletion reduces peak RNA levels and tim RNA cycling amplitude to about 15% of normal, and there are similar effects on tim protein (TIM). The cycling amplitude of other clock genes is also strongly reduced, in these cases due to increases in trough levels. These data underscore the important contribution of the upstream E-box enhancer region to tim expression and of TIM to clock gene transcriptional repression in fly heads. Surprisingly, tim expression in clock neurons is only modestly affected by the biggest upstream deletion and is similarly affected by a deletion of the intronic E-box region. This distinction between clock neurons and glia is paralleled by a dramatically enhanced accessibility of the intronic enhancer region within clock neurons. This distinctive feature of tim chromatin was revealed by ATAC-seq (assay for transposase-accessible chromatin with sequencing) assays of purified neurons and glia as well as of fly heads. The enhanced cell type-specific accessibility of the intronic enhancer region explains the resilience of clock neuron tim expression and circadian behavior to deletion of the otherwise more prominent upstream tim E-box region.