Accessible Chromatin Regions Research Articles

Dissociation of Human and Mouse Tumor Tissue Samples for Single-cell RNA Sequencing.

Human tumor samples hold a plethora of information about their microenvironment and immune repertoire. Effective dissociation of human tissue samples into viable cell suspensions is a required input for the single-cell RNA sequencing (scRNAseq) pipeline. Unlike bulk RNA sequencing approaches, scRNAseq enables us to infer the transcriptional heterogeneity in tumor specimens at the single-cell level. Incorporating this approach in recent years has led to many discoveries, such as identifying immune and tumor cellular states and programs associated with clinical responses to immunotherapies and other types of treatments. Moreover, single-cell technologies applied to dissociated tissues can be used to identify accessible chromatin regions T and B cell receptor repertoire, and the expression of proteins, using DNA barcoded antibodies (CITEseq). The viability and quality of the dissociated sample are critical variables when using these technologies, as these can dramatically affect the cross-contamination of single cells with ambient RNA, the quality of the data, and interpretation. Moreover, long dissociation protocols can lead to the elimination of sensitive cell populations and the upregulation of a stress response gene signature. To overcome these limitations, we devised a rapid universal dissociation protocol, which has been validated on multiple types of human and murine tumors. The process begins with mechanical and enzymatic dissociation, followed by filtration, red blood lysis, and live dead enrichment, suitable for samples with a low input of cells (e.g., needle core biopsies). This protocol ensures a clean and viable single-cell suspension paramount to the successful generation of Gel Bead-In Emulsions (GEMs), barcoding, and sequencing.

Single-nucleus multiome analysis of human cerebellum in Alzheimer's disease-related dementia.

Although human cerebellum is known to be neuropathologically impaired in Alzheimer's disease (AD) and AD-related dementias (ADRD), the cell type-specific transcriptional and epigenomic changes that contribute to this pathology are not well understood. Here, we report single-nucleus multiome (snRNA-seq and snATAC-seq) analysis of 103,861 nuclei isolated from cerebellum from 9 human cases of AD/ADRD and 8 controls, and with frontal cortex of 6 AD donors for additional comparison. Using peak-to-gene linkage analysis, we identified 431,834 significant linkages between gene expression and cell subtype-specific chromatin accessibility regions enriched for candidate cis-regulatory elements (cCREs). These cCREs were associated with AD/ADRD-specific transcriptomic changes and disease-related gene regulatory networks, especially for RAR Related Orphan Receptor A (RORA) and E74 Like ETS Transcription Factor 1 (ELF1) in cerebellar Purkinje cells and granule cells, respectively. Trajectory analysis of granule cell populations further identified disease-relevant transcription factors, such as RORA, and their regulatory targets. Finally, we prioritized two likely causal genes, including Seizure Related 6 Homolog Like 2 (SEZ6L2) in Purkinje cells and KAT8 Regulatory NSL Complex Subunit 1 (KANSL1) in granule cells, through integrative analysis of cCREs derived from snATAC-seq, genome-wide AD/ADRD loci, and Hi-C looping data. This first cell subtype-specific regulatory landscape in the human cerebellum identified here offer novel genomic and epigenomic insights into the neuropathology and pathobiology of AD/ADRD and other neurological disorders if broadly applied.

Central role of SUMOylation in the regulation of chromatin interactions and transcriptional outputs of the androgen receptor in prostate cancer cells.

The androgen receptor (AR) is pivotal in prostate cancer (PCa) progression and represents a critical therapeutic target. AR-mediated gene regulation involves intricate interactions with nuclear proteins, with many mediating and undergoing post-translational modifications that present alternative therapeutic avenues. Through chromatin proteomics in PCa cells, we identified SUMO ligases together with nuclear receptor coregulators and pioneer transcription factors within the AR's protein network. Intriguingly, this network displayed a significant association with SUMO2/3. To elucidate the influence of SUMOylation on AR chromatin interactions and subsequent gene regulation, we inhibited SUMOylation using ML-792 (SUMOi). While androgens generally facilitated the co-occupancy of SUMO2/3 and AR on chromatin, SUMOi induced divergent effects dependent on the type of AR-binding site (ARB). SUMOi augmented AR's pioneer-like binding on inaccessible chromatin regions abundant in androgen response elements (AREs) and diminished its interaction with accessible chromatin regions sparse in AREs yet rich in pioneer transcription factor motifs. The SUMOi-impacted ARBs divergently influenced AR-regulated genes; those associated with AR-mediated activation played roles in negative regulation of cell proliferation, while those with AR-mediated repression were involved in pattern formation. In conclusion, our findings underscore the pervasive influence of SUMOylation in shaping AR's role in PCa cells, potentially unveiling new therapeutic strategies.

ScaDA: A novel statistical method for differential analysis of single-cell chromatin accessibility sequencing data.

Single-cell ATAC-seq sequencing data (scATAC-seq) has been widely used to investigate chromatin accessibility on the single-cell level. One important application of scATAC-seq data analysis is differential chromatin accessibility (DA) analysis. However, the data characteristics of scATAC-seq such as excessive zeros and large variability of chromatin accessibility across cells impose a unique challenge for DA analysis. Existing statistical methods focus on detecting the mean difference of the chromatin accessible regions while overlooking the distribution difference. Motivated by real data exploration that distribution difference exists among cell types, we introduce a novel composite statistical test named "scaDA", which is based on zero-inflated negative binomial model (ZINB), for performing differential distribution analysis of chromatin accessibility by jointly testing the abundance, prevalence and dispersion simultaneously. Benefiting from both dispersion shrinkage and iterative refinement of mean and prevalence parameter estimates, scaDA demonstrates its superiority to both ZINB-based likelihood ratio tests and published methods by achieving the highest power and best FDR control in a comprehensive simulation study. In addition to demonstrating the highest power in three real sc-multiome data analyses, scaDA successfully identifies differentially accessible regions in microglia from sc-multiome data for an Alzheimer's disease (AD) study that are most enriched in GO terms related to neurogenesis and the clinical phenotype of AD, and AD-associated GWAS SNPs.

Current genomic deep learning models display decreased performance in cell type-specific accessible regions

BackgroundA number of deep learning models have been developed to predict epigenetic features such as chromatin accessibility from DNA sequence. Model evaluations commonly report performance genome-wide; however, cis regulatory elements (CREs), which play critical roles in gene regulation, make up only a small fraction of the genome. Furthermore, cell type-specific CREs contain a large proportion of complex disease heritability.ResultsWe evaluate genomic deep learning models in chromatin accessibility regions with varying degrees of cell type specificity. We assess two modeling directions in the field: general purpose models trained across thousands of outputs (cell types and epigenetic marks) and models tailored to specific tissues and tasks. We find that the accuracy of genomic deep learning models, including two state-of-the-art general purpose models―Enformer and Sei―varies across the genome and is reduced in cell type-specific accessible regions. Using accessibility models trained on cell types from specific tissues, we find that increasing model capacity to learn cell type-specific regulatory syntax―through single-task learning or high capacity multi-task models―can improve performance in cell type-specific accessible regions. We also observe that improving reference sequence predictions does not consistently improve variant effect predictions, indicating that novel strategies are needed to improve performance on variants.ConclusionsOur results provide a new perspective on the performance of genomic deep learning models, showing that performance varies across the genome and is particularly reduced in cell type-specific accessible regions. We also identify strategies to maximize performance in cell type-specific accessible regions.

Super-enhancer-driven IRF2BP2 enhances ALK activity and promotes neuroblastoma cell proliferation.

Super-enhancers (SEs) typically govern the expression of critical oncogenes and play a fundamental role in the initiation and progression of cancer. Focusing on genes that are abnormally regulated by SE in cancer may be a new strategy for understanding pathogenesis. In the context of this investigation, we have identified a previously unreported SE-driven gene IRF2BP2 in neuroblastoma (NB). The expression and prognostic value of IRF2BP2 were detected in public databases and clinical samples. The effect of IRF2BP2 on NB cell growth and apoptosis was evaluated through in vivo and in vitro functional loss experiments. The molecular mechanism of IRF2BP2 was investigated by the study of chromatin regulatory regions and transcriptome sequencing. The sustained high expression of IRF2BP2 results from the activation of a novel SE established by NB master transcription factors MYCN, MEIS2, and HAND2, and they form a new complex that regulates the gene network associated with the proliferation of NB cell populations. We also observed a significant enrichment of the AP-1 family at the binding sites of IRF2BP2. Remarkably, within NB cells, AP-1 plays a pivotal role in shaping the chromatin accessibility landscape, thereby exposing the binding site for IRF2BP2. This orchestrated action enables AP-1 and IRF2BP2 to collaboratively stimulate the expression of the NB susceptibility gene ALK, thereby upholding the highly proliferative phenotype characteristic of NB. Our findings indicate that SE-driven IRF2BP2 can bind to AP-1 to maintain the survival of tumor cells via regulating chromatin accessibility of the NB susceptibility gene ALK.

Epigenomic identification of vernalization cis-regulatory elements in winter wheat

BackgroundWinter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown.ResultsIn this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues—leaf, axillary bud, and shoot apex—during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering.ConclusionsOur study offers critical insights into wheat vernalization’s epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

INO80 regulates chromatin accessibility to facilitate suppression of sex-linked gene expression during mouse spermatogenesis.

Chromatin accessibility and DNA repair factor localization at the sex chromosomes are facilitated by INO80, which regulates sex-linked gene silencing during meiotic progression in spermatocytes.

Terminal deoxynucleotidyl transferase and CD84 identify human multi-potent lymphoid progenitors

Lymphoid specification in human hematopoietic progenitors is not fully understood. To better associate lymphoid identity with protein-level cell features, we conduct a highly multiplexed single-cell proteomic screen on human bone marrow progenitors. This screen identifies terminal deoxynucleotidyl transferase (TdT), a specialized DNA polymerase intrinsic to VDJ recombination, broadly expressed within CD34+ progenitors prior to B/T cell emergence. While these TdT+ cells coincide with granulocyte-monocyte progenitor (GMP) immunophenotype, their accessible chromatin regions show enrichment for lymphoid-associated transcription factor (TF) motifs. TdT expression on GMPs is inversely related to the SLAM family member CD84. Prospective isolation of CD84lo GMPs demonstrates robust lymphoid potentials ex vivo, while still retaining significant myeloid differentiation capacity, akin to LMPPs. This multi-omic study identifies human bone marrow lymphoid-primed progenitors, further defining the lympho-myeloid axis in human hematopoiesis.

The transcriptional co-repressor Runx1t1 is essential for MYCN-driven neuroblastoma tumorigenesis

MYCN oncogene amplification is frequently observed in aggressive childhood neuroblastoma. Using an unbiased large-scale mutagenesis screen in neuroblastoma-prone transgenic mice, we identify a single germline point mutation in the transcriptional corepressor Runx1t1, which abolishes MYCN-driven tumorigenesis. This loss-of-function mutation disrupts a highly conserved zinc finger domain within Runx1t1. Deletion of one Runx1t1 allele in an independent Runx1t1 knockout mouse model is also sufficient to prevent MYCN-driven neuroblastoma development, and reverse ganglia hyperplasia, a known pre-requisite for tumorigenesis. Silencing RUNX1T1 in human neuroblastoma cells decreases colony formation in vitro, and inhibits tumor growth in vivo. Moreover, RUNX1T1 knockdown inhibits the viability of PAX3-FOXO1 fusion-driven rhabdomyosarcoma and MYC-driven small cell lung cancer cells. Despite the role of Runx1t1 in MYCN-driven tumorigenesis neither gene directly regulates the other. We show RUNX1T1 forms part of a transcriptional LSD1-CoREST3-HDAC repressive complex recruited by HAND2 to enhancer regions to regulate chromatin accessibility and cell-fate pathway genes.

Current genomic deep learning models display decreased performance in cell type specific accessible regions.

A number of deep learning models have been developed to predict epigenetic features such as chromatin accessibility from DNA sequence. Model evaluations commonly report performance genome-wide; however, cis regulatory elements (CREs), which play critical roles in gene regulation, make up only a small fraction of the genome. Furthermore, cell type specific CREs contain a large proportion of complex disease heritability. We evaluate genomic deep learning models in chromatin accessibility regions with varying degrees of cell type specificity. We assess two modeling directions in the field: general purpose models trained across thousands of outputs (cell types and epigenetic marks), and models tailored to specific tissues and tasks. We find that the accuracy of genomic deep learning models, including two state-of-the-art general purpose models - Enformer and Sei - varies across the genome and is reduced in cell type specific accessible regions. Using accessibility models trained on cell types from specific tissues, we find that increasing model capacity to learn cell type specific regulatory syntax - through single-task learning or high capacity multi-task models - can improve performance in cell type specific accessible regions. We also observe that improving reference sequence predictions does not consistently improve variant effect predictions, indicating that novel strategies are needed to improve performance on variants. Our results provide a new perspective on the performance of genomic deep learning models, showing that performance varies across the genome and is particularly reduced in cell type specific accessible regions. We also identify strategies to maximize performance in cell type specific accessible regions.

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice.

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

Cellular zinc status alters chromatin accessibility and binding of p53 to DNA.

Zn2+ is an essential metal required by approximately 850 human transcription factors. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Using ATAC-seq to profile regions of accessible chromatin coupled with transcription factor enrichment analysis, we examined how increases and decreases in the labile zinc pool affect chromatin accessibility and transcription factor enrichment. We found 685 transcription factor motifs were differentially enriched, corresponding to 507 unique transcription factors. The pattern of perturbation and the types of transcription factors were notably different at promoters versus intergenic regions, with zinc-finger transcription factors strongly enriched in intergenic regions in elevated Zn2+ To test whether ATAC-seq and transcription factor enrichment analysis predictions correlate with changes in transcription factor binding, we used ChIP-qPCR to profile six p53 binding sites. We found that for five of the six targets, p53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc alter chromatin accessibility and transcription factor binding to DNA.

A spatially resolved multiomic single-cell atlas of soybean development.

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

Defining the Regulatory Logic of Breast Cancer Using Single-Cell Epigenetic and Transcriptome Profiling.

Annotation of the cis-regulatory elements that drive transcriptional dysregulation in cancer cells is critical to improving our understanding of tumor biology. Herein, we present a compendium of matched chromatin accessibility (scATAC-seq) and transcriptome (scRNA-seq) profiles at single-cell resolution from human breast tumors and healthy mammary tissues processed immediately following surgical resection. We identify the most likely cell-of-origin for luminal breast tumors and basal breast tumors and then introduce a novel methodology that implements linear mixed-effects models to systematically quantify associations between regions of chromatin accessibility (i.e. regulatory elements) and gene expression in malignant cells versus normal mammary epithelial cells. These data unveil regulatory elements with that switch from silencers of gene expression in normal cells to enhancers of gene expression in cancer cells, leading to the upregulation of clinically relevant oncogenes. To translate the utility of this dataset into tractable models, we generated matched scATAC-seq and scRNA-seq profiles for breast cancer cell lines, revealing, for each subtype, a conserved oncogenic gene expression program between in vitro and in vivo cells. Together, this work highlights the importance of non-coding regulatory mechanisms that underlie oncogenic processes and the ability of single-cell multi-omics to define the regulatory logic of BC cells at single-cell resolution.

Single-nucleus chromatin landscape dataset of mouse brain development and aging

The development and aging of the brain constitute a lifelong dynamic process, marked by structural and functional changes that entail highly coordinated cellular differentiation and epigenetic regulatory mechanisms. Chromatin accessibility serves as the foundational basis for genetic activity. However, the holistic and dynamic chromatin landscape that spans various brain regions throughout development and ageing remains predominantly unexplored. In this study, we employed single-nucleus ATAC-seq to generate comprehensive chromatin accessibility maps, incorporating data from 69,178 cells obtained from four distinct brain regions – namely, the olfactory bulb (OB), cerebellum (CB), prefrontal cortex (PFC), and hippocampus (HP) – across key developmental time points at 7 P, 3 M, 12 M, and 18 M. We delineated the distribution of cell types across different age stages and brain regions, providing insight into chromatin accessible regions and key transcription factors specific to different cell types. Our data contribute to understanding the epigenetic basis of the formation of different brain regions, providing a dynamic landscape and comprehensive resource for revealing gene regulatory programs during brain development and aging.

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars.

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates.

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

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.