Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position

Genetic and Epigenetic Characterization of Sarcoma Stem Cells Across Subtypes Identifies EZH2 as a Therapeutic Target

High-grade complex karyotype soft tissue sarcomas (STS) are a heterogeneous and aggressive set of cancers that share a common treatment strategy. Disease progression and failure to respond to anthracycline based chemotherapy, standard first-line treatment, is associated with poor patient outcomes. To address this, we investigated the contribution of STS cancer stem cells (STS-CSCs) to doxorubicin resistance. We identified a positive correlation between CSC abundance and doxorubicin IC 50 in resistant cell lines. We investigated if a common genetic signature across STS-CSCs could be targeted. Utilizing patient derived samples from five sarcoma subtypes we identified Enhancer of Zeste homolog 2 (EZH2), a member of the polycomb repressive complex 2 (PRC2) responsible for H3K27 methylation as being enriched in the CSC population. EZH2 activity and a shared epigenetic profile was observed across subtypes. Targeting of EZH2 using Tazemetostat, an FDA approved inhibitor specifically ablated the STS-CSC population. Treatment of doxorubicin resistant cell lines with tazemetostat resulted in a decrease in the STS-CSC population. Further, co-treatment was not only synergistic in the parent cell lines, but restored chemosensitivity in doxorubicin resistant lines. These data confirm the presence of shared genetic programs across distinct subtypes of CSC-STS that can be therapeutically targeted.

Clinical drug screening reveals clofazimine potentiates the efficacy while reducing the toxicity of anti-PD-1 and CTLA-4 immunotherapy

Emerging as the most potent and durable combinational immunotherapy, dual anti-PD-1 and CTLA-4 immune checkpoint blockade (ICB) therapy notoriously increases grade 3-5 immune-related adverse events (irAEs) in patients. Accordingly, attempts to improve the antitumor potency of anti-PD-1+CTLA-4 ICB by including additional therapeutics have been largely discouraged due to concerns of further increasing fatal toxicity. Here, we screened ∼3,000 Food and Drug Administration (FDA)-approved drugs and identified clofazimine as a potential third agent to optimize anti-PD-1+CTLA-4 ICB. Remarkably, clofazimine outperforms ICB dose reduction or steroid treatment in reversing lethality of irAEs, but unlike the detrimental effect of steroids on antitumor efficacy, clofazimine potentiates curative responses in anti-PD-1+CTLA-4 ICB. Mechanistically, clofazimine promotes E2F1 activation in CD8+ T cells to overcome resistance and counteracts pathogenic Th17 cells to abolish irAEs. Collectively, clofazimine potentiates the antitumor efficacy of anti-PD-1+CTLA-4 ICB, curbs intractable irAEs, and may fill a desperate clinical need to improve patient survival.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Emerging Approaches to Profile Accessible Chromatin from Formalin-Fixed Paraffin-Embedded Sections

Nucleosomes are non-uniformly distributed across eukaryotic genomes, with stretches of 'open' chromatin strongly associated with transcriptionally active promoters and enhancers. Understanding chromatin accessibility patterns in normal tissue and how they are altered in pathologies can provide critical insights to development and disease. With the advent of high-throughput sequencing, a variety of strategies have been devised to identify open regions across the genome, including DNase-seq, MNase-seq, FAIRE-seq, ATAC-seq, and NicE-seq. However, the broad application of such methods to FFPE (formalin-fixed paraffin-embedded) tissues has been curtailed by the major technical challenges imposed by highly fixed and often damaged genomic material. Here, we review the most common approaches for mapping open chromatin regions, recent optimizations to overcome the challenges of working with FFPE tissue, and a brief overview of a typical data pipeline with analysis considerations.

Epigenetic Regulation of Cardiomyocyte Maturation by Arginine Methyltransferase CARM1

BACKGROUND

During the neonatal stage, the cardiomyocyte undergoes a constellation of molecular, cytoarchitectural, and functional changes known collectively as cardiomyocyte maturation to increase myocardial contractility and cardiac output. Despite the importance of cardiomyocyte maturation, the molecular mechanisms governing this critical process remain largely unexplored.

METHODS

We leveraged an in vivo mosaic knockout system to characterize the role of Carm1, the founding member of protein arginine methyltransferase, in cardiomyocyte maturation. Using a battery of assays, including immunohistochemistry, immuno-electron microscopy imaging, and action potential recording, we assessed the effect of loss of Carm1 function on cardiomyocyte cell growth, myofibril expansion, T-tubule formation, and electrophysiological maturation. Genome-wide transcriptome profiling, H3R17me2a chromatin immunoprecipitation followed by sequencing, and assay for transposase-accessible chromatin with high-throughput sequencing were used to investigate the mechanisms by which CARM1 (coactivator-associated arginine methyltransferase 1) regulates cardiomyocyte maturation. Finally, we interrogated the human syntenic region to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks for single-nucleotide polymorphisms associated with human heart diseases.

RESULTS

We report that mosaic ablation of Carm1 disrupts multiple aspects of cardiomyocyte maturation cell autonomously, leading to reduced cardiomyocyte size and sarcomere thickness, severe loss and disorganization of T tubules, and compromised electrophysiological maturation. Genomics study demonstrates that CARM1 directly activates genes that underlie cardiomyocyte cytoarchitectural and electrophysiological maturation. Moreover, our study reveals significant enrichment of human heart disease-associated single-nucleotide polymorphisms in the human genomic region syntenic to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks.

CONCLUSIONS

This study establishes a critical and multifaceted role for CARM1 in regulating cardiomyocyte maturation and demonstrates that deregulation of CARM1-dependent cardiomyocyte maturation gene expression may contribute to human heart diseases.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

BACKGROUND

As genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype.

RESULTS

In this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation .

CONCLUSIONS

GaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

BACKGROUND

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

RESULTS

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

CONCLUSIONS

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

CD-Loop: a chromatin loop detection method based on the diffusion model

Motivation

In recent years, there have been significant advances in various chromatin conformation capture techniques, and annotating the topological structure from Hi-C contact maps has become crucial for studying the three-dimensional structure of chromosomes. However, the structure and function of chromatin loops are highly dynamic and diverse, influenced by multiple factors. Therefore, obtaining the three-dimensional structure of the genome remains a challenging task. Among many chromatin loop prediction methods, it is difficult to fully extract features from the contact map and make accurate predictions at low sequencing depths.

Results

In this study, we put forward a deep learning framework based on the diffusion model called CD-Loop for predicting accurate chromatin loops. First, by pre-training the input data, we obtain prior probabilities for predicting the classification of the Hi-C contact map. Then, by combining the denoising process based on the diffusion model and the prior probability obtained by pre-training, candidate loops were predicted from the input Hi-C contact map. Finally, CD-Loop uses a density-based clustering algorithm to cluster the candidate chromatin loops and predict the final chromatin loops. We compared CD-Loop with the currently popular methods, such as Peakachu, Chromosight, and Mustache, and found that in different cell types, species, and sequencing depths, CD-Loop outperforms other methods in loop annotation. We conclude that CD-Loop can accurately predict chromatin loops and reveal cell-type specificity. The code is available at https://github.com/wangyang199897/CD-Loop.

VGLL2 and TEAD1 fusion proteins drive YAP/TAZ-independent transcription and tumorigenesis by engaging p300

Studies on Hippo pathway regulation of tumorigenesis largely center on YAP and TAZ, the transcriptional co-regulators of TEAD. Here, we present an oncogenic mechanism involving VGLL and TEAD fusions that is Hippo pathway-related but YAP/TAZ-independent. We characterize two recurrent fusions, VGLL2-NCOA2 and TEAD1-NCOA2, recently identified in spindle cell rhabdomyosarcoma. We demonstrate that, in contrast to VGLL2 and TEAD1, the fusion proteins are strong activators of TEAD-dependent transcription, and their function does not require YAP/TAZ. Furthermore, we identify that VGLL2 and TEAD1 fusions engage specific epigenetic regulation by recruiting histone acetyltransferase p300 to control TEAD-mediated transcriptional and epigenetic landscapes. We showed that small molecule p300 inhibition can suppress fusion proteins-induced oncogenic transformation both in vitro and in vivo. Overall, our study reveals a molecular basis for VGLL involvement in cancer and provides a framework for targeting tumors carrying VGLL, TEAD, or NCOA translocations.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Nuclear RNA homeostasis promotes systems-level coordination of cell fate and senescence

Understanding cellular coordination remains a challenge despite knowledge of individual pathways. The RNA exosome, targeting a wide range of RNA substrates, is often downregulated in cellular senescence. Utilizing an auxin-inducible system, we observed that RNA exosome depletion in embryonic stem cells significantly affects the transcriptome and proteome, causing pluripotency loss and pre-senescence onset. Mechanistically, exosome depletion triggers acute nuclear RNA aggregation, disrupting nuclear RNA-protein equilibrium. This disturbance limits nuclear protein availability and hinders polymerase initiation and engagement, reducing gene transcription. Concurrently, it promptly disrupts nucleolar transcription, ribosomal processes, and nuclear exporting, resulting in a translational shutdown. Prolonged exosome depletion induces nuclear structural changes resembling senescent cells, including aberrant chromatin compaction, chromocenter disassembly, and intensified heterochromatic foci. These effects suggest that the dynamic turnover of nuclear RNA orchestrates crosstalk between essential processes to optimize cellular function. Disruptions in nuclear RNA homeostasis result in systemic functional decline, altering the cell state and promoting senescence.

A data-adaptive methods in detecting exogenous methyltransferase accessible chromatin in human genome using nanopore sequencing

MOTIVATION

Identifying chromatin accessibility is one of the key steps in studying the regulation of eukaryotic genomes. The combination of exogenous methyltransferase and nanopore sequencing provides an strategy to identify open chromatin over long genomic ranges at the single-molecule scale. However, endogenous methylation, non-open-chromatin-specific exogenous methylation and base-calling errors limit the accuracy and hinders its application to complex genomes.

RESULTS

We systematically evaluated the impact of these three influence factors, and developed a model-based computational method, methyltransferase accessible genome region finder (MAGNIFIER), to address the issues. By incorporating control data, MAGNIFIER attenuates the three influence factors with data-adaptive comparison strategy. We demonstrate that MAGNIFIER is not only sensitive to identify the open chromatin with much improved accuracy, but also able to detect the chromatin accessibility of repetitive regions that are missed by NGS-based methods. By incorporating long-read RNA-seq data, we revealed the association between the accessible Alu elements and non-classic gene isoforms.

AVAILABILITY AND IMPLEMENTATION

Freely available on web at https://github.com/Goatofmountain/MAGNIFIER.

LATS1 controls CTCF chromatin occupancy and hormonal response of 3D-grown breast cancer cells

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

Integration Analysis of Single-Cell Multi-Omics Reveals Prostate Cancer Heterogeneity

Prostate cancer (PCa) is an extensive heterogeneous disease with a complex cellular ecosystem in the tumor microenvironment (TME). However, the manner in which heterogeneity is shaped by tumors and stromal cells, or vice versa, remains poorly understood. In this study, single-cell RNA sequencing, spatial transcriptomics, and bulk ATAC-sequence are integrated from a series of patients with PCa and healthy controls. A stemness subset of club cells marked with SOX9highARlow expression is identified, which is markedly enriched after neoadjuvant androgen-deprivation therapy (ADT). Furthermore, a subset of CD8+CXCR6+ T cells that function as effector T cells is markedly reduced in patients with malignant PCa. For spatial transcriptome analysis, machine learning and computational intelligence are comprehensively utilized to identify the cellular diversity of prostate cancer cells and cell-cell communication in situ. Macrophage and neutrophil state transitions along the trajectory of cancer progression are also examined. Finally, the immunosuppressive microenvironment in advanced PCa is found to be associated with the infiltration of regulatory T cells (Tregs), potentially induced by an FAP+ fibroblast subset. In summary, the cellular heterogeneity is delineated in the stage-specific PCa microenvironment at single-cell resolution, uncovering their reciprocal crosstalk with disease progression, which can be helpful in promoting PCa diagnosis and therapy.

Blimp-1 and c-Maf regulate immune gene networks to protect against distinct pathways of pathobiont-induced colitis

Intestinal immune responses to microbes are controlled by the cytokine IL-10 to avoid immune pathology. Here, we use single-cell RNA sequencing of colon lamina propria leukocytes (LPLs) along with RNA-seq and ATAC-seq of purified CD4+ T cells to show that the transcription factors Blimp-1 (encoded by Prdm1) and c-Maf co-dominantly regulate Il10 while negatively regulating proinflammatory cytokines in effector T cells. Double-deficient Prdm1fl/flMaffl/flCd4Cre mice infected with Helicobacter hepaticus developed severe colitis with an increase in TH1/NK/ILC1 effector genes in LPLs, while Prdm1fl/flCd4Cre and Maffl/flCd4Cre mice exhibited moderate pathology and a less-marked type 1 effector response. LPLs from infected Maffl/flCd4Cre mice had increased type 17 responses with increased Il17a and Il22 expression and an increase in granulocytes and myeloid cell numbers, resulting in increased T cell-myeloid-neutrophil interactions. Genes over-expressed in human inflammatory bowel disease showed differential expression in LPLs from infected mice in the absence of Prdm1 or Maf, revealing potential mechanisms of human disease.

YY1 binding is a gene-intrinsic barrier to Xist-mediated gene silencing

X chromosome inactivation (XCI) in mammals is mediated by Xist RNA which functions in cis to silence genes on a single X chromosome in XX female cells, thereby equalising levels of X-linked gene expression relative to XY males. XCI progresses over a period of several days, with some X-linked genes silencing faster than others. The chromosomal location of a gene is an important determinant of silencing rate, but uncharacterised gene-intrinsic features also mediate resistance or susceptibility to silencing. In this study, we examine mouse embryonic stem cell lines with an inducible Xist allele (iXist-ChrX mESCs) and integrate allele-specific data of gene silencing and decreasing inactive X (Xi) chromatin accessibility over time courses of Xist induction with cellular differentiation. Our analysis reveals that motifs bound by the transcription factor YY1 are associated with persistently accessible regulatory elements, including many promoters and enhancers of slow-silencing genes. We further show that YY1 is evicted relatively slowly from target sites on Xi, and that silencing of X-linked genes is increased upon YY1 degradation. Together our results suggest that YY1 acts as a barrier to Xist-mediated silencing until the late stages of the XCI process.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Epigenetic control of microglial immune responses

Microglia, the major population of brain-resident macrophages, are now recognized as a heterogeneous population comprising several cell subtypes with different (so far mostly supposed) functions in health and disease. A number of studies have performed molecular characterization of these different microglial activation states over the last years making use of "omics" technologies, that is transcriptomics, proteomics and, less frequently, epigenomics profiling. These approaches offer the possibility to identify disease mechanisms, discover novel diagnostic biomarkers, and develop new therapeutic strategies. Here, we focus on epigenetic profiling as a means to understand microglial immune responses beyond what other omics methods can offer, that is, revealing past and present molecular responses, gene regulatory networks and potential future response trajectories, and defining cell subtype-specific disease relevance through mapping non-coding genetic variants. We review the current knowledge in the field regarding epigenetic regulation of microglial identity and function, provide an exemplary analysis that demonstrates the advantages of performing joint transcriptomic and epigenomic profiling of single microglial cells and discuss how comprehensive epigenetic analyses may enhance our understanding of microglial pathophysiology.

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.

Integrative ATAC-seq and RNA-seq analysis of myogenic differentiation of ovine skeletal muscle satellite cell

Skeletal muscle satellite cells (SMSCs) play an important role in regulating muscle growth and regeneration. Chromatin accessibility allows physical interactions that synergistically regulate gene expression through enhancers, promoters, insulators, and chromatin binding factors. However, the chromatin accessibility altas and its regulatory role in ovine myoblast differentiation is still unclear. Therefore, ATAC-seq and RNA-seq analysis were performed on ovine SMSCs at the proliferation stage (SCG) and differentiation stage (SCD). 17,460 DARs (differential accessibility regions) and 3732 DEGs (differentially expressed genes) were identified. Based on joint analysis of ATAC-seq and RNA-seq, we revealed that PI3K-Akt, TGF-β and other signaling pathways regulated SMSCs differentiation. We identified two novel candidate genes, FZD5 and MAP2K6, which may affect the proliferation and differentiation of SMSCs. Our data identify potential cis regulatory elements of ovine SMSCs. This study can provide a reference for exploring the mechanisms of the differentiation and regeneration of SMSCs in the future.

3D epigenomics and 3D epigenopathies

Mammalian genomes are intricately compacted to form sophisticated 3-dimensional structures within the tiny nucleus, so called 3D genome folding. Despite their shapes reminiscent of an entangled yarn, the rapid development of molecular and next-generation sequencing technologies (NGS) has revealed that mammalian genomes are highly organized in a hierarchical order that delicately affects transcription activities. An increasing amount of evidence suggests that 3D genome folding is implicated in diseases, giving us a clue on how to identify novel therapeutic approaches. In this review, we will study what 3D genome folding means in epigenetics, what types of 3D genome structures there are, how they are formed, and how the technologies have developed to explore them. We will also discuss the pathological implications of 3D genome folding. Finally, we will discuss how to leverage 3D genome folding and engineering for future studies. [BMB Reports 2024; 57(5): 216-231].

The accessible chromatin landscape of lipopolysaccharide-induced systemic inflammatory response identifying epigenome signatures and transcription regulatory networks in chickens

Lipopolysaccharide (LPS) can induce systemic inflammatory response (SIR) in animals. Understanding the regulatory mechanism of SIR and therapies to ensure healthy growth is urgently needed. Chromatin remodeling plays a crucial role in the expression of genes involved in immune diseases. In the present study, the ATAC-seq analysis revealed 3491 differential open chromatin sites in the spleen of chicks with SIR induced by LPS challenge, and we presented the motifs on these sites and the associated transcription factors. The regulatory network was presented by combining the differential open chromatin data with the mRNAs and exploded cytokines. Interestingly, the LPS challenge could regulate the mRNA expression of 202 genes through chromatin reprogramming, including critical genes such as TLE1 and JUN, which regulate signaling pathways such as I-κB kinase/NF-κB, Toll-like receptor, and downstream cytokine genes. Furthermore, dietary daidzein could inhibit DNA topoisomerase II, which reprograms the spatial conformation of chromatin in the inflammatory response and attenuates SIR. In conclusion, we successfully identified key genes directly regulated by chromatin reprogramming in SIR and demonstrated the chromatin epigenome signatures and transcriptional regulatory network, which provides an important reference for further research on avian epigenetics. There is great potential for alleviating SIR using dietary daidzein.

Clock-dependent chromatin accessibility rhythms regulate circadian transcription

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

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Reduction of ETV1 is Identified as a Prominent Feature of Age-Related Cataract

PURPOSE

To identify the inactive genes in cataract lenses and explore their function in lens epithelial cells (LECs).

METHODS

Lens epithelium samples obtained from both age-related cataract (ARC) patients and normal donors were subjected to two forms of histone H3 immunoprecipitation: H3K9ac and H3K27me3 chromatin immunoprecipitation (ChIP), followed by ChIP-seq. The intersection set of "active genes in normal controls" and "repressed genes in cataract lenses" was identified. To validate the role of a specific gene, ETV1, within this set, quantitative polymerase chain reaction (qPCR), western blot, and immunofluorescence were performed using clinical lens epithelium samples. Small interference RNA (siRNA) was utilized to reduce the mRNA level of ETV1 in cultured LECs. Following this, transwell assay and western blot was performed to examine the migration ability of the cells. Furthermore, RNA-seq analysis was conducted on both cell samples with ETV1 knockdown and control cells. Additionally, the expression level of ETV1 in LECs was examined using qPCR under H2O2 treatment.

RESULTS

Six genes were identified in the intersection set of "active genes in normal controls" and "repressed genes in ARC lenses". Among these genes, ETV1 showed the most significant fold-change decrease in the cataract samples compared to the control samples. After ETV1 knockdown by siRNA in cultured LECs, reduced cell migration was observed, along with a decrease in the expression of β-Catenin and Vimentin, two specific genes associated with cell migration. In addition, under the oxidative stress induced by H2O2 treatment, the expression level of ETV1 in LECs significantly decreased.

CONCLUSIONS

Based on the findings of this study, it can be concluded that ETV1 is significantly reduced in human ARC lenses. The repression of ETV1 in ARC lenses appears to contribute to the disrupted differentiation of lens epithelium, which is likely caused by the inhibition of both cell differentiation and migration processes.

Multi-omics analysis-based macrophage differentiation-associated papillary thyroid cancer patient classifier

BACKGROUND

The reclassification of Papillary Thyroid Carcinoma (PTC) is an area of research that warrants attention. The connection between thyroid cancer, inflammation, and immune responses necessitates considering the mechanisms of differential prognosis of thyroid tumors from an immunological perspective. Given the high adaptability of macrophages to environmental stimuli, focusing on the differentiation characteristics of macrophages might offer a novel approach to address the issues related to PTC subtyping.

METHODS

Single-cell RNA sequencing data of medullary cells infiltrated by papillary thyroid carcinoma obtained from public databases was subjected to dimensionality reduction clustering analysis. The RunUMAP and FindAllMarkers functions were utilized to identify the gene expression matrix of different clusters. Cell differentiation trajectory analysis was conducted using the Monocle R package. A complex regulatory network for the classification of Immune status and Macrophage differentiation-associated Papillary Thyroid Cancer Classification (IMPTCC) was constructed through quantitative multi-omics analysis. Immunohistochemistry (IHC) staining was utilized for pathological histology validation.

RESULTS

Through the integration of single-cell RNA and bulk sequencing data combined with multi-omics analysis, we identified crucial transcription factors, immune cells/immune functions, and signaling pathways. Based on this, regulatory networks for three IMPTCC clusters were established.

CONCLUSION

Based on the co-expression network analysis results, we identified three subtypes of IMPTCC: Immune-Suppressive Macrophage differentiation-associated Papillary Thyroid Carcinoma Classification (ISMPTCC), Immune-Neutral Macrophage differentiation-associated Papillary Thyroid Carcinoma Classification (INMPTCC), and Immune-Activated Macrophage differentiation-associated Papillary Thyroid Carcinoma Classification (IAMPTCC). Each subtype exhibits distinct metabolic, immune, and regulatory characteristics corresponding to different states of macrophage differentiation.

PU.1 and BCL11B sequentially cooperate with RUNX1 to anchor mSWI/SNF to poise the T cell effector landscape

Adaptive immunity relies on specialized effector functions elicited by lymphocytes, yet how antigen recognition activates appropriate effector responses through nonspecific signaling intermediates is unclear. Here we examined the role of chromatin priming in specifying the functional outputs of effector T cells and found that most of the cis-regulatory landscape active in effector T cells was poised early in development before the expression of the T cell antigen receptor. We identified two principal mechanisms underpinning this poised landscape: the recruitment of the nucleosome remodeler mammalian SWItch/Sucrose Non-Fermentable (mSWI/SNF) by the transcription factors RUNX1 and PU.1 to establish chromatin accessibility at T effector loci; and a 'relay' whereby the transcription factor BCL11B succeeded PU.1 to maintain occupancy of the chromatin remodeling complex mSWI/SNF together with RUNX1, after PU.1 silencing during lineage commitment. These mechanisms define modes by which T cells acquire the potential to elicit specialized effector functions early in their ontogeny and underscore the importance of integrating extrinsic cues to the developmentally specified intrinsic program.

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.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Determinant of m6A regional preference by transcriptional dynamics

N6-Methyladenosine (m6A) is the most abundant chemical modification occurring on eukaryotic mRNAs, and has been reported to be involved in almost all stages of mRNA metabolism. The distribution of m6A sites is notably asymmetric along mRNAs, with a strong preference toward the 3' terminus of the transcript. How m6A regional preference is shaped remains incompletely understood. In this study, by performing m6A-seq on chromatin-associated RNAs, we found that m6A regional preference arises during transcription. Nucleosome occupancy is remarkedly increased in the region downstream of m6A sites, suggesting an intricate interplay between m6A methylation and nucleosome-mediated transcriptional dynamics. Notably, we found a remarkable slowdown of Pol-II movement around m6A sites. In addition, inhibiting Pol-II movement increases nearby m6A methylation levels. By analyzing massively parallel assays for m6A, we found that RNA secondary structures inhibit m6A methylation. Remarkably, the m6A sites associated with Pol-II pausing tend to be embedded within RNA secondary structures. These results suggest that Pol-II pausing could affect the accessibility of m6A motifs to the methyltransferase complex and subsequent m6A methylation by mediating RNA secondary structure. Overall, our study reveals a crucial role of transcriptional dynamics in the formation of m6A regional preference.

Class I HDAC inhibitors enhance antitumor efficacy and persistence of CAR-T cells by activation of the Wnt pathway

Epigenetic modification shapes differentiation trajectory and regulates the exhaustion state of chimeric antigen receptor T (CAR-T) cells. Limited efficacy induced by terminal exhaustion closely ties with intrinsic transcriptional regulation. However, the comprehensive regulatory mechanisms remain largely elusive. Here, we identify class I histone deacetylase inhibitors (HDACi) as boosters of CAR-T cell function by high-throughput screening of chromatin-modifying drugs, in which M344 and chidamide enhance memory maintenance and resistance to exhaustion of CAR-T cells that induce sustained antitumor efficacy both in vitro and in vivo. Mechanistically, HDACi decrease HDAC1 expression and enhance H3K27ac activity. Multi-omics analyses from RNA-seq, ATAC-seq, and H3K27ac CUT&Tag-seq show that HDACi upregulate expression of TCF4, LEF1, and CTNNB1, which subsequently activate the canonical Wnt/β-catenin pathway. Collectively, our findings elucidate the functional roles of class I HDACi in enhancing CAR-T cell function, which provides the basis and therapeutic targets for synergic combination of CAR-T cell therapy and HDACi treatment.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

PeNAC67-PeKAN2-PeSCL23 and B-class MADS-box transcription factors synergistically regulate the specialization process from petal to lip in Phalaenopsis equestris

The molecular basis of orchid flower development involves a specific regulatory program in which MADS-box transcription factors play a central role. The recent 'perianth code' model hypothesizes that two types of higher-order heterotetrameric complexes, namely SP complex and L complex, play pivotal roles in the orchid perianth organ formation. Therefore, we explored their roles and searched for other components of the regulatory network.Through the combined analysis for transposase-accessible chromatin with high-throughput sequencing and RNA sequencing of the lip-like petal and lip from Phalaenopsis equestris var.trilip, transcription factor-(TF) genes involved in lip development were revealed. PeNAC67 encoding a NAC-type TF and PeSCL23 encoding a GRAS-type TF were differentially expressed between the lip-like petal and the lip. PeNAC67 interacted with and stabilized PeMADS3, which positively regulated the development of lip-like petal to lip. PeSCL23 and PeNAC67 competitively bound with PeKAN2 and positively regulated the development of lip-like petal to petal by affecting the level of PeMADS3. PeKAN2 as an important TF that interacts with PeMADS3 and PeMADS9 can promote lip development. These results extend the 'perianth code' model and shed light on the complex regulation of orchid flower development.

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.

Signaling Transduction Pathways and G-Protein-Coupled Receptors in Different Stages of the Embryonic Diapause Termination Process in Artemia

Artemia is a widely distributed small aquatic crustacean, renowned for its ability to enter a state of embryonic diapause. The embryonic diapause termination (EDT) is closely linked to environmental cues, but the precise underlying mechanisms remain elusive. In this study, ATAC-seq and RNA-seq sequencing techniques were employed to explore the gene expression profiles in Artemia cysts 30 min after EDT. These profiles were compared with those during diapause and 5 h after EDT. The regulatory mechanisms governing the EDT process were analyzed through Gene Ontology (GO) enrichment analysis of differentially expressed genes. Furthermore, the active G-protein-coupled receptors (GPCRs) were identified through structural analysis. The results unveiled that the signaling transduction during EDT primarily hinges on GPCRs and the cell surface receptor signaling pathway, but distinct genes are involved across different stages. Hormone-mediated signaling pathways and the tachykinin receptor signaling pathway exhibited heightened activity in the '0-30 min' group, whereas the Wnt signaling pathway manifested its function solely in the '30 min-5 h' group. These results imply a complete divergence in the mechanisms of signal regulation during these two stages. Moreover, through structural analysis, five GPCRs operating at different stages of EDT were identified. These findings provide valuable insights into the signal regulation mechanisms governing Artemia diapause.

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.

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.

Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data

Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N = 29635, FDR < 0.01 and fold change > 50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

scifi-ATAC-seq: massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called scifi-ATAC-seq, single-cell combinatorial fluidic indexing ATAC-sequencing, which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using the 10X Genomics platform. With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples can be indexed in a single emulsion reaction, representing an approximately 20-fold increase in throughput compared to the standard 10X Genomics workflow.

The Role of Immune Cells in DKD: Mechanisms and Targeted Therapies

Diabetic kidney disease (DKD), is a common microvascular complication and a major cause of death in patients with diabetes. Disorders of immune cells and immune cytokines can accelerate DKD development of in a number of ways. As the kidney is composed of complex and highly differentiated cells, the interactions among different cell types and immune cells play important regulatory roles in disease development. Here, we summarize the latest research into the molecular mechanisms underlying the interactions among various immune and renal cells in DKD. In addition, we discuss the most recent studies related to single cell technology and bioinformatics analysis in the field of DKD. The aims of our review were to explore immune cells as potential therapeutic targets in DKD and provide some guidance for future clinical treatments.

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.

Fibroblast regulates angiogenesis in assembled oral cancer organoid: A possible role of NNMT

OBJECTIVE

Tumour angiogenesis is affected by various cell types in the tumour microenvironment (TME), including cancer cells and cancer-associated fibroblasts (CAFs). Here, an assembled organoid model was generated to investigate the mechanism by which the TME regulates angiogenesis in oral squamous cell carcinoma (OSCC).

MATERIALS AND METHODS

Secretion of vascular endothelial growth factor-A (VEGFA) was analysed to compare the proangiogenic properties of OSCC cells and corresponding CAFs. Cell aggregates consisting of endothelial cells (ECs), CAFs and cancer cells were generated to construct assembled organoids. Nicotinamide N-methyltransferase (NNMT) was pharmacologically or genetically inhibited to block the activation of CAFs. ATAC-seq was employed to test the transcriptional network of fibroblasts overexpressing NNMT.

RESULTS

Compared with cancer cells, CAFs secreted more VEGFA. Coculture with CAFs more effectively promoted the sprouting of ECs. Blockade of CAF activation via inhibition of NNMT drastically reduced the expression of CD31 in the assembled organoids. Overexpression of NNMT enhanced the transcription of genes related to angiogenesis in fibroblasts. Specifically, NNMT orchestrated the enrichment of the transcription factor JUNB at the promoter of VEGFA.

CONCLUSIONS

We clarify that stromal NNMT enables the steady reproduction of angiogenesis in assembled oral cancer organoids, providing a novel target for exploiting antiangiogenic therapy.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

DNA methylation and chromatin accessibility predict age in the domestic dog

Across mammals, the epigenome is highly predictive of chronological age. These "epigenetic clocks," most of which have been built using DNA methylation (DNAm) profiles, have gained traction as biomarkers of aging and organismal health. While the ability of DNAm to predict chronological age has been repeatedly demonstrated, the ability of other epigenetic features to predict age remains unclear. Here, we use two types of epigenetic information-DNAm, and chromatin accessibility as measured by ATAC-seq-to develop age predictors in peripheral blood mononuclear cells sampled from a population of domesticated dogs. We measured DNAm and ATAC-seq profiles for 71 dogs, building separate predictive clocks from each, as well as the combined dataset. We also use fluorescence-assisted cell sorting to quantify major lymphoid populations for each sample. We found that chromatin accessibility can accurately predict chronological age (R2 ATAC = 26%), though less accurately than the DNAm clock (R2 DNAm = 33%), and the clock built from the combined datasets was comparable to both (R2 combined = 29%). We also observed various populations of CD62L+ T cells significantly correlated with dog age. Finally, we found that all three clocks selected features that were in or near at least two protein-coding genes: BAIAP2 and SCARF2, both previously implicated in processes related to cognitive or neurological impairment. Taken together, these results highlight the potential of chromatin accessibility as a complementary epigenetic resource for modeling and investigating biologic age.

Targeting the chromatin structural changes of antitumor immunity

Epigenomic imbalance drives abnormal transcriptional processes, promoting the onset and progression of cancer. Although defective gene regulation generally affects carcinogenesis and tumor suppression networks, tumor immunogenicity and immune cells involved in antitumor responses may also be affected by epigenomic changes, which may have significant implications for the development and application of epigenetic therapy, cancer immunotherapy, and their combinations. Herein, we focus on the impact of epigenetic regulation on tumor immune cell function and the role of key abnormal epigenetic processes, DNA methylation, histone post-translational modification, and chromatin structure in tumor immunogenicity, and introduce these epigenetic research methods. We emphasize the value of small-molecule inhibitors of epigenetic modulators in enhancing antitumor immune responses and discuss the challenges of developing treatment plans that combine epigenetic therapy and immunotherapy through the complex interaction between cancer epigenetics and cancer immunology.

Applications of single‑cell omics and spatial transcriptomics technologies in gastric cancer (Review)

Gastric cancer (GC) is a prominent contributor to global cancer-related mortalities, and a deeper understanding of its molecular characteristics and tumor heterogeneity is required. Single-cell omics and spatial transcriptomics (ST) technologies have revolutionized cancer research by enabling the exploration of cellular heterogeneity and molecular landscapes at the single-cell level. In the present review, an overview of the advancements in single-cell omics and ST technologies and their applications in GC research is provided. Firstly, multiple single-cell omics and ST methods are discussed, highlighting their ability to offer unique insights into gene expression, genetic alterations, epigenomic modifications, protein expression patterns and cellular location in tissues. Furthermore, a summary is provided of key findings from previous research on single-cell omics and ST methods used in GC, which have provided valuable insights into genetic alterations, tumor diagnosis and prognosis, tumor microenvironment analysis, and treatment response. In summary, the application of single-cell omics and ST technologies has revealed the levels of cellular heterogeneity and the molecular characteristics of GC, and holds promise for improving diagnostics, personalized treatments and patient outcomes in GC.

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.