Cell Type-specific Transcription Factors Research Articles

Context transcription factors establish cooperative environments and mediate enhancer communication.

Many enhancers control gene expression by assembling regulatory factor clusters, also referred to as condensates. This process is vital for facilitating enhancer communication and establishing cellular identity. However, how DNA sequence and transcription factor (TF) binding instruct the formation of high regulatory factor environments remains poorly understood. Here we developed a new approach leveraging enhancer-centric chromatin accessibility quantitative trait loci (caQTLs) to nominate regulatory factor clusters genome-wide. By analyzing TF-binding signatures within the context of caQTLs and comparing episomal versus endogenous enhancer activities, we discovered a class of regulators, 'context-only' TFs, that amplify the activity of cell type-specific caQTL-binding TFs, that is, 'context-initiator' TFs. Similar to super-enhancers, enhancers enriched for context-only TF-binding sites display high coactivator binding and sensitivity to bromodomain-inhibiting molecules. We further show that binding sites for context-only and context-initiator TFs underlie enhancer coordination, providing a mechanistic rationale for how a loose TF syntax confers regulatory specificity.

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

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

Evaluating Methods for the Prediction of Cell Type-Specific Enhancers in the Mammalian Cortex.

Identifying cell type-specific enhancers in the brain is critical to building genetic tools for investigating the mammalian brain. Computational methods for functional enhancer prediction have been proposed and validated in the fruit fly and not yet the mammalian brain. We organized the 'Brain Initiative Cell Census Network (BICCN) Challenge: Predicting Functional Cell Type-Specific Enhancers from Cross-Species Multi-Omics' to assess machine learning and feature-based methods designed to nominate enhancer DNA sequences to target cell types in the mouse cortex. Methods were evaluated based on in vivo validation data from hundreds of cortical cell type-specific enhancers that were previously packaged into individual AAV vectors and retro-orbitally injected into mice. We find that open chromatin was a key predictor of functional enhancers, and sequence models improved prediction of non-functional enhancers that can be deprioritized as opposed to pursued for in vivo testing. Sequence models also identified cell type-specific transcription factor codes that can guide designs of in silico enhancers. This community challenge establishes a benchmark for enhancer prioritization algorithms and reveals computational approaches and molecular information that are crucial for the identification of functional enhancers for mammalian cortical cell types. The results of this challenge bring us closer to understanding the complex gene regulatory landscape of the mammalian brain and help us design more efficient genetic tools and potential gene therapies for human neurological diseases.

CDK9 phosphorylates RUNX1 to promote megakaryocytic fate in megakaryocytic-erythroid progenitors

AbstractThe specification of megakaryocytic (Mk) or erythroid (E) lineages from primary human megakaryocytic-erythroid progenitors (MEPs) is crucial for hematopoietic homeostasis, yet the underlying mechanisms regulating fate specification remain elusive. In this study, we identify RUNX1 as a key modulator of gene expression during MEP fate specification. Overexpression of RUNX1 in primary human MEPs promotes Mk specification, whereas pan-RUNX inhibition favors E specification. Although total RUNX1 levels do not differ between Mk progenitors (MkPs) and E progenitors (ErPs), there are higher levels of serine-phosphorylated RUNX1 in MkPs than ErPs, and mutant RUNX1 with phosphorylated-serine/threonine mimetic mutations (RUNX1-4D) significantly enhances the functional efficacy of RUNX1. To model the effects of RUNX1 variants, we use human erythroleukemia (HEL) cell lines expressing wild-type (WT), phosphomimetic (RUNX1-4D), and nonphosphorylatable (RUNX1-4A) mutants showing that the 3 forms of RUNX1 differentially regulate expression of 2625 genes. Both WT and RUNX1-4D variants increase expression in 40%, and decrease expression in another 40%, with lesser effects of RUNX1-4A. We find a significant overlap between the upregulated genes in WT and RUNX1-4D–expressing HEL cells and those upregulated in primary human MkPs vs MEPs. Although inhibition of known RUNX1 serine/threonine kinases does not affect phosphoserine RUNX1 levels in primary MEPs, specific inhibition of cyclin dependent kinase 9 (CDK9) in MEPs leads to both decreased RUNX1 phosphorylation and increased E commitment. Collectively, our findings show that serine/threonine phosphorylation of RUNX1 promotes Mk fate specification and introduce a novel kinase for RUNX1 linking the fundamental transcriptional machinery with activation of a cell type–specific transcription factor.

Identifying cell type-specific transcription factor-mediated activity immune modules reveal implications for immunotherapy and molecular classification of pan-cancer.

Systematic investigation of tumor-infiltrating immune (TII) cells is important to the development of immunotherapies, and the clinical response prediction in cancers. There exists complex transcriptional regulation within TII cells, and different immune cell types display specific regulation patterns. To dissect transcriptional regulation in TII cells, we first integrated the gene expression profiles from single-cell datasets, and proposed a computational pipeline to identify TII cell type-specific transcription factor (TF) mediated activity immune modules (TF-AIMs). Our analysis revealed key TFs, such as BACH2 and NFKB1 play important roles in B and NK cells, respectively. We also found some of these TF-AIMs may contribute to tumor pathogenesis. Based on TII cell type-specific TF-AIMs, we identified eight CD8+ T cell subtypes. In particular, we found the PD1 + CD8+ T cell subset and its specific TF-AIMs associated with immunotherapy response. Furthermore, the TII cell type-specific TF-AIMs displayed the potential to be used as predictive markers for immunotherapy response of cancer patients. At the pan-cancer level, we also identified and characterized six molecular subtypes across 9680 samples based on the activation status of TII cell type-specific TF-AIMs. Finally, we constructed a user-friendly web interface CellTF-AIMs (http://bio-bigdata.hrbmu.edu.cn/CellTF-AIMs/) for exploring transcriptional regulatory pattern in various TII cell types. Our study provides valuable implications and a rich resource for understanding the mechanisms involved in cancer microenvironment and immunotherapy.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BackgroundInteractions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.ResultsWe comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.ConclusionsOur study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.

GET: a foundation model of transcription across human cell types.

Transcriptional regulation, involving the complex interplay between regulatory sequences and proteins, directs all biological processes. Computational models of transcription lack generalizability to accurately extrapolate in unseen cell types and conditions. Here, we introduce GET, an interpretable foundation model designed to uncover regulatory grammars across 213 human fetal and adult cell types. Relying exclusively on chromatin accessibility data and sequence information, GET achieves experimental-level accuracy in predicting gene expression even in previously unseen cell types. GET showcases remarkable adaptability across new sequencing platforms and assays, enabling regulatory inference across a broad range of cell types and conditions, and uncovering universal and cell type specific transcription factor interaction networks. We evaluated its performance on prediction of regulatory activity, inference of regulatory elements and regulators, and identification of physical interactions between transcription factors. Specifically, we show GET outperforms current models in predicting lentivirus-based massive parallel reporter assay readout with reduced input data. In fetal erythroblasts, we identify distal (>1Mbp) regulatory regions that were missed by previous models. In B cells, we identified a lymphocyte-specific transcription factor-transcription factor interaction that explains the functional significance of a leukemia-risk predisposing germline mutation. In sum, we provide a generalizable and accurate model for transcription together with catalogs of gene regulation and transcription factor interactions, all with cell type specificity.

Mapping Single-Cell Transcriptomes of Endometrium Reveals Potential Biomarkers in Endometrial Cancer.

The heterogeneity and dynamic changes of endometrial cells have a significant impact on health as they determine the normal function of the endometrium during the menstrual cycle. Dysfunction of the endometrium can lead to the occurrence of various gynecological diseases. Therefore, deconvolution of immune microenvironment that drives transcriptional programs throughout the menstrual cycle is key to understand regulatory biology of endometrium. Herein, we comprehensively analyzed single-cell transcriptome of 59,397 cells across ten human endometrium samples and revealed the dynamic cellular heterogeneity throughout the menstrual cycle. We identified two perivascular cell subtypes, four epithelial subtypes and four fibroblast cell types in endometrium. Moreover, we inferred the cell type-specific transcription factor (TF) activities and linked critical TFs to transcriptional output of diverse immune cell types, highlighting the importance of transcriptional regulation in endometrium. Dynamic interactions between various types of cells in endometrium contribute to a range of biological pathways regulating differentiation of secretory. Integration of the molecular biomarkers identified in endometrium and bulk transcriptome of 535 endometrial cancers (EC), we revealed five RNA-based molecular subtypes of EC with highly intratumoral heterogeneity and different clinical manifestations. Mechanism analysis uncovered clinically relevant pathways for pathogenesis of EC. In summary, our results revealed the dynamic immune microenvironment of endometrium and provided novel insights into future development of RNA-based treatments for endometriosis and endometrial carcinoma.

Single-cell insights: pioneering an integrated atlas of chromatin accessibility and transcriptomic landscapes in diabetic cardiomyopathy

BackgroundDiabetic cardiomyopathy (DCM) poses a growing health threat, elevating heart failure risk in diabetic individuals. Understanding DCM is crucial, with fibroblasts and endothelial cells playing pivotal roles in driving myocardial fibrosis and contributing to cardiac dysfunction. Advances in Multimodal single-cell profiling, such as scRNA-seq and scATAC-seq, provide deeper insights into DCM’s unique cell states and molecular landscape for targeted therapeutic interventions.MethodsSingle-cell RNA and ATAC data from 10x Multiome libraries were processed using Cell Ranger ARC v2.0.1. Gene expression and ATAC data underwent Seurat and Signac filtration. Differential gene expression and accessible chromatin regions were identified. Transcription factor activity was estimated with chromVAR, and Cis-coaccessibility networks were calculated using Cicero. Coaccessibility connections were compared to the GeneHancer database. Gene Ontology analysis, biological process scoring, cell-cell communication analysis, and gene-motif correlation was performed to reveal intricate molecular changes. Immunofluorescent staining utilized various antibodies on paraffin-embedded tissues to verify the findings.ResultsThis study integrated scRNA-seq and scATAC-seq data obtained from hearts of WT and DCM mice, elucidating molecular changes at the single-cell level throughout the diabetic cardiomyopathy progression. Robust and accurate clustering analysis of the integrated data revealed altered cell proportions, showcasing decreased endothelial cells and macrophages, coupled with increased fibroblasts and myocardial cells in the DCM group, indicating enhanced fibrosis and endothelial damage. Chromatin accessibility analysis unveiled unique patterns in cell types, with heightened transcriptional activity in myocardial cells. Subpopulation analysis highlighted distinct changes in cardiomyocytes and fibroblasts, emphasizing pathways related to fatty acid metabolism and cardiac contraction. Fibroblast-centered communication analysis identified interactions with endothelial cells, implicating VEGF receptors. Endothelial cell subpopulations exhibited altered gene expressions, emphasizing contraction and growth-related pathways. Candidate regulators, including Tcf21, Arnt, Stat5a, and Stat5b, were identified, suggesting their pivotal roles in DCM development. Immunofluorescence staining validated marker genes of cell subpopulations, confirming PDK4, PPARγ and Tpm1 as markers for metabolic pattern-altered cardiomyocytes, activated fibroblasts and endothelial cells with compromised proliferation.ConclusionOur integrated scRNA-seq and scATAC-seq analysis unveils intricate cell states and molecular alterations in diabetic cardiomyopathy. Identified cell type-specific changes, transcription factors, and marker genes offer valuable insights. The study sheds light on potential therapeutic targets for DCM.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5′ to 3′). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Abstract 1686: Dissecting the 3D enhancer-promoter hubs in prostate tumorigenesis

Abstract The transcription process is controlled by regulatory elements that are often located in non-coding regions. Among regulatory elements, activities of enhancers that are bound by cell-type specific transcription factors (TFs) are tightly linked to cell identity. Enhancers loop to interact with the promoters of target genes at a distance to drive different cell fates. Evidently, transcriptional dysregulation linked to enhancers that are differentially activated in cancer is an underlying factor in tumorigenesis. To identify enhancers linked to prostate tumorigenesis, we analyzed hundreds of in-house and public H3K27ac ChIP-seq data generated from normal prostate and prostate cancer cell lines and tissues. We identified thousands of prostate cancer-specific enhancers (PSEs) that were activated exclusively in cancer conditions. To further identify the target genes of PSEs, we performed 3C-derived assays, such as Hi-C and Micro-C, and mapped genome-wide interactions. However, they were limited by the resolution to detect enhancer-promoter interactions. The recent development of Region Capture Micro-C (RCMC), a hybridization-based capture Micro-C method, allowed us to deeply map the highly nested interactions between regulatory elements. By utilizing RCMC, we generated high-resolution, focal contact maps that revealed comprehensive enhancer-promoter interactions near PSEs in prostate cancer cells. Interestingly, we found that PSEs loop to other PSEs clustered in the proximal genomic region, in addition to promoters of significantly upregulated genes in cancer compared to normal cells. We termed these regions as multi-connected "3D enhancer-promoter hubs”. To mechanistically dissect these hubs, we deleted 5+ PSEs one by one using the CRISPR/Cas9 system and performed multi-omics analyses, including RNA-seq and ChIP-seq. We found that single deletions of a subset of PSEs led to the downregulation of multiple genes looped to the deleted PSE as well as to other PSEs. Moreover, we found that deletion reduced the activity of other connected PSEs. When we further evaluated the chromatin contacts of regulatory elements in the PSE-deleted cells, we observed that many chromatin contacts remained unaffected. Finally, we identified TFs that bind to PSEs, and we are currently in the process of investigating enhancer-promoter hub connectivity by performing knockdown experiments on these TFs. Furthermore, we are using molecular assays to identify PSEs that could exert a phenotypic shift in cancer cell identity to drive prostate tumorigenesis. Overall, our study not only provides mechanistic insights into the spatial regulation of 3D enhancer-promoter hubs but also identifies new non-coding drivers in prostate tumorigenesis. Citation Format: Huan Cao, Zexun Wu, Leonardo Gonzalez-Smith, Seolyn Yang, Suhn K. Rhie. Dissecting the 3D enhancer-promoter hubs in prostate tumorigenesis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1686.

Deciphering the spatiotemporal transcriptional and chromatin accessibility of human retinal organoid development at the single-cell level

Molecular information on the early stages of human retinal development remains scarce due to limitations in obtaining early human eye samples. Pluripotent stem cell-derived retinal organoids (ROs) provide an unprecedented opportunity for studying early retinogenesis. Using a combination of single cell RNA-seq and spatial transcriptomics we present for the first-time a single cell spatiotemporal transcriptome of RO development. Our data demonstrate that ROs recapitulate key events of retinogenesis including optic vesicle/cup formation, presence of a putative ciliary margin zone, emergence of retinal progenitor cells and their orderly differentiation to retinal neurons. Combining the scRNA- with scATAC-seq data, we were able to reveal cell-type specific transcription factor binding motifs on accessible chromatin at each stage of organoid development, and to show that chromatin accessibility is highly correlated to the developing human retina, but with some differences in the temporal emergence and abundance of some of the retinal neurons.

Exploring intrinsic variability between cultured nasal and bronchial epithelia in cystic fibrosis

The nasal and bronchial epithelium are unified parts of the respiratory tract that are affected in the monogenic disorder cystic fibrosis (CF). Recent studies have uncovered that nasal and bronchial tissues exhibit intrinsic variability, including differences in mucociliary cell composition and expression of unique transcriptional regulatory proteins which relate to germ layer origin. In the present study, we explored whether intrinsic differences between nasal and bronchial epithelial cells persist in cell cultures and affect epithelial cell functioning in CF. Comparison of air–liquid interface (ALI) differentiated epithelial cells from subjects with CF revealed distinct mucociliary differentiation states of nasal and bronchial cultures. Moreover, using RNA sequencing we identified cell type-specific signature transcription factors in differentiated nasal and bronchial epithelial cells, some of which were already poised for expression in basal progenitor cells as evidenced by ATAC sequencing. Analysis of differentiated nasal and bronchial epithelial 3D organoids revealed distinct capacities for fluid secretion, which was linked to differences in ciliated cell differentiation. In conclusion, we show that unique phenotypical and functional features of nasal and bronchial epithelial cells persist in cell culture models, which can be further used to investigate the effects of tissue-specific features on upper and lower respiratory disease development in CF.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease.

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Analysis of potential biomarkers for diabetic kidney disease based on single-cell RNA-sequencing integrated with a single-cell sequencing assay for transposase-accessible chromatin.

Diabetic kidney disease (DKD) is a renal microvascular disease caused by hyperglycemia that involves metabolic remodeling, oxidative stress, inflammation, and other factors. The mechanism is complex and not fully unraveled. We performed an integrated single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) and single-cell RNA-sequencing (scRNA-seq) analyses of kidneys from db/db and db/m mice to identify differential open chromatin regions and gene expression, particularly in genes related to proximal tubular reabsorption and secretion. We identified 9,776 differentially expressed genes (DEGs) and 884 cell type-specific transcription factors (TFs) across 15 cell types. Glucose and lipid transporters, and TFs related to the circadian rhythm in the proximal tubules had significantly higher expression in db/db mice than in db/m mice (P<0.01). Crosstalk between podocytes and tubular cells in the proximal tubules was enhanced, and renal inflammation, oxidative stress, and fibrosis pathways were activated in db/db mice. Western blotting and immunohistochemical staining results showed that Wfdc2 expression in the urine and kidneys of DKD patients was higher than that in non-diabetic kidney disease (NDKD) controls. The revealed landscape of chromatin accessibility and transcriptional profiles in db/db mice provide insights into the pathological mechanism of DKD.

Single-cell ATAC sequencing illuminates the cis-regulatory differentiation of taxol biosynthesis between leaf mesophyll and leaf epidermal cells in Taxus mairei

Taxol is a widely used anticancer drug; however, its accumulation pattern can vary largely among leaf tissues. Taxol biosynthesis is driven by a series of transcription factors (TFs), but the extent of chromatin changes that functionally regulate gene expression associated with taxol biosynthesis remain poorly understood. To assess the effects of tissue heterogeneity on taxol biosynthesis, we applied single-cell transposase-accessible chromatin sequencing to Taxus mairei leaves and identified differentially accessible regions. We captured 9488 independent cells with 42,192 fragment reads per cell. In total, 66.64% were confidently mapped onto the T. mairei genome. Projection analysis divided all the cells into nine clusters, which suggests a high degree of cell heterogeneity in T. mairei leaves. Based on the cell type markers, Clusters 2 and 7 were annotated as leaf mesophyll cells, and Clusters 5 and 6 were annotated as leaf epidermal cells. A differential accessibility analysis identified 9600 and 8538 leaf mesophyll and epidermal, respectively, cell-specific open chromatin regions. We predicted the genomic features of the genes involved in different stages of taxol biosynthesis, and there was a close correlation between chromatin open peaks and the expression levels of several taxol biosynthesis-related genes. We further identified a number of cell type-specific TF motifs and reference TFs that appeared to be involved in taxol biosynthesis. The study of T. mairei leaves at a single-cell resolution provides a valuable resource for determining the basic principles of the cell type-specific regulation of taxol biosynthesis.

Smarca4 deficiency induces Pttg1 oncogene upregulation and hyperproliferation of tubular and interstitial cells during kidney development.

Kidney formation and nephrogenesis are controlled by precise spatiotemporal gene expression programs, which are coordinately regulated by cell-cycle, cell type-specific transcription factors and epigenetic/chromatin regulators. However, the roles of epigenetic/chromatin regulators in kidney development and disease remain poorly understood. In this study, we investigated the impact of deleting the chromatin remodeling factor Smarca4 (Brg1), a human Wilms tumor-associated gene, in Wnt4-expressing cells. Smarca4 deficiency led to severe tubular defects and a shortened medulla. Through unbiased single-cell RNA sequencing analyses, we identified multiple types of Wnt4 Cre-labeled interstitial cells, along with nephron-related cells. Smarca4 deficiency increased interstitial cells but markedly reduced tubular cells, resulting in cells with mixed identity and elevated expression of cell-cycle regulators and genes associated with extracellular matrix and epithelial-to-mesenchymal transition/fibrosis. We found that Smarca4 loss induced a significant upregulation of the oncogene Pttg1 and hyperproliferation of Wnt4 Cre-labeled cells. These changes in the cellular state could hinder the cellular transition into characteristic tubular structures, eventually leading to fibrosis. In conclusion, our findings shed light on novel cell types and genes associated with Wnt4 Cre-labeled cells and highlight the critical role of Smarca4 in regulating tubular cell differentiation and the expression of the cancer-causing gene Pttg1 in the kidney. These findings may provide valuable insights into potential therapeutic strategies for renal cell carcinoma resulting from SMARCA4 deficiency.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3.

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

METB-03. THE DISTINCT GENOMICS OF HIGH-GRADE INFANT GLIOMAS AMONG PEDIATRIC GLIOMAS HARBORING RECEPTOR TYROSINE KINASE FUSIONS

Abstract Gliomas account for about 60% of pediatric brain tumors, about half of which are classified as high-grade. Recent studies have shown that infants with high-grade gliomas (HGG) harboring receptor tyrosine kinase (RTK) fusions have a more favorable prognosis than pediatric HGG patients. The molecular mechanisms for this favorable, but from the histology divergent clinical outcome, are yet unknown. We analyzed n=18 brain tumor samples (n=7 infant, n=11 pediatric) of patients with low-grade glioma (LGGs) or HGG presenting with RTK fusions using single-cell RNA sequencing. Considering the different fusion genes (ALK n=6, NTRK n=8, and MET n=4), age groups, and glioma grades, we evaluated the metaprograms of tumor cells and their proportions using non-negative matrix factorization. Preliminary results show similar metaprograms comparing infant and pediatric LGG. Cells detected as MES-like did not occur in LGG but were found in several HGG samples. However, in HGG, the MES-like cells show different expression patterns for infant and pediatric samples. In addition, compared with pediatric HGG, infant HGG has a higher proportion of programs that are predominantly present in LGG, such as AC-like and NPC-like. These initial findings serve as a starting point for cell type-specific transcription factor regulatory networks and are an important step toward understanding the underlying biology of this distinct pediatric brain tumor entity.

Single-Cell RNA Sequencing of Sox17-Expressing Lineages Reveals Distinct Gene Regulatory Networks and Dynamic Developmental Trajectories.

During early embryogenesis, the transcription factor SOX17 contributes to hepato-pancreato-biliary system formation and vascular-hematopoietic emergence. To better understand Sox17 function in the developing endoderm and endothelium, we developed a dual-color temporal lineage-tracing strategy in mice combined with single-cell RNA sequencing to analyze 6934 cells from Sox17-expressing lineages at embryonic days 9.0-9.5. Our analyses showed 19 distinct cellular clusters combined from all 3 germ layers. Differential gene expression, trajectory and RNA-velocity analyses of endothelial cells revealed a heterogenous population of uncommitted and specialized endothelial subtypes, including 2 hemogenic populations that arise from different origins. Similarly, analyses of posterior foregut endoderm revealed subsets of hepatic, pancreatic, and biliary progenitors with overlapping developmental potency. Calculated gene-regulatory networks predict gene regulons that are dominated by cell type-specific transcription factors unique to each lineage. Vastly different Sox17 regulons found in endoderm versus endothelial cells support the differential interactions of SOX17 with other regulatory factors thereby enabling lineage-specific regulatory actions.