From reads to insight: a hitchhiker's guide to ATAC-seq data analysis

Sample Preparation and Integrative Data Analysis of a Droplet-based Single-Cell ATAC-sequencing Using Murine Thymic Epithelial Cells

Accessible chromatin regions modulate gene expression by acting as cis-regulatory elements. Understanding the epigenetic landscape by mapping accessible regions of DNA is therefore imperative to decipher mechanisms of gene regulation under specific biological contexts of interest. The assay for transposase-accessible chromatin sequencing (ATAC-seq) has been widely used to detect accessible chromatin and the recent introduction of single-cell technology has increased resolution to the single-cell level. In a recent study, we used droplet-based, single-cell ATAC-seq technology (scATAC-seq) to reveal the epigenetic profile of the transit-amplifying subset of thymic epithelial cells (TECs), which was identified previously using single-cell RNA-sequencing technology (scRNA-seq). This protocol allows the preparation of nuclei from TECs in order to perform droplet-based scATAC-seq and its integrative analysis with scRNA-seq data obtained from the same cell population. Integrative analysis has the advantage of identifying cell types in scATAC-seq data based on cell cluster annotations in scRNA-seq analysis.

In Situ Hybridization as a Method to Examine Gene Regulatory Activity In Vivo

Transcription factor-enhancer binding events are among the most well-studied protein-DNA interactions, allowing researchers to determine mechanisms of transcriptional activation or repression during development. While large-scale ChIP-sequence datasets, together with computational predictions and chromatin accessibility data, yield information on potential transcription factor binding activities, reporter gene assays provide measurable information on whether these binding activities are functional in particular cell types during development. Here, we present a detailed protocol to examine enhancer activity in Drosophila embryos using cloning, transgenesis, and in situ hybridization.

Progress in molecular biology and translational science: Epigenetics in cardiovascular health and disease

Conrad Waddington's epigenetics landscape has provided a metaphorical framework for how cells progress from undifferentiated states to one of several discrete, distinct, differentiated cell fates. The understanding of epigenetics has evolved over time, with DNA methylation being the most studied epigenetic modification, followed by histone modifications and non-coding RNA. Cardiovascular diseases (CVD) are leading contributors to death worldwide, with the prevalence of CVDs increasing across the last couple of decades. Significant amount of resources being poured into researching key mechanisms and underpinnings of the various CVDs. These molecular studies looked at the genetics, epigenetics as well as the transcriptomics of various cardiovascular conditions, aiming to provide mechanistic insights. It has paved the way for therapeutics to be developed and in recent years, epi-drugs for the treatment of CVDs. This chapter aims to cover the various roles of epigenetics in the context of cardiovascular health and disease. The following will be examined in detail: the developments in basic experimental techniques used to study epigenetics, the role of epigenetics in various CVDs (hypertension, atrial fibrillation, atherosclerosis, and heart failure), and current advances in epi-therapeutics, providing a holistic view of the current concerted efforts in advancing the field of epigenetics in CVDs.

Epigenomic Profiling of Dorsal Root Ganglia upon Regenerative and Non-regenerative Axonal Injury

RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and assay for transposase-accessible chromatin sequencing (ATAC-seq) are genome-wide techniques that provide information relative to gene expression, chromatin binding sites, and chromatin accessibility, respectively. Here we describe RNA-seq, H3K9ac, H3K27ac and H3K27me3 ChIP-seq, and ATAC-seq in dorsal root ganglia (DRG) after sciatic nerve or dorsal column axotomy, to characterize the transcriptional and epigenetic signatures of DRG upon regenerative vs non-regenerative axonal lesion.

GAGAM v1.2: An Improvement on Peak Labeling and Genomic Annotated Gene Activity Matrix Construction

Single-cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) is rapidly becoming a powerful technology for assessing the epigenetic landscape of thousands of cells. However, the sparsity of the resulting data poses significant challenges to their interpretability and informativeness. Different computational methods are available, proposing ways to generate significant features from accessibility data and process them to obtain meaningful results. Foremost among them is the peak calling, which interprets the raw scATAC-seq data generating the peaks as features. However, scATAC-seq data are not trivially comparable with single-cell RNA sequencing (scRNA-seq) data, an increasingly pressing challenge since the necessity of multimodal experiments integration. For this reason, this study wants to improve the concept of the Gene Activity Matrix (GAM), which links the accessibility data to the genes, by proposing an improved version of the Genomic-Annotated Gene Activity Matrix (GAGAM) concept. Specifically, this paper presents GAGAM v1.2, a new and better version of GAGAM v1.0. GAGAM aims to label the peaks and link them to the genes through functional annotation of the whole genome. Using genes as features in scATAC-seq datasets makes different datasets comparable and allows linking gene accessibility and expression. This link is crucial for gene regulation understanding and fundamental for the increasing impact of multi-omics data. Results confirm that our method performs better than the previous GAMs and shows a preliminary comparison with scRNA-seq data.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Integrated chromatin accessibility and DNA methylation analysis to reveal the critical epigenetic modification and regulatory mechanism in gonadal differentiation of the sequentially hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic and economically farmed fish that undergoes sex reversal from ovary to testis via ovotestis during gonadal development. The epigenetic changes that are associated with gonadal development in this species remain unclear.

METHODS

We produced DNA methylome, transcriptome, and chromatin accessibility maps of the key stages of gonad development: ovary, ovotestis, and testis. The expression of the key candidate genes was detected using qRT-PCR and in situ hybridization and the methylation levels were analysed using bisulphite sequencing PCR. Promoter activity and regulation were assessed using dual-luciferase reporter assays.

RESULTS

Gonadal development exhibits highly dynamic transcriptomic, DNA methylation, and chromatin accessibility changes. We found that DNA methylation status, especially of the transcription start site, was significantly negatively correlated with gene expression while chromatin accessibility exhibited no correlation with gene expression during gonadal development. The epigenetic signatures revealed many novel regulatory elements and genes involved in sex reversal, which were validated. DNA methylation detection and site mutation of plastin-2 promoter, as a candidate gene, revealed that DNA methylation could impact the binding of transcription factor dmrt1 and foxl2 through methylation and demethylation to regulate plastin-2 expression during gonadal development.

CONCLUSIONS

These data provide novel insights into epigenetic modification and help elucidate the potential molecular mechanism by which dynamic modification of DNA methylation plays a crucial role in gonadal development.

Using "Galaxy-rCASC": A Public Galaxy Instance for Single-Cell RNA-Seq Data Analysis

rCASC is a modular workflow providing an integrated environment for single-cell RNA-seq (scRNA-Seq) data analysis exploiting Docker containers to achieve functional and computational reproducibility. It was initially developed as an R package usable also through a Java GUI. However, the Java frontend cannot be employed when running rCASC on a remote server, a typical setup due to the significant computational resources commonly needed to analyze scRNA-Seq data.To allow the use of rCASC through a graphical user interface on the client side and to harness the many advantages provided by the Galaxy platform, we have made rCASC available as a Galaxy set of tools, also providing a dedicated public instance of Galaxy named "Galaxy-rCASC." To integrate rCASC into Galaxy, all its functions, originally implemented as a set of Docker containers to maximize reproducibility, have been extensively reworked to become independent from the R package functions that launch them in the original implementation. Furthermore, suitable Galaxy wrappers have been developed for most functions of rCASC. We provide a detailed reference document to the use of Galaxy-rCASC with insights and explanations on the platform functionalities, parameters, and output while guiding the reader through the typical rCASC analysis workflow of a scRNA-Seq dataset.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

The early function of cortisol in liver during Aeromonas hydrophila infection: Dynamics of the transcriptome and accessible chromatin landscapes

In China, channel catfish (Ictalurus punctatus) is an important aquaculture species; however, haemorrhagic disease (Aeromonas hydrophila induced disease) in these fish has caused tremendous economic loss due to high morbidity and mass mortality in the breeding industry. The role of cortisol in bacterial diseases, particularly in the acute phase, remains unclear. In this study, liver transcriptome (RNA-seq) and chromatin accessibility (ATAC-seq) analyses were employed to investigate the early functional role of cortisol in Aeromonas hydrophila-stimulated responses. Our experiments confirmed that A. hydrophila infection can initially significantly increase serum cortisol levels at 1 h after infection. At this time point, the increased serum cortisol levels can significantly regulate A. hydrophila-regulated genes by affecting both transcriptome and chromatin accessibility. Cross-analysis of RNA-seq and ATAC-seq revealed that a certain gene group (92 target_DEGs) was regulated at an early time point by cortisol. KEGG enrichment analysis revealed that the top three pathways according to target_DEGs were cancer, glutathione metabolism, and the Notch signalling pathway. The protein-protein interaction analysis of target_DEGs revealed that they may be primarily involved in cell proliferation, CD8⁺ T cell function, glutathione synthesis, and activation of the NF-κB signalling pathway. This suggests that after the emergence of immune stress, the early regulation of cortisol is positive against the immune response. It is possible that in this situation, the animal is attempting to avoid dangerous situations and risks and then cope with the imbalance produced by the stressor to ultimately restore homeostasis. Our results will contribute to future research on fish and provide valuable insight regarding the mechanism of immune regulation by cortisol and the study of bacterial haemorrhagic disease in channel catfish.

Dynamic chromatin accessibility deploys heterotypic cis/trans-acting factors driving stomatal cell-fate commitment

Chromatin architecture and transcription factor (TF) binding underpin cell-fate specification during development, but their mutual regulatory relationships remain unclear. Here we report an atlas of dynamic chromatin landscapes during stomatal cell-lineage progression, in which sequential cell-state transitions are governed by lineage-specific bHLH TFs. Major reprogramming of chromatin accessibility occurs at the proliferation-to-differentiation transition. We discover novel co-cis regulatory elements (CREs) signifying the early precursor stage, BBR/BPC (GAGA) and bHLH (E-box) motifs, where master-regulatory bHLH TFs, SPEECHLESS and MUTE, consecutively bind to initiate and terminate the proliferative state, respectively. BPC TFs complex with MUTE to repress SPEECHLESS expression through a local deposition of repressive histone marks. We elucidate the mechanism by which cell-state-specific heterotypic TF complexes facilitate cell-fate commitment by recruiting chromatin modifiers via key co-CREs.

Cross-regulome profiling of RNA polymerases highlights the regulatory role of polymerase III on mRNA transcription by maintaining local chromatin architecture

BACKGROUND

Mammalian cells have three types of RNA polymerases (Pols), Pol I, II, and III. However, the extent to which these polymerases are cross-regulated and the underlying mechanisms remain unclear.

RESULTS

We employ genome-wide profiling after acute depletion of Pol I, Pol II, or Pol III to assess cross-regulatory effects between these Pols. We find that these enzymes mainly affect the transcription of their own target genes, while certain genes are transcribed by the other polymerases. Importantly, the most active type of crosstalk is exemplified by the fact that Pol III depletion affects Pol II transcription. Pol II genes with transcription changes upon Pol III depletion are enriched in diverse cellular functions, and Pol III binding sites are found near their promoters. However, these Pol III binding sites do not correspond to transfer RNAs. Moreover, we demonstrate that Pol III regulates Pol II transcription and chromatin binding of the facilitates chromatin transcription (FACT) complex to alter local chromatin structures, which in turn affects the Pol II transcription rate.

CONCLUSIONS

Our results support a model suggesting that RNA polymerases show cross-regulatory effects: Pol III affects local chromatin structures and the FACT-Pol II axis to regulate the Pol II transcription rate at certain gene loci. This study provides a new perspective for understanding the dysregulation of Pol III in various tissues affected by developmental diseases.

Omics technologies to understand drug toxicity mechanisms

Drug side effects are an important study subject in pharmacology. Recent omics technologies provide a range of omics data and help to understand the biological mechanisms involved in drug effects. These modern technologies provide significant support to all biological disciplines, including drug toxicology. In this review, we provide an overview the use of omics applications to understand drug side effects at the molecular level. We discuss by available omics technologies, their possible uses, as well as their advantages and limitations.

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

Little skate genome provides insights into genetic programs essential for limb-based locomotion

The little skate Leucoraja erinacea, a cartilaginous fish, displays pelvic fin driven walking-like behavior using genetic programs and neuronal subtypes similar to those of land vertebrates. However, mechanistic studies on little skate motor circuit development have been limited, due to a lack of high-quality reference genome. Here, we generated an assembly of the little skate genome, with precise gene annotation and structures, which allowed post-genome analysis of spinal motor neurons (MNs) essential for locomotion. Through interspecies comparison of mouse, skate and chicken MN transcriptomes, shared and divergent gene expression profiles were identified. Comparison of accessible chromatin regions between mouse and skate MNs predicted shared transcription factor (TF) motifs with divergent ones, which could be used for achieving differential regulation of MN-expressed genes. A greater number of TF motif predictions were observed in MN-expressed genes in mouse than in little skate. These findings suggest conserved and divergent molecular mechanisms controlling MN development of vertebrates during evolution, which might contribute to intricate gene regulatory networks in the emergence of a more sophisticated motor system in tetrapods.

Integrated multi-omics approach revealed cellular senescence landscape

Cellular senescence is a complex multifactorial biological phenomenon that plays essential roles in aging, and aging-related diseases. During this process, the senescent cells undergo gene expression altering and chromatin structure remodeling. However, studies on the epigenetic landscape of senescence using integrated multi-omics approaches are limited. In this research, we performed ATAC-seq, RNA-seq and ChIP-seq on different senescent types to reveal the landscape of senescence and identify the prime regulatory elements. We also obtained 34 key genes and deduced that NAT1, PBX1 and RRM2, which interacted with each other, could be the potential markers of aging and aging-related diseases. In summary, our work provides the landscape to study accessibility dynamics and transcriptional regulations in cellular senescence. The application of this technique in different types of senescence allows us to identify the regulatory elements responsible for the substantial regulation of transcription, providing the insights into molecular mechanisms of senescence.

Approaches to benchmark and characterize in vitro human model systems

In vitro human models, such as gastruloids and organoids, are complex three-dimensional (3D) structures often consist of cells from multiple germ layers that possess some attributes of a developing embryo or organ. To use these models to interrogate human development and organogenesis, these 3D models must accurately recapitulate aspects of their in vivo counterparts. Recent advances in single-cell technologies, including sequencing and spatial approaches, have enabled efforts to better understand and directly compare organoids with native tissues. For example, single-cell genomic efforts have created cell and organ atlases that enable benchmarking of in vitro models and can also be leveraged to gain novel biological insights that can be used to further improve in vitro models. This Spotlight discusses the state of current in vitro model systems, the efforts to create large publicly available atlases of the developing human and how these data are being used to improve organoids. Limitations and perspectives on future efforts are also discussed.

Omics technologies in allergy and asthma research: An EAACI position paper

Allergic diseases and asthma are heterogenous chronic inflammatory conditions with several distinct complex endotypes. Both environmental and genetic factors can influence the development and progression of allergy. Complex pathogenetic pathways observed in allergic disorders present a challenge in patient management and successful targeted treatment strategies. The increasing availability of high-throughput omics technologies, such as genomics, epigenomics, transcriptomics, proteomics, and metabolomics allows studying biochemical systems and pathophysiological processes underlying allergic responses. Additionally, omics techniques present clinical applicability by functional identification and validation of biomarkers. Therefore, finding molecules or patterns characteristic for distinct immune-inflammatory endotypes, can subsequently influence its development, progression, and treatment. There is a great potential to further increase the effectiveness of single omics approaches by integrating them with other omics, and nonomics data. Systems biology aims to simultaneously and longitudinally understand multiple layers of a complex and multifactorial disease, such as allergy, or asthma by integrating several, separated data sets and generating a complete molecular profile of the condition. With the use of sophisticated biostatistics and machine learning techniques, these approaches provide in-depth insight into individual biological systems and will allow efficient and customized healthcare approaches, called precision medicine. In this EAACI Position Paper, the Task Force "Omics technologies in allergic research" broadly reviewed current advances and applicability of omics techniques in allergic diseases and asthma research, with a focus on methodology and data analysis, aiming to provide researchers (basic and clinical) with a desk reference in the field. The potential of omics strategies in understanding disease pathophysiology and key tools to reach unmet needs in allergy precision medicine, such as successful patients' stratification, accurate disease prognosis, and prediction of treatment efficacy and successful prevention measures are highlighted.

Multiomic Analysis of Neurons with Divergent Projection Patterns Identifies Novel Regulators of Axon Pathfinding

Axon pathfinding is a key step in neural circuits formation. However, the transcriptional mechanisms regulating its progression remain poorly understood. The binary decision of crossing or avoiding the midline taken by some neuronal axons during development represents a robust model to investigate the mechanisms that control the selection of axonal trajectories. Here, to identify novel regulators of axon guidance, this work compares the transcriptome and chromatin occupancy profiles of two neuronal subpopulations, ipsilateral (iRGC) and contralateral retinal ganglion cells (cRGC), with similar functions but divergent axon trajectories. These analyses retrieved a number of genes encoding for proteins not previously implicated in axon pathfinding. In vivo functional experiments confirm the implication of some of these candidates in axonal navigation. Among the candidate genes, γ-synuclein is identified as essential for inducing midline crossing. Footprint and luciferase assays demonstrate that this small-sized protein is regulated by the transcription factor (TF) Pou4f1 in cRGCs. It is also shown that Lhx2/9 are specifically expressed in iRGCs and control a program that partially overlaps with that regulated by Zic2, previously described as essential for iRGC specification. Overall, the analyses identify dozens of new molecules potentially involved in axon guidance and reveal the regulatory logic behind the selection of axonal trajectories.

Intrinsic bias estimation for improved analysis of bulk and single-cell chromatin accessibility profiles using SELMA

Genome-wide profiling of chromatin accessibility by DNase-seq or ATAC-seq has been widely used to identify regulatory DNA elements and transcription factor binding sites. However, enzymatic DNA cleavage exhibits intrinsic sequence biases that confound chromatin accessibility profiling data analysis. Existing computational tools are limited in their ability to account for such intrinsic biases and not designed for analyzing single-cell data. Here, we present Simplex Encoded Linear Model for Accessible Chromatin (SELMA), a computational method for systematic estimation of intrinsic cleavage biases from genomic chromatin accessibility profiling data. We demonstrate that SELMA yields accurate and robust bias estimation from both bulk and single-cell DNase-seq and ATAC-seq data. SELMA can utilize internal mitochondrial DNA data to improve bias estimation. We show that transcription factor binding inference from DNase footprints can be improved by incorporating estimated biases using SELMA. Furthermore, we show strong effects of intrinsic biases in single-cell ATAC-seq data, and develop the first single-cell ATAC-seq intrinsic bias correction model to improve cell clustering. SELMA can enhance the performance of existing bioinformatics tools and improve the analysis of both bulk and single-cell chromatin accessibility sequencing data.

Genome-wide profiling of chromatin accessibility by DNase-seq or ATAC-seq has been widely used to identify regulatory DNA elements and transcription factor binding sites. Here the authors develop a computational model, SELMA, to estimate and correct enzymatic cleavage biases in chromatin accessibility profiling data.

The Landscape of Accessible Chromatin during Yak Adipocyte Differentiation

Although significant advancement has been made in the study of adipogenesis, knowledge about how chromatin accessibility regulates yak adipogenesis is lacking. We here described genome-wide dynamic chromatin accessibility in preadipocytes and adipocytes by using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), and thus revealed the unique characteristics of open chromatin during yak adipocyte differentiation. The chromatin accessibility of preadipocytes and adipocytes exhibited a similar genomic distribution, displaying a preferential location within the intergenic region, intron, and promoter. The pathway enrichment analysis identified that genes with differential chromatin accessibility were involved in adipogenic metabolism regulation pathways, such as the peroxisome proliferator activated receptor-γ (PPAR) signaling pathway, wingless-type MMTV integration site (Wnt) signaling pathway, and extracellular matrix-receptor (ECM–receptor) interaction. Integration of ATAC-seq and mRNA-seq revealed that genes with a high expression were associated with high levels of chromatin accessibility, especially within 1 kb upstream and downstream of the transcription start site. In addition, we identified a series of transcription factors (TFs) related to adipogenesis and created the TF regulatory network, providing the possible interactions between TFs during yak adipogenesis. This study is crucial for advancing the understanding of transcriptional regulatory mechanisms of adipogenesis and provides valuable information for understanding the adaptation of plateau species to high-altitude environments by maintaining whole body homeostasis through fat metabolism.

Decoding mechanism of action and sensitivity to drug candidates from integrated transcriptome and chromatin state

Omics-based technologies are driving major advances in precision medicine, but efforts are still required to consolidate their use in drug discovery. In this work, we exemplify the use of multi-omics to support the development of 3-chloropiperidines, a new class of candidate anticancer agents. Combined analyses of transcriptome and chromatin accessibility elucidated the mechanisms underlying sensitivity to test agents. Furthermore, we implemented a new versatile strategy for the integration of RNA- and ATAC-seq (Assay for Transposase-Accessible Chromatin) data, able to accelerate and extend the standalone analyses of distinct omic layers. This platform guided the construction of a perturbation-informed basal signature predicting cancer cell lines' sensitivity and to further direct compound development against specific tumor types. Overall, this approach offers a scalable pipeline to support the early phases of drug discovery, understanding of mechanisms, and potentially inform the positioning of therapeutics in the clinic.

WhichTF is functionally important in your open chromatin data?

We present WhichTF, a computational method to identify functionally important transcription factors (TFs) from chromatin accessibility measurements. To rank TFs, WhichTF applies an ontology-guided functional approach to compute novel enrichment by integrating accessibility measurements, high-confidence pre-computed conservation-aware TF binding sites, and putative gene-regulatory models. Comparison with prior sheer abundance-based methods reveals the unique ability of WhichTF to identify context-specific TFs with functional relevance, including NF-κB family members in lymphocytes and GATA factors in cardiac cells. To distinguish the transcriptional regulatory landscape in closely related samples, we apply differential analysis and demonstrate its utility in lymphocyte, mesoderm developmental, and disease cells. We find suggestive, under-characterized TFs, such as RUNX3 in mesoderm development and GLI1 in systemic lupus erythematosus. We also find TFs known for stress response, suggesting routine experimental caveats that warrant careful consideration. WhichTF yields biological insight into known and novel molecular mechanisms of TF-mediated transcriptional regulation in diverse contexts, including human and mouse cell types, cell fate trajectories, and disease-associated cells.

Transcription factors (TFs), a class of DNA binding proteins, regulate tissue- and cell-type-specific expression of genes. Identifying the critical TFs in a given cellular context leads to investigating molecular regulatory mechanisms in development, differentiation, and disease. Because there are more than 1,500 human TFs, experimental measurements of genome-wide occupancy across all TFs have been challenging. While computational approaches play pivotal roles, most existing methods rely on statistical enrichment, focusing either on sequence motif similarity recognized by TFs or the similarity of the genomic region of interest with the previously characterized TF occupancy profile. Here we propose WhichTF as an alternative, incorporating curated biomedical knowledge from ontology and integrating it with the high-confidence prediction of conserved TF binding sites in user-provided genomic regions of interest. We develop a new WhichTF score to rank TFs and demonstrate its applicability across human and mouse cell types, cellular differentiation trajectories, and disease-associated cells.

Integrated Analysis of Transcriptome, microRNAs, and Chromatin Accessibility Revealed Potential Early B-Cell Factor1-Regulated Transcriptional Networks during the Early Development of Fetal Brown Adipose Tissues in Rabbits

In domestic mammals, cold stress decreases the survival rate of newborns and increases the cost of management. Brown adipose tissue (BAT) is the main thermogenic organ for cubs, and well-developed fetal BAT (FBAT) is beneficial for newborns to maintain core temperatures during the first several days of life. However, our knowledge of the epigenetic mechanisms during the early development of FBAT remains largely unknown. Rabbits (Oryctolagus cuniculus) are economically important domestic animals. In this study, a histological analysis shows that the tissue content, thermogenic capacity, and lipid content of FBAT dramatically increases from gestational day 21 (G21) to gestational day 24 (G24) in rabbits. RNA-seq, microRNA-seq (miRNA-seq), and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) show that many genes, miRNAs, and chromatin-accessible regions (referred to as peaks) were identified as significantly changed from G21 to G24, respectively. The upregulated genes from G21 to G24 were significantly enriched in the mitochondrial metabolism and thermogenesis-related signal pathways. The integrated analysis of transcriptome and chromatin accessibility reveals that the peaks in the promoters have a more regulatory effect than peaks in other genomic elements on the expression of their nearby genes in FBATs. The upregulated genes that are associated with increased chromatin accessibility in the promoter regions are involved in the energy metabolism-related signaling pathways. The genes that have a greater tendency to be regulated by miRNAs than the chromatin accessibility in gene promoters are involved in the apelin, insulin, and endocytosis signaling pathways. Furthermore, genome-wide transcription factor (TF) footprinting analysis identifies early B-cell factor1 (EBF1) as playing a key role during early FBAT development. The carbon metabolism, citrate cycle, and PPAR signaling pathways are significantly enriched by the predicted EBF1-regulated cascade TF-network. In conclusion, our work provides a framework for understanding epigenetics regulatory mechanisms underlying the early development of FBAT and identifies potential TF involved in the early development of FBAT in rabbits.

ATAC-STARR-seq reveals transcription factor-bound activators and silencers within chromatin-accessible regions of the human genome

Massively parallel reporter assays (MPRAs) test the capacity of putative gene regulatory elements to drive transcription on a genome-wide scale. Most gene regulatory activity occurs within accessible chromatin, and recently described methods have combined assays that capture these regions-such as assay for transposase-accessible chromatin using sequencing (ATAC-seq)-with self-transcribing active regulatory region sequencing (STARR-seq) to selectively assay the regulatory potential of accessible DNA (ATAC-STARR-seq). Here, we report an integrated approach that quantifies activating and silencing regulatory activity, chromatin accessibility, and transcription factor (TF) occupancy with one assay using ATAC-STARR-seq. Our strategy, including important updates to the ATAC-STARR-seq assay and workflow, enabled high-resolution testing of ∼50 million unique DNA fragments tiling ∼101,000 accessible chromatin regions in human lymphoblastoid cells. We discovered that 30% of all accessible regions contain an activator, a silencer, or both. Although few MPRA studies have explored silencing activity, we demonstrate that silencers occur at similar frequencies to activators, and they represent a distinct functional group enriched for unique TF motifs and repressive histone modifications. We further show that Tn5 cut-site frequencies are retained in the ATAC-STARR plasmid library compared to standard ATAC-seq, enabling TF occupancy to be ascertained from ATAC-STARR data. With this approach, we found that activators and silencers cluster by distinct TF footprint combinations, and these groups of activity represent different gene regulatory networks of immune cell function. Altogether, these data highlight the multilayered capabilities of ATAC-STARR-seq to comprehensively investigate the regulatory landscape of the human genome all from a single DNA fragment source.

Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells

Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.

Epigenetic Application of ATAC-Seq Based on Tn5 Transposase Purification Technology

Background

Assays of transposase accessible chromatin sequencing (ATAC-seq) is an efficient assay to investigate chromatin accessibility, which depends on the activity of a robust Tn5 transposase to fragment the genome while cutting in the sequencing adapters.

Methods

We propose reliable approaches for purifying hyperactive Tn5 transposase by chitin magnetic bead sorting. Double-stranded DNA of J76 cells and 293T cells were digested and subjected to tagmentation as test samples with Tn5 transposase, and libraries were established and sequenced. Sequencing data was then analyzed for peak calling, GO enrichment, and motif analysis.

Results

We report a set of rapid, efficient, and low-cost methods for ATAC-seq library construction and data analysis, through large-scale and rapid sequencing. These methods can provide a reference for the study of epigenetic regulation of gene expression.

Banp regulates DNA damage response and chromosome segregation during the cell cycle in zebrafish retina

Btg3-associated nuclear protein (Banp) was originally identified as a nuclear matrix-associated region (MAR)-binding protein and it functions as a tumor suppressor. At the molecular level, Banp regulates transcription of metabolic genes via a CGCG-containing motif called the Banp motif. However, its physiological roles in embryonic development are unknown. Here, we report that Banp is indispensable for the DNA damage response and chromosome segregation during mitosis. Zebrafish banp mutants show mitotic cell accumulation and apoptosis in developing retina. We found that DNA replication stress and tp53-dependent DNA damage responses were activated to induce apoptosis in banp mutants, suggesting that Banp is required for regulation of DNA replication and DNA damage repair. Furthermore, consistent with mitotic cell accumulation, chromosome segregation was not smoothly processed from prometaphase to anaphase in banp morphants, leading to a prolonged M-phase. Our RNA- and ATAC-sequencing identified 31 candidates for direct Banp target genes that carry the Banp motif. Interestingly, a DNA replication fork regulator, wrnip1, and two chromosome segregation regulators, cenpt and ncapg, are included in this list. Thus, Banp directly regulates transcription of wrnip1 for recovery from DNA replication stress, and cenpt and ncapg for chromosome segregation during mitosis. Our findings provide the first in vivo evidence that Banp is required for cell-cycle progression and cell survival by regulating DNA damage responses and chromosome segregation during mitosis.

In order for a cell to divide, it must progress through a series of carefully controlled steps known as the cell cycle. First, the cell replicates its DNA and both copies get segregated to opposite ends. The cell then splits into two and each new cell receives a copy of the duplicated genetic material. If any of the stages in the cell cycle become disrupted or mis-regulated this can lead to uncontrolled divisions that may result in cancer.

Researchers have often used a structure within the eye known as the retina to study the cell cycle in zebrafish and other animals as cells in the retina rapidly divide in a highly controlled manner. A protein called Banp is known to help stop tumors from growing in humans and mice, but its normal role in the body, particularly the cell cycle, has remained unclear. To investigate, Babu et al. studied the retina of mutant zebrafish that were unable to make the Banp protein.

The experiments revealed that two stress responses indicating DNA damage or defects in copying DNA were active in the retinal cells of the mutant zebrafish. This suggested that Banp allows cell to progress through the cell cycle by repairing any DNA damage that may arise during replication. Banp does this by activating the gene for another protein called Wrnip1. Babu et al. also found that Banp helps segregate the two copies of DNA during cell division by promoting the activation of two other proteins called Cenpt and Ncapg. Further experiments identified 31 genes that were directly regulated by Banp.

These findings demonstrate that Banp is required for zebrafish cells to be able to accurately copy their DNA and divide in to two new cells. In the future, the work of Babu et al. will provide a useful resource to investigate how tumors grow and spread around the body, and may contribute to the development of new treatments for cancer.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

Association of a common genetic variant with Parkinson's disease is mediated by microglia

Studies of multiple neurodegenerative disorders have identified many genetic variants that are associated with risk of disease throughout a lifetime. For example, Parkinson's disease (PD) risk is attributed in part to both coding mutations in the leucine-rich repeat kinase 2 (LRRK2) gene and to a common noncoding variation in the 5' region of the LRRK2 locus, as identified by genome-wide association studies (GWAS). However, the mechanisms linking GWAS variants to pathogenicity are largely unknown. Here, we found that the influence of PD-associated noncoding variation on LRRK2 expression is specifically propagated through microglia and not by other cell types that express LRRK2 in the human brain. We find microglia-specific regulatory chromatin regions that modulate the LRRK2 expression in human frontal cortex and substantia nigra and confirm these results in a human-induced pluripotent stem cell-derived microglia model. We showed, using a large-scale clustered regularly interspaced short palindromic repeats interference (CRISPRi) screen, that a regulatory DNA element containing the single-nucleotide variant rs6581593 influences the LRRK2 expression in microglia. Our study demonstrates that cell type should be considered when evaluating the role of noncoding variation in disease pathogenesis and sheds light on the mechanism underlying the association of the 5' region of LRRK2 with PD risk.

Loss of Monoallelic Expression of IGF2 in the Adult Liver Via Alternative Promoter Usage and Chromatin Reorganization

In mammals, genomic imprinting operates via gene silencing mechanisms. Although conservation of the imprinting mechanism at the H19/IGF2 locus has been generally described in pigs, tissue-specific imprinting at the transcript level, monoallelic-to-biallelic conversion, and spatio-temporal chromatin reorganization remain largely uninvestigated. Here, we delineate spatially regulated imprinting of IGF2 transcripts, age-dependent hepatic mono- to biallelic conversion, and reorganization of topologically associating domains at the porcine H19/IGF2 locus for better translation to human and animal research. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of normal and parthenogenetic porcine embryos revealed the paternally hypermethylated H19 differentially methylated region and paternal expression of IGF2. Using a polymorphism-based approach and omics datasets from chromatin immunoprecipitation sequencing (ChIP–seq), whole-genome sequencing (WGS), RNA-seq, and Hi-C, regulation of IGF2 during development was analyzed. Regulatory elements in the liver were distinguished from those in the muscle where the porcine IGF2 transcript was monoallelically expressed. The IGF2 transcript from the liver was biallelically expressed at later developmental stages in both pigs and humans. Chromatin interaction was less frequent in the adult liver compared to the fetal liver and skeletal muscle. The duration of genomic imprinting effects within the H19/IGF2 locus might be reduced in the liver with biallelic conversion through alternative promoter usage and chromatin remodeling. Our integrative omics analyses of genome, epigenome, and transcriptome provided a comprehensive view of imprinting status at the H19/IGF2 cluster.

Inhibition of Cancer Cell Migration and Glycolysis by Terahertz Wave Modulation via Altered Chromatin Accessibility

Metastasis and metabolic disorders contribute to most cancer deaths and are potential drug targets in cancer treatment. However, corresponding drugs inevitably induce myeloid suppression and gastrointestinal toxicity. Here, we report a nonpharmaceutical and noninvasive electromagnetic intervention technique that exhibited long-term inhibition of cancer cells. Firstly, we revealed that optical radiation at the specific wavelength of 3.6 μm (i.e., 83 THz) significantly increased binding affinity between DNA and histone via molecular dynamics simulations, providing a theoretical possibility for THz modulation- (THM-) based cancer cell intervention. Subsequent cell functional assays demonstrated that low-power 3.6 μm THz wave could successfully inhibit cancer cell migration by 50% and reduce glycolysis by 60%. Then, mRNA sequencing and assays for transposase-accessible chromatin using sequencing (ATAC-seq) indicated that low-power THM at 3.6 μm suppressed the genes associated with glycolysis and migration by reducing the chromatin accessibility of certain gene loci. Furthermore, THM at 3.6 μm on HCT-116 cancer cells reduced the liver metastasis by 60% in a metastatic xenograft mouse model by splenic injection, successfully validated the inhibition of cancer cell migration by THM in vivo. Together, this work provides a new paradigm for electromagnetic irradiation-induced epigenetic changes and represents a theoretical basis for possible innovative therapeutic applications of THM as the future of cancer treatments.

Epigenetic Memories in Hematopoietic Stem and Progenitor Cells

The recent development of next-generation sequencing (NGS) technologies has contributed to research into various biological processes. These novel NGS technologies have revealed the involvement of epigenetic memories in trained immunity, which are responses to transient stimulation and result in better responses to secondary challenges. Not only innate system cells, such as macrophages, monocytes, and natural killer cells, but also bone marrow hematopoietic stem cells (HSCs) have been found to gain memories upon transient stimulation, leading to the enhancement of responses to secondary challenges. Various stimuli, including microbial infection, can induce the epigenetic reprogramming of innate immune cells and HSCs, which can result in an augmented response to secondary stimulation. In this review, we introduce novel NGS technologies and their application to unraveling epigenetic memories that are key in trained immunity and summarize the recent findings in trained immunity. We also discuss our most recent finding regarding epigenetic memory in aged HSCs, which may be associated with the exposure of HSCs to aging-related stresses.

Genome-Wide Analysis of Yeast Metabolic Cycle through Metabolic Network Models Reveals Superiority of Integrated ATAC-seq Data over RNA-seq Data

Saccharomyces cerevisiae undergoes robust oscillations to regulate its physiology for adaptation and survival under nutrient-limited conditions. Environmental cues can induce rhythmic metabolic alterations in order to facilitate the coordination of dynamic metabolic behaviors. Of such metabolic processes, the yeast metabolic cycle enables adaptation of the cells to varying nutritional status through oscillations in gene expression and metabolite production levels. In this process, yeast metabolism is altered between diverse cellular states based on changing oxygen consumption levels: quiescent (reductive charging [RC]), growth (oxidative [OX]), and proliferation (reductive building [RB]) phases. We characterized metabolic alterations during the yeast metabolic cycle using a variety of approaches. Gene expression levels are widely used for condition-specific metabolic simulations, whereas the use of epigenetic information in metabolic modeling is still limited despite the clear relationship between epigenetics and metabolism. This prompted us to investigate the contribution of epigenomic information to metabolic predictions for progression of the yeast metabolic cycle. In this regard, we determined altered pathways through the prediction of regulated reactions and corresponding model genes relying on differential chromatin accessibility levels. The predicted metabolic alterations were confirmed via data analysis and literature. We subsequently utilized RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data sets in the contextualization of the yeast model. The use of ATAC-seq data considerably enhanced the predictive capability of the model. To the best of our knowledge, this is the first attempt to use genome-wide chromatin accessibility data in metabolic modeling. The preliminary results showed that epigenomic data sets can pave the way for more accurate metabolic simulations.

IMPORTANCE Dynamic chromatin organization mediates the emergence of condition-specific phenotypes in eukaryotic organisms. Saccharomyces cerevisiae can alter its metabolic profile via regulation of genome accessibility and robust transcriptional oscillations under nutrient-limited conditions. Thus, both epigenetic information and transcriptomic information are crucial in the understanding of condition-specific metabolic behavior in this organism. Based on genome-wide alterations in chromatin accessibility and transcription, we investigated the yeast metabolic cycle, which is a remarkable example of coordinated and dynamic yeast behavior. In this regard, we assessed the use of ATAC-seq and RNA-seq data sets in condition-specific metabolic modeling. To our knowledge, this is the first attempt to use chromatin accessibility data in the reconstruction of context-specific metabolic models, despite the extensive use of transcriptomic data. As a result of comparative analyses, we propose that the incorporation of epigenetic information is a promising approach in the accurate prediction of metabolic dynamics.

Bridging Glycomics and Genomics: New Uses of Functional Genetics in the Study of Cellular Glycosylation

All living cells are coated with a diverse collection of carbohydrate molecules called glycans. Glycans are key regulators of cell behavior and important therapeutic targets for human disease. Unlike proteins, glycans are not directly templated by discrete genes. Instead, they are produced through multi-gene pathways that generate a heterogenous array of glycoprotein and glycolipid antigens on the cell surface. This genetic complexity has sometimes made it challenging to understand how glycosylation is regulated and how it becomes altered in disease. Recent years, however, have seen the emergence of powerful new functional genomics technologies that allow high-throughput characterization of genetically complex cellular phenotypes. In this review, we discuss how these techniques are now being applied to achieve a deeper understanding of glyco-genomic regulation. We highlight specifically how methods like ChIP-seq, RNA-seq, CRISPR genomic screening and scRNA-seq are being used to map the genomic basis for various cell-surface glycosylation states in normal and diseased cell types. We also offer a perspective on how emerging functional genomics technologies are likely to create further opportunities for studying cellular glycobiology in the future. Taken together, we hope this review serves as a primer to recent developments at the glycomics-genomics interface.

Maternal Factor Dppa3 Activates 2C-Like Genes and Depresses DNA Methylation in Mouse Embryonic Stem Cells

Mouse embryonic stem cells (ESCs) contain a rare cell population of “two-cell embryonic like” cells (2CLCs) that display similar features to those found in the two-cell (2C) embryo and thus represent an in vitro model for studying the progress of zygotic genome activation (ZGA). However, the positive regulator determinants of the 2CLCs’ conversion and ZGA have not been completely elucidated. Here, we identify a new regulator promoting 2CLCs and ZGA transcripts. Through a combination of overexpression (OE), knockdown (KD), together with transcriptional analysis and methylome analysis, we find that Dppa3 regulates the 2CLC-associated transcripts, DNA methylation, and 2CLC population in ESCs. The differentially methylated regions (DMRs) analysis identified 6,920 (98.2%) hypomethylated, whilst only 129 (1.8%) hypermethylated, regions in Dppa3 OE ESCs, suggesting that Dppa3 facilitates 2CLCs reprogramming. The conversion to 2CLCs by overexpression of Dppa3 is also associated with DNA damage response. Dppa3 knockdown manifest impairs transition into the 2C-like state. Global DNA methylome and chromatin state analysis of Dppa3 OE ESCs reveal that Dppa3 facilitates the chromatin configuration to 2CLCs reversion. Our finding for the first time elucidates a novel role of Dppa3 in mediating the 2CLC conversion, and suggests that Dppa3 is a new regulator for ZGA progress.

Identification and classification of cis-regulatory elements in the amphipod crustacean Parhyale hawaiensis

Emerging research organisms enable the study of biology that cannot be addressed using classical 'model' organisms. New data resources can accelerate research in such animals. Here, we present new functional genomic resources for the amphipod crustacean Parhyale hawaiensis, facilitating the exploration of gene regulatory evolution using this emerging research organism. We use Omni-ATAC-seq to identify accessible chromatin genome-wide across a broad time course of Parhyale embryonic development. This time course encompasses many major morphological events, including segmentation, body regionalization, gut morphogenesis and limb development. In addition, we use short- and long-read RNA-seq to generate an improved Parhyale genome annotation, enabling deeper classification of identified regulatory elements. We discover differential accessibility, predict nucleosome positioning, infer transcription factor binding, cluster peaks based on accessibility dynamics, classify biological functions and correlate gene expression with accessibility. Using a Minos transposase reporter system, we demonstrate the potential to identify novel regulatory elements using this approach. This work provides a platform for the identification of novel developmental regulatory elements in Parhyale, and offers a framework for performing such experiments in other emerging research organisms.

FTO mediates LINE1 m6A demethylation and chromatin regulation in mESCs and mouse development

N6-methyladenosine (m6A) is the most abundant internal modification on mammalian messenger RNA. It is installed by a writer complex and can be reversed by erasers such as the fat mass and obesity-associated protein FTO. Despite extensive research, the primary physiological substrates of FTO in mammalian tissues and development remain elusive. Here, we show that FTO mediates m6A demethylation of long-interspersed element-1 (LINE1) RNA in mouse embryonic stem cells (mESCs), regulating LINE1 RNA abundance and the local chromatin state, which in turn modulates the transcription of LINE1-containing genes. FTO-mediated LINE1 RNA m6A demethylation also plays regulatory roles in shaping chromatin state and gene expression during mouse oocyte and embryonic development. Our results suggest broad effects of LINE1 RNA m6A demethylation by FTO in mammals.

Investigating chromatin accessibility during development and differentiation by ATAC-sequencing to guide the identification of cis-regulatory elements

Mapping accessible chromatin across time scales can give insights into its dynamic nature, for example during cellular differentiation and tissue or organism development. Analysis of such data can be utilised to identify functional cis-regulatory elements (CRE) and transcription factor binding sites and, when combined with transcriptomics, can reveal gene regulatory networks (GRNs) of expressed genes. Chromatin accessibility mapping is a powerful approach and can be performed using ATAC-sequencing (ATAC-seq), whereby Tn5 transposase inserts sequencing adaptors into genomic DNA to identify differentially accessible regions of chromatin in different cell populations. It requires low sample input and can be performed and analysed relatively quickly compared with other methods. The data generated from ATAC-seq, along with other genomic approaches, can help uncover chromatin packaging and potential cis-regulatory elements that may be responsible for gene expression. Here, we describe the ATAC-seq approach and give examples from mainly vertebrate embryonic development, where such datasets have identified the highly dynamic nature of chromatin, with differing landscapes between cellular precursors for different lineages.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

Translator: A Transfer Learning Approach to Facilitate Single-Cell ATAC-Seq Data Analysis from Reference Dataset

Recent advances in single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) have allowed simultaneous epigenetic profiling over thousands of individual cells to dissect the cellular heterogeneity and elucidate regulatory mechanisms at the finest possible resolution. However, scATAC-seq is challenging to model computationally due to the ultra-high dimensionality, low signal-to-noise ratio, complex feature interactions, and high vulnerability to various confounding factors. In this study, we present Translator, an efficient transfer learning approach to capture generalizable chromatin interactions from high-quality (HQ) reference scATAC-seq data to obtain robust cell representations in low-to-moderate quality target scATAC-seq data. We applied Translator on various simulated and real scATAC-seq datasets and demonstrated that Translator could learn more biologically meaningful cell representations than other methods by incorporating information learned from the reference data, thus facilitating various downstream analyses such as clustering and motif enrichment measurements. Moreover, Translator's block-wise deep learning framework can handle nonlinear relationships with restricted connections using fewer parameters to boost computational efficiency through Graphics Processing Unit (GPU) parallelism. Finally, we have implemented Translator as a free software package available for the community to leverage large-scale, HQ reference data to study target scATAC-seq data.

Twist2-driven chromatin remodeling governs the postnatal maturation of dermal fibroblasts

Dermal fibroblasts lose stem cell potency after birth, which prevents regenerative healing. However, the underlying intracellular mechanisms are largely unknown. We uncover the postnatal maturation of papillary fibroblasts (PFs) driven by the extensive Twist2-mediated remodeling of chromatin accessibility. A loss of the regenerative ability of postnatal PFs occurs with decreased H3K27ac levels. Single-cell transcriptomics, assay for transposase-accessible chromatin sequencing (ATAC-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) reveal the postnatal maturation trajectory associated with the loss of the regenerative trajectory in PFs, which is characterized by a marked decrease in chromatin accessibility and H3K27ac modifications. Histone deacetylase inhibition delays spontaneous chromatin remodeling, thus maintaining the regenerative ability of postnatal PFs. Genomic analysis identifies Twist2 as a major regulator within chromatin regions with decreased accessibility during the postnatal period. When Twist2 is genetically deleted in dermal fibroblasts, the intracellular cascade of postnatal maturation is significantly delayed. Our findings reveal the comprehensive intracellular mechanisms underlying intrinsic postnatal changes in dermal fibroblasts.

High-throughput single-сell sequencing in cancer research

With advances in sequencing and instrument technology, bioinformatics analysis is being applied to batches of massive cells at single-cell resolution. High-throughput single-cell sequencing can be utilized for multi-omics characterization of tumor cells, stromal cells or infiltrated immune cells to evaluate tumor progression, responses to environmental perturbations, heterogeneous composition of the tumor microenvironment, and complex intercellular interactions between these factors. Particularly, single-cell sequencing of T cell receptors, alone or in combination with single-cell RNA sequencing, is useful in the fields of tumor immunology and immunotherapy. Clinical insights obtained from single-cell analysis are critically important for exploring the biomarkers of disease progression or antitumor treatment, as well as for guiding precise clinical decision-making for patients with malignant tumors. In this review, we summarize the clinical applications of single-cell sequencing in the fields of tumor cell evolution, tumor immunology, and tumor immunotherapy. Additionally, we analyze the tumor cell response to antitumor treatment, heterogeneity of the tumor microenvironment, and response or resistance to immune checkpoint immunotherapy. The limitations of single-cell analysis in cancer research are also discussed.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

Chromatin Structure and Dynamics: Focus on Neuronal Differentiation and Pathological Implication

Chromatin structure is an essential regulator of gene expression. Its state of compaction contributes to the regulation of genetic programs, in particular during differentiation. Epigenetic processes, which include post-translational modifications of histones, DNA methylation and implication of non-coding RNA, are powerful regulators of gene expression. Neurogenesis and neuronal differentiation are spatio-temporally regulated events that allow the formation of the central nervous system components. Here, we review the chromatin structure and post-translational histone modifications associated with neuronal differentiation. Studying the impact of histone modifications on neuronal differentiation improves our understanding of the pathophysiological mechanisms of chromatinopathies and opens up new therapeutic avenues. In addition, we will discuss techniques for the analysis of histone modifications on a genome-wide scale and the pathologies associated with the dysregulation of the epigenetic machinery.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

IL-17 Producing Lymphocytes Cause Dry Eye and Corneal Disease With Aging in RXRα Mutant Mouse

Purpose

To investigate IL-17 related mechanisms for developing dry eye disease in the Pinkie mouse strain with a loss of function RXRα mutation.

Methods

Measures of dry eye disease were assessed in the cornea and conjunctiva. Expression profiling was performed by single-cell RNA sequencing (scRNA-seq) to compare gene expression in conjunctival immune cells. Conjunctival immune cells were immunophenotyped by flow cytometry and confocal microscopy. The activity of RXRα ligand 9-cis retinoic acid (RA) was evaluated in cultured monocytes and γδ T cells.

Results

Compared to wild type (WT) C57BL/6, Pinkie has increased signs of dry eye disease, including decreased tear volume, corneal barrier disruption, corneal/conjunctival cornification and goblet cell loss, and corneal vascularization, opacification, and ulceration with aging. ScRNA-seq of conjunctival immune cells identified γδ T cells as the predominant IL-17 expressing population in both strains and there is a 4-fold increased percentage of γδ T cells in Pinkie. Compared to WT, IL-17a, and IL-17f significantly increased in Pinkie with conventional T cells and γδ T cells as the major producers. Flow cytometry revealed an increased number of IL-17⁺ γδ T cells in Pinkie. Tear concentration of the IL-17 inducer IL-23 is significantly higher in Pinkie. 9-cis RA treatment suppresses stimulated IL-17 production by γδ T and stimulatory activity of monocyte supernatant on γδ T cell IL-17 production. Compared to WT bone marrow chimeras, Pinkie chimeras have increased IL-17⁺ γδ T cells in the conjunctiva after desiccating stress and anti-IL-17 treatment suppresses dry eye induced corneal MMP-9 production/activity and conjunctival goblet cell loss.

Conclusion

These findings indicate that RXRα suppresses generation of dry eye disease-inducing IL-17 producing lymphocytes s in the conjunctiva and identifies RXRα as a potential therapeutic target in dry eye.

Chromatin accessibility and transcriptome integrative analysis revealed AP-1-mediated genes potentially modulate histopathology features in psoriasis

Background

Psoriasis is a chronic and hyperproliferative skin disease featured by hyperkeratosis with parakeratosis, Munro micro-abscess, elongation of rete pegs, granulosa thinning, and lymphocyte infiltration. We previously profiled gene expression and chromatin accessibility of psoriatic skins by transcriptome sequencing and ATAC-seq. However, integrating both of these datasets to unravel gene expression regulation is lacking. Here, we integrated transcriptome and ATAC-seq of the same psoriatic and normal skin tissues, trying to leverage the potential role of chromatin accessibility and their function in histopathology features.

Results

By inducing binding and expression target analysis (BETA) algorithms, we explored the target prediction of transcription factors binding in 15 psoriatic and 19 control skins. BETA identified 408 upregulated genes (rank product < 0.01) and 133 downregulated genes linked with chromatin accessibility. We noticed that cumulative fraction of genes in upregulation group was statistically higher than background, while that of genes in downregulation group was not significant. KEGG pathway analysis showed that the upregulated 408 genes were enriched in TNF, NOD, and IL-17 signaling pathways. In addition, the motif module in BETA suggested the 57 upregulated genes are targeted by transcription factor AP-1, indicating that increased chromatin accessibility facilitated the binding of AP-1 to the target regions and further induced expression of relevant genes. Among these genes, SQLE, STRN, EIF4, and MYO1B expression was increased in patients with hyperkeratosis, parakeratosis, and acanthosis thickening.

Conclusions

In summary, with the advantage of BETA, we identified a series of genes that contribute to the disease pathogenesis, especially in modulating histopathology features, providing us with new clues in treating psoriasis.

Bibliometric review of ATAC-Seq and its application in gene expression

With recent advances in high-throughput next-generation sequencing, it is possible to describe the regulation and expression of genes at multiple levels. An assay for transposase-accessible chromatin using sequencing (ATAC-seq), which uses Tn5 transposase to sequence protein-free binding regions of the genome, can be combined with chromatin immunoprecipitation coupled with deep sequencing (ChIP-seq) and ribonucleic acid sequencing (RNA-seq) to provide a detailed description of gene expression. Here, we reviewed the literature on ATAC-seq and described the characteristics of ATAC-seq publications. We then briefly introduced the principles of RNA-seq, ChIP-seq and ATAC-seq, focusing on the main features of the techniques. We built a phylogenetic tree from species that had been previously studied by using ATAC-seq. Studies of Mus musculus and Homo sapiens account for approximately 90% of the total ATAC-seq data, while other species are still in the process of accumulating data. We summarized the findings from human diseases and other species, illustrating the cutting-edge discoveries and the role of multi-omics data analysis in current research. Moreover, we collected and compared ATAC-seq analysis pipelines, which allowed biological researchers who lack programming skills to better analyze and explore ATAC-seq data. Through this review, it is clear that multi-omics analysis and single-cell sequencing technology will become the mainstream approach in future research.

Proteomic analysis of alcohol-associated hepatitis reveals glycoprotein NMB (GPNMB) as a novel hepatic and serum biomarker

Alcohol consumption remains a leading cause of liver disease worldwide, resulting in a complex array of hepatic pathologies, including steatosis, steatohepatitis, and cirrhosis. Individuals who progress to a rarer form of alcohol-associated liver disease (ALD), alcohol-associated hepatitis (AH), require immediate life-saving intervention in the form of liver transplantation. Rapid onset of AH is poorly understood and the metabolic mechanisms contributing to the progression to liver failure remain undetermined. While multiple mechanisms have been identified that contribute to ALD, no cures exist and mortality from AH remains high. To identify novel pathways associated with AH, our group utilized proteomics to investigate AH-specific biomarkers in liver explant tissues. The goal of the present study was to determine changes in the proteome as well as epigenetic changes occurring in AH. Protein abundance and acetylomic analyses were performed utilizing nHPLC-MS/MS, revealing significant changes to proteins associated with metabolic and inflammatory fibrosis pathways. Here, we describe a novel hepatic and serum biomarker of AH, glycoprotein NMB (GPNMB). The anti-inflammatory protein GPNMB was significantly increased in AH explant liver and serum compared to healthy donors by 50-fold and 6.5-fold, respectively. Further, bioinformatics analyses identified an AH-dependent decrease in protein abundance across fatty acid degradation, biosynthesis of amino acids, and carbon metabolism. The greatest increases in protein abundance were observed in pathways for focal adhesion, lysosome, phagosome, and actin cytoskeleton. In contrast with the hyperacetylation observed in murine models of ALD, protein acetylation was decreased in AH compared to normal liver across fatty acid degradation, biosynthesis of amino acids, and carbon metabolism. Interestingly, immunoblot analysis found epigenetic marks were significantly increased in AH explants, including Histone H3K9 and H2BK5 acetylation. The increased acetylation of histones likely plays a role in the altered proteomic profile observed, including increases in GPNMB. Indeed, our results reveal that the AH proteome is dramatically impacted through unanticipated and unknown mechanisms. Understanding the origin and consequences of these changes will yield new mechanistic insight for ALD as well as identify novel hepatic and serum biomarkers, such as GPNMB.