Advances of DNase-seq for mapping active gene regulatory elements across the genome in animals

Massively parallel characterization of transcriptional regulatory elements

The human genome contains millions of candidate cis-regulatory elements (cCREs) with cell-type-specific activities that shape both health and many disease states1. However, we lack a functional understanding of the sequence features that control the activity and cell-type-specific features of these cCREs. Here we used lentivirus-based massively parallel reporter assays (lentiMPRAs) to test the regulatory activity of more than 680,000 sequences, representing an extensive set of annotated cCREs among three cell types (HepG2, K562 and WTC11), and found that 41.7% of these sequences were active. By testing sequences in both orientations, we find promoters to have strand-orientation biases and their 200-nucleotide cores to function as non-cell-type-specific 'on switches' that provide similar expression levels to their associated gene. By contrast, enhancers have weaker orientation biases, but increased tissue-specific characteristics. Utilizing our lentiMPRA data, we develop sequence-based models to predict cCRE function and variant effects with high accuracy, delineate regulatory motifs and model their combinatorial effects. Testing a lentiMPRA library encompassing 60,000 cCREs in all three cell types further identified factors that determine cell-type specificity. Collectively, our work provides an extensive catalogue of functional CREs in three widely used cell lines and showcases how large-scale functional measurements can be used to dissect regulatory grammar.

Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair

Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.

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.

Microsatellites used in forensics are in regions enriched for trait-associated variants

The 20 short tandem repeat (STR) loci of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS loci are thought to contain little information about ancestry or traits. However, in the past 20 years, a growing field has identified hundreds of thousands of genotype-trait associations. Here, we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. Although this study cannot establish or quantify associations between CODIS genotypes and phenotypes, we find that the regions around the CODIS loci are enriched for both known pathogenic variants (> 90th percentile) and for trait-associated SNPs identified in genome-wide association studies (GWAS) (≥ 95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs.

Exploring the Role of Enhancer-Mediated Transcriptional Regulation in Precision Biology

The emergence of precision biology has been driven by the development of advanced technologies and techniques in high-resolution biological research systems. Enhancer-mediated transcriptional regulation, a complex network of gene expression and regulation in eukaryotes, has attracted significant attention as a promising avenue for investigating the underlying mechanisms of biological processes and diseases. To address biological problems with precision, large amounts of data, functional information, and research on the mechanisms of action of biological molecules is required to address biological problems with precision. Enhancers, including typical enhancers and super enhancers, play a crucial role in gene expression and regulation within this network. The identification and targeting of disease-associated enhancers hold the potential to advance precision medicine. In this review, we present the concepts, progress, importance, and challenges in precision biology, transcription regulation, and enhancers. Furthermore, we propose a model of transcriptional regulation for multi-enhancers and provide examples of their mechanisms in mammalian cells, thereby enhancing our understanding of how enhancers achieve precise regulation of gene expression in life processes. Precision biology holds promise in providing new tools and platforms for discovering insights into gene expression and disease occurrence, ultimately benefiting individuals and society as a whole.

Gene expression in organoids: an expanding horizon

Recent development of human three-dimensional organoid cultures has opened new doors and opportunities ranging from modelling human development in vitro to personalised cancer therapies. These new in vitro systems are opening new horizons to the classic understanding of human development and disease. However, the complexity and heterogeneity of these models requires cutting-edge techniques to capture and trace global changes in gene expression to enable identification of key players and uncover the underlying molecular mechanisms. Rapid development of sequencing approaches made possible global transcriptome analyses and epigenetic profiling. Despite challenges in organoid culture and handling, these techniques are now being adapted to embrace organoids derived from a wide range of human tissues. Here, we review current state-of-the-art multi-omics technologies, such as single-cell transcriptomics and chromatin accessibility assays, employed to study organoids as a model for development and a platform for precision medicine.

Microsatellites used in forensics are located in regions unusually rich in trait-associated variants

The 20 short tandem repeat (STR) markers of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS markers are thought to contain information relevant to identification only (such as a human fingerprint would), with little information about ancestry or traits. However, in the past 20 years, a quickly growing field has identified hundreds of thousands of genotype-trait associations. Here we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. We find that the regions around the CODIS markers are enriched for both known pathogenic variants (>90th percentile) and for SNPs identified as trait-associated in genome-wide association studies (GWAS) (≥95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs. Although it is not obvious how much phenotypic information CODIS would need to convey to strain the "DNA fingerprint" analogy, the CODIS markers, considered as a set, are in regions unusually dense with variants with known phenotypic associations.

CFA: An explainable deep learning model for annotating the transcriptional roles of cis-regulatory modules based on epigenetic codes

Metazoa gene expression is controlled by modular DNA segments called cis-regulatory modules (CRMs). CRMs can convey promoter/enhancer/insulator roles, generating additional regulation layers in transcription. Experiments for understanding CRM roles are low-throughput and costly. Large-scale CRM function investigation still depends on computational methods. However, existing in silico tools only recognize enhancers or promoters exclusively, thus accumulating errors when considering CRM promoter/enhancer/insulator roles altogether. Currently, no algorithm can concurrently consider these CRM roles. In this research, we developed the CRM Function Annotator (CFA) model. CFA provides complete CRM transcriptional role labeling based on epigenetic profiling interpretation. We demonstrated that CFA achieves high performance (test macro auROC/auPRC = 94.1%/90.3%) and outperforms existing tools in promoter/enhancer/insulator identification. CFA is also inspected to recognize explainable epigenetic codes consistent with previous findings when labeling CRM roles. By considering the higher-order combinations of the epigenetic codes, CFA significantly reduces false-positive rates in CRM transcriptional role annotation. CFA is available at https://github.com/cobisLab/CFA/.

regCNN: identifying Drosophila genome-wide cis-regulatory modules via integrating the local patterns in epigenetic marks and transcription factor binding motifs

Transcription regulation in metazoa is controlled by the binding events of transcription factors (TFs) or regulatory proteins on specific modular DNA regulatory sequences called cis-regulatory modules (CRMs). Understanding the distributions of CRMs on a genomic scale is essential for constructing the metazoan transcriptional regulatory networks that help diagnose genetic disorders. While traditional reporter-assay CRM identification approaches can provide an in-depth understanding of functions of some CRM, these methods are usually cost-inefficient and low-throughput. It is generally believed that by integrating diverse genomic data, reliable CRM predictions can be made. Hence, researchers often first resort to computational algorithms for genome-wide CRM screening before specific experiments. However, current existing in silico methods for searching potential CRMs were restricted by low sensitivity, poor prediction accuracy, or high computation time from TFBS composition combinatorial complexity. To overcome these obstacles, we designed a novel CRM identification pipeline called regCNN by considering the base-by-base local patterns in TF binding motifs and epigenetic profiles. On the test set, regCNN shows an accuracy/auROC of 84.5%/92.5% in CRM identification. And by further considering local patterns in epigenetic profiles and TF binding motifs, it can accomplish 4.7% (92.5%–87.8%) improvement in the auROC value over the average value-based pure multi-layer perceptron model. We also demonstrated that regCNN outperforms all currently available tools by at least 11.3% in auROC values. Finally, regCNN is verified to be robust against its resizing window hyperparameter in dealing with the variable lengths of CRMs. The model of regCNN can be downloaded athttp://cobisHSS0.im.nuk.edu.tw/regCNN/.

Epitome: predicting epigenetic events in novel cell types with multi-cell deep ensemble learning

The accumulation of large epigenomics data consortiums provides us with the opportunity to extrapolate existing knowledge to new cell types and conditions. We propose Epitome, a deep neural network that learns similarities of chromatin accessibility between well characterized reference cell types and a query cellular context, and copies over signal of transcription factor binding and modification of histones from reference cell types when chromatin profiles are similar to the query. Epitome achieves state-of-the-art accuracy when predicting transcription factor binding sites on novel cellular contexts and can further improve predictions as more epigenetic signals are collected from both reference cell types and the query cellular context of interest.

SWI/SNF (BAF) complexes: From framework to a functional role in endothelial mechanotransduction

Endothelial cells (ECs) are constantly subjected to an array of mechanical cues, especially shear stress, due to their luminal placement in the blood vessels. Blood flow can regulate various aspects of endothelial biology and pathophysiology by regulating the endothelial processes at the transcriptomic, proteomic, miRNomic, metabolomics, and epigenomic levels. ECs sense, respond, and adapt to altered blood flow patterns and shear profiles by specialized mechanisms of mechanosensing and mechanotransduction, resulting in qualitative and quantitative differences in their gene expression. Chromatin-regulatory proteins can regulate transcriptional activation by modifying the organization of nucleosomes at promoters, enhancers, silencers, insulators, and locus control regions. Recent research efforts have illustrated that SWI/SNF (SWItch/Sucrose Non-Fermentable) or BRG1/BRM-associated factor (BAF) complex regulates DNA accessibility and chromatin structure. Since the discovery, the gene-regulatory mechanisms of the BAF complex associated with chromatin remodeling have been intensively studied to investigate its role in diverse disease phenotypes. Thus far, it is evident that (1) the SWI/SNF complex broadly regulates the activity of transcriptional enhancers to control lineage-specific differentiation and (2) mutations in the BAF complex proteins lead to developmental disorders and cancers. It is unclear if blood flow can modulate the activity of SWI/SNF complex to regulate EC differentiation and reprogramming. This review emphasizes the integrative role of SWI/SNF complex from a structural and functional standpoint with a special reference to cardiovascular diseases (CVDs). The review also highlights how regulation of this complex by blood flow can lead to the discovery of new therapeutic interventions for the treatment of endothelial dysfunction in vascular diseases.

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.

Tissue and sex-specific programming of DNA methylation by perinatal lead exposure: implications for environmental epigenetics studies

Early developmental environment can influence long-term health through reprogramming of the epigenome. Human environmental epigenetics studies rely on surrogate tissues, such as blood, to assess the effects of environment on disease-relevant but inaccessible target tissues. However, the extent to which environment-induced epigenetic changes are conserved between these tissues is unclear. A better understanding of this conservation is imperative for effective design and interpretation of human environmental epigenetics studies. The Toxicant Exposures and Responses by Genomic and Epigenomic Regulators of Transcription (TaRGET II) consortium was established by the National Institute of Environmental Health Sciences to address the utility of surrogate tissues as proxies for toxicant-induced epigenetic changes in target tissues. We and others have recently reported that perinatal exposure to lead (Pb) is associated with adverse metabolic outcomes. Here, we investigated the sex-specific effects of perinatal exposure to a human environmentally relevant level of Pb on DNA methylation in paired liver and blood samples from adult mice using enhanced reduced-representation bisulphite sequencing. Although Pb exposure ceased at 3 weeks of age, we observed thousands of sex-specific differentially methylated cytosines in the blood and liver of Pb-exposed animals at 5 months of age, including 44 genomically imprinted loci. We observed significant tissue overlap in the genes mapping to differentially methylated cytosines. A small but significant subset of Pb-altered genes exhibit basal sex differences in gene expression in the mouse liver. Collectively, these data identify potential molecular targets for Pb-induced metabolic diseases, and inform the design of more robust human environmental epigenomics studies.

Retinoic Acid-Regulated Target Genes During Development: Integrative Genomics Analysis

Retinoic acid (RA), a major natural active metabolite of vitamin A (VA) is well known to play critical roles in embryonic development. The effects of RA are mediated by nuclear receptors (RARs), which regulate the expression of gene batteries involved in cell growth and differentiation. Since the early 1990s several laboratories have focused on understanding how RA-regulated genes and RAR binding sites operate by studying the differentiation of embryonal carcinoma cells and embryonic stem cells. The development of hybridization-based microarray technology and high performance software analysis programs has allowed the characterization of thousands of RA-regulated genes. During the two last decades, publication of the genome sequence of various organisms has allowed advances in massive parallel sequencing and bioinformatics analysis of genome-wide data sets. These new generation sequencing (NGS) technologies have revolutionized the field by providing a global integrated picture of RA-regulated gene networks and the regulatory programs involved in cell fate decisions during embryonal carcinoma and embryonic stem cells differentiation. Now the challenge is to reconstruct the RA-regulated gene networks at the single cell level during the development of specialized embryonic tissues.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

Detect accessible chromatin using ATAC-sequencing, from principle to applications

Background

Chromatin accessibility is crucial for gene expression regulation in specific cells and in multiple biological processes. Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) is an effective way to reveal chromatin accessibility at a genome-wide level. Through ATAC-seq, produced reads from a small number of cells reflect accessible regions that correspond to nucleosome positioning and transcription factor binding sites, due to probing hyperactive Tn5 transposase to DNA sequence.

Conclusion

In this review, we summarize both principle and features of ATAC-seq, highlight its applications in basic and clinical research. ATAC-seq has generated comprehensive chromatin accessible maps, and is becoming a powerful tool to understand dynamic gene expression regulation in stem cells, early embryos and tumors.

The chromatin remodeling protein Lsh alters nucleosome occupancy at putative enhancers and modulates binding of lineage specific transcription factors

Dynamic regulation of chromatin accessibility is a key feature of cellular differentiation during embryogenesis, but the precise factors that control access to chromatin remain largely unknown. Lsh/HELLS is critical for normal development and mutations of Lsh in human cause the ICF (Immune deficiency, Centromeric instability, Facial anomalies) syndrome, a severe immune disorder with multiple organ deficiencies. We report here that Lsh, previously known to regulate DNA methylation level, has a genome wide chromatin remodeling function. Using micrococcal nuclease (MNase)-seq analysis, we demonstrate that Lsh protects MNase accessibility at transcriptional regulatory regions characterized by DNase I hypersensitivity and certain histone 3 (H3) tail modifications associated with enhancers. Using an auxin-inducible degron system, allowing proteolytical degradation of Lsh, we show that Lsh mediated changes in nucleosome occupancy are independent of DNA methylation level and are characterized by reduced H3 occupancy. While Lsh mediated nucleosome occupancy prevents binding sites for transcription factors in wild type cells, depletion of Lsh leads to an increase in binding of ectopically expressed tissue specific transcription factors to their respective binding sites. Our data suggests that Lsh mediated chromatin remodeling can modulate nucleosome positioning at a subset of putative enhancers contributing to the preservation of cellular identity through regulation of accessibility.