Chromatin accessibility: a window into the genome

How chromatin senses plant hormones

Plant hormones activate receptors, initiating intracellular signaling pathways. Eventually, hormone-specific transcription factors become active in the nucleus, facilitating hormone-induced transcriptional regulation. Chromatin plays a fundamental role in the regulation of transcription, the process by which genetic information encoded in DNA is converted into RNA. The structure of chromatin, a complex of DNA and proteins, directly influences the accessibility of genes to the transcriptional machinery. The different signaling pathways and transcription factors involved in the transmission of information from the receptors to the nucleus have been readily explored, but not so much for the specific mechanisms employed by the cell to ultimately instruct the chromatin changes necessary for a fast and robust transcription activation, specifically for plant hormone responses. In this review, we will focus on the advancements in understanding how chromatin receives plant hormones, facilitating the changes necessary for fast, robust, and specific transcriptional regulation.

The importance of functional genomics studies in precision rheumatology

Rheumatic diseases, those that affect the musculoskeletal system, cause significant morbidity. Among risk factors of these diseases is a significant genetic component. Recent advances in high-throughput omics techniques now allow a comprehensive profiling of patients at a genetic level through genome-wide association studies. Without functional interpretation of variants identified through these studies, clinical insight remains limited. Strategies include statistical fine-mapping that refine the list of variants in loci associated with disease, whilst colocalization techniques attempt to attribute function to variants that overlap a genetically active chromatin annotation. Functional validation using genome editing techniques can be used to further refine genetic signals and identify key pathways in cell types relevant to rheumatic disease biology. Insight gained from the combination of genetic studies and functional validation can be used to improve precision medicine in rheumatic diseases by allowing risk prediction and drug repositioning.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

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

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

Century-old chromatin architecture revealed in formalin-fixed vertebrates

Gene expression is regulated by changes in chromatin architecture intrinsic to cellular differentiation and as an active response to environmental stimuli. Chromatin dynamics are a major driver of phenotypic diversity, regulation of development, and manifestation of disease. Remarkably, we know little about the evolutionary dynamics of chromatin reorganisation through time, data essential to characterise the impact of environmental stress during the ongoing biodiversity extinction crisis (20th-21st century). Linking the disparate fields of chromatin biology and museum science through their common use of the preservative formaldehyde (a constituent of formalin), we have generated historical chromatin profiles in museum specimens up to 117 years old. Historical chromatin profiles are reproducible, tissue-specific, sex-specific, and environmental condition-dependent in vertebrate specimens. Additionally, we show that over-fixation modulates differential chromatin accessibility to enable semi-quantitative estimates of relative gene expression in vertebrates and a yeast model. Our approach transforms formalin-fixed biological collections into an accurate, comprehensive, and global record of environmental impact on gene expression and phenotype.

Assessing next-generation sequencing-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators (TRs) with query gene sets. Identification of TRs is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Transcriptomic Approaches to Cardiomyocyte-Biomaterial Interactions: A Review

Biomaterials, essential for supporting, enhancing, and repairing damaged tissues, play a critical role in various medical applications. This Review focuses on the interaction of biomaterials and cardiomyocytes, emphasizing the unique significance of transcriptomic approaches in understanding their interactions, which are pivotal in cardiac bioengineering and regenerative medicine. Transcriptomic approaches serve as powerful tools to investigate how cardiomyocytes respond to biomaterials, shedding light on the gene expression patterns, regulatory pathways, and cellular processes involved in these interactions. Emerging technologies such as bulk RNA-seq, single-cell RNA-seq, single-nucleus RNA-seq, and spatial transcriptomics offer promising avenues for more precise and in-depth investigations. Longitudinal studies, pathway analyses, and machine learning techniques further improve the ability to explore the complex regulatory mechanisms involved. This review also discusses the challenges and opportunities of utilizing transcriptomic techniques in cardiomyocyte-biomaterial research. Although there are ongoing challenges such as costs, cell size limitation, sample differences, and complex analytical process, there exist exciting prospects in comprehensive gene expression analyses, biomaterial design, cardiac disease treatment, and drug testing. These multimodal methodologies have the capacity to deepen our understanding of the intricate interaction network between cardiomyocytes and biomaterials, potentially revolutionizing cardiac research with the aim of promoting heart health, and they are also promising for studying interactions between biomaterials and other cell types.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Spatially resolved epigenome sequencing via Tn5 transposition and deterministic DNA barcoding in tissue

Spatial epigenetic mapping of tissues enables the study of gene regulation programs and cellular functions with the dependency on their local tissue environment. Here we outline a complete procedure for two spatial epigenomic profiling methods: spatially resolved genome-wide profiling of histone modifications using in situ cleavage under targets and tagmentation (CUT&Tag) chemistry (spatial-CUT&Tag) and transposase-accessible chromatin sequencing (spatial-ATAC-sequencing) for chromatin accessibility. Both assays utilize in-tissue Tn5 transposition to recognize genomic DNA loci followed by microfluidic deterministic barcoding to incorporate spatial address codes. Furthermore, these two methods do not necessitate prior knowledge of the transcription or epigenetic markers for a given tissue or cell type but permit genome-wide unbiased profiling pixel-by-pixel at the 10 μm pixel size level and single-base resolution. To support the widespread adaptation of these methods, details are provided in five general steps: (1) sample preparation; (2) Tn5 transposition in spatial-ATAC-sequencing or antibody-controlled pA-Tn5 tagmentation in CUT&Tag; (3) library preparation; (4) next-generation sequencing; and (5) data analysis using our customed pipelines available at: https://github.com/dyxmvp/Spatial_ATAC-seq and https://github.com/dyxmvp/spatial-CUT-Tag . The whole procedure can be completed on four samples in 2-3 days. Familiarity with basic molecular biology and bioinformatics skills with access to a high-performance computing environment are required. A rudimentary understanding of pathology and specimen sectioning, as well as deterministic barcoding in tissue-specific skills (e.g., design of a multiparameter barcode panel and creation of microfluidic devices), are also advantageous. In this protocol, we mainly focus on spatial profiling of tissue region-specific epigenetic landscapes in mouse embryos and mouse brains using spatial-ATAC-sequencing and spatial-CUT&Tag, but these methods can be used for other species with no need for species-specific probe design.

NFIC mediates m6A mRNA methylation to orchestrate transcriptional and post-transcriptional regulation to represses malignant phenotype of non-small cell lung cancer cells

BACKGROUND

Multiple genetic and epigenetic regulatory mechanisms are crucial in the development and tumorigenesis process. Transcriptional regulation often involves intricate relationships and networks with post-transcriptional regulatory molecules, impacting the spatial and temporal expression of genes. However, the synergistic relationship between transcription factors and N6-methyladenosine (m6A) modification in regulating gene expression, as well as their influence on the mechanisms underlying the occurrence and progression of non-small cell lung cancer (NSCLC), requires further investigation. The present study aimed to investigate the synergistic relationship between transcription factors and m6A modification on NSCLC.

METHODS

The transcription factor NFIC and its potential genes was screened by analyzing publicly available datasets (ATAC-seq, DNase-seq, and RNA-seq). The association of NFIC and its potential target genes were validated through ChIP-qPCR and dual-luciferase reporter assays. Additionally, the roles of NFIC and its potential genes in NSCLC were detected in vitro and in vivo through silencing and overexpression assays.

RESULTS

Based on multi-omics data, the transcription factor NFIC was identified as a potential tumor suppressor of NSCLC. NFIC was significantly downregulated in both NSCLC tissues and cells, and when NFIC was overexpressed, the malignant phenotype and total m6A content of NSCLC cells was suppressed, while the PI3K/AKT pathway was inactivated. Additionally, we discovered that NFIC inhibits the expression of METTL3 by directly binding to its promoter region, and METTL3 regulates the expression of KAT2A, a histone acetyltransferase, by methylating the m6A site in the 3'UTR of KAT2A mRNA in NSCLC cells. Intriguingly, NFIC was also found to negatively regulate the expression of KAT2A by directly binding to its promoter region.

CONCLUSIONS

Our findings demonstrated that NFIC suppresses the malignant phenotype of NSCLC cells by regulating gene expression at both the transcriptional and post-transcriptional levels. A deeper comprehension of the genetic and epigenetic regulatory mechanisms in tumorigenesis would be beneficial for the development of personalized treatment strategies.

Strigolactone induces D14-dependent large-scale changes in gene expression requiring SWI/SNF chromatin remodellers

Strigolactones (SL) function as plant hormones in control of multiple aspects of plant development, mostly via the regulation of gene expression. Immediate early-gene regulation by SL remains unexplored due to difficulty in dissecting early from late gene expression responses to SL. We used synthetic SL, rac-GR24 treatment of protoplasts and RNA-seq to explore early SL-induced changes in gene expression over time (5-180 minutes) and discovered rapid, dynamic and SL receptor D14-dependent regulation of gene expression in response to rac-GR24. Importantly, we discovered a significant dependence of SL signalling on chromatin remodelling processes, as the induction of a key SL-induced transcription factor BRANCHED1 requires the SWI/SNF chromatin remodelling ATPase SPLAYED (SYD) and leads to upregulation of a homologue SWI/SNF ATPase BRAHMA. ATAC-seq profiling of genome-wide changes in chromatin accessibility in response to rac-GR24 identified large-scale changes, with over 1400 differentially accessible regions. These changes in chromatin accessibility often precede transcriptional changes and are likely to harbour SL cis-regulatory elements. Importantly, we discovered that this early and extensive modification of the chromatin landscape also requires SYD. This study, therefore, provides evidence that SL signalling requires regulation of chromatin accessibility, and it identifies genomic locations harbouring likely SL cis-regulatory sequences.

OliTag-seq enhances in cellulo detection of CRISPR-Cas9 off-targets

The potential for off-target mutations is a critical concern for the therapeutic application of CRISPR-Cas9 gene editing. Current detection methodologies, such as GUIDE-seq, exhibit limitations in oligonucleotide integration efficiency and sensitivity, which could hinder their utility in clinical settings. To address these issues, we introduce OliTag-seq, an in-cellulo assay specifically engineered to enhance the detection of off-target events. OliTag-seq employs a stable oligonucleotide for precise break tagging and an innovative triple-priming amplification strategy, significantly improving the scope and accuracy of off-target site identification. This method surpasses traditional assays by providing comprehensive coverage across various sgRNAs and genomic targets. Our research particularly highlights the superior sensitivity of induced pluripotent stem cells (iPSCs) in detecting off-target mutations, advocating for using patient-derived iPSCs for refined off-target analysis in therapeutic gene editing. Furthermore, we provide evidence that prolonged Cas9 expression and transient HDAC inhibitor treatments enhance the assay's ability to uncover off-target events. OliTag-seq merges the high sensitivity typical of in vitro assays with the practical application of cellular contexts. This approach significantly improves the safety and efficacy profiles of CRISPR-Cas9 interventions in research and clinical environments, positioning it as an essential tool for the precise assessment and refinement of genome editing applications.

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.

Nucleosome-binding by TP53, TP63, and TP73 is determined by the composition, accessibility, and helical orientation of their binding sites

The p53 family of transcription factors plays key roles in driving development and combating cancer by regulating gene expression. TP53, TP63, and TP73-the three members of the p53 family-regulate gene expression by binding to their DNA binding sites, many of which are situated within nucleosomes. To thoroughly examine the nucleosome-binding abilities of the p53 family, we used Pioneer-seq, a technique that assesses a transcription factor's binding affinity to its DNA binding sites at all possible positions within the nucleosome core particle. Using Pioneer-seq, we analyzed the binding affinity of TP53, TP63, and TP73 to 10 p53-family binding sites across the nucleosome core particle. We found that the affinity of TP53, TP63, and TP73 for nucleosomes was largely determined by the positioning of p53-family binding sites within nucleosomes; p53-family members bind strongly to the more accessible edges of nucleosomes but weakly to the less accessible centers of nucleosomes. We also found that the DNA-helical orientation of p53-family binding sites within nucleosomal DNA impacted the nucleosome-binding affinity of p53-family members. The composition of their binding sites also impacted each p53-family member's nucleosome-binding affinities only when the binding site was located in an accessible location. Taken together, our results show that the accessibility, composition, and helical orientation of p53-family binding sites collectively determine the nucleosome-binding affinities of TP53, TP63, and TP73. These findings help explain the rules underlying p53-family-nucleosome binding and thus provide requisite insight into how we may better control gene-expression changes involved in development and tumor suppression.

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.

Chromatin accessibility profiling reveals that human fibroblasts respond to mechanical stimulation in a cell-specific manner

Fibroblasts in the skin are highly heterogeneous, both in vivo and in vitro. One difference between follicular (dermal papilla fibroblasts [DP]) and interfollicular fibroblasts (papillary fibroblasts [PFi]) in vitro is their ability to differentiate in response to osteogenic media (OM), or mechanical stimulation. Here, we asked whether differences in the ability of DP and PFi to respond to differentiation stimuli are due to differences in chromatin accessibility. We performed chromatin accessibility and transcriptional profiling of DP and PFi in human skin, which arise from a common progenitor during development, yet display distinct characteristics in adult tissue and in vitro. We found that cells cultured in growth media had unique chromatin accessibility profiles; however, these profiles control similar functional networks. Upon introduction of a chemical perturbation (OM) to promote differentiation, we observed a divergence not only in the accessible chromatin signatures but also in the functional networks controlled by these signatures. The biggest divergence between DP and PFi was observed when we applied 2 perturbations to cells: growth in OM and mechanical stimulation (a shock wave [OMSW]). DP readily differentiate into bone in OMSW conditions, while PFi lack differentiation capability in vitro. In the DP we found a number of uniquely accessible promoters that controlled osteogenic interaction networks associated with bone and differentiation functions. Using ATAC-seq and RNA-seq we found that the combination of 2 stimuli (OMSW) could result in significant changes in chromatin accessibility associated with osteogenic differentiation, but only within the DP (capable of osteogenic differentiation). De novo motif analysis identified enrichment of motifs bound by the TEA domain (TEAD) family of transcription factors, and inter-cell comparisons (UpSet analysis) displayed large groups of genes to be unique to single cell types and conditions. Our results suggest that these 2 stimuli (OMSW) elicit cell-specific responses by modifying chromatin accessibility of osteogenic-related gene promoters.

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.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Integrated ATAC-seq and RNA-seq data analysis identifies transcription factors related to rice stripe virus infection in Oryza sativa

Animal studies have shown that virus infection causes changes in host chromatin accessibility, but little is known about changes in chromatin accessibility of plants infected by viruses and its potential impact. Here, rice infected by rice stripe virus (RSV) was used to investigate virus-induced changes in chromatin accessibility. Our analysis identified a total of 6462 open- and 3587 closed-differentially accessible chromatin regions (DACRs) in rice under RSV infection by ATAC-seq. Additionally, by integrating ATAC-seq and RNA-seq, 349 up-regulated genes in open-DACRs and 126 down-regulated genes in closed-DACRs were identified, of which 34 transcription factors (TFs) were further identified by search of upstream motifs. Transcription levels of eight of these TFs were validated by reverse transcription-PCR. Importantly, four of these TFs (OsWRKY77, OsWRKY28, OsZFP12 and OsERF91) interacted with RSV proteins and are therefore predicted to play important roles in RSV infection. This is the first application of ATAC-seq and RNA-seq techniques to analyse changes in rice chromatin accessibility caused by RSV infection. Integrating ATAC-seq and RNA-seq provides a new approach to select candidate TFs in response to virus infection.

Epigenetic Effects of Resveratrol on Oncogenic Signaling in Breast Cancer

The crosstalk between oncogenic signaling pathways plays a crucial role in driving cancer development. We previously demonstrated that dietary polyphenols, specifically resveratrol (RSV) and other stilbenoids, epigenetically target oncogenes for silencing via DNA hypermethylation in breast cancer. In the present study, we identify signal transduction regulators among RSV-hypermethylated targets and investigate the functional role of RSV-mediated DNA hypermethylation in the regulation of Hedgehog and Wnt signaling. Non-invasive ER-positive MCF-7 and highly invasive triple-negative MCF10CA1a human breast cancer cell lines were used as experimental models. Upon 9-day exposure to 15 µM RSV, pyrosequencing and qRT-PCR were performed to assess DNA methylation and expression of GLI2 and WNT4, which are upstream regulators of the Hedgehog and Wnt pathways, respectively. Our results showed that RSV led to a DNA methylation increase within GLI2 and WNT4 enhancers, which was accompanied by decreases in gene expression. Consistently, we observed the downregulation of genes downstream of the Hedgehog and Wnt signaling, including common targets shared by both pathways, CCND1 and CYR61. Further analysis using chromatin immunoprecipitation identified increased H3K27 trimethylation and decreased H3K9 and H3K27 acetylation, along with abolishing OCT1 transcription factor binding. Those changes indicate a transcriptionally silent chromatin state at GLI2 and WNT4 enhancers. The inhibition of the Wnt signal transduction was confirmed using a phospho-antibody array that demonstrated suppression of positive and stimulation of negative Wnt regulators. In conclusion, our results provide scientific evidence for dietary polyphenols as epigenetics-modulating agents that act to re-methylate and silence oncogenes, reducing the oncogenic signal transduction. Targeting such an action could be an effective strategy in breast cancer prevention and/or adjuvant therapy.

Experimental evidence for cancer resistance in a bat species

Mammals exhibit different rates of cancer, with long-lived species generally showing greater resistance. Although bats have been suggested to be resistant to cancer due to their longevity, this has yet to be systematically examined. Here, we investigate cancer resistance across seven bat species by activating oncogenic genes in their primary cells. Both in vitro and in vivo experiments suggest that Myotis pilosus (MPI) is particularly resistant to cancer. The transcriptomic and functional analyses reveal that the downregulation of three genes (HIF1A, COPS5, and RPS3) largely contributes to cancer resistance in MPI. Further, we identify the loss of a potential enhancer containing the HIF1A binding site upstream of COPS5 in MPI, resulting in the downregulation of COPS5. These findings not only provide direct experimental evidence for cancer resistance in a bat species but also offer insights into the natural mechanisms of cancer resistance in mammals.

Eomes-dependent mitochondrial regulation promotes survival of pathogenic CD4+ T cells during inflammation

The mechanisms whereby Eomes controls tissue accumulation of T cells and strengthens inflammation remain ill-defined. Here, we show that Eomes deletion in antigen-specific CD4+ T cells is sufficient to protect against central nervous system (CNS) inflammation. While Eomes is dispensable for the initial priming of CD4+ T cells, it is required for long-term maintenance of CNS-infiltrating CD4+ T cells. We reveal that the impact of Eomes on effector CD4+ T cell longevity is associated with sustained expression of multiple genes involved in mitochondrial organization and functions. Accordingly, epigenetic studies demonstrate that Eomes supports mitochondrial function by direct binding to either metabolism-associated genes or mitochondrial transcriptional modulators. Besides, the significance of these findings was confirmed in CD4+ T cells from healthy donors and multiple sclerosis patients. Together, our data reveal a new mechanism by which Eomes promotes severity and chronicity of inflammation via the enhancement of CD4+ T cell mitochondrial functions and resistance to stress-induced cell death.

Decoding the Epigenetics and Chromatin Loop Dynamics of Androgen Receptor-Mediated Transcription

Androgen receptor (AR)-mediated transcription plays a critical role in normal prostate development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigated the AR-driven epigenetic and chromosomal chromatin looping changes. We collected a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we found that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we showed that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then proposed and experimentally validated an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide new insight into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Chromatin accessibility mediated transcriptome changes contribute to flavor substance alterations and jasmonic acid hyperaccumulation during oolong tea withering process

During the oolong tea withering process, abiotic stresses induce significant changes in the content of various flavor substances and jasmonic acid (JA). However, the changes in chromatin accessibility during withering and their potential impact remain poorly understood. By integrating ATAC-seq, RNA-seq, metabolite, and hormone assays, we characterized the withering treatment-induced changes in chromatin accessibility, gene expression levels, important metabolite contents, and JA and JA-ILE contents. Additionally, we analyzed the effects of chromatin accessibility alterations on gene expression changes, content changes of important flavor substances, and JA hyperaccumulation. Our analysis identified a total of 3451 open- and 13 426 close-differentially accessible chromatin regions (DACRs) under withering treatment. Our findings indicate that close-DACRs-mediated down-regulated differentially expressed genes (DEGs) resulted in the reduced accumulation of multiple catechins during withering, whereas open-DACRs-mediated up-regulated DEGs contributed to the increased accumulation of important terpenoids, JA, JA-ILE and short-chain C5/C6 volatiles. We further highlighted important DACRs-mediated DEGs associated with the synthesis of catechins, terpenoids, JA and JA and short-chain C5/C6 volatiles and confirmed the broad effect of close-DACRs on catechin synthesis involving almost all enzymes in the pathway during withering. Importantly, we identified a novel MYB transcription factor (CsMYB83) regulating catechin synthesis and verified the binding of CsMYB83 in the promoter-DACRs regions of key catechin synthesis genes using DAP-seq. Overall, our results not only revealed a landscape of chromatin alters-mediated transcription, flavor substance and hormone changes under oolong tea withering, but also provided target genes for flavor improvement breeding in tea plant.

Endothelial cells adopt a pro-reparative immune responsive signature during cardiac injury

Modulation of the heart's immune microenvironment is crucial for recovery after ischemic events such as myocardial infarction (MI). Endothelial cells (ECs) can have immune regulatory functions; however, interactions between ECs and the immune environment in the heart after MI remain poorly understood. We identified an EC-specific IFN responsive and immune regulatory gene signature in adult and pediatric heart failure (HF) tissues. Single-cell transcriptomic analysis of murine hearts subjected to MI uncovered an EC population (IFN-ECs) with immunologic gene signatures similar to those in human HF. IFN-ECs were enriched in regenerative-stage mouse hearts and expressed genes encoding immune responsive transcription factors (Irf7, Batf2, and Stat1). Single-cell chromatin accessibility studies revealed an enrichment of these TF motifs at IFN-EC signature genes. Expression of immune regulatory ligand genes by IFN-ECs suggests bidirectional signaling between IFN-ECs and macrophages in regenerative-stage hearts. Our data suggest that ECs may adopt immune regulatory signatures after cardiac injury to accompany the reparative response. The presence of these signatures in human HF and murine MI models suggests a potential role for EC-mediated immune regulation in responding to stress induced by acute injury in MI and chronic adverse remodeling in HF.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

WRN loss accelerates abnormal adipocyte metabolism in Werner syndrome

BACKGROUND

Metabolic dysfunction is one of the main symptoms of Werner syndrome (WS); however, the underlying mechanisms remain unclear. Here, we report that loss of WRN accelerates adipogenesis at an early stage both in vitro (stem cells) and in vivo (zebrafish). Moreover, WRN depletion causes a transient upregulation of late-stage of adipocyte-specific genes at an early stage.

METHODS

In an in vivo study, we generated wrn-/- mutant zebrafish and performed histological stain and Oil Red O staining to assess the fat metabolism. In an in vitro study, we used RNA-seq and ATAC-seq to profile the transcriptional features and chromatin accessibility in WRN depleted adipocytes. Moreover, we performed ChIP-seq to further study the regulatory mechanisms of metabolic dysfunction in WS.

RESULTS

Our findings show that mechanistically WRN deficiency causes SMARCA5 upregulation. SMARCA5 is crucial in chromatin remodeling and gene regulation. Additionally, rescuing WRN could normalize SMARCA5 expression and adipocyte differentiation. Moreover, we find that nicotinamide riboside (NR) supplementation restores adipocyte metabolism in both stem cells and zebrafish models.

CONCLUSIONS

Our findings unravel a new mechanism for the influence of WRN in the early stage of adipogenesis and provide a possible treatment for metabolic dysfunction in WS. These data provide promising insights into potential therapeutics for ageing and ageing-related diseases.

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

When Dad's Stress Gets under Kid's Skin-Impacts of Stress on Germline Cargo and Embryonic Development

Multiple lines of evidence suggest that paternal psychological stress contributes to an increased prevalence of neuropsychiatric and metabolic diseases in the progeny. While altered paternal care certainly plays a role in such transmitted disease risk, molecular factors in the germline might additionally be at play in humans. This is supported by findings on changes to the molecular make up of germ cells and suggests an epigenetic component in transmission. Several rodent studies demonstrate the correlation between paternal stress induced changes in epigenetic modifications and offspring phenotypic alterations, yet some intriguing cases also start to show mechanistic links in between sperm and the early embryo. In this review, we summarise efforts to understand the mechanism of intergenerational transmission from sperm to the early embryo. In particular, we highlight how stress alters epigenetic modifications in sperm and discuss the potential for these modifications to propagate modified molecular trajectories in the early embryo to give rise to aberrant phenotypes in adult offspring.

High-capacity sample multiplexing for single cell chromatin accessibility profiling

Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

Dose rate dependent reduction in chromatin accessibility at transcriptional start sites long time after exposure to gamma radiation

Ionizing radiation (IR) impact cellular and molecular processes that require chromatin remodelling relevant for cellular integrity. However, the cellular implications of ionizing radiation (IR) delivered per time unit (dose rate) are still debated. This study investigates whether the dose rate is relevant for inflicting changes to the epigenome, represented by chromatin accessibility, or whether it is the total dose that is decisive. CBA/CaOlaHsd mice were whole-body exposed to either chronic low dose rate (2.5 mGy/h for 54 d) or the higher dose rates (10 mGy/h for 14 d and 100 mGy/h for 30 h) of gamma radiation (60Co, total dose: 3 Gy). Chromatin accessibility was analysed in liver tissue samples using Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-Seq), both one day after and over three months post-radiation (>100 d). The results show that the dose rate contributes to radiation-induced epigenomic changes in the liver at both sampling timepoints. Interestingly, chronic low dose rate exposure to a high total dose (3 Gy) did not inflict long-term changes to the epigenome. In contrast to the acute high dose rate given to the same total dose, reduced accessibility at transcriptional start sites (TSS) was identified in genes relevant for the DNA damage response and transcriptional activity. Our findings link dose rate to essential biological mechanisms that could be relevant for understanding long-term changes after ionizing radiation exposure. However, future studies are needed to comprehend the biological consequence of these findings.

Integration of RNA-seq and ATAC-seq analyzes the effect of low dose neutron-γ radiation on gene expression of lymphocytes from oilfield logging workers

Objective: Although radiation workers are exposed to much lower doses of neutron-γ rays than those suffered in nuclear explosions and accidents, it does not mean that their health is not affected by radiation. Lower doses of radiation do not always cause morphological aberrations in chromosomes, so more sophisticated tests must be sought to specific alterations in the exposed cells. Our goal was to characterize the specific gene expression in lymphocytes from logging workers who were continuously exposed to low doses of neutron-γ radiation. We hypothesized that the combination of cell type-specific transcriptomes and open chromatin profiles would identify lymphocyte-specific gene alterations induced by long-term radiation with low-dose neutron-γ-rays and discover new regulatory pathways and transcriptional regulatory elements. Methods: Lymphocytes were extracted from workers who have been occupationally exposed to neutron-γ and workers unexposed to radiation in the same company. mRNA-seq and ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) were performed, followed integrative analysis to identify specific gene regulatory regions induced by neutron-γ radiation. A qPCR assay was then performed to verify the downregulation of RNA coding for ribosomal proteins and flow cytometry was used to detect ribosomal protein expression and cell cycle alterations. Results: We identified transcripts that were specifically induced by neutron-γ radiation and discovered differential open chromatin regions that correlated with these gene activation patterns. Notably, we observed a downward trend in the expression of both differentially expressed genes and open chromatin peaks. Our most significant finding was that the differential peak upregulated in ATAC-seq, while the differential gene was downregulated in the ribosome pathway. We confirmed that neutron-γ radiation leads to transcriptional inhibition by analyzing the most enriched promoters, examining RPS18 and RPS27A expression by qPCR, and analyzing protein-protein interactions of the differential genes. Ribosomal protein expression and cell cycle were also affected by neutron-γ as detected by flow cytometry. Conclusion: We have comprehensively analyzed the genetic landscape of human lymphocytes based on chromatin accessibility and transcript levels, enabling the identification of novel neutron-γ induced signature genes not previously known. By comparing fine-mapping of open chromatin and RNA reads, we have determined that neutron-γ specifically leads to downregulation of genes in the ribosome pathway, with pseudogenes potentially playing a crucial role.

Prospects and perspectives: inferring physiological and regulatory targets for CAM from molecular and modelling approaches

BACKGROUND AND SCOPE

This review summarizes recent advances in our understanding of Crassulacean Acid Metabolism (CAM) by integrating evolutionary, ecological, physiological, metabolic and molecular perspectives. A number of key control loops which moderate the expression of CAM phases, and their metabolic and molecular control, are explored. These include nocturnal stomatal opening, activation of phosphoenolpyruvate carboxylase by a specific protein kinase, interactions with circadian clock control, as well as daytime decarboxylation and activation of Rubisco. The vacuolar storage and release of malic acid and the interplay between the supply and demand for carbohydrate reserves are also key metabolic control points.

FUTURE OPPORTUNITIES

We identify open questions and opportunities, with experimentation informed by top-down molecular modelling approaches allied with bottom-up mechanistic modelling systems. For example, mining transcriptomic datasets using high-speed systems approaches will help to identify targets for future genetic manipulation experiments to define the regulation of CAM (whether circadian or metabolic control). We emphasize that inferences arising from computational approaches or advanced nuclear sequencing techniques can identify potential genes and transcription factors as regulatory targets. However, these outputs then require systematic evaluation, using genetic manipulation in key model organisms over a developmental progression, combining gene silencing and metabolic flux analysis and modelling to define functionality across the CAM day-night cycle. From an evolutionary perspective, the origins and function of CAM succulents and responses to water deficits are set against the mesophyll and hydraulic limitations imposed by cell and tissue succulence in contrasting morphological lineages. We highlight the interplay between traits across shoots (3D vein density, mesophyll conductance and cell shrinkage) and roots (xylem embolism and segmentation). Thus, molecular, biophysical and biochemical processes help to curtail water losses and exploit rapid rehydration during restorative rain events. In the face of a changing climate, we hope such approaches will stimulate opportunities for future research.

Transcriptional responses of cancer cells to heat shock-inducing stimuli involve amplification of robust HSF1 binding

Responses of cells to stimuli are increasingly discovered to involve the binding of sequence-specific transcription factors outside of known target genes. We wanted to determine to what extent the genome-wide binding and function of a transcription factor are shaped by the cell type versus the stimulus. To do so, we induced the Heat Shock Response pathway in two different cancer cell lines with two different stimuli and related the binding of its master regulator HSF1 to nascent RNA and chromatin accessibility. Here, we show that HSF1 binding patterns retain their identity between basal conditions and under different magnitudes of activation, so that common HSF1 binding is globally associated with distinct transcription outcomes. HSF1-induced increase in DNA accessibility was modest in scale, but occurred predominantly at remote genomic sites. Apart from regulating transcription at existing elements including promoters and enhancers, HSF1 binding amplified during responses to stimuli may engage inactive chromatin.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

cisDynet: An integrated platform for modeling gene-regulatory dynamics and networks

Chromatin accessibility sequencing has been widely used for uncovering genetic regulatory mechanisms and inferring gene regulatory networks. However, effectively integrating large-scale chromatin accessibility datasets has posed a significant challenge. This is due to the lack of a comprehensive end-to-end solution, as many existing tools primarily emphasize data preprocessing and overlook downstream analyses. To bridge this gap, we have introduced cisDynet, a holistic solution that combines streamlined data preprocessing using Snakemake and R functions with advanced downstream analysis capabilities. cisDynet excels in conventional data analyses, encompassing peak statistics, peak annotation, differential analysis, motif enrichment analysis, and more. Additionally, it allows to perform sophisticated data exploration, such as tissue-specific peak identification, time course data modeling, integration of RNA-seq data to establish peak-to-gene associations, constructing regulatory networks, and conducting enrichment analysis of genome-wide association study (GWAS) variants. As a proof of concept, we applied cisDynet to reanalyze comprehensive ATAC-seq datasets across various tissues from the Encyclopedia of DNA Elements (ENCODE) project. The analysis successfully delineated tissue-specific open chromatin regions (OCRs), established connections between OCRs and target genes, and effectively linked these discoveries with 1861 GWAS variants. Furthermore, cisDynet was instrumental in dissecting the time course open chromatin data of mouse embryonic development, revealing the dynamic behavior of OCRs over developmental stages and identifying key transcription factors governing differentiation trajectories. In summary, cisDynet offers researchers a user-friendly solution that minimizes the need for extensive coding, ensures the reproducibility of results, and greatly simplifies the exploration of epigenomic data.

Prevention, Diagnosis and Eradication of Mycoplasma Contamination in Cell Culture

Mycoplasmas, which are the smallest and simplest prokaryotes, lack a cell wall but possess the ability to undergo self-replication. Mycoplasma contamination is a common problem for laboratories engaging in cell culture. Due to their small size, Mycoplasmas can easily permeate filters designed to prevent bacterial and fungal contamination in cell culture. Although Mycoplasma contamination usually does not result in cell death, it can significantly affect cell proliferation, metabolism, and cause chromosomal aberrations. Therefore, it is crucial to detect and eliminate Mycoplasma contamination in cell culture. This step-by-step protocol presents a comprehensive approach to prevent Mycoplasma contamination in cell culture, as well as to detect and eradicate Mycoplasma to ensure accurate experimental and sequencing results.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

Chromatin accessibility memory of donor cells disrupts bovine somatic cell nuclear transfer blastocysts development

The post-transfer developmental capacity of bovine somatic cell nuclear transfer (SCNT) blastocysts is reduced, implying that abnormalities in gene expression regulation are present at blastocyst stage. Chromatin accessibility, as an indicator for transcriptional regulatory elements mediating gene transcription activity, has heretofore been largely unexplored in SCNT embryos, especially at blastocyst stage. In the present study, single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) of in vivo and SCNT blastocysts were conducted to segregate lineages and demonstrate the aberrant chromatin accessibility of transcription factors (TFs) related to inner cell mass (ICM) development in SCNT blastocysts. Pseudotime analysis of lineage segregation further reflected dysregulated chromatin accessibility dynamics of TFs in the ICM of SCNT blastocysts compared to their in vivo counterparts. ATAC- and ChIP-seq results of SCNT donor cells revealed that the aberrant chromatin accessibility in the ICM of SCNT blastocysts was due to the persistence of chromatin accessibility memory at corresponding loci in the donor cells, with strong enrichment of trimethylation of histone H3 at lysine 4 (H3K4me3) at these loci. Correction of the aberrant chromatin accessibility through demethylation of H3K4me3 by KDM5B diminished the expression of related genes (e.g., BCL11B) and significantly improved the ICM proliferation in SCNT blastocysts. This effect was confirmed by knocking down BCL11B in SCNT embryos to down-regulate p21 and alleviate the inhibition of ICM proliferation. These findings expand our understanding of the chromatin accessibility abnormalities in SCNT blastocysts and BCL11B may be a potential target to improve SCNT efficiency.

Prediction accuracy of regulatory elements from sequence varies by functional sequencing technique

Introduction

Various sequencing based approaches are used to identify and characterize the activities of cis-regulatory elements in a genome-wide fashion. Some of these techniques rely on indirect markers such as histone modifications (ChIP-seq with histone antibodies) or chromatin accessibility (ATAC-seq, DNase-seq, FAIRE-seq), while other techniques use direct measures such as episomal assays measuring the enhancer properties of DNA sequences (STARR-seq) and direct measurement of the binding of transcription factors (ChIP-seq with transcription factor-specific antibodies). The activities of cis-regulatory elements such as enhancers, promoters, and repressors are determined by their sequence and secondary processes such as chromatin accessibility, DNA methylation, and bound histone markers.

Methods

Here, machine learning models are employed to evaluate the accuracy with which cis-regulatory elements identified by various commonly used sequencing techniques can be predicted by their underlying sequence alone to distinguish between cis-regulatory activity that is reflective of sequence content versus secondary processes.

Results and discussion

Models trained and evaluated on D. melanogaster sequences identified through DNase-seq and STARR-seq are significantly more accurate than models trained on sequences identified by H3K4me1, H3K4me3, and H3K27ac ChIP-seq, FAIRE-seq, and ATAC-seq. These results suggest that the activity detected by DNase-seq and STARR-seq can be largely explained by underlying DNA sequence, independent of secondary processes. Experimentally, a subset of DNase-seq and H3K4me1 ChIP-seq sequences were tested for enhancer activity using luciferase assays and compared with previous tests performed on STARR-seq sequences. The experimental data indicated that STARR-seq sequences are substantially enriched for enhancer-specific activity, while the DNase-seq and H3K4me1 ChIP-seq sequences are not. Taken together, these results indicate that the DNase-seq approach identifies a broad class of regulatory elements of which enhancers are a subset and the associated data are appropriate for training models for detecting regulatory activity from sequence alone, STARR-seq data are best for training enhancer-specific sequence models, and H3K4me1 ChIP-seq data are not well suited for training and evaluating sequence-based models for cis-regulatory element prediction.

Regulation of chromatin accessibility and gene expression in the developing hippocampal primordium by LIM-HD transcription factor LHX2

In the mammalian cerebral cortex, the hippocampal primordium (Hcp) occupies a discrete position in the dorsal telencephalic neuroepithelium adjacent to the neocortical primordium (Ncp). We examined transcriptomic and chromatin-level features that distinguish the Hcp from the Ncp in the mouse during the early neurogenic period, embryonic day (E)12.5. ATAC-seq revealed that the Hcp was more accessible than the Ncp at this stage. Motif analysis of the differentially accessible loci in these tissues revealed LHX2 as a candidate transcription factor for modulating gene regulatory networks (GRNs). We analyzed LHX2 occupancy profiles and compared these with transcriptomic data from control and Lhx2 mutant Hcp and Ncp at E12.5. Our results revealed that LHX2 directly regulates distinct genes in the Hcp and Ncp within a set of common pathways that control fundamental aspects of development namely pluripotency, axon pathfinding, Wnt, and Hippo signaling. Loss of Lhx2 caused a decrease in accessibility, specifically in hippocampal chromatin, suggesting that this factor may play a unique role in hippocampal development. We identified 14 genes that were preferentially enriched in the Hcp, for which LHX2 regulates both chromatin accessibility and mRNA expression, which have not thus far been examined in hippocampal development. Together, these results provide mechanistic insight into how LHX2 function in the Hcp may contribute to the process by which the hippocampus acquires features distinct from the neocortex.

Loss of epigenetic suppression of retrotransposons with oncogenic potential in aging mammary luminal epithelial cells

A primary function of DNA methylation in mammalian genomes is to repress transposable elements (TEs). The widespread methylation loss that is commonly observed in cancer cells results in the loss of epigenetic repression of TEs. The aging process is similarly characterized by changes to the methylome. However, the impact of these epigenomic alterations on TE silencing and the functional consequences of this have remained unclear. To assess the epigenetic regulation of TEs in aging, we profiled DNA methylation in human mammary luminal epithelial cells (LEps)-a key cell lineage implicated in age-related breast cancers-from younger and older women. We report here that several TE subfamilies function as regulatory elements in normal LEps, and a subset of these display consistent methylation changes with age. Methylation changes at these TEs occurred at lineage-specific transcription factor binding sites, consistent with loss of lineage specificity. Whereas TEs mainly showed methylation loss, CpG islands (CGIs) that are targets of the Polycomb repressive complex 2 (PRC2) show a gain of methylation in aging cells. Many TEs with methylation loss in aging LEps have evidence of regulatory activity in breast cancer samples. We furthermore show that methylation changes at TEs impact the regulation of genes associated with luminal breast cancers. These results indicate that aging leads to DNA methylation changes at TEs that undermine the maintenance of lineage specificity, potentially increasing susceptibility to breast cancer.

Comparative valuation of the chlorpyrifos, acetamiprid, and lambda-cyhalothrin toxicity and their hematological and histopathological consequences in pigeons

Globally agrochemicals are widely used in the agricultural sectors, posing potential eco-toxicological risks and disrupting various lifeforms including birds. Thus, the current work was conducted to compare the acute toxic impacts of pesticides (e.g., chlorpyrifos, acetamiprid, and lambda-cyhalothrin) on the pigeon's health. In total 50 adult pigeons were purchased from a local market where these pigeons were fed on pollution-free food. Post adaptation period (15 days), the pigeons were arbitrarily separated into five distinct groups after having been identified in this manner by chance (each group containing 10 pigeons). Control group (group 1) was not treated with any pesticide while the remaining groups (groups 2, 3, and 4) were treated with 0.25-mg/kg body weight of chlorpyrifos, acetamiprid, lambda-cyhalothrin, and a mixture of all three pesticides (group 5), respectively. After 36 days of exposure, the groups that had been exposed to the pesticide showed a significant (p < 0.05) increase in both the total number of platelets and the number of white blood cells (WBCs), in comparison to the control group. On the other hand, the groups that were exposed to the insecticides had significantly lower levels of red blood cells (RBCs), hemoglobin (Hb), and packed cell volume (PCV) (p < 0.05). The value of mean corpuscular volume (MCV) was significantly (p < 0.05) reduced in acetamiprid-exposed group, while a significant increase was observed in other pesticide-exposed groups. Obvious histopathological changes were observed in the tissues of control group and no such changes were reported by control group. Necrosis, pyknosis, lymphocyte infiltration, congestion of blood, dissolution of plasma membrane, and vacuolation were observed in the livers of pesticide-treated pigeons. The intestinal study showed the formation of goblet cells, villi rupturing, degeneration of serosa, necrosis, and pyknosis in treated groups. Renal alterations, dilation of renal tubules, reduction of glomerulus tissue, and edema were observed. This study manifests that the uncontrolled use of pesticides impairs ecosystems and poses a substantial health risk to wildlife and ultimately to human.

Predicting chromosomal compartments directly from the nucleotide sequence with DNA-DDA

Three-dimensional (3D) genome architecture is characterized by multi-scale patterns and plays an essential role in gene regulation. Chromatin conformation capturing experiments have revealed many properties underlying 3D genome architecture, such as the compartmentalization of chromatin based on transcriptional states. However, they are complex, costly and time consuming, and therefore only a limited number of cell types have been examined using these techniques. Increasing effort is being directed towards deriving computational methods that can predict chromatin conformation and associated structures. Here we present DNA-delay differential analysis (DDA), a purely sequence-based method based on chaos theory to predict genome-wide A and B compartments. We show that DNA-DDA models derived from a 20 Mb sequence are sufficient to predict genome wide compartmentalization at the scale of 100 kb in four different cell types. Although this is a proof-of-concept study, our method shows promise in elucidating the mechanisms responsible for genome folding as well as modeling the impact of genetic variation on 3D genome architecture and the processes regulated thereby.

C/EBPα Confers Dependence to Fatty Acid Anabolic Pathways and Vulnerability to Lipid Oxidative Stress-Induced Ferroptosis in FLT3-Mutant Leukemia

Although transcription factor CCAAT-enhancer binding protein α (C/EBPα) is critical for normal and leukemic differentiation, its role in cell and metabolic homeostasis is largely unknown in cancer. Here, multiomics analyses uncovered a coordinated activation of C/EBPα and Fms-like tyrosine kinase 3 (FLT3) that increased lipid anabolism in vivo and in patients with FLT3-mutant acute myeloid leukemia (AML). Mechanistically, C/EBPα regulated the fatty acid synthase (FASN)-stearoyl-CoA desaturase (SCD) axis to promote fatty acid (FA) biosynthesis and desaturation. We further demonstrated that FLT3 or C/EBPα inactivation decreased monounsaturated FA incorporation to membrane phospholipids through SCD downregulation. Consequently, SCD inhibition enhanced susceptibility to lipid redox stress that was exploited by combining FLT3 and glutathione peroxidase 4 inhibition to trigger lipid oxidative stress, enhancing ferroptotic death of FLT3-mutant AML cells. Altogether, our study reveals a C/EBPα function in lipid homeostasis and adaptation to redox stress, and a previously unreported vulnerability of FLT3-mutant AML to ferroptosis with promising therapeutic application.

SIGNIFICANCE

FLT3 mutations are found in 30% of AML cases and are actionable by tyrosine kinase inhibitors. Here, we discovered that C/EBPα regulates FA biosynthesis and protection from lipid redox stress downstream mutant-FLT3 signaling, which confers a vulnerability to ferroptosis upon FLT3 inhibition with therapeutic potential in AML. This article is highlighted in the In This Issue feature, p. 1501.

Simultaneous Mapping of DNA Binding and Nucleosome Positioning with SpLiT-ChEC

The organization of chromatin - including the positions of nucleosomes and the binding of other proteins to DNA - helps define transcriptional profiles in eukaryotic organisms. While techniques like ChIP-Seq and MNase-Seq can map protein-DNA and nucleosome localization separately, assays designed to simultaneously capture nucleosome positions and protein-DNA interactions can produce a detailed picture of the chromatin landscape. Most assays that monitor chromatin organization and protein binding rely on antibodies, which often exhibit nonspecific binding, and/or the addition of bulky adducts to the DNA-binding protein being studied, which can affect their expression and activity. Here, we describe SpyCatcher Linked Targeting of Chromatin Endogenous Cleavage (SpLiT-ChEC), where a 13-amino acid SpyTag peptide, appended to a protein of interest, serves as a highly-specific targeting moiety for in situ enzymatic digestion. The SpyTag/SpyCatcher system forms a covalent bond, linking the target protein and a co-expressed MNase-SpyCatcher fusion construct. SpyTagged proteins are expressed from endogenous loci, whereas MNase-SpyCatcher expression is induced immediately before harvesting cultures. MNase is activated with high concentrations of calcium, which primarily digests DNA near target protein binding sites. By sequencing the DNA fragments released by targeted MNase digestion, we found that this method recovers information on protein binding and proximal nucleosome positioning. SpLiT-ChEC provides precise temporal control that we anticipate can be used to monitor chromatin under various conditions and at distinct points in the cell cycle.

Characterization and Optimization of Multiomic Single-Cell Epigenomic Profiling

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n = 18) and a solid tumor type, ovarian cancer (OC, n = 18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell type identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n = 6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

Transposons (TEs) account for more than 80% of the wheat genome, the highest among all known crop species. They play an important role in shaping the elaborate genomic landscape, which is the key to the speciation of wheat. In this study, we analyzed the association between TEs, chromatin states, and chromatin accessibility in Aegilops tauschii, the D genome donor of bread wheat. We found that TEs contributed to the complex but orderly epigenetic landscape as chromatin states showed diverse distributions on TEs of different orders or superfamilies. TEs also contributed to the chromatin state and openness of potential regulatory elements, affecting the expression of TE-related genes. Some TE superfamilies, such as hAT-Ac, carry active/open chromatin regions. In addition, the histone mark H3K9ac was found to be associated with the accessibility shaped by TEs. These results suggest the role of diversiform TEs in shaping the epigenetic landscape and in gene expression regulation in Aegilops tauschii. This has positive implications for understanding the transposon roles in Aegilops tauschii or the wheat D genome.

Integration of ATAC-Seq and RNA-Seq reveals FOSL2 drives human liver progenitor-like cell aging by regulating inflammatory factors

BACKGROUND

Human primary hepatocytes (PHCs) are considered to be the best cell source for cell-based therapies for the treatment of end-stage liver disease and acute liver failure. To obtain sufficient and high-quality functional human hepatocytes, we have established a strategy to dedifferentiate human PHCs into expandable hepatocyte-derived liver progenitor-like cells (HepLPCs) through in vitro chemical reprogramming. However, the reduced proliferative capacity of HepLPCs after long-term culture still limits their utility. Therefore, in this study, we attempted to explore the potential mechanism related to the proliferative ability of HepLPCs in vitro culture.

RESULTS

In this study, analysis of assay for transposase accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were performed for PHCs, proliferative HepLPCs (pro-HepLPCs) and late-passage HepLPCs (lp-HepLPCs). Genome-wide transcriptional and chromatin accessibility changes during the conversion and long-term culture of HepLPCs were studied. We found that lp-HepLPCs exhibited an aged phenotype characterized by the activation of inflammatory factors. Epigenetic changes were found to be consistent with our gene expression findings, with promoter and distal regions of many inflammatory-related genes showing increased accessibility in the lp-HepLPCs. FOSL2, a member of the AP-1 family, was found to be highly enriched in the distal regions with increased accessibility in lp-HepLPCs. Its depletion attenuated the expression of aging- and senescence-associated secretory phenotype (SASP)-related genes and resulted in a partial improvement of the aging phenotype in lp-HepLPCs.

CONCLUSIONS

FOSL2 may drive the aging of HepLPCs by regulating inflammatory factors and its depletion may attenuate this phenotypic shift. This study provides a novel and promising approach for the long-term in vitro culture of HepLPCs.