ATAC-seq Data Research Articles

SIRT7 knockdown promotes gemcitabine sensitivity of pancreatic cancer cell via upregulation of GLUT3 expression

Gemcitabine serves as a first-line chemotherapeutic treatment for pancreatic cancer (PC), but it is prone to rapid drug resistance. Increasing the sensitivity of PC to gemcitabine has long been a focus of research. Fasting interventions may augment the effects of chemotherapy and present new options. SIRT7 is known to link metabolism with various cellular processes through post-translational modifications. We found upregulation of SIRT7 in PC cells is associated with poor prognosis and gemcitabine resistance. Cross-analysis of RNA-seq and ATAC-seq data suggested that GLUT3 might be a downstream target gene of SIRT7. Subsequent investigations demonstrated that SIRT7 directly interacts with the enhancer region of GLUT3 to desuccinylate H3K122. Our group's another study revealed that GLUT3 can transport gemcitabine in breast cancer cells. Here, we found GLUT3 KD reduces the sensitivity of PC cells to gemcitabine, and SIRT7 KD-associated gemcitabine-sensitizing could be reversed by GLUT3 KD. While fasting mimicking induced upregulation of SIRT7 expression in PC cells, knocking down SIRT7 enhanced sensitivity to gemcitabine through upregulating GLUT3 expression. We further confirmed the effect of SIRT7 deficiency on the sensitivity of gemcitabine under fasting conditions using a mouse xenograft model. In summary, our study demonstrates that SIRT7 can regulate GLUT3 expression by binding to its enhancer and altering H3K122 succinylation levels, thus affecting gemcitabine sensitivity in PC cells. Additionally, combining SIRT7 knockdown with fasting may improve the efficacy of gemcitabine. This unveils a novel mechanism by which SIRT7 influences gemcitabine sensitivity in PC and offer innovative strategies for clinical combination therapy with gemcitabine.

Transcriptomic and epigenomic landscapes of muscle growth during the postnatal period of broilers

BackgroundBroilers stand out as one of the fastest-growing livestock globally, making a substantial contribution to animal meat production. However, the molecular and epigenetic mechanisms underlying the rapid growth and development of broiler chickens are still unclear. This study aims to explore muscle development patterns and regulatory networks during the postnatal rapid growth phase of fast-growing broilers. We measured the growth performance of Cornish (CC) and White Plymouth Rock (RR) over a 42-d period. Pectoral muscle samples from both CC and RR were randomly collected at day 21 after hatching (D21) and D42 for RNA-seq and ATAC-seq library construction.ResultsThe consistent increase in body weight and pectoral muscle weight across both breeds was observed as they matured, with CC outpacing RR in terms of weight at each stage of development. Differential expression analysis identified 398 and 1,129 genes in the two dimensions of breeds and ages, respectively. A total of 75,149 ATAC-seq peaks were annotated in promoter, exon, intron and intergenic regions, with a higher number of peaks in the promoter and intronic regions. The age-biased genes and breed-biased genes of RNA-seq were combined with the ATAC-seq data for subsequent analysis. The results spotlighted the upregulation of ACTC1 and FDPS at D21, which were primarily associated with muscle structure development by gene cluster enrichment. Additionally, a noteworthy upregulation of MUSTN1, FOS and TGFB3 was spotted in broiler chickens at D42, which were involved in cell differentiation and muscle regeneration after injury, suggesting a regulatory role of muscle growth and repair.ConclusionsThis work provided a regulatory network of postnatal broiler chickens and revealed ACTC1 and MUSTN1 as the key responsible for muscle development and regeneration. Our findings highlight that rapid growth in broiler chickens triggers ongoing muscle damage and subsequent regeneration. These findings provide a foundation for future research to investigate the functional aspects of muscle development.

P-045 Deciphering meiotic chromatin organization by SCP3

Abstract Study question What’s the pattern of genome-wide chromatin occupancy of SCP3 during meiotic prophase, and how its interaction with chromatin is regulated? Summary answer We document SCP3 enriched at active chromatin regions with largely inherited from leptotene to pachytene stage, which facilitated by transcription and condensate assembly in spermatocytes. What is known already Meiosis is a critical step during the process of gametogenesis to form haploid cells which finally develop into functional gametes. The synaptonemal complex (SC) is a meiosis-specific protein assembly which bridges homologous chromosomes to establish the architecture of meiotic chromosome. Genetic ablation of genes encoding SC components impedes meiosis and gametogenesis in mice and mutation of these genes is frequently found in infertile patients. SCP3 localization changes with chromosome morphology in autosomes and is involved in maintaining the association of sex chromosomes. Male mice with null mutation of SCP3 failed to form SCs, resulting in developmental arrest during meiotic prophase. Study design, size, duration This study included male mice at post-natal day 10 (P10), 4 weeks or adult (8 weeks or older) were maintained on C57BL/6J background. Spermatocytes were performed modified CUT&Tag using antibodies. ChIP-seq data was analyzed together with WGBS and ATAC-seq data. Participants/materials, setting, methods Spermatocytes were obtained by a discontinuous BSA density gradient. Cell pellet was resuspended in culture medium (MEM with 10% FBS) at a concentration of 105 cells/ml and cultured at 37 °C in a 5% CO2 incubator. Spermatocytes were cultured with or without 1,6-hexanediol or α-amanitin. Spermatocytes with or without drug treatment were collected at indicated time points for further examination and then modified CUT&Tag and co-IP using anti-SCP3 and anti-POL II antibodies. Main results and the role of chance Chromatin occupancy of SCP3 in leptotene and pachytene spermatocytes, SCP3-Core sites retained upon treatment of 1,6-hexanediol which dissolves SCP3 condensate, and SCP1 occupancy at meiotic chromatin in pachytene spermatocytes. When we determined inter-peak distances and width of peaks, we observe that: (1) SCP3 occupancy in leptotene stage had greater inter-peak distances and narrower peak width than SCP3 occupancy in pachytene stage, demonstrating gradual accumulated SCP3 at chromatin during meiotic progression; (2) SCP3-Core sites have broader peak width than SCP3-NonCore sites, in agreement with stability of SCP3 occupancy at SCP3-Core sites; (3) SCP1 occupancy had greater inter-peak distances and narrower peak width than SCP3 occupancy in pachytene stage, supporting that SCP1 peaks are a subset of SCP3 peaks and SCP1’s major function in connecting LE and CE respectively. Generally, at pachytene-stage meiosis, SCP3 mainly occupies chromatin regions including transcriptionally active sites, Deu and Alu elements, while SCP1 mainly co-localizes with SCP3 at a subset of regions with open chromatin and high cohesin enrichment to connect meiotic chromosome axes with CE. In this way, a local active chromatin environment is created at meiotic chromosome axes, facilitating homologous chromosome recombination. Limitations, reasons for caution Our test was not successful to examine SCP3 occupancy with as few as less than 1000 spermatocytes using our protocol. It is improvement of efficiency of tagmentation or development of more efficient antibodies against SCP3 will be helpful for clinical sample detection. Wider implications of the findings SCP3 occupancy was largely inherited from leptotene to pachytene stage, facilitated by transcription and condensate assembly, and enriched at specific SINE repeats. Further exploration using knockout mouse models of meiotic axial proteins and epigenetic modifiers will provide valuable information on meiotic chromatin architecture and crosstalk among meiotic chromatin events. Trial registration number N/A

Fibroblasts as an in vitro model of circadian genetic and genomic studies

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 h period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

ScATAcat: cell-type annotation for scATAC-seq data.

Cells whose accessibility landscape has been profiled with scATAC-seq cannot readily be annotated to a particular cell type. In fact, annotating cell-types in scATAC-seq data is a challenging task since, unlike in scRNA-seq data, we lack knowledge of 'marker regions' which could be used for cell-type annotation. Current annotation methods typically translate accessibility to expression space and rely on gene expression patterns. We propose a novel approach, scATAcat, that leverages characterized bulk ATAC-seq data as prototypes to annotate scATAC-seq data. To mitigate the inherent sparsity of single-cell data, we aggregate cells that belong to the same cluster and create pseudobulk. To demonstrate the feasibility of our approach we collected a number of datasets with respective annotations to quantify the results and evaluate performance for scATAcat. scATAcat is available as a python package at https://github.com/aybugealtay/scATAcat.

Immunological Roles ofCCL18 inPan‑Cancer andIts Potential Value inEndometrial Cancer.

Endometrial cancer (EC) is one of the most prevalent malignancies in the female reproductive system. However, the potential functions and mechanisms of immune-related genes in the onset and progression of EC remain unclear. The immune-related gene CCL18 has been implicated in apoptosis, proliferation, invasion, metastasis, and drug resistance in various types of tumors. Nevertheless, its role in pan-cancer has been poorly investigated, and its expression value and prognostic significance in endometrial cancer (EC) have not been explored. Therefore, the objective of this study was to identify potential immune-related prognostic biomarkers for EC by utilizing the cancer genome atlas (TCGA), immunology database and analysis portal (ImmPort) database, and Gene Expression Omnibus (GEO). Immunohistochemistry staining results from EC tissue chips demonstrated elevated expression levels of inflammatory chemokine protein 18 (CCL18) in EC compared to normal endometrium. This study offers a potential therapeutic strategy for EC treatment by identifying regulatory targets through microRNA sequencing data. Additionally, drug prediction was based on CCL18 targets. Furthermore, an analysis of CCL18 expression in pan-cancer was conducted, and the results revealed its high expression in various types of cancer, including EC and bladder cancer. Through analysis of the ATAC-seq data, we found that SIX1, SOX3, and TWIST2 may regulate CCL18 transcription by binding to the gene promoter of CCL18 in EC. This study indicated that CCL18 could be a potential biomarker in pan-cancer and EC.

Her2 amplification, Rel-A, and Bach1 can influence APOBEC3A expression in breast cancer cells.

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells.

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.

Transcriptomic and epigenomic insights into pectoral muscle fiber formation at the late embryonic development in pure chicken lines

Long-term intensive genetic selection has led to significant differences between broiler and layer chickens, which are evident during the embryonic period. Despite this, there is a paucity of research on the genetic regulation of the initial formation of muscle fiber morphology in chick embryos. Embryonic d 17 (E17) is the key time point for myoblast fusion completion and muscle fiber morphology formation in chickens. This study aimed to explore the genetic regulatory mechanisms underlying the early muscle fiber morphology establishment in broiler chickens of Cornish (CC) and White Plymouth Rock (RR) and layer chickens of White Leghorn (WW) at E17 using the transcriptomic and chromatin accessibility sequencing of pectoral major muscles. The results showed that broiler chickens exhibited significant higher embryo weight and pectoral major muscle weight at E17 compared to layer chickens (P = 0.000). A total of 1,278, 1,248, and 892 differentially expressed genes (DEGs) of RNA-seq data were identified between CC vs. WW, RR vs. WW, and CC vs. RR, separately. All DEGs were combined for cluster analysis and they were divided into 6 clusters, including cluster 1 with higher expression in broilers and cluster 6 with higher expression in layers. DEGs in cluster 1 were enriched in terms related to macrophage activation (P = 0.002) and defense response to bacteria (P = 0.002), while DEGs in cluster 6 showed enrichment in protein-DNA complex (P = 0.003) and monooxygenase activity (P = 0.000). ATAC-seq data analysis identified a total of 38,603 peaks, with 13,051 peaks for CC, 18,780 peaks for RR, and 6,772 peaks for WW. Integrative analysis of transcriptomic and chromatin accessibility data revealed GOLM1, ISLR2, and TOPAZ1 were commonly upregulated genes in CC and RR. Furthermore, screening of all upregulated DEGs in cluster 1 from CC and RR identified GOLM1, ISLR2, and HNMT genes associated with neuroimmune functions and MYOM3 linked to muscle morphology development, showing significantly elevated expression in broiler chickens compared to layer chickens. These findings suggest active neural system connectivity during the initial formation of muscle fiber morphology in embryonic period, highlighting the early interaction between muscle fiber formation morphology and the nervous system. This study provides novel insights into late chick embryo development and lays a deeper foundation for further research.

Abstract A013: Functional and clinical consequences of hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer

Abstract Recent whole-genome sequencing studies have unveiled intriguing hotspot mutations within non-coding genomic regions, yet their implications in cancer remain unknown. Here, we characterized the role of twin hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer: chr8: 29,952,919 G>A and chr8: 29,952,919 C>T. Either mutation is found in ~20% of bladder tumors, and the mutations do not typically co-occur. Through a comprehensive analysis integrating ChiP-seq, ATAC-seq, and Hi-C data from bladder cancer cell lines, we observed that these mutations co-localize with enhancer elements that interact with multiple nearby genes. In vitro dual luciferase reporter assays showed that both mutations are functional and likely affect gene expression in a genomic context. CRISPR-Cas9-mediated knocking out of the putative enhancer region identified LEPROTL1, DCTN6, and SARAF as target genes that are likely co-regulated, which was further supported by strong positive correlations of their mRNA expression in TCGA data. Isogenic cell clones generated by CRISPR-mediated base-editing showed that the mutations result in simultaneous upregulation of the three target genes, which were also expressed at higher levels in TCGA tumors harboring the non-coding mutations. Little is known about the function of LEPROTL1, DCTN6, and SARAF. Utilizing siRNA knockdown models complemented by live-cell imaging, viability assays, and cell cycle analysis, we found that they promote cell proliferation and cell cycle progression. Furthermore, our base-edited isogenic cell lines showed more rapid proliferation and cell cycle progression, as well as resistance to cisplatin. The mutations were also associated with shorter disease-specific survival in a TCGA cohort, and DCTN6 showed a negative association with overall survival. These results underscore the clinical significance of the non-coding mutations. Currently, we are assessing the effect of these mutations on tumor growth and cisplatin response in vivo using xenograft mouse models and evaluating the mechanisms by which LEPROTL1, DCTN6, and SARAF modulate cell cycle progression. Taken together, our data suggest that non-coding enhancer mutations near LEPROTL1 contribute to bladder carcinogenesis by increasing the expression of three co-regulated novel oncogenes, all of which were previously uncharacterized in bladder cancer. Since we found the mutations to be clinically significant, our study also raises the possibility of using the non-coding mutations and their target genes as clinical biomarkers and potential therapeutic drug targets. Citation Format: Kelly Butler, Bilal Lone, Alexandra Dobbins, Arup Chakraborty, Salah Boudajadi, Andrea Apolo, Rouf Banday. Functional and clinical consequences of hotspot non-coding mutations upstream of the LEPROTL1 gene in bladder cancer [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr A013.

Abstract PR009: Exploring the functional consequences of APOBEC3-induced non-coding hotspot mutations in bladder cancer using massively parallel reporter assays and CRISPR-mediated base editing

Abstract Cancer research has traditionally focused on protein-coding genes, but approximately 98% of the human genome is non-coding, comprising gene regulatory elements such as enhancers, silencers, and promoters. Mutations within these regions can foster tumorigenesis by modulating the expression of oncogenes. Recent whole-genome sequencing studies have revealed several recurrently mutated non-coding genomic regions, including in bladder cancer, but their role in tumorigenesis remains unknown. In this study, we delved into whole-genome sequencing data from the bladder cancer TCGA cohort and an independent set of bladder tumors, uncovering hundreds of recurrent mutations. Mutational signature analysis identified over 65% of these mutations as APOBEC3-induced. We substantiated the involvement of APOBEC3 enzymes by using conventional deamination assays and novel approach dubbed as APOBEC3-mediated mutation capture assay. Predominantly, these mutations occurred in the hairpin loop, a favored motif for APOBEC3 enzymes. To systematically assess which of the mutations are functional and prioritize them for biological characterization, we developed Massively Parallel Reporter Assays (MPRA). Two types of MPRA libraries were constructed, one for evaluating promoter activity and the other for assessing enhancer activity of the mutations. Each library construction entailed 36,000 oligonucleotides, comprehensively representing all mutation hotspots and their corresponding wild-type alleles. The application of MPRA resulted in the identification of a substantial number of functionally relevant mutations. Zooming in on the top 20 genomic regions, we employed CRISPR-based approaches to pinpoint target genes and elucidate their functional consequences. Most of these regions co-localized with enhancers, as determined by analysis of ATAC-seq, ChIP-seq, and Hi-C data from bladder cancer cell lines. Perturbation of hotspots regions resulted in dysregulation of multiple genes. The co-regulation patterns of genes were further supported by strong correlation of their mRNA expression in RNA-seq datasets (TCGA, SRA, etc.). Notable examples included twin hotspot mutations upstream of the LEPROTL1 gene, elevating the expression of LEPROTL1, DCTN6 and SARAF, resulting in aggressive tumor growth in in vitro functional experiments. These three genes emerged as novel oncogenes in cancer. Another functional hotspot influenced the expression of two known long non-coding RNAs, NEAT1 and MALAT1, previously characterized as cancer genes. Additionally, some recurrently mutated regions or their target genes were associated with worse overall survival in the TCGA cohort. In conclusion, our findings propose that non-coding mutations likely drive tumorigenesis, contributing to tumor heterogeneity. Employing a combination of functional, multi-omics, and in vivo approaches, we aim to comprehensively characterize the role of hotspot non-coding mutations in bladder cancer and assess their clinical significance. Citation Format: Rouf Banday, Bilal Lone, Kelly Butler, Dhanusha Yesudhas, Alexandra Dobbins, Arup Chakraborty. Exploring the functional consequences of APOBEC3-induced non-coding hotspot mutations in bladder cancer using massively parallel reporter assays and CRISPR-mediated base editing [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr PR009.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales.

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Cactus: A user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner.

Decoding cell replicational age from single-cell ATAC-seq data

Decoding cell replicational age from single-cell ATAC-seq data

Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Disease related changes in ATAC-seq of iPSC-derived motor neuron lines from ALS patients and controls

Amyotrophic Lateral Sclerosis (ALS), like many other neurodegenerative diseases, is highly heritable, but with only a small fraction of cases explained by monogenic disease alleles. To better understand sporadic ALS, we report epigenomic profiles, as measured by ATAC-seq, of motor neuron cultures derived from a diverse group of 380 ALS patients and 80 healthy controls. We find that chromatin accessibility is heavily influenced by sex, the iPSC cell type of origin, ancestry, and the inherent variance arising from sequencing. Once these covariates are corrected for, we are able to identify ALS-specific signals in the data. Additionally, we find that the ATAC-seq data is able to predict ALS disease progression rates with similar accuracy to methods based on biomarkers and clinical status. These results suggest that iPSC-derived motor neurons recapitulate important disease-relevant epigenomic changes.

Abstract 3114: Time Course Multi-omics Approach Investigating Gene Expression Regulation During Vascular Smooth Muscle Cell Phenotypic Transitioning

Coronary artery disease (CAD) is the leading cause of death in the United States. There is a need for novel therapeutic targets that prevent life-threatening consequences directly associated with plaque accumulation in the vascular wall—the underlying cause of CAD. During the progression of atherosclerosis, vascular smooth muscle cells (VSMCs) transdifferentiate from a contractile to a synthetic phenotype that can either promote fibrous cap stability (fibroblast-like) or contribute to inflammation and plaque accumulation (macrophage or osteoblast-like).The molecular mechanisms driving the functional phenotypic switching of VSMCs are not fully understood and represent a unique opportunity for discovering new treatments directly targeting plaque formation in the vascular wall. We measured gene expression regulation during VSMC dedifferentiation using time-course chromatin remodeling (ATACseq) and nascent transcription data (PROseq). We stimulated VSMCs from three healthy heart transplant donors with PDGF-BB for 0, 4, 8, 12, and 16 hours and performed sequencing in duplicate for each donor at each time point for a total of 30 PROseq and 30 ATACseq libraries. We identified 231,315 regions of accessible chromatin (peaks) in the ATACseq data and 12,954 genes in the PROseq data across the five time points. Likelihood ratio testing detected 13,285 peaks and 7,438 genes with time-dynamic changes in read counts representing six patterns of variable accessibility and nascent transcription following PDGF-BB stimulation: gradual increase/decrease (I/D), immediate I/D, and transient I/D. Transcription factor binding site (TFBS) motif enrichment of dynamic peaks revealed 31 TFs unique to one time dynamic. We identified 157 genes directly downstream of peaks enriched for the 31 TFBSs sharing similar time-specific fold changes. THAP1 binding site motifs had the highest number of cis-regulatory connections (596) and downstream genes (80). Our results capture the first ever time-course ATACseq and PROseq analysis conducted in human VSMCs and provide the ability to predict regulatory edges between TFs and genes in order to map the regulatory cascading events occurring during VSMC switching from a contractile to a synthetic state.

The upregulation and transcriptional regulatory mechanisms of Extra spindle pole bodies like 1 in bladder cancer: An immunohistochemistry and high-throughput screening Evaluation

BackgroundThis study aimed to explore the expression level and transcriptional regulation mechanism of Extra Spindle Pole Bodies Like 1 (ESPL1) in bladder cancer (BC). MethodsA multicentre database of samples (n = 1391) was assayed for ESPL1 mRNA expression in BC and validated at the protein level by immunohistochemical (IHC) staining of in-house samples (n = 202). Single-cell sequencing (scRNA-seq) analysis and enrichment analysis explored ESPL1 distribution and their accompanying molecular mechanisms. ATAC-seq, ChIP-seq and Hi-C data from multiple platforms were used to investigate ESPL1 upstream transcription factors (TFs) and potential epigenetic regulatory mechanisms. Immune-related analysis, drug sensitivity and molecular docking of ESPL1 were also calculated. Furthermore, upstream microRNAs and the binding sites of ESPL1 were predicted. The expression level and early screening efficacy of miR-299-5p in blood (n = 6625) and tissues (n = 537) were examined. ResultsESPL1 was significantly overexpressed at the mRNA level (p < 0.05, SMD = 0.75; 95 % CI = 0.09, 1.40), and IHC staining of in-house samples verified this finding (p < 0.0001). ESPL1 was predominantly distributed in BC epithelial cells. Coexpressed genes of ESPL1 were enriched in cell cycle–related signalling pathways, and ESPL1 might be involved in the communication between epithelial and residual cells in the Hippo, ErbB, PI3K-Akt and Ras signalling pathways. Three TFs (H2AZ, IRF5 and HIF1A) were detected upstream of ESPL1 and presence of promoter-super enhancer and promoter-typical enhancer loops. ESPL1 expression was correlated with various immune cell infiltration levels. ESPL1 expression might promote BC growth and affect the sensitivity and therapeutic efficacy of paclitaxel and gemcitabine in BC patients. As an upstream regulator of ESPL1, miR-299-5p expression was downregulated in both the blood and tissues, possessing great potential for early screening. ConclusionsESPL1 expression was upregulated in BC and was mainly distributed in epithelial cells. Elevated ESPL1 expression was associated with TFs at the upstream transcription start site (TSS) and distant chromatin loops of regulatory elements. ESPL1 might be an immune-related predictive and diagnostic marker for BC, and the overexpression of ESPL1 played a cancer-promoting role and affected BC patients’ sensitivity to drug therapy. miR-299-5p was downregulated in BC blood and tissues and was also expected to be a novel marker for early screening.

Overexpressed KCNK1 regulates potassium channels affecting molecular mechanisms and biological pathways in bladder cancer

BackgroundThis study aimed to explore the expression, molecular mechanism and its biological function of potassium two pore domain channel subfamily K member 1 (KCNK1) in bladder cancer (BC).MethodsWe integrated large numbers of external samples (n = 1486) to assess KCNK1 mRNA expression levels and collected in-house samples (n = 245) for immunohistochemistry (IHC) experiments to validate at the KCNK1 protein level. Single-cell RNA sequencing (scRNA-seq) analysis was performed to further assess KCNK1 expression and cellular communication. The transcriptional regulatory mechanisms of KCNK1 expression were explored by ChIP-seq, ATAC-seq and ChIA-PET data. Highly expressed co-expressed genes (HECEGs) of KCNK1 were used to explore potential signalling pathways. Furthermore, the immunoassay, clinical significance and molecular docking of KCNK1 were calculated.ResultsKCNK1 mRNA was significantly overexpressed in BC (SMD = 0.58, 95% CI [0.05; 1.11]), validated at the protein level (p < 0.0001). Upregulated KCNK1 mRNA exhibited highly distinguishing ability between BC and control samples (AUC = 0.82 [0.78–0.85]). Further, scRNA-seq analysis revealed that KCNK1 expression was predominantly clustered in BC epithelial cells and tended to increase with cellular differentiation. BC epithelial cells were involved in cellular communication mainly through the MK signalling pathway. Secondly, the KCNK1 transcription start site (TSS) showed promoter-enhancer interactions in three-dimensional space, while being transcriptionally regulated by GRHL2 and FOXA1. Most of the KCNK1 HECEGs were enriched in cell cycle–related signalling pathways. KCNK1 was mainly involved in cellular metabolism–related pathways and regulated cell membrane potassium channel activity. KCNK1 expression was associated with the level of infiltration of various immune cells. Immunotherapy and chemotherapy (docetaxel, paclitaxel and vinblastine) were more effective in BC patients in the high KCNK1 expression group. KCNK1 expression correlated with age, pathology grade and pathologic_M in BC patients.ConclusionsKCNK1 was significantly overexpressed in BC. A complex and sophisticated three-dimensional spatial transcriptional regulatory network existed in the KCNK1 TSS and promoted the upregulated of KCNK1 expression. The high expression of KCNK1 might be involved in the cell cycle, cellular metabolism, and tumour microenvironment through the regulation of potassium channels, and ultimately contributed to the deterioration of BC.

Alternative splicing and its regulation in the malaria vector Anopheles gambiae

IntroductionAlternative splicing (AS) is a highly conserved mechanism that allows for the expansion of the coding capacity of the genome, through modifications of the way that multiple isoforms are expressed or used to generate different phenotypes. Despite its importance in physiology and disease, genome-wide studies of AS are lacking in most insects, including mosquitoes. Even for model organisms, chromatin associated processes involved in the regulation AS are poorly known.MethodsIn this study, we investigated AS in the mosquito Anopheles gambiae in the context of tissue-specific gene expression and mosquito responses to a Plasmodium falciparum infection, as well as the relationship between patterns of differential isoform expression and usage with chromatin accessibility changes. For this, we combined RNA-seq and ATAC-seq data from A. gambiae midguts and salivary glands, infected and non-infected.ResultsWe report differences between tissues in the expression of 392 isoforms and in the use of 247 isoforms. Secondly, we find a clear and significant association between chromatin accessibility states and tissue-specific patterns of AS. The analysis of differential accessible regions located at splicing sites led to the identification of several motifs resembling the binding sites of Drosophila transcription factors. Finally, the genome-wide analysis of tissue-dependent enhancer activity revealed that approximately 20% of A. gambiae transcriptional enhancers annotate to a differentially expressed or used isoform, and that their activation status is linked to AS differences between tissues.ConclusionThis research elucidates the role of AS in mosquito vector gene expression and identifies regulatory regions potentially involved in AS regulation, which could be important in the development of novel strategies for vector control.