Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity

Enhancers in Plant Development, Adaptation and Evolution

Understanding plant responses to developmental and environmental cues is crucial for studying morphological divergence and local adaptation. Gene expression changes, governed by cis-regulatory modules (CRMs) including enhancers, are a major source of plant phenotypic variation. However, while genome-wide approaches have revealed thousands of putative enhancers in mammals, far fewer have been identified and functionally characterized in plants. This review provides an overview of how enhancers function to control gene regulation, methods to predict DNA sequences that may have enhancer activity, methods utilized to functionally validate enhancers and the current knowledge of enhancers in plants, including how they impact plant development, response to environment and evolutionary adaptation.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

The combinatorial binding syntax of transcription factors in forebrain-specific enhancers

Tissue-specific gene regulation in mammals involves the coordinated binding of multiple transcription factors (TFs). Using the forebrain as a model, we investigated the syntax of TF occupancy to determine tissue-specific enhancer regions. We analyzed forebrain-exclusive enhancers from the VISTA Enhancer Browser and a curated set of 23 TFs relevant to forebrain development and disease. Our findings revealed multiple distinct patterns of combinatorial TF binding, with the HES5-FOXP2-GATA3 triad being the most frequent in forebrain-specific enhancers. This syntactic structure was detected in 2614 enhancers from a genome-wide catalog of 25,000 predicted human forebrain enhancers. Notably, this catalog represents a computationally predicted dataset, distinct from the in vivo validated set of enhancers obtained from the VISTA Enhancer Browser. The shortlisted 2614 enhancers were further analyzed using genome-wide epigenetic data and evaluated for evolutionary conservation and disease relevance. Our findings highlight the value of these 2614 enhancers in forebrain-specific gene regulation and provide a framework for discovering tissue-specific enhancers, enhancing the understanding of enhancer function.

Cross-cellular analysis of chromatin accessibility markers H3K4me3 and DNase in the context of detecting cell-identity genes: An "all-or-nothing" approach

Cell identity is often associated to a subset of highly-expressed genes that define the cell processes, as opposed to essential genes that are always active. Cell-specific genes may be defined in opposition to essential genes, or via experimental means. Detection of said cell-specific genes is often a primary goal in the study of novel biosamples. Chromatin accessibility markers (such as DNase and H3K4me3) help identify actively transcribed genes, but data can be difficult to come by for entirely novel biosamples. In this study, we investigate the possibility of associating the cell-specificity status of genes with chromatin accessibility markers from different cell lines, and we suggest that the number of cell lines in which a gene is found to be marked by DNase/H3K4me3 is predictive of the essentiality status itself. We define a measure called the Cross-cellular Chromatin Openness (CCO) level, and show that it is associated with the essentiality status using two differentiation experiments. We then compare the CCO-level predictive power to existing scRNA-Seq and bulk RNA-Seq methods, showing it has good concordance when applicable.

Liver-specific gene PGRMC1 blocks c-Myc-induced hepatocarcinogenesis through ER stress-independent PERK activation

Roles of liver-specific genes (LSGs) in tumor initiation and progression are rarely explored in hepatocellular carcinoma (HCC). Here we show that LSGs are generally downregulated in HCC tumor tissues compared to non-HCC liver tissues, and low-LSG HCCs show poor prognosis and the activated c-Myc pathway. Among the c-Myc- and patient prognosis-associated LSGs, PGRMC1 significantly blocks c-Myc-induced orthotopic HCC formation. The role of PGRMC1 depends on its localization to the endoplasmic reticulum (ER) membrane, where PGRMC1 interacts with PERK through their ER luminal domains. This interaction in turn activates PERK in an ER stress-independent manner, which phosphorylates eIF2α and consequently inhibits c-Myc protein translation. In HCC patients, PGRMC1 level is significantly reduced in tumor tissues and negatively associated with the c-Myc signature. Patients with low-PGRMC1 in their tumors have poor prognosis. Collectively, deregulated LSGs in HCC are associated with the c-Myc pathway activation and PGRMC1 blocks c-Myc-induced hepatic carcinogenesis through promoting ER stress-independent PERK activation.

Integrating multi-omics features enables non-invasive early diagnosis and treatment response prediction of diffuse large B-cell lymphoma

BACKGROUND

Multi-omics features of cell-free DNA (cfDNA) can effectively improve the performance of non-invasive early diagnosis and prognosis of cancer. However, multimodal characterization of cfDNA remains technically challenging.

METHODS

We developed a comprehensive multi-omics solution (COMOS) to specifically obtain an extensive fragmentomics landscape, presented by breakpoint characteristics of nucleosomes, CpG islands, DNase clusters and enhancers, besides typical methylation, copy number alteration of cfDNA. The COMOS was tested on 214 plasma samples of diffuse large B-cell lymphoma (DLBCL) and matched healthy controls.

RESULTS

For early diagnosis, COMOS improved the area under the curve (AUC) value to .993 compared with the individual omics model, with a sensitivity of 95% at 98% specificity. Detection sensitivity achieved 91% at 99% specificity in early-stage patients, while the AUC values of the individual omics model were 0.942, 0.968, 0.989, 0.935, 0.921, 0.781 and 0.917, respectively, with lower sensitivity and specificity. In the treatment response cohort, COMOS yielded a superior sensitivity of 88% at 86% specificity (AUC, 0.903). COMOS has achieved excellent performance in early diagnosis and treatment response prediction.

CONCLUSIONS

Our study provides an effectively improved approach with high accuracy for the diagnosis and prognosis of DLBCL, showing great potential for future clinical application.

KEY POINTS

A comprehensive multi-omics solution to specifically obtain an extensive fragmentomics landscape, presented by breakpoint characteristics of nucleosomes, CpG islands, DNase clusters and enhancers, besides typical methylation, copy number alteration of cfDNA. Integrated model of cfDNA multi-omics could be used for non-invasive early diagnosis of DLBCL. Integrated model of cfDNA multi-omics could effectively evaluate the efficacy of R-CHOP before DLBCL treatment.

Role of the Egr2 Promoter Antisense RNA in Modulating the Schwann Cell Chromatin Landscape

Background: Schwann cells (SCs) and their plasticity contribute to the peripheral nervous system's capacity for nerve regeneration after injury. The Egr2/Krox20 promoter antisense RNA (Egr2-AS) recruits chromatin remodeling complexes to inhibit Egr2 transcription following peripheral nerve injury. Methods: RNA-seq and ATAC-seq were performed on control cells, Lenti-GFP-transduced cells, and cells overexpressing Egr2-AS (Lenti-AS). Egr2 AS-RNA was cloned into the pLVX-DsRed-Express2-N1 lentiviral expression vector (Clontech, Mountain View, CA, USA), and the levels of AS-RNA expression were determined. Ezh2 and Wdr5 were immunoprecipitated from rat SCs and RT-qPCR was performed against AS-Egr2 RNA. ChIP followed by DNA purification columns was used to perform qPCR for relevant promoters. Hi-C, HiC-DC+, R, Bioconductor, and TOBIAS were used for significant and differential loop analysis, identifications of COREs and CORE-promotor loops, comparisons of TF activity at promoter sites, and identification of site-specific TF footprints. OnTAD was used to detect TADs, and Juicer was used to identify A/B compartments. Results: Here we show that a Neuregulin-ErbB2/3 signaling axis mediates binding of the Egr2-AS to YY1Ser184 and regulates its expression. Egr2-AS modulates the chromatin accessibility of Schwann cells and interacts with two distinct histone modification complexes. It binds to EZH2 and WDR5 and enables targeting of H3K27me3 and H3K4me3 to promoters of Egr2 and C-JUN, respectively. Expression of the Egr2-AS results in reorganization of the global chromatin landscape and quantitative changes in the loop formation and contact frequency at domain boundaries exhibiting enrichment for AP-1 genes. In addition, the Egr2-AS induces changes in the hierarchical TADs and increases transcription factor binding scores on an inter-TAD loop between a super-enhancer regulatory hub and the promoter of mTOR. Conclusions: Our results show that Neuregulin-ErbB2/3-YY1 regulates the expression of Egr2-AS, which mediates remodeling of the chromatin landscape in Schwann cells.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Construction of transcript regulation mechanism prediction models based on binding motif environment of transcription factor AoXlnR in Aspergillus oryzae

DNA-binding transcription factors (TFs) play a central role in transcriptional regulation mechanisms, mainly through their specific binding to target sites on the genome and regulation of the expression of downstream genes. Therefore, a comprehensive analysis of the function of these TFs will lead to the understanding of various biological mechanisms. However, the functions of TFs in vivo are diverse and complicated, and the identified binding sites on the genome are not necessarily involved in the regulation of downstream gene expression. In this study, we investigated whether DNA structural information around the binding site of TFs can be used to predict the involvement of the binding site in the regulation of the expression of genes located downstream of the binding site. Specifically, we calculated the structural parameters based on the DNA shape around the DNA binding motif located upstream of the gene whose expression is directly regulated by one TF AoXlnR from Aspergillus oryzae, and showed that the presence or absence of expression regulation can be predicted from the sequence information with high accuracy ([Formula: see text]-1.0) by machine learning incorporating these parameters.

Predmoter-cross-species prediction of plant promoter and enhancer regions

Motivation

Identifying cis-regulatory elements (CREs) is crucial for analyzing gene regulatory networks. Next generation sequencing methods were developed to identify CREs but represent a considerable expenditure for targeted analysis of few genomic loci. Thus, predicting the outputs of these methods would significantly cut costs and time investment.

Results

We present Predmoter, a deep neural network that predicts base-wise Assay for Transposase Accessible Chromatin using sequencing (ATAC-seq) and histone Chromatin immunoprecipitation DNA-sequencing (ChIP-seq) read coverage for plant genomes. Predmoter uses only the DNA sequence as input. We trained our final model on 21 species for 13 of which ATAC-seq data and for 17 of which ChIP-seq data was publicly available. We evaluated our models on Arabidopsis thaliana and Oryza sativa. Our best models showed accurate predictions in peak position and pattern for ATAC- and histone ChIP-seq. Annotating putatively accessible chromatin regions provides valuable input for the identification of CREs. In conjunction with other in silico data, this can significantly reduce the search space for experimentally verifiable DNA-protein interaction pairs.

Availability and implementation

The source code for Predmoter is available at: https://github.com/weberlab-hhu/Predmoter. Predmoter takes a fasta file as input and outputs h5, and optionally bigWig and bedGraph files.

Genome-wide screening and functional validation of methylation barriers near promoters

CpG islands near promoters are normally unmethylated despite being surrounded by densely methylated regions. Aberrant hypermethylation of these CpG islands has been associated with the development of various human diseases. Although local genetic elements have been speculated to play a role in protecting promoters from methylation, only a limited number of methylation barriers have been identified. In this study, we conducted an integrated computational and experimental investigation of colorectal cancer methylomes. Our study revealed 610 genes with disrupted methylation barriers. Genomic sequences of these barriers shared a common 41-bp sequence motif (MB-41) that displayed homology to the chicken HS4 methylation barrier. Using the CDKN2A (P16) tumor suppressor gene promoter, we validated the protective function of MB-41 and showed that loss of such protection led to aberrant hypermethylation. Our findings highlight a novel sequence signature of cis-acting methylation barriers in the human genome that safeguard promoters from silencing.

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.

Cell-free DNA fragmentomics and second malignant neoplasm risk in patients with PTEN hamartoma tumor syndrome

Individuals with PTEN hamartoma tumor syndrome (PHTS) harbor pathogenic germline PTEN variants that confer a significantly increased lifetime risk of various organ-specific cancers including second primary malignant neoplasms (SMNs). Currently, there are no reliable biomarkers that can predict individual-level cancer risk. Despite the highly promising value of cell-free DNA (cfDNA) as a biomarker for underlying sporadic cancers, the utility of cfDNA in individuals with known cancer-associated germline variants and subclinical cancers remains poorly understood. We perform ultra-low-pass whole-genome sequencing (ULP-WGS) of cfDNA from plasma samples from patients with PHTS and cancer as well as those without cancer. Analysis of cfDNA reveals that patients with PHTS and SMNs have distinct cfDNA size distribution, aberrant genome-wide fragmentation, and differential fragment end motif frequencies. Our work provides evidence that cfDNA profiles may be used as a marker for SMN risk in patients with PHTS.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

ROCCO: a robust method for detection of open chromatin via convex optimization

MOTIVATION

Analysis of open chromatin regions across multiple samples from two or more distinct conditions can determine altered gene regulatory patterns associated with biological phenotypes and complex traits. The ATAC-seq assay allows for tractable genome-wide open chromatin profiling of large numbers of samples. Stable, broadly applicable genomic annotations of open chromatin regions are not available. Thus, most studies first identify open regions using peak calling methods for each sample independently. These are then heuristically combined to obtain a consensus peak set. Reconciling sample-specific peak results post hoc from larger cohorts is particularly challenging, and informative spatial features specific to open chromatin signals are not leveraged effectively.

RESULTS

We propose a novel method, ROCCO, that determines consensus open chromatin regions across multiple samples simultaneously. ROCCO employs robust summary statistics and solves a constrained optimization problem formulated to account for both enrichment and spatial dependence of open chromatin signal data. We show this formulation admits attractive theoretical and conceptual properties as well as superior empirical performance compared to current methodology.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, and usage demos for ROCCO are available on GitHub at: https://github.com/nolan-h-hamilton/ROCCO. ROCCO can also be installed as a stand-alone binary utility using pip/PyPI.

Integrative approaches based on genomic techniques in the functional studies on enhancers

With the development of sequencing technology and the dramatic drop in sequencing cost, the functions of noncoding genes are being characterized in a wide variety of fields (e.g. biomedicine). Enhancers are noncoding DNA elements with vital transcription regulation functions. Tens of thousands of enhancers have been identified in the human genome; however, the location, function, target genes and regulatory mechanisms of most enhancers have not been elucidated thus far. As high-throughput sequencing techniques have leapt forwards, omics approaches have been extensively employed in enhancer research. Multidimensional genomic data integration enables the full exploration of the data and provides novel perspectives for screening, identification and characterization of the function and regulatory mechanisms of unknown enhancers. However, multidimensional genomic data are still difficult to integrate genome wide due to complex varieties, massive amounts, high rarity, etc. To facilitate the appropriate methods for studying enhancers with high efficacy, we delineate the principles, data processing modes and progress of various omics approaches to study enhancers and summarize the applications of traditional machine learning and deep learning in multi-omics integration in the enhancer field. In addition, the challenges encountered during the integration of multiple omics data are addressed. Overall, this review provides a comprehensive foundation for enhancer analysis.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

Current Status of Next-Generation Sequencing in Bone Genetic Diseases

The development of next-generation sequencing (NGS) has dramatically increased the speed and volume of genetic analysis. Furthermore, the range of applications of NGS is rapidly expanding to include genome, epigenome (such as DNA methylation), metagenome, and transcriptome analyses (such as RNA sequencing and single-cell RNA sequencing). NGS enables genetic research by offering various sequencing methods as well as combinations of methods. Bone tissue is the most important unit supporting the body and is a reservoir of calcium and phosphate ions, which are important for physical activity. Many genetic diseases affect bone tissues, possibly because metabolic mechanisms in bone tissue are complex. For instance, the presence of specialized immune cells called osteoclasts in the bone tissue, which absorb bone tissue and interact with osteoblasts in complex ways to support normal vital functions. Moreover, the many cell types in bones exhibit cell-specific proteins for their respective activities. Mutations in the genes encoding these proteins cause a variety of genetic disorders. The relationship between age-related bone tissue fragility (also called frailty) and genetic factors has recently attracted attention. Herein, we discuss the use of genomic, epigenomic, transcriptomic, and metagenomic analyses in bone genetic disorders.

Whole genome deconvolution unveils Alzheimer's resilient epigenetic signature

Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) accurately depicts the chromatin regulatory state and altered mechanisms guiding gene expression in disease. However, bulk sequencing entangles information from different cell types and obscures cellular heterogeneity. To address this, we developed Cellformer, a deep learning method that deconvolutes bulk ATAC-seq into cell type-specific expression across the whole genome. Cellformer enables cost-effective cell type-specific open chromatin profiling in large cohorts. Applied to 191 bulk samples from 3 brain regions, Cellformer identifies cell type-specific gene regulatory mechanisms involved in resilience to Alzheimer's disease, an uncommon group of cognitively healthy individuals that harbor a high pathological load of Alzheimer's disease. Cell type-resolved chromatin profiling unveils cell type-specific pathways and nominates potential epigenetic mediators underlying resilience that may illuminate therapeutic opportunities to limit the cognitive impact of the disease. Cellformer is freely available to facilitate future investigations using high-throughput bulk ATAC-seq data.

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.

Cell-type specific transcriptional adaptations of nucleus accumbens interneurons to amphetamine

Parvalbumin-expressing (PV+) interneurons of the nucleus accumbens (NAc) play an essential role in the addictive-like behaviors induced by psychostimulant exposure. To identify molecular mechanisms of PV+ neuron plasticity, we isolated interneuron nuclei from the NAc of male and female mice following acute or repeated exposure to amphetamine (AMPH) and sequenced for cell type-specific RNA expression and chromatin accessibility. AMPH regulated the transcription of hundreds of genes in PV+ interneurons, and this program was largely distinct from that regulated in other NAc GABAergic neurons. Chromatin accessibility at enhancers predicted cell-type specific gene regulation, identifying transcriptional mechanisms of differential AMPH responses. Finally, we assessed expression of PV-enriched, AMPH-regulated genes in an Mecp2 mutant mouse strain that shows heightened behavioral sensitivity to psychostimulants to explore the functional importance of this transcriptional program. Together these data provide novel insight into the cell-type specific programs of transcriptional plasticity in NAc neurons that underlie addictive-like behaviors.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

Translating non-coding genetic associations into a better understanding of immune-mediated disease

Genome-wide association studies have identified hundreds of genetic loci that are associated with immune-mediated diseases. Most disease-associated variants are non-coding, and a large proportion of these variants lie within enhancers. As a result, there is a pressing need to understand how common genetic variation might affect enhancer function and thereby contribute to immune-mediated (and other) diseases. In this Review, we first describe statistical and experimental methods to identify causal genetic variants that modulate gene expression, including statistical fine-mapping and massively parallel reporter assays. We then discuss approaches to characterise the mechanisms by which these variants modulate immune function, such as clustered regularly interspaced short palindromic repeats (CRISPR)-based screens. We highlight examples of studies that, by elucidating the effects of disease variants within enhancers, have provided important insights into immune function and uncovered key pathways of disease.

FAIRE-MS reveals mitotic retention of transcriptional regulators on a proteome-wide scale

Mitosis entails global and dramatic alterations, such as higher-order chromatin organization disruption, concomitant with global transcription downregulation. Cells reliably re-establishing gene expression patterns upon mitotic exit and maintaining cellular identities remain poorly understood. Previous studies indicated that certain transcription factors (TFs) remain associated with individual loci during mitosis and serve as mitotic bookmarkers. However, it is unclear which regulatory factors remain bound to the compacted mitotic chromosomes. We developed formaldehyde-assisted isolation of regulatory elements-coupled mass spectrometry (FAIRE-MS) that combines FAIRE-based open chromatin-associated protein pull-down and mass spectrometry (MS) to quantify the open chromatin-associated proteome during the interphase and mitosis. We identified 189 interphase and mitosis maintained (IM) regulatory factors using FAIRE-MS and found intrinsically disordered proteins and regions (IDP(R)s) are highly enriched, which plays a crucial role in liquid-liquid phase separation (LLPS) and chromatin organization during the cell cycle. Notably, in these IDP(R)s, we identified mitotic bookmarkers, such as CEBPB, HMGB1, and TFAP2A, and several factors, including MAX, HMGB3, hnRNP A2/B1, FUS, hnRNP D, and TIAL1, which are at least partially bound to the mitotic chromosome. Furthermore, it will be essential to study whether these IDP(R)s through LLPS helps cells transit from mitosis to the G1 phase during the cell cycle.

The RNA editing landscape in acute myeloid leukemia reveals associations with disease mutations and clinical outcome

Several studies have documented aberrant RNA editing patterns across multiple tumors across large patient cohorts from The Cancer Genome Atlas (TCGA). However, studies on understanding the role of RNA editing in acute myeloid leukemia (AML) have been limited to smaller sample sizes. Using high throughput transcriptomic data from the TCGA, we demonstrated higher levels of editing as a predictor of poor outcome within the AML patient samples. Moreover, differential editing patterns were observed across individual AML genotypes. We also could demonstrate a negative association between the degree of editing and mRNA abundance for some transcripts, identifying the potential regulatory potential of RNA-editing in altering gene expression in AML. Further edQTL analysis suggests potential cis-regulatory mechanisms in RNA editing variation. Our work suggests a functional and regulatory role of RNA editing in the pathogenesis of AML and we extended our analysis to gain insight into the factors influencing altered levels of editing.

Alternative transcription start sites contribute to acute-stress-induced transcriptome response in human skeletal muscle

Background

More than half of human protein-coding genes have an alternative transcription start site (TSS). We aimed to investigate the contribution of alternative TSSs to the acute-stress-induced transcriptome response in human tissue (skeletal muscle) using the cap analysis of gene expression approach. TSSs were examined at baseline and during recovery after acute stress (a cycling exercise).

Results

We identified 44,680 CAGE TSS clusters (including 3764 first defined) belonging to 12,268 genes and annotated for the first time 290 TSSs belonging to 163 genes. The transcriptome dynamically changes during the first hours after acute stress; the change in the expression of 10% of genes was associated with the activation of alternative TSSs, indicating differential TSSs usage. The majority of the alternative TSSs do not increase proteome complexity suggesting that the function of thousands of alternative TSSs is associated with the fine regulation of mRNA isoform expression from a gene due to the transcription factor-specific activation of various alternative TSSs. We identified individual muscle promoter regions for each TSS using muscle open chromatin data (ATAC-seq and DNase-seq). Then, using the positional weight matrix approach we predicted time course activation of “classic” transcription factors involved in response of skeletal muscle to contractile activity, as well as diversity of less/un-investigated factors.

Conclusions

Transcriptome response induced by acute stress related to activation of the alternative TSSs indicates that differential TSSs usage is an essential mechanism of fine regulation of gene response to stress stimulus. A comprehensive resource of accurate TSSs and individual promoter regions for each TSS in muscle was created. This resource together with the positional weight matrix approach can be used to accurate prediction of TFs in any gene(s) of interest involved in the response to various stimuli, interventions or pathological conditions in human skeletal muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00399-8.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Non-coding genetic variation in regulatory elements determines thrombosis and hemostasis phenotypes

Since the early inception of genome-wide association studies (GWAS), it became clear that, in all diseases or traits studied, most genetic variants are likely to exert their effect on gene expression mainly by altering the function of regulatory elements. At the same time, the regulation of the gene expression field broadened its boundaries, from the univocal relationship between regulatory elements and genes to include genome organization, long-range DNA interactions, and epigenetics. Next-generation sequencing has introduced genome-wide approaches that have greatly improved our understanding of the general principles of gene expression. However, elucidating how these apply in every single genomic locus still requires painstaking experimental work, in which several independent lines of evidence are required, and often this is helped by rare genetic variants in individuals with rare diseases. This review will focus on the non-coding features of the genome involved in transcriptional regulation, that when altered, leads to known cases of inherited (familial) thrombotic and hemostatic phenotypes, emphasizing the role of enhancers and super-enhancers.

Chromatin Dynamics in Digestive System Cancer: Commander and Regulator

Digestive system tumors have a poor prognosis due to complex anatomy, insidious onset, challenges in early diagnosis, and chemoresistance. Epidemiological statistics has verified that digestive system tumors rank first in tumor-related death. Although a great number of studies are devoted to the molecular biological mechanism, early diagnostic markers, and application of new targeted drugs in digestive system tumors, the therapeutic effect is still not satisfactory. Epigenomic alterations including histone modification and chromatin remodeling are present in human cancers and are now known to cooperate with genetic changes to drive the cancer phenotype. Chromatin is the carrier of genetic information and consists of DNA, histones, non-histone proteins, and a small amount of RNA. Chromatin and nucleosomes control the stability of the eukaryotic genome and regulate DNA processes such as transcription, replication, and repair. The dynamic structure of chromatin plays a key role in this regulatory function. Structural fluctuations expose internal DNA and thus provide access to the nuclear machinery. The dynamic changes are affected by various complexes and epigenetic modifications. Variation of chromatin dynamics produces early and superior regulation of the expression of related genes and downstream pathways, thereby controlling tumor development. Intervention at the chromatin level can change the process of cancer earlier and is a feasible option for future tumor diagnosis and treatment. In this review, we introduced chromatin dynamics including chromatin remodeling, histone modifications, and chromatin accessibility, and current research on chromatin regulation in digestive system tumors was also summarized.

DiMeLo-seq: a long-read, single-molecule method for mapping protein-DNA interactions genome wide

Studies of genome regulation routinely use high-throughput DNA sequencing approaches to determine where specific proteins interact with DNA, and they rely on DNA amplification and short-read sequencing, limiting their quantitative application in complex genomic regions. To address these limitations, we developed directed methylation with long-read sequencing (DiMeLo-seq), which uses antibody-tethered enzymes to methylate DNA near a target protein's binding sites in situ. These exogenous methylation marks are then detected simultaneously with endogenous CpG methylation on unamplified DNA using long-read, single-molecule sequencing technologies. We optimized and benchmarked DiMeLo-seq by mapping chromatin-binding proteins and histone modifications across the human genome. Furthermore, we identified where centromere protein A localizes within highly repetitive regions that were unmappable with short sequencing reads, and we estimated the density of centromere protein A molecules along single chromatin fibers. DiMeLo-seq is a versatile method that provides multimodal, genome-wide information for investigating protein-DNA interactions.

MethReg: estimating the regulatory potential of DNA methylation in gene transcription

Epigenome-wide association studies often detect many differentially methylated sites, and many are located in distal regulatory regions. To further prioritize these significant sites, there is a critical need to better understand the functional impact of CpG methylation. Recent studies demonstrated that CpG methylation-dependent transcriptional regulation is a widespread phenomenon. Here, we present MethReg, an R/Bioconductor package that analyzes matched DNA methylation and gene expression data, along with external transcription factor (TF) binding information, to evaluate, prioritize and annotate CpG sites with high regulatory potential. At these CpG sites, TF-target gene associations are often only present in a subset of samples with high (or low) methylation levels, so they can be missed by analyses that use all samples. Using colorectal cancer and Alzheimer's disease datasets, we show MethReg significantly enhances our understanding of the regulatory roles of DNA methylation in complex diseases.

Comprehensive epigenomic profiling of human alveolar epithelial differentiation identifies key epigenetic states and transcription factor co-regulatory networks for maintenance of distal lung identity

BACKGROUND

Disruption of alveolar epithelial cell (AEC) differentiation is implicated in distal lung diseases such as chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, and lung adenocarcinoma that impact morbidity and mortality worldwide. Elucidating underlying disease pathogenesis requires a mechanistic molecular understanding of AEC differentiation. Previous studies have focused on changes of individual transcription factors, and to date no study has comprehensively characterized the dynamic, global epigenomic alterations that facilitate this critical differentiation process in humans.

RESULTS

We comprehensively profiled the epigenomic states of human AECs during type 2 to type 1-like cell differentiation, including the methylome and chromatin functional domains, and integrated this with transcriptome-wide RNA expression data. Enhancer regions were drastically altered during AEC differentiation. Transcription factor binding analysis within enhancer regions revealed diverse interactive networks with enrichment for many transcription factors, including NKX2-1 and FOXA family members, as well as transcription factors with less well characterized roles in AEC differentiation, such as members of the MEF2, TEAD, and AP1 families. Additionally, associations among transcription factors changed during differentiation, implicating a complex network of heterotrimeric complex switching in driving differentiation. Integration of AEC enhancer states with the catalog of enhancer elements in the Roadmap Epigenomics Mapping Consortium and Encyclopedia of DNA Elements (ENCODE) revealed that AECs have similar epigenomic structures to other profiled epithelial cell types, including human mammary epithelial cells (HMECs), with NKX2-1 serving as a distinguishing feature of distal lung differentiation.

CONCLUSIONS

Enhancer regions are hotspots of epigenomic alteration that regulate AEC differentiation. Furthermore, the differentiation process is regulated by dynamic networks of transcription factors acting in concert, rather than individually. These findings provide a roadmap for understanding the relationship between disruption of the epigenetic state during AEC differentiation and development of lung diseases that may be therapeutically amenable.

Methods of massive parallel reporter assays for investigation of enhancers

The correct deployment of genetic programs for development and differentiation relies on finely coordinated regulation of specific gene sets. Genomic regulatory elements play an exceptional role in this process. There are few types of gene regulatory elements, including promoters, enhancers, insulators and silencers. Alterations of gene regulatory elements may cause various pathologies, including cancer, congenital disorders and autoimmune diseases. The development of high-throughput genomic assays has made it possible to significantly accelerate the accumulation of information about the characteristic epigenetic properties of regulatory elements. In combination with high-throughput studies focused on the genome-wide distribution of epigenetic marks, regulatory proteins and the spatial structure of chromatin, this significantly expands the understanding of the principles of epigenetic regulation of genes and allows potential regulatory elements to be searched for in silico. However, common experimental approaches used to study the local characteristics of chromatin have a number of technical limitations that may reduce the reliability of computational identification of genomic regulatory sequences. Taking into account the variability of the functions of epigenetic determinants and complex multicomponent regulation of genomic elements activity, their functional verification is often required. A plethora of methods have been developed to study the functional role of regulatory elements on the genome scale. Common experimental approaches for in silico identification of regulatory elements and their inherent technical limitations will be described. The present review is focused on original high-throughput methods of enhancer activity reporter analysis that are currently used to validate predicted regulatory elements and to perform de novo searches. The methods described allow assessing the functional role of the nucleotide sequence of a regulatory element, to determine its exact boundaries and to assess the influence of the local state of chromatin on the activity of enhancers and gene expression. These approaches have contributed substantially to the understanding of the fundamental principles of gene regulation.

Programmed genomic instability regulates neural transdifferentiation of human brain microvascular pericytes

BACKGROUND

Transdifferentiation describes transformation in vivo of specialized cells from one lineage into another. While there is extensive literature on forced induction of lineage reprogramming in vitro, endogenous mechanisms that govern transdifferentiation remain largely unknown. The observation that human microvascular pericytes transdifferentiate into neurons provided an opportunity to explore the endogenous molecular basis for lineage reprogramming.

RESULTS

We show that abrupt destabilization of the higher-order chromatin topology that chaperones lineage memory of pericytes is driven by transient global transcriptional arrest. This leads within minutes to localized decompression of the repressed competing higher-order chromatin topology and expression of pro-neural genes. Transition to neural lineage is completed by probabilistic induction of R-loops in key myogenic loci upon re-initiation of RNA polymerase activity, leading to depletion of the myogenic transcriptome and emergence of the neurogenic transcriptome.

CONCLUSIONS

These findings suggest that the global transcriptional landscape not only shapes the functional cellular identity of pericytes, but also stabilizes lineage memory by silencing the competing neural program within a repressed chromatin state.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

A systematic genome-wide mapping of oncogenic mutation selection during CRISPR-Cas9 genome editing

Recent studies have reported that genome editing by CRISPR–Cas9 induces a DNA damage response mediated by p53 in primary cells hampering their growth. This could lead to a selection of cells with pre-existing p53 mutations. In this study, employing an integrated computational and experimental framework, we systematically investigated the possibility of selection of additional cancer driver mutations during CRISPR-Cas9 gene editing. We first confirm the previous findings of the selection for pre-existing p53 mutations by CRISPR-Cas9. We next demonstrate that similar to p53, wildtype KRAS may also hamper the growth of Cas9-edited cells, potentially conferring a selective advantage to pre-existing KRAS-mutant cells. These selective effects are widespread, extending across cell-types and methods of CRISPR-Cas9 delivery and the strength of selection depends on the sgRNA sequence and the gene being edited. The selection for pre-existing p53 or KRAS mutations may confound CRISPR-Cas9 screens in cancer cells and more importantly, calls for monitoring patients undergoing CRISPR-Cas9-based editing for clinical therapeutics for pre-existing p53 and KRAS mutations.

CRISPR-Cas9 gene editing can induce a p53 mediated damage response. Here the authors investigate the possibility of selection of pre-existing cancer driver mutations during CRISPR-Cas9 knockout based gene editing and identify KRAS mutants that may confer a selected advantage to edited cells.

A mouse-specific retrotransposon drives a conserved Cdk2ap1 isoform essential for development

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

Decoding the Equine Genome: Lessons from ENCODE

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

Tree-Based Co-Clustering Identifies Chromatin Accessibility Patterns Associated With Hematopoietic Lineage Structure

Chromatin accessibility, as measured by ATACseq, varies between hematopoietic cell types in different lineages of the hematopoietic differentiation tree, e.g. T cells vs. B cells, but methods that associate variation in chromatin accessibility to the lineage structure of the differentiation tree are lacking. Using an ATACseq dataset recently published by the ImmGen consortium, we construct associations between chromatin accessibility and hematopoietic cell types using a novel co-clustering approach that accounts for the structure of the hematopoietic, differentiation tree. Under a model in which all loci and cell types within a co-cluster have a shared accessibility state, we show that roughly 80% of cell type associated accessibility variation can be captured through 12 cell type clusters and 20 genomic locus clusters, with the cell type clusters reflecting coherent components of the differentiation tree. Using publicly available ChIPseq datasets, we show that our clustering reflects transcription factor binding patterns with implications for regulation across cell types. We show that traditional methods such as hierarchical and kmeans clusterings lead to cell type clusters that are more dispersed on the tree than our tree-based algorithm. We provide a python package, chromcocluster, that implements the algorithms presented.

Brain regionalization by Polycomb-group proteins and chromatin accessibility

During brain development, neural precursor cells (NPCs) in different brain regions produce different types of neurons, and each of these regions plays a different role in the adult brain. Therefore, precise regionalization is essential in the early stages of brain development, and irregular regionalization has been proposed as the cause of neurodevelopmental disorders. The mechanisms underlying brain regionalization have been well studied in terms of morphogen-induced expression of critical transcription factors for regionalization. NPC potential in different brain regions is defined by chromatin structures that regulate the plasticity of gene expression. Herein, we present recent findings on the importance of chromatin structure in brain regionalization, particularly with respect to its regulation by Polycomb-group proteins and chromatin accessibility.

Sites of active gene regulation in the prenatal frontal cortex and their role in neuropsychiatric disorders

Common genetic variation appears to largely influence risk for neuropsychiatric disorders through effects on gene regulation. It is therefore possible to shed light on the biology of these conditions by testing for enrichment of associated genetic variation within regulatory genomic regions operating in specific tissues or cell types. Here, we have used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) to map open chromatin (an index of active regulatory genomic regions) in bulk tissue, NeuN+ and NeuN- nuclei from the prenatal human frontal cortex, and tested enrichment of single-nucleotide polymorphism (SNP) heritability for five neuropsychiatric disorders (autism spectrum disorder, attention deficit hyperactivity disorder [ADHD], bipolar disorder, major depressive disorder, and schizophrenia) within these regions. We observed significant enrichment of SNP heritability for ADHD, major depressive disorder, and schizophrenia within open chromatin regions (OCRs) mapped in bulk fetal frontal cortex, and for all five tested neuropsychiatric conditions when we restricted these sites to those overlapping histone modifications indicative of enhancers (H3K4me1) or promoters (H3K4me3) in fetal brain. SNP heritability for neuropsychiatric disorders was significantly enriched in OCRs identified in fetal frontal cortex NeuN- as well as NeuN+ nuclei overlapping fetal brain H3K4me1 or H3K4me3 sites. We additionally demonstrate the utility of our mapped OCRs for prioritizing potentially functional SNPs at genome-wide significant risk loci for neuropsychiatric disorders. Our data provide evidence for an early neurodevelopmental component to a range of neuropsychiatric conditions and highlight an important role for regulatory genomic regions active within both NeuN+ and NeuN- cells of the prenatal brain.

Sex dependent glial-specific changes in the chromatin accessibility landscape in late-onset Alzheimer's disease brains

BACKGROUND

In the post-GWAS era, there is an unmet need to decode the underpinning genetic etiologies of late-onset Alzheimer's disease (LOAD) and translate the associations to causation.

METHODS

We conducted ATAC-seq profiling using NeuN sorted-nuclei from 40 frozen brain tissues to determine LOAD-specific changes in chromatin accessibility landscape in a cell-type specific manner.

RESULTS

We identified 211 LOAD-specific differential chromatin accessibility sites in neuronal-nuclei, four of which overlapped with LOAD-GWAS regions (±100 kb of SNP). While the non-neuronal nuclei did not show LOAD-specific differences, stratification by sex identified 842 LOAD-specific chromatin accessibility sites in females. Seven of these sex-dependent sites in the non-neuronal samples overlapped LOAD-GWAS regions including APOE. LOAD loci were functionally validated using single-nuclei RNA-seq datasets.

CONCLUSIONS

Using brain sorted-nuclei enabled the identification of sex-dependent cell type-specific LOAD alterations in chromatin structure. These findings enhance the interpretation of LOAD-GWAS discoveries, provide potential pathomechanisms, and suggest novel LOAD-loci.

Interferon stimulated binding of ISRE is cell type specific and is predicted by homeostatic chromatin state

•

IFN stimulated binding of ISRE by ISGF3 is cell specific, particularly for ISRE in enhancer regions.

•

IFN stimulated binding of ISRE in enhancer regions associates with differential expression.

•

The homeostatic, chromatin state of an ISRE is predictive of IFN stimulated binding.

The type I interferon (IFN) signaling pathway involves binding of the transcription factor ISGF3 to IFN-stimulated response elements, ISREs. Gene expression under IFN stimulation is known to vary across cell types, but variation in ISGF3 binding to ISRE across cell types has not been characterized. We examined ISRE binding patterns under IFN stimulation across six cell types using existing ChIPseq datasets. We find that ISRE binding is largely cell specific for ISREs distal to transcription start sites (TSS) and largely conserved across cell types for ISREs proximal to TSS. We show that bound ISRE distal to TSS associate with differential expression of ISGs, although at weaker levels than bound ISRE proximal to TSS. Using existing ATACseq and ChIPseq datasets, we show that the chromatin state of ISRE at homeostasis is cell type specific and is predictive of cell specific, ISRE binding under IFN stimulation. Our results support a model in which the chromatin state of ISRE in enhancer elements is modulated in a cell type specific manner at homeostasis, leading to cell type specific differences in ISRE binding patterns under IFN stimulation.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Genome-Wide Identification of DNase I Hypersensitive Sites in Plants

DNase I hypersensitive site (DHS) mapping combined with high-throughput sequencing (DNase-seq) enables the identification of cis-regulatory DNA elements (CREs) genome wide. However, despite the wide applications of DNase-seq in plants, its application to the highly repetitive genomes of plants has lagged. Here, we describe a modified DNase-seq method, making it more practical for application to plants with genomes enriched with repetitive DNA. This approach adopts a double-hit-based strategy, in which small (<250-bp) DNA fragments digested by DNase I are selected and used for sequencing library construction. Using these protocols, we have conducted DNase-seq in plants with high content of repetitive DNA, including maize, sugarcane, and tetraploid cotton. Genome-wide maps of DHS and CREs have been created using these DNase-seq datasets. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Nuclei isolation Basic Protocol 2: DNase I digestion Basic Protocol 3: Target DNA isolation Basic Protocol 4: Library construction and validation.

Undisclosed, unmet and neglected challenges in multi-omics studies

Multi-omics approaches have become a reality in both large genomics projects and small laboratories. However, the multi-omics research community still faces a number of issues that have either not been sufficiently discussed or for which current solutions are still limited. In this Perspective, we elaborate on these limitations and suggest points of attention for future research. We finally discuss new opportunities and challenges brought to the field by the rapid development of single-cell high-throughput molecular technologies.

Genome-wide R-loop Landscapes during Cell Differentiation and Reprogramming

DNA:RNA hybrids play key roles in both physiological and disease states by regulating chromatin and genome organization. Their homeostasis during cell differentiation and cell plasticity remains elusive. Using an isogenic human stem cell platform, we systematically characterize R-loops, DNA methylation, histone modifications, and chromatin accessibility in pluripotent cells and their lineage-differentiated derivatives. We confirm that a portion of R-loops formed co-transcriptionally at pluripotency genes in pluripotent stem cells and at lineage-controlling genes in differentiated lineages. Notably, a subset of R-loops maintained after differentiation are associated with repressive chromatin marks on silent pluripotency genes and undesired lineage genes. Moreover, in reprogrammed pluripotent cells, cell-of-origin-specific R-loops are initially present but are resolved with serial passaging. Our analysis suggests a multifaceted role of R-loops in cell fate determination that may serve as an additional layer of modulation on cell fate memory and cell plasticity.

Prediction of disease-associated functional variants in noncoding regions through a comprehensive analysis by integrating datasets and features

One of the greatest challenges in human genetics is deciphering the link between functional variants in noncoding sequences and the pathophysiology of complex diseases. To address this issue, many methods have been developed to sort functional single-nucleotide variants (SNVs) for neutral SNVs in noncoding regions. In this study, we integrated well-established features and commonly used datasets and merged them into large-scale datasets based on a random forest model, which yielded promising performance and outperformed some cutting-edge approaches. Our analyses of feature importance and data coverage also provide certain clues for future research in enhancing the prediction of functional noncoding SNVs.