ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation

Identifying a gene signature of metastatic potential by linking pre-metastatic state to ultimate metastatic fate

Identifying the key molecular pathways that enable metastasis by analyzing the eventual metastatic tumor is challenging because the state of the founder subclone likely changes following metastatic colonization. To address this challenge, we labeled primary mouse pancreatic ductal adenocarcinoma (PDAC) subclones with DNA barcodes to characterize their pre-metastatic state using ATAC-seq and RNA-seq and determine their relative in vivo metastatic potential prospectively. We identified a gene signature separating metastasis-high and metastasis-low subclones orthogonal to the normal-to-PDAC and classical-to-basal axes. The metastasis-high subclones feature activation of IL-1 pathway genes and high NF-κB and Zeb/Snail family activity and the metastasis-low subclones feature activation of neuroendocrine, motility, and Wnt pathway genes and high CDX2 and HOXA13 activity. In a functional screen, we validated novel mediators of PDAC metastasis in the IL-1 pathway, including the NF-κB targets Fos and Il23a, and beyond the IL-1 pathway including Myo1b and Tmem40. We scored human PDAC tumors for our signature of metastatic potential from mouse and found that metastases have higher scores than primary tumors. Moreover, primary tumors with higher scores are associated with worse prognosis. We also found that our metastatic potential signature is enriched in other human carcinomas, suggesting that it is conserved across epithelial malignancies. This work establishes a strategy for linking cancer cell state to future behavior, reveals novel functional regulators of PDAC metastasis, and establishes a method for scoring human carcinomas based on metastatic potential.

Integrative multi-omics increase resolution of the sea urchin posterior gut gene regulatory network at single-cell level

Drafting gene regulatory networks (GRNs) requires embryological knowledge pertaining to the cell type families, information on the regulatory genes, causal data from gene knockdown experiments and validations of the identified interactions by cis-regulatory analysis. We use multi-omics involving next-generation sequencing to obtain the necessary information for drafting the Strongylocentrotus purpuratus (Sp) posterior gut GRN. Here, we present an update to the GRN using: (1) a single-cell RNA-sequencing-derived cell atlas highlighting the 2 day-post-fertilization (dpf) sea urchin gastrula cell type families, as well as the genes expressed at the single-cell level; (2) a set of putative cis-regulatory modules and transcription factor-binding sites obtained from chromatin accessibility ATAC-seq data; and (3) interactions directionality obtained from differential bulk RNA sequencing following knockdown of the transcription factor Sp-Pdx1, a key regulator of gut patterning in sea urchins. Combining these datasets, we draft the GRN for the hindgut Sp-Pdx1-positive cells in the 2 dpf gastrula embryo. Overall, our data suggest the complex connectivity of the posterior gut GRN and increase the resolution of gene regulatory cascades operating within it.

Multimodal single-cell analyses reveal mechanisms of perianal fistula in diverse patients with Crohn's disease

BACKGROUND

Crohn's disease complicated by perianal fistulae is more prevalent and severe in patients of African ancestry.

METHODS

We profiled single cells from diverse patients with Crohn's disease with perianal fistula from colorectal mucosa and fistulous tracts. Immunofluorescence was performed to validate predicted cell-cell interactions. Unstimulated monocytes were chronically cultured in diverse cohorts. A subset was analyzed by single-nucleus RNA + ATAC sequencing.

FINDINGS

Fistulous tract cells from complete proctectomies demonstrated enrichment of myeloid cells compared to paired rectal tissues. Ligand-receptor analysis highlights myeloid-stromal cross-talk and cellular senescence, with cellular co-localization validated by immunofluorescence. Chitinase-3 like-protein-1 (CHI3L1) is a top upregulated gene in stromal cells from fistulae expressing both destructive and fibrotic gene signatures. Monocyte cultures from patients of African ancestry and controls demonstrated differences in CHI3L1 and oncostatin M (OSM) expression upon differentiation compared to individuals of European ancestry. Activating protein-1 footprints are present in ATAC-seq peaks in stress response genes, including CHI3L1 and OSM; genome-wide chromatin accessibility including JUN footprints was observed, consistent with reported mechanisms of inflammatory memory. Regulon analyses confirm known cell-specific transcription factor regulation and implicate novel ones in fibroblast subsets. All pseudo-bulked clusters demonstrate enrichment of genetic loci, establishing multicellular contributions. In the most significant African American Crohn's genetic locus, upstream of prostaglandin E receptor 4, lymphoid-predominant ATAC-seq peaks were observed, with predicted RORC footprints.

CONCLUSIONS

Population differences in myeloid-stromal cross-talk implicate fibrotic and destructive fibroblasts, senescence, epigenetic memory, and cell-specific enhancers in perianal fistula pathogenesis. The transcriptomic and epigenetic data provided here may guide optimization of promising mesenchymal stem cell therapies for perianal fistula.

FUNDING

This work was supported by grants U01DK062422, U24DK062429, and R01DK123758.

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Comprehensive mapping and modelling of the rice regulome landscape unveils the regulatory architecture underlying complex traits

Unraveling the regulatory mechanisms that govern complex traits is pivotal for advancing crop improvement. Here we present a comprehensive regulome atlas for rice (Oryza sativa), charting the chromatin accessibility across 23 distinct tissues from three representative varieties. Our study uncovers 117,176 unique open chromatin regions (OCRs), accounting for ~15% of the rice genome, a notably higher proportion compared to previous reports in plants. Integrating RNA-seq data from matched tissues, we confidently predict 59,075 OCR-to-gene links, with enhancers constituting 69.54% of these associations, including many known enhancer-to-gene links. Leveraging this resource, we re-evaluate genome-wide association study results and discover a previously unknown function of OsbZIP06 in seed germination, which we subsequently confirm through experimental validation. We optimize deep learning models to decode regulatory grammar, achieving robust modeling of tissue-specific chromatin accessibility. This approach allows to predict cross-variety regulatory dynamics from genomic sequences, shedding light on the genetic underpinnings of cis-regulatory divergence and morphological disparities between varieties. Overall, our study establishes a foundational resource for rice functional genomics and precision molecular breeding, providing valuable insights into regulatory mechanisms governing complex traits.

A single-cell multimodal view on gene regulatory network inference from transcriptomics and chromatin accessibility data

Eukaryotic gene regulation is a combinatorial, dynamic, and quantitative process that plays a vital role in development and disease and can be modeled at a systems level in gene regulatory networks (GRNs). The wealth of multi-omics data measured on the same samples and even on the same cells has lifted the field of GRN inference to the next stage. Combinations of (single-cell) transcriptomics and chromatin accessibility allow the prediction of fine-grained regulatory programs that go beyond mere correlation of transcription factor and target gene expression, with enhancer GRNs (eGRNs) modeling molecular interactions between transcription factors, regulatory elements, and target genes. In this review, we highlight the key components for successful (e)GRN inference from (sc)RNA-seq and (sc)ATAC-seq data exemplified by state-of-the-art methods as well as open challenges and future developments. Moreover, we address preprocessing strategies, metacell generation and computational omics pairing, transcription factor binding site detection, and linear and three-dimensional approaches to identify chromatin interactions as well as dynamic and causal eGRN inference. We believe that the integration of transcriptomics together with epigenomics data at a single-cell level is the new standard for mechanistic network inference, and that it can be further advanced with integrating additional omics layers and spatiotemporal data, as well as with shifting the focus towards more quantitative and causal modeling strategies.

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

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

Brain and cancer associated binding domain mutations provide insight into CTCF's relationship with chromatin and its ability to act as a chromatin organizer

Although only a fraction of CTCF motifs are bound in any cell type, and approximately half of the occupied sites overlap cohesin, the mechanisms underlying cell-type specific attachment and ability to function as a chromatin organizer remain unknown. To investigate the relationship between CTCF and chromatin we applied a combination of imaging, structural and molecular approaches, using a series of brain and cancer associated CTCF mutations that act as CTCF perturbations. We demonstrate that binding and the functional impact of WT and mutant CTCF depend not only on the unique properties of each protein, but also on the genomic context of bound sites. Our studies also highlight the reciprocal relationship between CTCF and chromatin, demonstrating that the unique binding properties of WT and mutant proteins have a distinct impact on accessibility, TF binding, cohesin overlap, chromatin interactivity and gene expression programs, providing insight into their cancer and brain related effects.

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully functional new organs, including new brains, de novo. The regeneration of a new brain requires the formation of diverse neural cell types and their assembly into an organized structure with correctly wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neural subpopulations, however, how these transcriptional programs are initiated in response to injury remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia, to study wound-induced transcriptional regulatory events that lead to the production of neurons and subsequently a functional brain. Footprinting analysis using chromatin accessibility data on a chromosome-scale genome assembly revealed that binding sites for the Nuclear Factor Y (NFY) transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal function. Single-cell transcriptome analysis combined with functional studies identified soxC+ stem cells as a putative progenitor population for multiple neural subtypes. Further, we found that wound-induced soxC expression is likely under direct transcriptional control by NFY, uncovering a mechanism for the initiation of a neural differentiation pathway by early wound-induced binding of a transcriptional regulator.

Intergenic risk variant rs56258221 skews the fate of naive CD4+ T cells via miR4464-BACH2 interplay in primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is an immune-mediated liver disease of unknown pathogenesis, with a high risk to develop cirrhosis and malignancies. Functional dysregulation of T cells and association with genetic polymorphisms in T cell-related genes were previously reported for PSC. Here, we genotyped a representative PSC cohort for several disease-associated risk loci and identified rs56258221 (BACH2/MIR4464) to correlate with not only the peripheral blood T cell immunophenotype but also the functional capacities of naive CD4+ T (CD4+ TN) cells in people with PSC. Mechanistically, rs56258221 leads to an increased expression of miR4464, in turn causing attenuated translation of BACH2, a major gatekeeper of T cell quiescence. Thereby, the fate of CD4+ TN is skewed toward polarization into pro-inflammatory subsets. Clinically, people with PSC carrying rs56258221 show signs of accelerated disease progression. The data presented here highlight the importance of assigning functional outcomes to disease-associated genetic polymorphisms as potential drivers of diseases.

PRDM6 promotes medulloblastoma by repressing chromatin accessibility and altering gene expression

SNCAIP duplication may promote Group 4 medulloblastoma via induction of PRDM6, a poorly characterized member of the PRDF1 and RIZ1 homology domain-containing (PRDM) family of transcription factors. Here, we investigated the function of PRDM6 in human hindbrain neuroepithelial stem cells and tested PRDM6 as a driver of Group 4 medulloblastoma. We report that human PRDM6 localizes predominantly to the nucleus, where it causes widespread repression of chromatin accessibility and complex alterations of gene expression patterns. Genome-wide mapping of PRDM6 binding reveals that PRDM6 binds to chromatin regions marked by histone H3 lysine 27 trimethylation that are located within, or proximal to, genes. Moreover, we show that PRDM6 expression in neuroepithelial stem cells promotes medulloblastoma. Surprisingly, medulloblastomas derived from PRDM6-expressing neuroepithelial stem cells match human Group 3, but not Group 4, medulloblastoma. We conclude that PRDM6 expression has oncogenic potential but is insufficient to drive Group 4 medulloblastoma from neuroepithelial stem cells. We propose that both PRDM6 and additional factors, such as specific cell-of-origin features, are required for Group 4 medulloblastoma. Given the lack of PRDM6 expression in normal tissues and its oncogenic potential shown here, we suggest that PRDM6 inhibition may have therapeutic value in PRDM6-expressing medulloblastomas.

Epigenomic landscape of the human dorsal root ganglion: sex differences and transcriptional regulation of nociceptive genes

Gene expression is influenced by chromatin architecture via controlled access of regulatory factors to DNA. To better understand gene regulation in the human dorsal root ganglion (hDRG) we used bulk and spatial transposase-accessible chromatin technology followed by sequencing (ATAC-seq). Using bulk ATAC-seq, we detected that in females diverse differentially accessible chromatin regions (DARs) mapped to the X chromosome and in males to autosomal genes. EGR1/3 and SP1/4 transcription factor binding motifs were abundant within DARs in females, and JUN, FOS and other AP-1 factors in males. To dissect the open chromatin profile in hDRG neurons, we used spatial ATAC-seq. The neuron cluster showed higher chromatin accessibility in GABAergic, glutamatergic, and interferon-related genes in females, and in Ca2+- signaling-related genes in males. Sex differences in transcription factor binding sites in neuron-proximal barcodes were consistent with the trends observed in bulk ATAC-seq data. We validated that EGR1 expression is biased to female hDRG compared to male. Strikingly, XIST, the long-noncoding RNA responsible for X inactivation, hybridization signal was found to be highly dispersed in the female neuronal but not non-neuronal nuclei suggesting weak X inactivation in female hDRG neurons. Our findings point to baseline epigenomic sex differences in the hDRG that likely underlie divergent transcriptional responses that determine mechanistic sex differences in pain.

Single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of single nucleus RNA sequencing (snRNA-seq) offers to revolutionize the study of Alzheimer's disease (AD). Integration with complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link cell subpopulations and molecular networks with a broader disease-relevant context. We report snRNA-seq profiles from superior frontal gyrus samples from 101 well characterized subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in hematological parameters. We observed an AD-associated CD83(+) microglial subtype with unique molecular networks and which is associated with immunoglobulin IgG4 production in the transverse colon. Our major observations were replicated in two additional, independent snRNA-seq data sets. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal disease biology.

Transposase-assisted tagmentation: an economical and scalable strategy for single-worm whole-genome sequencing

AlphaMissense identifies 23 million human missense variants as likely pathogenic, but only 0.1% have been clinically classified. To experimentally validate these predictions, chemical mutagenesis presents a rapid, cost-effective method to produce billions of mutations in model organisms. However, the prohibitive costs and limitations in the throughput of whole-genome sequencing (WGS) technologies, crucial for variant identification, constrain its widespread application. Here, we introduce a Tn5 transposase-assisted tagmentation technique for conducting WGS in Caenorhabditis elegans, Escherichia coli, Saccharomyces cerevisiae, and Chlamydomonas reinhardtii. This method, demands merely 20 min of hands-on time for a single-worm or single-cell clones and incurs a cost below 10 US dollars. It effectively pinpoints causal mutations in mutants defective in cilia or neurotransmitter secretion and in mutants synthetically sterile with a variant analogous to the B-Raf Proto-oncogene, Serine/Threonine Kinase (BRAF) V600E mutation. Integrated with chemical mutagenesis, our approach can generate and identify missense variants economically and efficiently, facilitating experimental investigations of missense variants in diverse species.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

TWEAK/Fn14 signalling driven super-enhancer reprogramming promotes pro-metastatic metabolic rewiring in triple-negative breast cancer

Triple Negative Breast Cancer (TNBC) is the most aggressive breast cancer subtype suffering from limited targeted treatment options. Following recent reports correlating Fibroblast growth factor-inducible 14 (Fn14) receptor overexpression in Estrogen Receptor (ER)-negative breast cancers with metastatic events, we show that Fn14 is specifically overexpressed in TNBC patients and associated with poor survival. We demonstrate that constitutive Fn14 signalling rewires the transcriptomic and epigenomic landscape of TNBC, leading to enhanced tumour growth and metastasis. We further illustrate that such mechanisms activate TNBC-specific super enhancers (SE) to drive the transcriptional activation of cancer dependency genes via chromatin looping. In particular, we uncover the SE-driven upregulation of Nicotinamide phosphoribosyltransferase (NAMPT), which promotes NAD+ and ATP metabolic reprogramming critical for filopodia formation and metastasis. Collectively, our study details the complex mechanistic link between TWEAK/Fn14 signalling and TNBC metastasis, which reveals several vulnerabilities which could be pursued for the targeted treatment of TNBC patients.

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.

Chromatin modifier Hmga2 promotes adult hematopoietic stem cell function and blood regeneration in stress conditions

The molecular mechanisms governing the response of hematopoietic stem cells (HSCs) to stress insults remain poorly defined. Here, we investigated effects of conditional knock-out or overexpression of Hmga2 (High mobility group AT-hook 2), a transcriptional activator of stem cell genes in fetal HSCs. While Hmga2 overexpression did not affect adult hematopoiesis under homeostasis, it accelerated HSC expansion in response to injection with 5-fluorouracil (5-FU) or in vitro treatment with TNF-α. In contrast, HSC and megakaryocyte progenitor cell numbers were decreased in Hmga2 KO animals. Transcription of inflammatory genes was repressed in Hmga2-overexpressing mice injected with 5-FU, and Hmga2 bound to distinct regions and chromatin accessibility was decreased in HSCs upon stress. Mechanistically, we found that casein kinase 2 (CK2) phosphorylates the Hmga2 acidic domain, promoting its access and binding to chromatin, transcription of anti-inflammatory target genes, and the expansion of HSCs under stress conditions. Notably, the identified stress-regulated Hmga2 gene signature is activated in hematopoietic stem progenitor cells of human myelodysplastic syndrome patients. In sum, these results reveal a TNF-α/CK2/phospho-Hmga2 axis controlling adult stress hematopoiesis.

REUNION: transcription factor binding prediction and regulatory association inference from single-cell multi-omics data

MOTIVATION

Profiling of gene expression and chromatin accessibility by single-cell multi-omics approaches can help to systematically decipher how transcription factors (TFs) regulate target gene expression via cis-region interactions. However, integrating information from different modalities to discover regulatory associations is challenging, in part because motif scanning approaches miss many likely TF binding sites.

RESULTS

We develop REUNION, a framework for predicting genome-wide TF binding and cis-region-TF-gene "triplet" regulatory associations using single-cell multi-omics data. The first component of REUNION, Unify, utilizes information theory-inspired complementary score functions that incorporate TF expression, chromatin accessibility, and target gene expression to identify regulatory associations. The second component, Rediscover, takes Unify estimates as input for pseudo semi-supervised learning to predict TF binding in accessible genomic regions that may or may not include detected TF motifs. Rediscover leverages latent chromatin accessibility and sequence feature spaces of the genomic regions, without requiring chromatin immunoprecipitation data for model training. Applied to peripheral blood mononuclear cell data, REUNION outperforms alternative methods in TF binding prediction on average performance. In particular, it recovers missing region-TF associations from regions lacking detected motifs, which circumvents the reliance on motif scanning and facilitates discovery of novel associations involving potential co-binding transcriptional regulators. Newly identified region-TF associations, even in regions lacking a detected motif, improve the prediction of target gene expression in regulatory triplets, and are thus likely to genuinely participate in the regulation.

AVAILABILITY AND IMPLEMENTATION

All source code is available at https://github.com/yangymargaret/REUNION.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development

STUDY QUESTION

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development?

SUMMARY ANSWER

Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce.

STUDY DESIGN, SIZE, DURATION

This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome.

MAIN RESULTS AND THE ROLE OF CHANCE

A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation.

LARGE SCALE DATA

All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501.

LIMITATIONS, REASONS FOR CAUTION

Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells.

WIDER IMPLICATIONS OF THE FINDINGS

The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

Mapping putative enhancers in mouse oocytes and early embryos reveals TCF3/12 as key folliculogenesis regulators

Dynamic epigenomic reprogramming occurs during mammalian oocyte maturation and early development. However, the underlying transcription circuitry remains poorly characterized. By mapping cis-regulatory elements using H3K27ac, we identified putative enhancers in mouse oocytes and early embryos distinct from those in adult tissues, enabling global transitions of regulatory landscapes around fertilization and implantation. Gene deserts harbour prevalent putative enhancers in fully grown oocytes linked to oocyte-specific genes and repeat activation. Embryo-specific enhancers are primed before zygotic genome activation and are restricted by oocyte-inherited H3K27me3. Putative enhancers in oocytes often manifest H3K4me3, bidirectional transcription, Pol II binding and can drive transcription in STARR-seq and a reporter assay. Finally, motif analysis of these elements identified crucial regulators of oogenesis, TCF3 and TCF12, the deficiency of which impairs activation of key oocyte genes and folliculogenesis. These data reveal distinctive regulatory landscapes and their interacting transcription factors that underpin the development of mammalian oocytes and early embryos.

PRMT1 promotes epigenetic reprogramming associated with acquired chemoresistance in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) carries a dismal prognosis due to therapeutic resistance. We show that PDAC cells undergo global epigenetic reprogramming to acquire chemoresistance, a process that is driven at least in part by protein arginine methyltransferase 1 (PRMT1). Genetic or pharmacological PRMT1 inhibition impairs adaptive epigenetic reprogramming and delays acquired resistance to gemcitabine and other common chemo drugs. Mechanistically, gemcitabine treatment induces translocation of PRMT1 into the nucleus, where its enzymatic activity limits the assembly of chromatin-bound MAFF/BACH1 transcriptional complexes. Cut&Tag chromatin profiling of H3K27Ac, MAFF, and BACH1 suggests a pivotal role for MAFF/BACH1 in global epigenetic response to gemcitabine, which is confirmed by genetically silencing MAFF. PRMT1 and MAFF/BACH1 signature genes identified by Cut&Tag analysis distinguish gemcitabine-resistant from gemcitabine-sensitive patient-derived xenografts of PDAC, supporting the PRMT1-MAFF/BACH1 epigenetic regulatory axis as a potential therapeutic avenue for improving the efficacy and durability of chemotherapies in patients of PDAC.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations

BACKGROUND

Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype.

RESULTS

In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases.

CONCLUSIONS

Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Molecular cascades and cell type-specific signatures in ASD revealed by single-cell genomics

Genomic profiling in postmortem brain from autistic individuals has consistently revealed convergent molecular changes. What drives these changes and how they relate to genetic susceptibility in this complex condition are not well understood. We performed deep single-nucleus RNA sequencing (snRNA-seq) to examine cell composition and transcriptomics, identifying dysregulation of cell type-specific gene regulatory networks (GRNs) in autism spectrum disorder (ASD), which we corroborated using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) and spatial transcriptomics. Transcriptomic changes were primarily cell type specific, involving multiple cell types, most prominently interhemispheric and callosal-projecting neurons, interneurons within superficial laminae, and distinct glial reactive states involving oligodendrocytes, microglia, and astrocytes. Autism-associated GRN drivers and their targets were enriched in rare and common genetic risk variants, connecting autism genetic susceptibility and cellular and circuit alterations in the human brain.

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Viral reprogramming of host transcription initiation

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Bone marrow stromal cells induce chromatin remodeling in multiple myeloma cells leading to transcriptional changes

The natural history of multiple myeloma is characterized by its localization to the bone marrow and its interaction with bone marrow stromal cells. The bone marrow stromal cells provide growth and survival signals, thereby promoting the development of drug resistance. Here, we show that the interaction between bone marrow stromal cells and myeloma cells (using human cell lines) induces chromatin remodeling of cis-regulatory elements and is associated with changes in the expression of genes involved in the cell migration and cytokine signaling. The expression of genes involved in these stromal interactions are observed in extramedullary disease in patients with myeloma and provides the rationale for survival of myeloma cells outside of the bone marrow microenvironment. Expression of these stromal interaction genes is also observed in a subset of patients with newly diagnosed myeloma and are akin to the transcriptional program of extramedullary disease. The presence of such adverse stromal interactions in newly diagnosed myeloma is associated with accelerated disease dissemination, predicts the early development of therapeutic resistance, and is of independent prognostic significance. These stromal cell induced transcriptomic and epigenomic changes both predict long-term outcomes and identify therapeutic targets in the tumor microenvironment for the development of novel therapeutic approaches.

JAK-STAT signaling maintains homeostasis in T cells and macrophages

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Uncovering uncharacterized binding of transcription factors from ATAC-seq footprinting data

Transcription factors (TFs) are crucial epigenetic regulators, which enable cells to dynamically adjust gene expression in response to environmental signals. Computational procedures like digital genomic footprinting on chromatin accessibility assays such as ATACseq can be used to identify bound TFs in a genome-wide scale. This method utilizes short regions of low accessibility signals due to steric hindrance of DNA bound proteins, called footprints (FPs), which are combined with motif databases for TF identification. However, while over 1600 TFs have been described in the human genome, only ~ 700 of these have a known binding motif. Thus, a substantial number of FPs without overlap to a known DNA motif are normally discarded from FP analysis. In addition, the FP method is restricted to organisms with a substantial number of known TF motifs. Here we present DENIS (DE Novo motIf diScovery), a framework to generate and systematically investigate the potential of de novo TF motif discovery from FPs. DENIS includes functionality (1) to isolate FPs without binding motifs, (2) to perform de novo motif generation and (3) to characterize novel motifs. Here, we show that the framework rediscovers artificially removed TF motifs, quantifies de novo motif usage during an early embryonic development example dataset, and is able to analyze and uncover TF activity in organisms lacking canonical motifs. The latter task is exemplified by an investigation of a scATAC-seq dataset in zebrafish which covers different cell types during hematopoiesis.

Detection of newly synthesized RNA reveals transcriptional reprogramming during ZGA and a role of Obox3 in totipotency acquisition

Zygotic genome activation (ZGA) after fertilization enables the maternal-to-zygotic transition. However, the global view of ZGA, particularly at initiation, is incompletely understood. Here, we develop a method to capture and sequence newly synthesized RNA in early mouse embryos, providing a view of transcriptional reprogramming during ZGA. Our data demonstrate that major ZGA gene activation begins earlier than previously thought. Furthermore, we identify a set of genes activated during minor ZGA, the promoters of which show enrichment of the Obox factor motif, and find that Obox3 or Obox5 overexpression in mouse embryonic stem cells activates ZGA genes. Notably, the expression of Obox factors is severely impaired in somatic cell nuclear transfer (SCNT) embryos, and restoration of Obox3 expression corrects the ZGA profile and greatly improves SCNT embryo development. Hence, our study reveals dynamic transcriptional reprogramming during ZGA and underscores the crucial role of Obox3 in facilitating totipotency acquisition.

The genome regulatory landscape of Atlantic salmon liver through smoltification

The anadromous Atlantic salmon undergo a preparatory physiological transformation before seawater entry, referred to as smoltification. Key molecular developmental processes involved in this life stage transition, such as remodeling of gill functions, are known to be synchronized and modulated by environmental cues like photoperiod. However, little is known about the photoperiod influence and genome regulatory processes driving other canonical aspects of smoltification such as the large-scale changes in lipid metabolism and energy homeostasis in the developing smolt liver. Here we generate transcriptome, DNA methylation, and chromatin accessibility data from salmon livers across smoltification under different photoperiod regimes. We find a systematic reduction of expression levels of genes with a metabolic function, such as lipid metabolism, and increased expression of energy related genes such as oxidative phosphorylation, during smolt development in freshwater. However, in contrast to similar studies of the gill, smolt liver gene expression prior to seawater transfer was not impacted by photoperiodic history. Integrated analyses of gene expression, chromatin accessibility, and transcription factor (TF) binding signatures highlight chromatin remodeling and TF dynamics underlying smolt gene regulatory changes. Differential peak accessibility patterns largely matched differential gene expression patterns during smoltification and we infer that ZNF682, KLFs, and NFY TFs are important in driving a liver metabolic shift from synthesis to break down of organic compounds in freshwater. Overall, chromatin accessibility and TFBS occupancy were highly correlated to changes in gene expression. On the other hand, we identified numerous differential methylation patterns across the genome, but associated genes were not functionally enriched or correlated to observed gene expression changes across smolt development. Taken together, this work highlights the relative importance of chromatin remodeling during smoltification and demonstrates that metabolic remodeling occurs as a preadaptation to life at sea that is not to a large extent driven by photoperiod history.

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

BACKGROUND

Motif finding in Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data is essential to reveal the intricacies of transcription factor binding sites (TFBSs) and their pivotal roles in gene regulation. Deep learning technologies including convolutional neural networks (CNNs) and graph neural networks (GNNs), have achieved success in finding ATAC-seq motifs. However, CNN-based methods are limited by the fixed width of the convolutional kernel, which makes it difficult to find multiple transcription factor binding sites with different lengths. GNN-based methods has the limitation of using the edge weight information directly, makes it difficult to aggregate the neighboring nodes' information more efficiently when representing node embedding.

RESULTS

To address this challenge, we developed a novel graph attention network framework named MMGAT, which employs an attention mechanism to adjust the attention coefficients among different nodes. And then MMGAT finds multiple ATAC-seq motifs based on the attention coefficients of sequence nodes and k-mer nodes as well as the coexisting probability of k-mers. Our approach achieved better performance on the human ATAC-seq datasets compared to existing tools, as evidenced the highest scores on the precision, recall, F1_score, ACC, AUC, and PRC metrics, as well as finding 389 higher quality motifs. To validate the performance of MMGAT in predicting TFBSs and finding motifs on more datasets, we enlarged the number of the human ATAC-seq datasets to 180 and newly integrated 80 mouse ATAC-seq datasets for multi-species experimental validation. Specifically on the mouse ATAC-seq dataset, MMGAT also achieved the highest scores on six metrics and found 356 higher-quality motifs. To facilitate researchers in utilizing MMGAT, we have also developed a user-friendly web server named MMGAT-S that hosts the MMGAT method and ATAC-seq motif finding results.

CONCLUSIONS

The advanced methodology MMGAT provides a robust tool for finding ATAC-seq motifs, and the comprehensive server MMGAT-S makes a significant contribution to genomics research. The open-source code of MMGAT can be found at https://github.com/xiaotianr/MMGAT , and MMGAT-S is freely available at https://www.mmgraphws.com/MMGAT-S/ .

Proteasome Inhibition Reprograms Chromatin Landscape in Breast Cancer

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNA polymerase II (RNAPII) transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed differentially open chromatin regions (DOCR). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic superenhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology.

SIGNIFICANCE

Our study provides a strong basis for understanding the mechanisms by which proteasome inhibitors exert anticancer effects. We find open chromatin regions that change during proteasome inhibition, are typically accessible in non-basal breast cancers.

Widespread impact of nucleosome remodelers on transcription at cis-regulatory elements

Nucleosome remodeling complexes and other regulatory factors work in concert to build a chromatin environment that directs the expression of a distinct set of genes in each cell using cis-regulatory elements (CREs), such as promoters and enhancers, that drive transcription of both mRNAs and CRE-associated non-coding RNAs (ncRNAs). Two classes of CRE-associated ncRNAs include upstream antisense RNAs (uaRNAs), which are transcribed divergently from a shared mRNA promoter, and enhancer RNAs (eRNAs), which are transcribed bidirectionally from active enhancers. The complicated network of CRE regulation by nucleosome remodelers remains only partially explored, with a focus on a select, limited number of remodelers. We endeavored to elucidate a remodeler-based regulatory network governing CRE-associated transcription (mRNA, eRNA, and uaRNA) in murine embryonic stem (ES) cells to test the hypothesis that many SNF2-family nucleosome remodelers collaborate to regulate the coding and non-coding transcriptome via alteration of underlying nucleosome architecture. Using depletion followed by transient transcriptome sequencing (TT-seq), we identified thousands of misregulated mRNAs and CRE-associated ncRNAs across the remodelers examined, identifying novel contributions by understudied remodelers in the regulation of coding and noncoding transcription. Our findings suggest that mRNA and eRNA transcription are coordinately co-regulated, while mRNA and uaRNAs sharing a common promoter are independently regulated. Subsequent mechanistic studies suggest that while remodelers SRCAP and CHD8 modulate transcription through classical mechanisms such as transcription factors and histone variants, a broad set of remodelers including SMARCAL1 indirectly contribute to transcriptional regulation through maintenance of genomic stability and proper Integrator complex localization. This study systematically examines the contribution of SNF2-remodelers to the CRE-associated transcriptome, identifying at least two classes for remodeler action.

Unravelling the genetic basis of Schizophrenia

Neuronal development is a highly regulated mechanism that is central to organismal function in animals. In humans, disruptions to this process can lead to a range of neurodevelopmental phenotypes, including Schizophrenia (SCZ). SCZ has a significant genetic component, whereby an individual with an SCZ affected family member is eight times more likely to develop the disease than someone with no family history of SCZ. By examining a combination of genomic, transcriptomic and epigenomic datasets, large-scale 'omics' studies aim to delineate the relationship between genetic variation and abnormal cellular activity in the SCZ brain. Herein, we provide a brief overview of some of the key omics methods currently being used in SCZ research, including RNA-seq, the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and high-throughput chromosome conformation capture (3C) approaches (e.g., Hi-C), as well as single-cell/nuclei iterations of these methods. We also discuss how these techniques are being employed to further our understanding of the genetic basis of SCZ, and to identify associated molecular pathways, biomarkers, and candidate drug targets.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Mutant FOXO1 controls an oncogenic network via enhancer accessibility

Transcriptional dysregulation is a hallmark of diffuse large B cell lymphoma (DLBCL), as transcriptional regulators are frequently mutated. However, our mechanistic understanding of how normal transcriptional programs are co-opted in DLBCL has been hindered by a lack of methodologies that provide the temporal resolution required to separate direct and indirect effects on transcriptional control. We applied a chemical-genetic approach to engineer the inducible degradation of the transcription factor FOXO1, which is recurrently mutated (mFOXO1) in DLBCL. The combination of rapid degradation of mFOXO1, nascent transcript detection, and assessment of chromatin accessibility allowed us to identify the direct targets of mFOXO1. mFOXO1 was required to maintain accessibility at specific enhancers associated with multiple oncogenes, and mFOXO1 degradation impaired RNA polymerase pause-release at some targets. Wild-type FOXO1 appeared to weakly regulate many of the same targets as mFOXO1 and was able to complement the degradation of mFOXO1 in the context of AKT inhibition.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

ARID1A orchestrates SWI/SNF-mediated sequential binding of transcription factors with ARID1A loss driving pre-memory B cell fate and lymphomagenesis

ARID1A, a subunit of the canonical BAF nucleosome remodeling complex, is commonly mutated in lymphomas. We show that ARID1A orchestrates B cell fate during the germinal center (GC) response, facilitating cooperative and sequential binding of PU.1 and NF-kB at crucial genes for cytokine and CD40 signaling. The absence of ARID1A tilts GC cell fate toward immature IgM+CD80-PD-L2- memory B cells, known for their potential to re-enter new GCs. When combined with BCL2 oncogene, ARID1A haploinsufficiency hastens the progression of aggressive follicular lymphomas (FLs) in mice. Patients with FL with ARID1A-inactivating mutations preferentially display an immature memory B cell-like state with increased transformation risk to aggressive disease. These observations offer mechanistic understanding into the emergence of both indolent and aggressive ARID1A-mutant lymphomas through the formation of immature memory-like clonal precursors. Lastly, we demonstrate that ARID1A mutation induces synthetic lethality to SMARCA2/4 inhibition, paving the way for potential precision therapy for high-risk patients.

SMARCA4 is a haploinsufficient B cell lymphoma tumor suppressor that fine-tunes centrocyte cell fate decisions

SMARCA4 encodes one of two mutually exclusive ATPase subunits in the BRG/BRM associated factor (BAF) complex that is recruited by transcription factors (TFs) to drive chromatin accessibility and transcriptional activation. SMARCA4 is among the most recurrently mutated genes in human cancer, including ∼30% of germinal center (GC)-derived Burkitt lymphomas. In mice, GC-specific Smarca4 haploinsufficiency cooperated with MYC over-expression to drive lymphomagenesis. Furthermore, monoallelic Smarca4 deletion drove GC hyperplasia with centroblast polarization via significantly increased rates of centrocyte recycling to the dark zone. Mechanistically, Smarca4 loss reduced the activity of TFs that are activated in centrocytes to drive GC-exit, including SPI1 (PU.1), IRF family, and NF-κB. Loss of activity for these factors phenocopied aberrant BCL6 activity within murine centrocytes and human Burkitt lymphoma cells. SMARCA4 therefore facilitates chromatin accessibility for TFs that shape centrocyte trajectories, and loss of fine-control of these programs biases toward centroblast cell-fate, GC hyperplasia and lymphoma.

Tissue-specific nonheritable influences drive endometrial immune system variation

Although human twin studies have revealed the combined contribution of heritable and environmental factors in shaping immune system variability in blood, the contribution of these factors to immune system variability in tissues remains unexplored. The human uterus undergoes constant regeneration and is exposed to distinct environmental factors. To assess uterine immune system variation, we performed a system-level analysis of endometrial and peripheral blood immune cells in monozygotic twins. Although most immune cell phenotypes in peripheral blood showed high genetic heritability, more variation was found in endometrial immune cells, indicating a stronger influence by environmental factors. Cytomegalovirus infection was identified to influence peripheral blood immune cell variability but had limited effect on endometrial immune cells. Instead, hormonal contraception shaped the local endometrial milieu and immune cell composition with minor influence on the systemic immune system. These results highlight that the magnitude of human immune system variation and factors influencing it can be tissue specific.

Epigenomic states contribute to coordinated allelic transcriptional bursting in iPSC reprogramming

Two alleles of a gene can be transcribed independently or coordinatedly, which can lead to temporal expression heterogeneity with potentially distinct impacts on cell fate. Here, we profiled genome-wide allelic transcriptional burst kinetics during the reprogramming of MEF to induced pluripotent stem cells. We show that the degree of coordination of allelic bursting differs among genes, and alleles of many reprogramming-related genes burst in a highly coordinated fashion. Notably, we show that the chromatin accessibility of the two alleles of highly coordinated genes is similar, unlike the semi-coordinated or independent genes, suggesting the degree of coordination of allelic bursting is linked to allelic chromatin accessibility. Consistently, we show that many transcription factors have differential binding affinity between alleles of semi-coordinated or independent genes. We show that highly coordinated genes are enriched with chromatin accessibility regulators such as H3K4me3, H3K4me1, H3K36me3, H3K27ac, histone variant H3.3, and BRD4. Finally, we demonstrate that enhancer elements are highly enriched in highly coordinated genes. Our study demonstrates that epigenomic states contribute to coordinated allelic bursting to fine-tune gene expression during induced pluripotent stem cell reprogramming.

MYC activity at enhancers drives prognostic transcriptional programs through an epigenetic switch

The transcription factor MYC is overexpressed in most cancers, where it drives multiple hallmarks of cancer progression. MYC is known to promote oncogenic transcription by binding to active promoters. In addition, MYC has also been shown to invade distal enhancers when expressed at oncogenic levels, but this enhancer binding has been proposed to have low gene-regulatory potential. Here, we demonstrate that MYC directly regulates enhancer activity to promote cancer type-specific gene programs predictive of poor patient prognosis. MYC induces transcription of enhancer RNA through recruitment of RNA polymerase II (RNAPII), rather than regulating RNAPII pause-release, as is the case at promoters. This process is mediated by MYC-induced H3K9 demethylation and acetylation by GCN5, leading to enhancer-specific BRD4 recruitment through its bromodomains, which facilitates RNAPII recruitment. We propose that MYC drives prognostic cancer type-specific gene programs through induction of an enhancer-specific epigenetic switch, which can be targeted by BET and GCN5 inhibitors.

Chromatin Remodeling in Patient-Derived Colorectal Cancer Models

Patient-Derived Organoids (PDO) and Xenografts (PDX) are the current gold standards for patient-derived models of cancer (PDMC). Nevertheless, how patient tumor cells evolve in these models and the impact on drug response remains unclear. Herein, the transcriptomic and chromatin accessibility landscapes of matched colorectal cancer (CRC) PDO, PDX, PDO-derived PDX (PDOX), and original patient tumors (PT) are compared. Two major remodeling axes are discovered. The first axis delineates PDMC from PT, and the second axis distinguishes PDX and PDO. PDOX are more similar to PDX than PDO, indicating the growth environment is a driving force for chromatin adaptation. Transcription factors (TF) that differentially bind to open chromatins between matched PDO and PDOX are identified. Among them, KLF14 and EGR2 footprints are enriched in PDOX relative to matched PDO, and silencing of KLF14 or EGR2 promoted tumor growth. Furthermore, EPHA4, a shared downstream target gene of KLF14 and EGR2, altered tumor sensitivity to MEK inhibitor treatment. Altogether, patient-derived CRC cells undergo both common and distinct chromatin remodeling in PDO and PDX/PDOX, driven largely by their respective microenvironments, which results in differences in growth and drug sensitivity and needs to be taken into consideration when interpreting their ability to predict clinical outcome.

GNNMF: a multi-view graph neural network for ATAC-seq motif finding

BACKGROUND

The Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) utilizes the Transposase Tn5 to probe open chromatic, which simultaneously reveals multiple transcription factor binding sites (TFBSs) compared to traditional technologies. Deep learning (DL) technology, including convolutional neural networks (CNNs), has successfully found motifs from ATAC-seq data. Due to the limitation of the width of convolutional kernels, the existing models only find motifs with fixed lengths. A Graph neural network (GNN) can work on non-Euclidean data, which has the potential to find ATAC-seq motifs with different lengths. However, the existing GNN models ignored the relationships among ATAC-seq sequences, and their parameter settings should be improved.

RESULTS

In this study, we proposed a novel GNN model named GNNMF to find ATAC-seq motifs via GNN and background coexisting probability. Our experiment has been conducted on 200 human datasets and 80 mouse datasets, demonstrated that GNNMF has improved the area of eight metrics radar scores of 4.92% and 6.81% respectively, and found more motifs than did the existing models.

CONCLUSIONS

In this study, we developed a novel model named GNNMF for finding multiple ATAC-seq motifs. GNNMF built a multi-view heterogeneous graph by using ATAC-seq sequences, and utilized background coexisting probability and the iterloss to find different lengths of ATAC-seq motifs and optimize the parameter sets. Compared to existing models, GNNMF achieved the best performance on TFBS prediction and ATAC-seq motif finding, which demonstrates that our improvement is available for ATAC-seq motif finding.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

Epigenetic signature and key transcriptional regulators of human antigen-specific type 1 regulatory T cells

Human adaptive immunity is orchestrated by effector and regulatory T (Treg) cells. Natural Tregs arise in the thymus where they are shaped to recognize self-antigens, while type 1 Tregs or Tr1 cells are induced from conventional peripheral CD4 + T cells in response to peripheral antigens, such as alloantigens and allergens. Tr1 cells have been developed as a potential therapy for inducing antigen-specific tolerance, because they can be rapidly differentiated in vitro in response to a target antigen. However, the epigenetic landscape and the identity of transcription factors (TFs) that regulate differentiation, phenotype, and functions of human antigen-specific Tr1 cells is largely unknown, hindering Tr1 research and broader clinical development. Here, we reveal the unique epigenetic signature of antigen-specific Tr1 cells, and TFs that regulate their differentiation, phenotype and function. We showed that in vitro induced antigen-specific Tr1 cells are distinct both clonally and transcriptionally from natural Tregs and other conventional CD4 + T cells on a single-cell level. An integrative analysis of Tr1 cell epigenome and transcriptome identified a TF signature unique to antigen-specific Tr1 cells, and predicted that IRF4, BATF, and MAF act as their transcriptional regulators. Using functional genomics, we showed that each of these TFs play a non-redundant role in regulating Tr1 cell differentiation, suppressive function, and expression of co-inhibitory and cytotoxic proteins. By using the Tr1-specific TF signature as a molecular fingerprint, we tracked Tr1 cells in peripheral blood of recipients of allogeneic hematopoietic stem cell transplantation treated with adoptive Tr1 cell therapy. Furthermore, the same signature identified Tr1 cells in resident CD4 + T cells in solid tumors. Altogether, these results reveal the epigenetic signature and the key transcriptional regulators of human Tr1 cells. These data will guide mechanistic studies of human Tr1 cell biology and the development and optimization of adoptive Tr1 cell therapies.

A screen for regeneration-associated silencer regulatory elements in zebrafish

Regeneration involves gene expression changes explained in part by context-dependent recruitment of transcriptional activators to distal enhancers. Silencers that engage repressive transcriptional complexes are less studied than enhancers and more technically challenging to validate, but they potentially have profound biological importance for regeneration. Here, we identified candidate silencers through a screening process that examined the ability of DNA sequences to limit injury-induced gene expression in larval zebrafish after fin amputation. A short sequence (s1) on chromosome 5 near several genes that reduce expression during adult fin regeneration could suppress promoter activity in stable transgenic lines and diminish nearby gene expression in knockin lines. High-resolution analysis of chromatin organization identified physical associations of s1 with gene promoters occurring preferentially during fin regeneration, and genomic deletion of s1 elevated the expression of these genes after fin amputation. Our study provides methods to identify "tissue regeneration silencer elements" (TRSEs) with the potential to reduce unnecessary or deleterious gene expression during regeneration.