Factors affecting the cleavage efficiency of the CRISPR-Cas9 system

The CRISPR-Cas system stands out as a promising genome editing tool due to its cost-effectiveness and time efficiency compared to other methods. This system has tremendous potential for treating various diseases, including genetic disorders and cancer, and promotes therapeutic research for a wide range of genetic diseases. Additionally, the CRISPR-Cas system simplifies the generation of animal models, offering a more accessible alternative to traditional methods. The CRISPR-Cas9 system can be used to cleave target DNA strands that need to be corrected, causing double-strand breaks (DSBs). DNA with DSBs can then be recovered by the DNA repair pathway that the CRISPR-Cas9 system uses to edit target gene sequences. High cleavage efficiency of the CRISPR-Cas9 system is thus imperative for effective gene editing. Herein, we explore several factors affecting the cleavage efficiency of the CRISPR-Cas9 system. These factors include the GC content of the protospacer-adjacent motif (PAM) proximal and distal regions, single-guide RNA (sgRNA) properties, and chromatin state. These considerations contribute to the efficiency of genome editing.

Environmental factor reversibly determines cellular identity through opposing Integrators that unify epigenetic and transcriptional pathways

Organisms must adapt to environmental stresses to ensure their survival and prosperity. Different types of stresses, including thermal, mechanical, and hypoxic stresses, can alter the cellular state that accompanies changes in gene expression but not the cellular identity determined by a chromatin state that remains stable throughout life. Some tissues, such as adipose tissue, demonstrate remarkable plasticity and adaptability in response to environmental cues, enabling reversible cellular identity changes; however, the mechanisms underlying these changes are not well understood. We hypothesized that positive and/or negative "Integrators" sense environmental cues and coordinate the epigenetic and transcriptional pathways required for changes in cellular identity. Adverse environmental factors such as pollution disrupt the coordinated control contributing to disease development. Further research based on this hypothesis will reveal how organisms adapt to fluctuating environmental conditions, such as temperature, extracellular matrix stiffness, oxygen, cytokines, and hormonal cues by changing their cellular identities.

Targeted recombination of homologous chromosomes using CRISPR-Cas9

CRISPR mutagenesis is an efficient way to disrupt specific target genes in many model organisms. We previously devised a targeted CRISPR recombination method to generate intragenic recombinants of alleles in Drosophila. Here, we assessed the applicability of CRISPR targeting-induced recombination to different genetic loci. We compared the ectopic recombination rates in the male germline by CRISPR targeting at two neighboring genetic loci within the genomic region that consists of the repressed chromatin domain of the Lobe gene, and the transcriptionally active domain of PRAS40. Targeting around the transcription initiation of PRAS40 resulted in higher recombination rates of homologous chromosomes than targeting at the Lobe intron. Based on the efficient homologous recombination by CRISPR targeting observed around transcriptionally active loci, we further investigated targeted recombination between P-elements that are inserted at different genomic locations. Male recombination by CRISPR targeting of P-elements located proximally and distally to the ebony gene produced recombinants deficient for the intervening region of ebony transcription. Taken together, we suggest that targeted homologous recombination by CRISPR targeting may have specific genetic applications, such as generation of allelic combinations or chromosomal variations.

Histone variants shape chromatin states in Arabidopsis

How different intrinsic sequence variations and regulatory modifications of histones combine in nucleosomes remain unclear. To test the importance of histone variants in the organization of chromatin we investigated how histone variants and histone modifications assemble in the Arabidopsis thaliana genome. We showed that a limited number of chromatin states divide euchromatin and heterochromatin into several subdomains. We found that histone variants are as significant as histone modifications in determining the composition of chromatin states. Particularly strong associations were observed between H2A variants and specific combinations of histone modifications. To study the role of H2A variants in organizing chromatin states we determined the role of the chromatin remodeler DECREASED IN DNA METHYLATION (DDM1) in the organization of chromatin states. We showed that the loss of DDM1 prevented the exchange of the histone variant H2A.Z to H2A.W in constitutive heterochromatin, resulting in significant effects on the definition and distribution of chromatin states in and outside of constitutive heterochromatin. We thus propose that dynamic exchanges of histone variants control the organization of histone modifications into chromatin states, acting as molecular landmarks.

Histone Maps in Gossypium darwinii Reveal Epigenetic Regulation Drives Subgenome Divergence and Cotton Domestication

The functional annotation of genomes, including chromatin modifications, is essential to understand the intricate architecture of chromatin and the consequential gene regulation. However, such an annotation remains limited for cotton genomes. Here, we conducted chromatin profiling in a wild allotetraploid cotton Gossypium darwinii (AD genome) by integrating the data of histone modification, transcriptome, and chromatin accessibility. We revealed that the A subgenome showed a higher level of active histone marks and lower level of repressive histone marks than the D subgenome, which was consistent with the expression bias between the two subgenomes. We show that the bias in transcription and histone modification between the A and D subgenomes may be caused by genes unique to the subgenome but not by homoeologous genes. Moreover, we integrate histone marks and open chromatin to define six chromatin states (S1-S6) across the cotton genome, which index different genomic elements including genes, promoters, and transposons, implying distinct biological functions. In comparison to the domesticated cotton species, we observed that 23.2% of genes in the genome exhibit a transition from one chromatin state to another at their promoter. Strikingly, the S2 (devoid of epigenetic marks) to S3 (enriched for the mark of open chromatin) was the largest transition group. These transitions occurred simultaneously with changes in gene expression, which were significantly associated with several domesticated traits in cotton. Collectively, our study provides a useful epigenetic resource for research on allopolyploid plants. The domestication-induced chromatin dynamics and associated genes identified here will aid epigenetic engineering, improving polyploid crops.

Chromatin-interacting RNA-binding proteins regulate transcription

RNA-binding proteins (RBPs) are essential regulators involved in the fate determination of diverse RNA species; however, emerging evidence indicates that a subset of RBPs may physically interact with chromatin and function at the transcriptional level. Here, we highlight the recently discovered mechanisms of chromatin-interacting RBPs (ChRBPs) in the regulation of chromatin/transcriptional activities.

Epigenetic Landscape Is Largely Shaped by Diversiform Transposons in Aegilops tauschii

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

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

The dynamics of chromatin states mediated by epigenetic modifications during somatic cell reprogramming

Somatic cell reprogramming (SCR) is the conversion of differentiated somatic cells into totipotent or pluripotent cells through a variety of methods. Somatic cell reprogramming also provides a platform to investigate the role of chromatin-based factors in establishing and maintaining totipotency or pluripotency, since high expression of totipotency- or pluripotency-related genes usually require an active chromatin state. Several studies in plants or mammals have recently shed light on the molecular mechanisms by which epigenetic modifications regulate the expression of totipotency or pluripotency genes by altering their chromatin states. In this review, we present a comprehensive overview of the dynamic changes in epigenetic modifications and chromatin states during reprogramming from somatic cells to totipotent or pluripotent cells. In addition, we illustrate the potential role of DNA methylation, histone modifications, histone variants, and chromatin remodeling during somatic cell reprogramming, which will pave the way to developing reliable strategies for efficient cellular reprogramming.

Bioinformatic Evaluation of Features on Cis-regulatory Elements at 6q25.1

Eukaryotic non-coding regulatory features contribute significantly to cellular plasticity which on aberration leads to cellular malignancy. Enhancers are cis-regulatory elements that contribute to the development of resistance to endocrine therapy in estrogen receptor (ER)-positive breast cancer leading to poor clinical outcome. ER is vital for therapeutic targets in ER-positive breast cancer. Here, we review and report the different regulatory features present on ER with the objective to delineate potential mechanisms which may contribute to development of resistance. The UCSC Genome Browser, data mining, and bioinformatics tools were used to review enhancers, transcription factors (TFs), histone marks, long non-coding RNAs (lncRNAs), and variants residing in the non-coding region of the ER gene. We report 7 enhancers, 3 of which were rich in TF-binding sites and histone marks in a cell line-specific manner. Furthermore, some enhancers contain estrogen resistance variants and sites for lncRNA. Our review speculates putative models suggesting potential aberrations in gene regulation and expression if these regulatory landscapes and assemblies are altered. This review gives an interesting perspective in designing integrated in vitro studies including non-coding elements to study development of endocrine resistance in ER-positive breast cancer.

DNA Sequence-Dependent Properties of Nucleosome Positioning in Regions of Distinct Chromatin States in Mouse Embryonic Stem Cells

Chromatin architecture is orchestrated, and plays crucial roles during the developmental process by regulating gene expression. In embryonic stem cells (ESCs), three types of chromatin states, including active, repressive and poised states, were previously identified and characterized with specific chromatin modification marks and different transcription activity, but it is largely unknown how nucleosomes are organized in these chromatin states. In this study, by using a DNA deformation energy model, we investigated the sequence-dependent nucleosome organization within the chromatin states in mouse ESCs. The results revealed that: (1) compared with poised genes, active genes are characterized with a higher level of nucleosome occupancy around their transcription start sites (TSS) and transcription termination sites (TTS), and both types of genes do not have a nucleosome-depleted region at their TTS, contrasting with the MNase-seq based result; (2) based on our previous DNA bending energy model, we developed an improved model capable of predicting both rotational positioning and nucleosome occupancy determined by a chemical mapping approach; (3) DNA bending-energy-based analyses demonstrated that the fragile nucleosomes positioned at both gene ends could be explained largely by enhanced rotational positioning signals encoded in DNA, but nucleosome phasing around the TSS of active genes was not determined by sequence preference; (4) the nucleosome occupancy landscape around the binding sites of some developmentally important transcription factors known to bind with different chromatin contexts, was also successfully predicted; (5) the difference of nucleosome occupancy around the TSS between CpG-rich and CpG-poor promoters was partly captured by our sequence-dependent model. Taken together, by developing an improved deformation-energy-based model, we revealed some sequence-dependent properties of the nucleosome arrangements in regions of distinct chromatin states in mouse ESCs.

Computational methods to explore chromatin state dynamics

The human genome is marked by several singular and combinatorial histone modifications that shape the different states of chromatin and its three-dimensional organization. Genome-wide mapping of these marks as well as histone variants and open chromatin regions is commonly carried out via profiling DNA-protein binding or via chromatin accessibility methods. After the generation of epigenomic datasets in a cell type, statistical models can be used to annotate the noncoding regions of DNA and infer the combinatorial histone marks or chromatin states (CS). These methods involve partitioning the genome and labeling individual segments based on their CS patterns. Chromatin labels enable the systematic discovery of genomic function and activity and can label the gene body, promoters or enhancers without using other genomic maps. CSs are dynamic and change under different cell conditions, such as in normal, preneoplastic or tumor cells. This review aims to explore the available computational tools that have been developed to capture CS alterations under two or more cellular conditions.

Differential condensation of sister chromatids acts with Cdc6 to ensure asynchronous S-phase entry in Drosophila male germline stem cell lineage

During Drosophila melanogaster male germline stem cell (GSC) asymmetric division, preexisting old versus newly synthesized histones H3 and H4 are asymmetrically inherited. However, the biological outcomes of this phenomenon have remained unclear. Here, we tracked old and new histones throughout the GSC cell cycle through the use of high spatial and temporal resolution microscopy. We found unique features that differ between old and new histone-enriched sister chromatids, including differences in nucleosome density, chromosomal condensation, and H3 Ser10 phosphorylation. These distinct chromosomal features lead to their differential association with Cdc6, a pre-replication complex component, and subsequent asynchronous DNA replication initiation in the resulting daughter cells. Disruption of asymmetric histone inheritance abolishes differential Cdc6 association and asynchronous S-phase entry, demonstrating that histone asymmetry acts upstream of these critical cell-cycle progression events. Furthermore, disruption of these GSC-specific chromatin features leads to GSC defects, indicating a connection between histone inheritance, cell-cycle progression, and cell fate determination.

Quantitative modelling of fine-scale variations in the Arabidopsis thaliana crossover landscape

In, essentially, all species where meiotic crossovers (COs) have been studied, they occur preferentially in open chromatin, typically near gene promoters and to a lesser extent, at the end of genes. Here, in the case of Arabidopsis thaliana, we unveil further trends arising when one considers contextual information, namely summarised epigenetic status, gene or intergenic region size, and degree of divergence between homologs. For instance, we find that intergenic recombination rate is reduced if those regions are less than 1.5 kb in size. Furthermore, we propose that the presence of single nucleotide polymorphisms enhances the rate of CO formation compared to when homologous sequences are identical, in agreement with previous works comparing rates in adjacent homozygous and heterozygous blocks. Lastly, by integrating these different effects, we produce a quantitative and predictive model of the recombination landscape that reproduces much of the experimental variation.

A unique histone 3 lysine 14 chromatin signature underlies tissue-specific gene regulation

Organismal development and cell differentiation critically depend on chromatin state transitions. However, certain developmentally regulated genes lack histone 3 lysine 9 and 27 acetylation (H3K9ac and H3K27ac, respectively) and histone 3 lysine 4 (H3K4) methylation, histone modifications common to most active genes. Here we describe a chromatin state featuring unique histone 3 lysine 14 acetylation (H3K14ac) peaks in key tissue-specific genes in Drosophila and human cells. Replacing H3K14 in Drosophila demonstrates that H3K14 is essential for expression of genes devoid of canonical histone modifications in the embryonic gut and larval wing imaginal disc, causing lethality and defective wing patterning. We find that the SWI/SNF protein Brahma (Brm) recognizes H3K14ac, that brm acts in the same genetic pathway as H3K14R, and that chromatin accessibility at H3K14ac-unique genes is decreased in H3K14R mutants. Our results show that acetylation of a single lysine is essential at genes devoid of canonical histone marks and uncover an important requirement for H3K14 in tissue-specific gene regulation.

Quantifying the Stability of Coupled Genetic and Epigenetic Switches With Variational Methods

The Waddington landscape provides an intuitive metaphor to view development as a ball rolling down the hill, with distinct phenotypes as basins and differentiation pathways as valleys. Since, at a molecular level, cell differentiation arises from interactions among the genes, a mathematical definition for the Waddington landscape can, in principle, be obtained by studying the gene regulatory networks. For eukaryotes, gene regulation is inextricably and intimately linked to histone modifications. However, the impact of such modifications on both landscape topography and stability of attractor states is not fully understood. In this work, we introduced a minimal kinetic model for gene regulation that combines the impact of both histone modifications and transcription factors. We further developed an approximation scheme based on variational principles to solve the corresponding master equation in a second quantized framework. By analyzing the steady-state solutions at various parameter regimes, we found that histone modification kinetics can significantly alter the behavior of a genetic network, resulting in qualitative changes in gene expression profiles. The emerging epigenetic landscape captures the delicate interplay between transcription factors and histone modifications in driving cell-fate decisions.

Why Do Some Sex Chromosomes Degenerate More Slowly Than Others? The Odd Case of Ratite Sex Chromosomes

The hallmark of sex chromosome evolution is the progressive suppression of recombination which leads to subsequent degeneration of the non-recombining chromosome. In birds, species belonging to the two major clades, Palaeognathae (including tinamous and flightless ratites) and Neognathae (all remaining birds), show distinctive patterns of sex chromosome degeneration. Birds are female heterogametic, in which females have a Z and a W chromosome. In Neognathae, the highly-degenerated W chromosome seems to have followed the expected trajectory of sex chromosome evolution. In contrast, among Palaeognathae, sex chromosomes of ratite birds are largely recombining. The underlying reason for maintenance of recombination between sex chromosomes in ratites is not clear. Degeneration of the W chromosome might have halted or slowed down due to a multitude of reasons ranging from selective processes, such as a less pronounced effect of sexually antagonistic selection, to neutral processes, such as a slower rate of molecular evolution in ratites. The production of genome assemblies and gene expression data for species of Palaeognathae has made it possible, during recent years, to have a closer look at their sex chromosome evolution. Here, we critically evaluate the understanding of the maintenance of recombination in ratites in light of the current data. We conclude by highlighting certain aspects of sex chromosome evolution in ratites that require further research and can potentially increase power for the inference of the unique history of sex chromosome evolution in this lineage of birds.

Changes in H3K27ac at Gene Regulatory Regions in Porcine Alveolar Macrophages Following LPS or PolyIC Exposure

Changes in chromatin structure, especially in histone modifications (HMs), linked with chromatin accessibility for transcription machinery, are considered to play significant roles in transcriptional regulation. Alveolar macrophages (AM) are important immune cells for protection against pulmonary pathogens, and must readily respond to bacteria and viruses that enter the airways. Mechanism(s) controlling AM innate response to different pathogen-associated molecular patterns (PAMPs) are not well defined in pigs. By combining RNA sequencing (RNA-seq) with chromatin immunoprecipitation and sequencing (ChIP-seq) for four histone marks (H3K4me3, H3K4me1, H3K27ac and H3K27me3), we established a chromatin state map for AM stimulated with two different PAMPs, lipopolysaccharide (LPS) and Poly(I:C), and investigated the potential effect of identified histone modifications on transcription factor binding motif (TFBM) prediction and RNA abundance changes in these AM. The integrative analysis suggests that the differential gene expression between non-stimulated and stimulated AM is significantly associated with changes in the H3K27ac level at active regulatory regions. Although global changes in chromatin states were minor after stimulation, we detected chromatin state changes for differentially expressed genes involved in the TLR4, TLR3 and RIG-I signaling pathways. We found that regions marked by H3K27ac genome-wide were enriched for TFBMs of TF that are involved in the inflammatory response. We further documented that TF whose expression was induced by these stimuli had TFBMs enriched within H3K27ac-marked regions whose chromatin state changed by these same stimuli. Given that the dramatic transcriptomic changes and minor chromatin state changes occurred in response to both stimuli, we conclude that regulatory elements (i.e. active promoters) that contain transcription factor binding motifs were already active/poised in AM for immediate inflammatory response to PAMPs. In summary, our data provides the first chromatin state map of porcine AM in response to bacterial and viral PAMPs, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the role of HMs, especially H3K27ac, in regulating transcription in AM in response to LPS and Poly(I:C).

Emergence of a High-Plasticity Cell State during Lung Cancer Evolution

Tumor evolution from a single cell into a malignant, heterogeneous tissue remains poorly understood. Here, we profile single-cell transcriptomes of genetically engineered mouse lung tumors at seven stages, from pre-neoplastic hyperplasia to adenocarcinoma. The diversity of transcriptional states increases over time and is reproducible across tumors and mice. Cancer cells progressively adopt alternate lineage identities, computationally predicted to be mediated through a common transitional, high-plasticity cell state (HPCS). Accordingly, HPCS cells prospectively isolated from mouse tumors and human patient-derived xenografts display high capacity for differentiation and proliferation. The HPCS program is associated with poor survival across human cancers and demonstrates chemoresistance in mice. Our study reveals a central principle underpinning intra-tumoral heterogeneity and motivates therapeutic targeting of the HPCS.

LncRNAs are regulated by chromatin states and affect the skeletal muscle cell differentiation

Objective

This study aims to clarify the mechanisms underlying transcriptional regulation and regulatory roles of lncRNAs in skeletal muscle cell differentiation.

Methods

We analysed the expression patterns of lncRNAs via time‐course RNA‐seq. Then, we further combined the ATAC‐seq and ChIP‐seq to investigate the governing mechanisms of transcriptional regulation of differentially expressed (DE) lncRNAs. Weighted correlation network analysis and GO analysis were conducted to identify the transcription factor (TF)‐lncRNA pairs related to skeletal muscle cell differentiation.

Results

We identified 385 DE lncRNAs during C2C12 differentiation, the transcription of which is determined by chromatin states around their transcriptional start sites. The TF‐lncRNA correlation network showed substantially concordant changes in DE lncRNAs between C2C12 differentiation and satellite cell rapid growth stages. Moreover, the up‐regulated lncRNAs showed a significant decrease following the differentiation capacity of satellite cells, which gradually declines during skeletal muscle development. Notably, inhibition of the lncRNA Atcayos and Trp53cor1 led to the delayed differentiation of satellite cells. Those lncRNAs were significantly up‐regulated during the rapid growth stage of satellite cells (4‐6 weeks) and down‐regulated with reduced differentiation capacity (8‐12 weeks). It confirms that these lncRNAs are positively associated with myogenic differentiation of satellite cells during skeletal muscle development.

Conclusions

This study extends the understanding of mechanisms governing transcriptional regulation of lncRNAs and provides a foundation for exploring their functions in skeletal muscle cell differentiation.

JARID2 and the PRC2 complex regulate the cell cycle in skeletal muscle

JARID2 is a noncatalytic member of the polycomb repressive complex 2 (PRC2) which methylates of histone 3 lysine 27 (H3K27). In this work, we show that JARID2 and the PRC2 complex regulate the cell cycle in skeletal muscle cells to control proliferation and mitotic exit. We found that the stable depletion of JARID2 leads to increased proliferation and cell accumulation in S phase. The regulation of the cell cycle by JARID2 is mediated by direct repression of both cyclin D1 and cyclin E1, both of which are targets of PRC2-mediated H3K27 methylation. Intriguingly, we also find that the retinoblastoma protein (RB1) is a direct target of JARID2 and the PRC2 complex. The depletion of JARID2 is not sufficient to activate RB1. However, the ectopic expression of RB1 can suppress cyclin D1 expression in JARID2-depleted cells. Transient depletion of JARID2 in skeletal muscle cells leads to a transient up-regulation of cyclin D1 that is quickly suppressed with no resulting effect on proliferation, Taken together, we show that JARID2 and the PRC2 complex regulate skeletal muscle proliferation in a precise manner that involves the repression of cyclin D1, thus restraining proliferation and repressing RB1, which is required for mitotic exit and terminal differentiation.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

Chromatin State-Based Analysis of Epigenetic H3K4me3 Marks of Arabidopsis in Response to Dark Stress

Light is essential to plant growth and development. Extended darkness causes dramatic gene expression changes, leading to leaf senescence, hypocotyl growth, petiole elongation, reduced leaf area, and early flowering, etc. However, the underlying mechanism of response to darkness at epigenetic levels remains largely unknown. In this study, we conducted ChIP-seq to generate global epigenomic profiles of H3K4me3 under 3-day extended darkness and normal light conditions in Arabidopsis. We applied chromatin state analysis together with self-organization mapping (SOM) to study the combination of epigenetic regulation under dark stress. The SOM map clusters the segments on the genome according to multiple diverse epigenomic datasets, which breaks the limit of dispersed distribution of epigenetic marks on the genome. Through SOM analysis, we also found that the signals of H3K4me3 were mainly increased after darkness. Analysis of H3K4me3-changed genes together with differentially expressed genes indicated that the genes showing dark-increased H3K4me3 were most involved in senescence and autophagy, and cross-talk existed between dark-induced and natural senescence. In summary, we studied the regulation of the epigenetic H3K4me3 marks of Arabidopsis in response to dark stress using chromatin state and SOM analyses. Our study revealed the regulatory mechanisms of the epigenome in response to dark stress, and SOM analysis based on chromatin states used in our study will also be helpful for other studies on dynamic changes of multiple epigenetic marks.

Decoding genetic and epigenetic information embedded in cell free DNA with adapted SALP-seq

Cell free DNA (cfDNA) in human plasma carries abundant information, which has therefore been the key sample for noninvasive prenatal testing (NIPT) and liquid biopsy. Especially by using the rapidly developed next-generation sequencing (NGS) techniques, the genetic and epigenetic information embedded in cfDNA has been effectively and extensively decoded. In this process, a key step is to construct the NGS library. Due to its high degradation, the single strand-based method was reported to be more qualified to construct the NGS library of cfDNA. In order to develop a new simple method for this application, this study adapted our recently developed single strand adaptor library preparation (SALP) method for constructing an NGS library of cfDNA. In the improved method, cfDNA was firstly denatured into single strands and then ligated with a single strand adaptor (SSA) that had a 3' overhang of 3 random bases by using T4 DNA ligase. The SSA-ligated DNA was converted into double-stranded DNA with an additional adenine at the other end by polymerizing with Taq polymerase. Next, a barcode T adaptor (BTA) was ligated to this end. Finally, the cfDNA ligated with two adaptors was amplified with the Illumina-compatible primers for NGS. Using the method, this study successfully sequenced 20 cfDNA samples from 16 esophageal cancer patients and 4 healthy people. By bioinformatics analysis, this study found the genetic and epigenetic difference between patients and healthy people and identified 23 epigenetic and 28 genetic altered esophageal cancer-specific genes.

A Pan-Cancer Analysis of Enhancer Expression in Nearly 9000 Patient Samples

The role of enhancers, a key class of non-coding regulatory DNA elements, in cancer development has increasingly been appreciated. Here, we present the detection and characterization of a large number of expressed enhancers in a genome-wide analysis of 8928 tumor samples across 33 cancer types using TCGA RNA-seq data. Compared with matched normal tissues, global enhancer activation was observed in most cancers. Across cancer types, global enhancer activity was positively associated with aneuploidy, but not mutation load, suggesting a hypothesis centered on "chromatin-state" to explain their interplay. Integrating eQTL, mRNA co-expression, and Hi-C data analysis, we developed a computational method to infer causal enhancer-gene interactions, revealing enhancers of clinically actionable genes. Having identified an enhancer ∼140 kb downstream of PD-L1, a major immunotherapy target, we validated it experimentally. This study provides a systematic view of enhancer activity in diverse tumor contexts and suggests the clinical implications of enhancers.

AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

As an economically important crop, apple is one of the most cultivated fruit trees in temperate regions worldwide. Recently, a large number of high-quality transcriptomic and epigenomic datasets for apple were made available to the public, which could be helpful in inferring gene regulatory relationships and thus predicting gene function at the genome level. Through integration of the available apple genomic, transcriptomic, and epigenomic datasets, we constructed co-expression networks, identified functional modules, and predicted chromatin states. A total of 112 RNA-seq datasets were integrated to construct a global network and a conditional network (tissue-preferential network). Furthermore, a total of 1,076 functional modules with closely related gene sets were identified to assess the modularity of biological networks and further subjected to functional enrichment analysis. The results showed that the function of many modules was related to development, secondary metabolism, hormone response, and transcriptional regulation. Transcriptional regulation is closely related to epigenetic marks on chromatin. A total of 20 epigenomic datasets, which included ChIP-seq, DNase-seq, and DNA methylation analysis datasets, were integrated and used to classify chromatin states. Based on the ChromHMM algorithm, the genome was divided into 620,122 fragments, which were classified into 24 states according to the combination of epigenetic marks and enriched-feature regions. Finally, through the collaborative analysis of different omics datasets, the online database AppleMDO (http://bioinformatics.cau.edu.cn/AppleMDO/) was established for cross-referencing and the exploration of possible novel functions of apple genes. In addition, gene annotation information and functional support toolkits were also provided. Our database might be convenient for researchers to develop insights into the function of genes related to important agronomic traits and might serve as a reference for other fruit trees.

Virus DNA Replication and the Host DNA Damage Response

Viral DNA genomes have limited coding capacity and therefore harness cellular factors to facilitate replication of their genomes and generate progeny virions. Studies of viruses and how they interact with cellular processes have historically provided seminal insights into basic biology and disease mechanisms. The replicative life cycles of many DNA viruses have been shown to engage components of the host DNA damage and repair machinery. Viruses have evolved numerous strategies to navigate the cellular DNA damage response. By hijacking and manipulating cellular replication and repair processes, DNA viruses can selectively harness or abrogate distinct components of the cellular machinery to complete their life cycles. Here, we highlight consequences for viral replication and host genome integrity during the dynamic interactions between virus and host.

Virus DNA Replication and the Host DNA Damage Response

Viral DNA genomes have limited coding capacity and therefore harness cellular factors to facilitate replication of their genomes and generate progeny virions. Studies of viruses and how they interact with cellular processes have historically provided seminal insights into basic biology and disease mechanisms. The replicative life cycles of many DNA viruses have been shown to engage components of the host DNA damage and repair machinery. Viruses have evolved numerous strategies to navigate the cellular DNA damage response. By hijacking and manipulating cellular replication and repair processes, DNA viruses can selectively harness or abrogate distinct components of the cellular machinery to complete their life cycles. Here, we highlight consequences for viral replication and host genome integrity during the dynamic interactions between virus and host.

Trithorax-group proteins ARABIDOPSIS TRITHORAX4 (ATX4) and ATX5 function in abscisic acid and dehydration stress responses

Trithorax-group proteins (TrxGs) play essential regulatory roles in chromatin modification to activate transcription. Although TrxGs have been shown to be extensively involved in the activation of developmental genes, how the specific TrxGs function in the dehydration and abscisic acid (ABA)-mediated modulation of downstream gene expression remains unknown. Here, we report that two evolutionarily conserved Arabidopsis thaliana TrxGs, ARABIDOPSIS TRITHORAX4 (ATX4) and ATX5, play essential roles in the drought stress response. atx4 and atx5 single loss-of-function mutants showed drought stress-tolerant and ABA-hypersensitive phenotypes during seed germination and seedling development, while the atx4 atx5 double mutant displayed further exacerbation of the phenotypes. Genome-wide RNA-sequencing analyses showed that ATX4 and ATX5 regulate the expression of genes functioning in dehydration stress. Intriguingly, ABA-HYPERSENSITIVE GERMINATION 3 (AHG3), an essential negative regulator of ABA signaling, acts genetically downstream of ATX4 and ATX5 in response to ABA. ATX4 and ATX5 directly bind to the AHG3 locus and trimethylate histone H3 of Lys 4 (H3K4). Moreover, ATX4 and ATX5 occupancies at AHG3 are dramatically increased under ABA treatment, and are also essential for RNA polymerase II (RNAPII) occupancies. Our findings reveal novel molecular functions of A. thaliana TrxGs in dehydration stress and ABA responses.

Systematic Epigenomic Analysis Reveals Chromatin States Associated with Melanoma Progression

The extent and nature of epigenomic changes associated with melanoma progression is poorly understood. Through systematic epigenomic profiling of 35 epigenetic modifications and transcriptomic analysis, we define chromatin state changes associated with melanomagenesis by using a cell phenotypic model of non-tumorigenic and tumorigenic states. Computation of specific chromatin state transitions showed loss of histone acetylations and H3K4me2/3 on regulatory regions proximal to specific cancer-regulatory genes in important melanoma-driving cell signaling pathways. Importantly, such acetylation changes were also observed between benign nevi and malignant melanoma human tissues. Intriguingly, only a small fraction of chromatin state transitions correlated with expected changes in gene expression patterns. Restoration of acetylation levels on deacetylated loci by histone deacetylase (HDAC) inhibitors selectively blocked excessive proliferation in tumorigenic cells and human melanoma cells, suggesting functional roles of observed chromatin state transitions in driving hyperproliferative phenotype. Through these results, we define functionally relevant chromatin states associated with melanoma progression.

Distal Limb Patterning Requires Modulation of cis-Regulatory Activities by HOX13

The combinatorial expression of Hox genes along the body axes is a major determinant of cell fate and plays a pivotal role in generating the animal body plan. Loss of HOXA13 and HOXD13 transcription factors (HOX13) leads to digit agenesis in mice, but how HOX13 proteins regulate transcriptional outcomes and confer identity to the distal-most limb cells has remained elusive. Here, we report on the genome-wide profiling of HOXA13 and HOXD13 in vivo binding and changes of the transcriptome and chromatin state in the transition from the early to the late-distal limb developmental program, as well as in Hoxa13^−/−; Hoxd13^−/−limbs. Our results show that proper termination of the early limb transcriptional program and activation of the late-distal limb program are coordinated by the dual action of HOX13 on cis-regulatory modules.

Local Chromatin Features Including PU.1 and IKAROS Binding and H3K4 Methylation Shape the Repertoire of Immunoglobulin Kappa Genes Chosen for V(D)J Recombination

V(D)J recombination is essential for the generation of diverse antigen receptor (AgR) repertoires. In B cells, immunoglobulin kappa (Igκ) light chain recombination follows immunoglobulin heavy chain (Igh) recombination. We recently developed the DNA-based VDJ-seq assay for the unbiased quantitation of Igh VH and DH repertoires. Integration of VDJ-seq data with genome-wide datasets revealed that two chromatin states at the recombination signal sequence (RSS) of VH genes are highly predictive of recombination in mouse pro-B cells. It is unknown whether local chromatin states contribute to Vκ gene choice during Igκ recombination. Here we adapt VDJ-seq to profile the Igκ VκJκ repertoire and present a comprehensive readout in mouse pre-B cells, revealing highly variable Vκ gene usage. Integration with genome-wide datasets for histone modifications, DNase hypersensitivity, transcription factor binding and germline transcription identified PU.1 binding at the RSS, which was unimportant for Igh, as highly predictive of whether a Vκ gene will recombine or not, suggesting that it plays a binary, all-or-nothing role, priming genes for recombination. Thereafter, the frequency with which these genes recombine was shaped both by the presence and level of enrichment of several other chromatin features, including H3K4 methylation and IKAROS binding. Moreover, in contrast to the Igh locus, the chromatin landscape of the promoter, as well as of the RSS, contributes to Vκ gene recombination. Thus, multiple facets of local chromatin features explain much of the variation in Vκ gene usage. Together, these findings reveal shared and divergent roles for epigenetic features and transcription factors in AgR V(D)J recombination and provide avenues for further investigation of chromatin signatures that may underpin V(D)J-mediated chromosomal translocations.

Defining Functional Genic Regions in the Human Genome through Integration of Biochemical, Evolutionary, and Genetic Evidence

The human genome is dominated by large tracts of DNA with extensive biochemical activity but no known function. In particular, it is well established that transcriptional activities are not restricted to known genes. However, whether this intergenic transcription represents activity with functional significance or noise is under debate, highlighting the need for an effective method of defining functional genomic regions. Moreover, these discoveries raise the question whether genomic regions can be defined as functional based solely on the presence of biochemical activities, without considering evolutionary (conservation) and genetic (effects of mutations) evidence. Here, computational models integrating genetic, evolutionary, and biochemical evidence are established that provide reliable predictions of human protein-coding and RNA genes. Importantly, in addition to sequence conservation, biochemical features allow accurate predictions of genic sequences with phenotypic evidence under strong purifying selection, suggesting that they can be used as an alternative measure of selection. Moreover, 18.5% of annotated noncoding RNAs exhibit higher degrees of similarity to phenotype genes and, thus, are likely functional. However, 64.5% of noncoding RNAs appear to belong to a sequence class of their own, and the remaining 17% are more similar to pseudogenes and random intergenic sequences that may represent noisy transcription.

Inference of cell type specific regulatory networks on mammalian lineages

Transcriptional regulatory networks are at the core of establishing cell type specific gene expression programs. In mammalian systems, such regulatory networks are determined by multiple levels of regulation, including by transcription factors, chromatin environment, and three-dimensional organization of the genome. Recent efforts to measure diverse regulatory genomic datasets across multiple cell types and tissues offer unprecedented opportunities to examine the context-specificity and dynamics of regulatory networks at a greater resolution and scale than before. In parallel, numerous computational approaches to analyze these data have emerged that serve as important tools for understanding mammalian cell type specific regulation. In this article, we review recent computational approaches to predict the expression and sequence-based regulators of a gene's expression level and examine long-range gene regulation. We highlight promising approaches, insights gained, and open challenges that need to be overcome to build a comprehensive picture of cell type specific transcriptional regulatory networks.

Chromatin signatures at Notch-regulated enhancers reveal large-scale changes in H3K56ac upon activation

The conserved Notch pathway functions in diverse developmental and disease-related processes, requiring mechanisms to ensure appropriate target selection and gene activation in each context. To investigate the influence of chromatin organisation and dynamics on the response to Notch signalling, we partitioned Drosophila chromatin using histone modifications and established the preferred chromatin conditions for binding of Su(H), the Notch pathway transcription factor. By manipulating activity of a co-operating factor, Lozenge/Runx, we showed that it can help facilitate these conditions. While many histone modifications were unchanged by Su(H) binding or Notch activation, we detected rapid changes in acetylation of H3K56 at Notch-regulated enhancers. This modification extended over large regions, required the histone acetyl-transferase CBP and was independent of transcription. Such rapid changes in H3K56 acetylation appear to be a conserved indicator of enhancer activation as they also occurred at the mammalian Notch-regulated Hey1 gene and at Drosophila ecdysone-regulated genes. This intriguing example of a core histone modification increasing over short timescales may therefore underpin changes in chromatin accessibility needed to promote transcription following signalling activation.

Functional chromatin features are associated with structural mutations in cancer

BACKGROUND

Structural mutations (SMs) play a major role in cancer development. In some cancers, such as breast and ovarian, DNA double-strand breaks (DSBs) occur more frequently in transcribed regions, while in other cancer types such as prostate, there is a consistent depletion of breakpoints in transcribed regions. Despite such regularity, little is understood about the mechanisms driving these effects. A few works have suggested that protein binding may be relevant, e.g. in studies of androgen receptor binding and active chromatin in specific cell types. We hypothesized that this behavior might be general, i.e. that correlation between protein-DNA binding (and open chromatin) and breakpoint locations is common across divergent cancers.

RESULTS

We investigated this hypothesis by comprehensively analyzing the relationship among 457 ENCODE protein binding ChIP-seq experiments, 125 DnaseI and 24 FAIRE experiments, and 14,600 SMs from 8 diverse cancer datasets covering 147 samples. In most cancers, including breast and ovarian, we found enrichment of protein binding and open chromatin in the vicinity of SM breakpoints at distances up to 200 kb. Furthermore, for all cancer types we observed an enhanced enrichment in regions distant from genes when compared to regions proximal to genes, suggesting that the SM-induction mechanism is independent from the bias of DSBs to occur near transcribed regions. We also observed a stronger effect for sites with more than one protein bound.

CONCLUSIONS

Protein binding and open chromatin state are associated with nearby SM breakpoints in many cancer datasets. These observations suggest a consistent mechanism underlying SM locations across different cancers.

Computational and experimental methods to decipher the epigenetic code

A multi-layered set of epigenetic marks, including post-translational modifications of histones and methylation of DNA, is finely tuned to define the epigenetic state of chromatin in any given cell type under specific conditions. Recently, the knowledge about the combinations of epigenetic marks occurring in the genome of different cell types under various conditions is rapidly increasing. Computational methods were developed for the identification of these states, unraveling the combinatorial nature of epigenetic marks and their association to genomic functional elements and transcriptional states. Nevertheless, the precise rules defining the interplay between all these marks remain poorly characterized. In this perspective we review the current state of this research field, illustrating the power and the limitations of current approaches. Finally, we sketch future avenues of research illustrating how the adoption of specific experimental designs coupled with available experimental approaches could be critical for a significant progress in this area.

Exploring the function of genetic variants in the non-coding genomic regions: approaches for identifying human regulatory variants affecting gene expression

Understanding the genetic basis of human traits/diseases and the underlying mechanisms of how these traits/diseases are affected by genetic variations is critical for public health. Current genome-wide functional genomics data uncovered a large number of functional elements in the noncoding regions of human genome, providing new opportunities to study regulatory variants (RVs). RVs play important roles in transcription factor bindings, chromatin states and epigenetic modifications. Here, we systematically review an array of methods currently used to map RVs as well as the computational approaches in annotating and interpreting their regulatory effects, with emphasis on regulatory single-nucleotide polymorphism. We also briefly introduce experimental methods to validate these functional RVs.

Regulation of cellular chromatin state: insights from quiescence and differentiation

The identity and functionality of eukaryotic cells is defined not just by their genomic sequence which remains constant between cell types, but by their gene expression profiles governed by epigenetic mechanisms. Epigenetic controls maintain and change the chromatin state throughout development, as exemplified by the setting up of cellular memory for the regulation and maintenance of homeotic genes in proliferating progenitors during embryonic development. Higher order chromatin structure in reversibly arrested adult stem cells also involves epigenetic regulation and in this review we highlight common trends governing chromatin states, focusing on quiescence and differentiation during myogenesis. Together, these diverse developmental modules reveal the dynamic nature of chromatin regulation providing fresh insights into the role of epigenetic mechanisms in potentiating development and differentiation.