Identifying gene regulatory elements by genome-wide recovery of DNase hypersensitive sites

Role of gene therapy in treatment of cancer with craniofacial regeneration-current molecular strategies, future perspectives, and challenges: a narrative review

Gene therapy involves the introduction of foreign genetic material into host tissue to alter the expression of genetic products. Gene therapy represents an opportunity to alter the course of various diseases. Hence, genetic products utilizing safe and reliable vectors with improved biotechnology will play a critical role in the treatment of various diseases in the future. This review summarizes various important vectors for gene therapy along with modern techniques for potential craniofacial regeneration using gene therapy. This review also explains current molecular approaches for the management and treatment of cancer using gene therapy. The existing literature was searched to find studies related to gene therapy and its role in craniofacial regeneration and cancer treatment. Various databases such as PubMed, Science Direct, Scopus, Web of Science, and Google Scholar were searched for English language articles using the keywords "gene therapy," "gene therapy in present scenario," "gene therapy in cancer," "gene therapy and vector," "gene therapy in diseases," and "gene therapy and molecular strategies."

Genome-wide mapping of DNase I hypersensitive sites in pineapple leaves

Pineapple [Ananas comosus (L.) Merr.] is the most economically important crop possessing crassulacean acid metabolism (CAM) photosynthesis which has a higher water use efficiency by control of nocturnal opening and diurnal closure of stomata. To provide novel insights into the diel regulatory landscape in pineapple leaves, we performed genome-wide mapping of DNase I hypersensitive sites (DHSs) in pineapple leaves at day (2a.m.) and night (10a.m.) using a simplified DNase-seq method. As a result, totally 33340 and 28753 DHSs were found in green-tip tissue, and 29597 and 40068 were identified in white-base tissue at 2a.m. and 10a.m., respectively. We observed that majority of the pineapple genes occupied less than two DHSs with length shorter than 1 kb, and the promotor DHSs showed a proximal trend to the transcription start site (>77% promotor DHSs within 1 kb). In addition, more intergenic DHSs were identified around transcription factors or transcription co-regulators (TFs/TCs) than other functional genes, indicating complex regulatory contexts around TFs/TCs. Through combined analysis of tissue preferential DHSs and genes, we respectively found 839 and 888 coordinately changed genes in green-tip at 2a.m. and 10a.m. (AcG2 and AcG10). Furthermore, AcG2-specific, AcG10-specific and common accessible DHSs were dissected from the total photosynthetic preferential DHSs, and the regulatory networks indicated dynamic regulations with multiple cis-regulatory elements occurred to genes preferentially expressed in photosynthetic tissues. Interestingly, binding motifs of several cycling TFs were identified in the DHSs of key CAM genes, revealing a circadian regulation to CAM coordinately diurnal expression. Our results provide a chromatin regulatory landscape in pineapple leaves during the day and night. This will provide important information to assist with deciphering the circadian regulation of CAM photosynthesis.

B lymphocytes in treatment-naive paediatric patients with lupus are epigenetically distinct from healthy children

BACKGROUND

SLE is likely triggered by gene-environment interactions. We have shown that most SLE-associated haplotypes encompass genomic regions enriched for epigenetic marks associated with enhancer function in lymphocytes, suggesting genetic risk is exerted through altered gene regulation. Data remain scarce on how epigenetic variance contributes to disease risk in paediatric SLE (pSLE). We aim to identify differences in epigenetically regulated chromatin architecture in treatment-naive patients with pSLE compared with healthy children.

METHODS

Using the assay for transposase-accessible chromatin with sequencing (ATACseq), we surveyed open chromatin in 10 treatment-naive patients with pSLE, with at least moderate disease severity, and 5 healthy children. We investigated whether regions of open chromatin unique to patients with pSLE demonstrate enrichment for specific transcriptional regulators, using standard computational approaches to identify unique peaks and a false discovery rate of <0.05. Further analyses for histone modification enrichment and variant calling were performed using bioinformatics packages in R and Linux.

RESULTS

We identified 30 139 differentially accessible regions (DAR) unique to pSLE B cells; 64.3% are more accessible in pSLE than healthy children. Many DAR are found in distal, intergenic regions and enriched for enhancer histone marks (p=0.027). B cells from adult patients with SLE contain more regions of inaccessible chromatin than those in pSLE. In pSLE B cells, 65.2% of the DAR are located within or near known SLE haplotypes. Further analysis revealed enrichment of transcription factor binding motifs within these DAR that may regulate genes involved in pro-inflammatory responses and cellular adhesion.

CONCLUSIONS

We demonstrate an epigenetically distinct profile in pSLE B cells when compared with healthy children and adults with lupus, indicating that pSLE B cells are predisposed for disease onset/development. Increased chromatin accessibility in non-coding genomic regions controlling activation of inflammation suggest that transcriptional dysregulation by regulatory elements controlling B cell activation plays an important role in pSLE pathogenesis.

Enhancer Clusters Drive Type I Interferon-Induced TRAIL Overexpression in Cancer, and Its Intracellular Protein Accumulation Fails to Induce Apoptosis

Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) is a cytokine produced and secreted by immune cells in response to an infection, often in response to interferon (IFN) stimulation. In cancer, it has also been shown that IFN stimulates the production of TRAIL, and it has been proposed that this TRAIL can induce apoptosis in an autocrine or paracrine manner in different cancer cells. Yet, the mechanism mediating TRAIL upregulation and the implications of TRAIL as an apoptotic molecule in cancer cells are still poorly understood. We show here that in certain cancer cells, TRAIL is upregulated by enhancer clusters, potent genomic regulatory regions containing densely packed enhancers that have combinatorial and additive activity and that are usually found to be associated with cancer-promoting genes. Moreover, we found that TRAIL upregulation by IFNα is mediated by these enhancer clusters in breast and lung cancer cells. Surprisingly, IFNα stimulation leads to the intracellular accumulation of TRAIL protein in these cancer cells. Consequently, this TRAIL is not capable of inducing apoptosis. Our study provides novel insights into the mechanism behind the interferon-mediated upregulation of TRAIL and its protein accumulation in cancer cells. Further investigation is required to understand the role of intracellular TRAIL or depict the mechanisms mediating its apoptosis impairment in cancer cells.

LangMoDHS: A deep learning language model for predicting DNase I hypersensitive sites in mouse genome

DNase I hypersensitive sites (DHSs) are a specific genomic region, which is critical to detect or understand cis-regulatory elements. Although there are many methods developed to detect DHSs, there is a big gap in practice. We presented a deep learning-based language model for predicting DHSs, named LangMoDHS. The LangMoDHS mainly comprised the convolutional neural network (CNN), the bi-directional long short-term memory (Bi-LSTM) and the feed-forward attention. The CNN and the Bi-LSTM were stacked in a parallel manner, which was helpful to accumulate multiple-view representations from primary DNA sequences. We conducted 5-fold cross-validations and independent tests over 14 tissues and 4 developmental stages. The empirical experiments showed that the LangMoDHS is competitive with or slightly better than the iDHS-Deep, which is the latest method for predicting DHSs. The empirical experiments also implied substantial contribution of the CNN, Bi-LSTM, and attention to DHSs prediction. We implemented the LangMoDHS as a user-friendly web server which is accessible at http:/www.biolscience.cn/LangMoDHS/. We used indices related to information entropy to explore the sequence motif of DHSs. The analysis provided a certain insight into the DHSs.

Enhancers as potential targets for engineering salinity stress tolerance in crop plants

Enhancers represent noncoding regulatory regions of the genome located distantly from their target genes. They regulate gene expression programs in a context-specific manner via interacting with promoters of one or more target genes and are generally associated with transcription factor binding sites and epi(genomic)/chromatin features, such as regions of chromatin accessibility and histone modifications. The enhancers are difficult to identify due to the modularity of their associated features. Although enhancers have been studied extensively in human and animals, only a handful of them has been identified in few plant species till date due to nonavailability of plant-specific experimental and computational approaches for their discovery. Being an important regulatory component of the genome, enhancers represent potential targets for engineering agronomic traits, including salinity stress tolerance in plants. Here, we provide a review of the available experimental and computational approaches along with the associated sequence and chromatin/epigenetic features for the discovery of enhancers in plants. In addition, we provide insights into the challenges and future prospects of enhancer research in plant biology with emphasis on potential applications in engineering salinity stress tolerance in crop plants.

RGD engineered dendrimer nanotherapeutic as an emerging targeted approach in cancer therapy

A bird's eye view is now demanded in the area of cancer research to suppress the suffering of cancer patient and mediate the lack of treatment related to chemotherapy. Chemotherapy is always preferred over surgery or radiation therapy, but they never met the patient's demand of safe medication. Targeted therapy has now been in research that could hinder the unnecessary effect of drug on normal cells but could affect the tumor cells in much efficient manner. Angiogenesis is process involved in development of new blood vessel that nourishes tumor growth. Integrin receptors are over expressed on cancer cells that play vital role in angiogenesis for growth and metastasis of tumor cell. A delivery of RGD based peptide to integrin targeted site could help in its successful binding and liberation of drug in tumor vasculature. Dendrimers, in addition to its excellent pharmacokinetic properties also helps to carry targeting ligand to site of tumor by successfully conjugating with them. The aim of this review is to bring light upon the role of integrin in cancer progression, interaction of RGD to integrin receptor and more importantly the RGD-dendrimer based targeted therapy for the treatment of various cancers.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

Genome-Wide Identification of DNase I Hypersensitive Sites in Plants

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

Simplified Epigenome Profiling Using Antibody-tethered Tagmentation

We previously introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling method in which antibody tethering of the Tn5 transposase to a chromatin epitope of interest maps specific chromatin features in small samples and single cells. With CUT&Tag, intact cells or nuclei are permeabilized, followed by successive addition of a primary antibody, a secondary antibody, and a chimeric Protein A-Transposase fusion protein that binds to the antibody. Addition of Mg⁺⁺ activates the transposase and inserts sequencing adapters into adjacent DNA in situ. We have since adapted CUT&Tag to also map chromatin accessibility by simply modifying the transposase activation conditions when using histone H3K4me2, H3K4me3, or Serine-5-phosphorylated RNA Polymerase II antibodies. Using these antibodies, we redirect the tagmentation of accessible DNA sites to produce chromatin accessibility maps with exceptionally high signal-to-noise and resolution. All steps from nuclei to amplified sequencing-ready libraries are performed in single PCR tubes using non-toxic reagents and inexpensive equipment, making our simplified strategy for simultaneous chromatin profiling and accessibility mapping suitable for the lab, home workbench, or classroom.

The era of gene therapy: From preclinical development to clinical application

Three decades of promise have culminated in the development of gene therapies that can be applied to a broad range of human diseases. After a brief history, we provide an overview of gene therapy types and delivery methods, gene editing technologies, regulatory affairs, clinical trials, approved products, ongoing challenges, and future goals. Information on clinical trials of candidates and on approved products for gene therapy developed between 1988 and 2020 is systematically collated. To obtain this global information, we scanned and reviewed more than 46,000 records of clinical trials from 17 clinical trial database providers. The medical benefits of transformative gene therapies are gradually being accepted by payors, and a significant increase in the number of gene therapy clinical trials and approved gene therapy products has resulted.

Multi-Dimensional Gene Regulation in Innate and Adaptive Lymphocytes: A View From Regulomes

The precise control of cytokine production by innate lymphoid cells (ILCs) and their T cell adaptive system counterparts is critical to mounting a proper host defense immune response without inducing collateral damage and autoimmunity. Unlike T cells that differentiate into functionally divergent subsets upon antigen recognition, ILCs are developmentally programmed to rapidly respond to environmental signals in a polarized manner, without the need of T cell receptor (TCR) signaling. The specification of cytokine production relies on dynamic regulation of cis-regulatory elements that involve multi-dimensional epigenetic mechanisms, including DNA methylation, transcription factor binding, histone modification and DNA-DNA interactions that form chromatin loops. How these different layers of gene regulation coordinate with each other to fine tune cytokine production, and whether ILCs and their T cell analogs utilize the same regulatory strategy, remain largely unknown. Herein, we review the molecular mechanisms that underlie cell identity and functionality of helper T cells and ILCs, focusing on networks of transcription factors and cis-regulatory elements. We discuss how higher-order chromatin architecture orchestrates these components to construct lineage- and state-specific regulomes that support ordered immunoregulation.

iDHS-Deep: an integrated tool for predicting DNase I hypersensitive sites by deep neural network

DNase I hypersensitive site (DHS) refers to the hypersensitive region of chromatin for the DNase I enzyme. It is an important part of the noncoding region and contains a variety of regulatory elements, such as promoter, enhancer, and transcription factor-binding site, etc. Moreover, the related locus of disease (or trait) are usually enriched in the DHS regions. Therefore, the detection of DHS region is of great significance. In this study, we develop a deep learning-based algorithm to identify whether an unknown sequence region would be potential DHS. The proposed method showed high prediction performance on both training datasets and independent datasets in different cell types and developmental stages, demonstrating that the method has excellent superiority in the identification of DHSs. Furthermore, for the convenience of related wet-experimental researchers, the user-friendly web-server iDHS-Deep was established at http://lin-group.cn/server/iDHS-Deep/, by which users can easily distinguish DHS and non-DHS and obtain the corresponding developmental stage ofDHS.

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA ('tagmentation'). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here, we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

Enhancer Predictions and Genome-Wide Regulatory Circuits

Spatiotemporal control of gene expression during development requires orchestrated activities of numerous enhancers, which are cis-regulatory DNA sequences that, when bound by transcription factors, support selective activation or repression of associated genes. Proper activation of enhancers is critical during embryonic development, adult tissue homeostasis, and regeneration, and inappropriate enhancer activity is often associated with pathological conditions such as cancer. Multiple consortia [e.g., the Encyclopedia of DNA Elements (ENCODE) Consortium and National Institutes of Health Roadmap Epigenomics Mapping Consortium] and independent investigators have mapped putative regulatory regions in a large number of cell types and tissues, but the sequence determinants of cell-specific enhancers are not yet fully understood. Machine learning approaches trained on large sets of these regulatory regions can identify core transcription factor binding sites and generate quantitative predictions of enhancer activity and the impact of sequence variants on activity. Here, we review these computational methods in the context of enhancer prediction and gene regulatory network models specifying cell fate.

Abiotic stress-mediated modulation of the chromatin landscape in Arabidopsis thaliana

Limited information is available on abiotic stress-mediated alterations of chromatin conformation influencing gene expression in plants. In order to characterize the effect of abiotic stresses on changes in chromatin conformation, we employed FAIRE-seq (formaldehyde-assisted isolation of regulatory element sequencing) and DNase-seq to isolate accessible regions of chromatin from Arabidopsis thaliana seedlings exposed to either heat, cold, salt, or drought stress. Approximately 25% of regions in the Arabidopsis genome were captured as open chromatin, the majority of which included promoters and exons. A large proportion of chromatin regions apparently did not change their conformation in response to any of the four stresses. Digital footprints present within these regions had differential enrichment of motifs for binding of 43 different transcription factors. Further, in contrast to drought and salt stress, both high and low temperature treatments resulted in increased accessibility of the chromatin. Also, pseudogenes attained increased chromatin accessibility in response to cold and drought stresses. The highly accessible and inaccessible chromatin regions of seedlings exposed to drought stress correlated with the Ser/Thr protein kinases (MLK1 and MLK2)-mediated reduction and increase in H3 phosphorylation (H3T3Ph), respectively. The presented results provide a deeper understanding of abiotic stress-mediated chromatin modulation in plants.

lentiMPRA and MPRAflow for high-throughput functional characterization of gene regulatory elements

Massively parallel reporter assays (MPRAs) can simultaneously measure the function of thousands of candidate regulatory sequences (CRSs) in a quantitative manner. In this method, CRSs are cloned upstream of a minimal promoter and reporter gene, alongside a unique barcode, and introduced into cells. If the CRS is a functional regulatory element, it will lead to the transcription of the barcode sequence, which is measured via RNA sequencing and normalized for cellular integration via DNA sequencing of the barcode. This technology has been used to test thousands of sequences and their variants for regulatory activity, to decipher the regulatory code and its evolution, and to develop genetic switches. Lentivirus-based MPRA (lentiMPRA) produces 'in-genome' readouts and enables the use of this technique in hard-to-transfect cells. Here, we provide a detailed protocol for lentiMPRA, along with a user-friendly Nextflow-based computational pipeline-MPRAflow-for quantifying CRS activity from different MPRA designs. The lentiMPRA protocol takes ~2 months, which includes sequencing turnaround time and data processing with MPRAflow.

Index and biological spectrum of human DNase I hypersensitive sites

DNase I hypersensitive sites (DHSs) are generic markers of regulatory DNA^1–5 and contain genetic variations associated with diseases and phenotypic traits^6–8. We created high-resolution maps of DHSs from 733 human biosamples encompassing 438 cell and tissue types and states, and integrated these to delineate and numerically index approximately 3.6 million DHSs within the human genome sequence, providing a common coordinate system for regulatory DNA. Here we show that these maps highly resolve the cis-regulatory compartment of the human genome, which encodes unexpectedly diverse cell- and tissue-selective regulatory programs at very high density. These programs can be captured comprehensively by a simple vocabulary that enables the assignment to each DHS of a regulatory barcode that encapsulates its tissue manifestations, and global annotation of protein-coding and non-coding RNA genes in a manner orthogonal to gene expression. Finally, we show that sharply resolved DHSs markedly enhance the genetic association and heritability signals of diseases and traits. Rather than being confined to a small number of distal elements or promoters, we find that genetic signals converge on congruently regulated sets of DHSs that decorate entire gene bodies. Together, our results create a universal, extensible coordinate system and vocabulary for human regulatory DNA marked by DHSs, and provide a new global perspective on the architecture of human gene regulation.

High-resolution maps of DNase I hypersensitive sites from 733 human biosamples are used to identify and index regulatory DNA within the human genome.

Plasma Membrane Proteomic Profile Discovers Macrophage-capping Protein Related to Latent HIV-1

BACKGROUND

Due to the persistence of latent HIV-infected cellular reservoirs, HIV virus can not be eradicated completely.

OBJECTIVE

To identify proteins related to HIV latency, we performed a subcellular proteomic study in HIV latent cell lines.

METHODS

An established HIV-1 latent cell model (J-Lat Tat-GFP Clone A7 cells, A7 cells) and its parental cell line (Jurkat cells) were used. The plasma membrane (PM) fraction from cultured cells was enriched through aqueous two-phase partition. PM proteins were extracted and then separated using two-dimensional electrophoresis (2DE). Differentially expressed proteins were identified by mass spectrometry, and verified by western blotting.

RESULTS

Thirteen non-redundant proteins were identified to be differentially expressed in the A7 PM fraction compared to those in the Jurkat PM. Eight had a PM location through Gene Ontology (GO) analysis. A differential protein network of CAPG-ACTR3-CD3D was detected to have interactions with HIV Vpr, Tat, gp160, etc. through STRING software analysis. One of the differential proteins (Macrophage-capping protein (CAPG)) was verified by western blotting to be down- regulated in two cell lines and HIV resting CD4+ T cells negatively selected from patients.

CONCLUSION

We identified 13 proteins in A7 compared to Jurkat cells. CAPG may be a potential biomarker related to HIV latency.

Genetic and epigenetic regulation of abdominal aortic aneurysms

Abdominal aortic aneurysms (AAAs) are focal dilations of the aorta that develop from degenerative changes in the media and adventitia of the vessel. Ruptured AAAs have a mortality of up to 85%, thus it is important to identify patients with AAA at increased risk for rupture who would benefit from increased surveillance and/or surgical repair. Although the exact genetic and epigenetic mechanisms regulating AAA formation are not completely understood, Mendelian cases of AAA, which result from pathologic variants in a single gene, have helped provide a basic understanding of AAA pathophysiology. More recently, genome wide associated studies (GWAS) have identified additional variants, termed single nucleotide polymorphisms, in humans that may be associated with AAAs. While some variants may be associated with AAAs and play causal roles in aneurysm pathogenesis, it should be emphasized that the majority of SNPs do not actually cause disease. In addition to GWAS, other studies have uncovered epigenetic causes of disease that regulate expression of genes known to be important in AAA pathogenesis. This review describes many of these genetic and epigenetic contributors of AAAs, which altogether provide a deeper insight into AAA pathogenesis.

Genetic Insights Into Smooth Muscle Cell Contributions to Coronary Artery Disease

Coronary artery disease is a complex cardiovascular disease involving an interplay of genetic and environmental influences over a lifetime. Although considerable progress has been made in understanding lifestyle risk factors, genetic factors identified from genome-wide association studies may capture additional hidden risk undetected by traditional clinical tests. These genetic discoveries have highlighted many candidate genes and pathways dysregulated in the vessel wall, including those involving smooth muscle cell phenotypic modulation and injury responses. Here, we summarize experimental evidence for a few genome-wide significant loci supporting their roles in smooth muscle cell biology and disease. We also discuss molecular quantitative trait locus mapping as a powerful discovery and fine-mapping approach applied to smooth muscle cell and coronary artery disease-relevant tissues. We emphasize the critical need for alternative genetic strategies, including cis/trans-regulatory network analysis, genome editing, and perturbations, as well as single-cell sequencing in smooth muscle cell tissues and model organisms, under both normal and disease states. By integrating multiple experimental and analytical modalities, these multidimensional datasets should improve the interpretation of coronary artery disease genome-wide association studies and molecular quantitative trait locus signals and inform candidate targets for therapeutic intervention or risk prediction.

Genome-Wide Identification of Regulatory DNA Elements in Crop Plants

cis-regulatory DNA elements (CREs) are noncoding but functional DNA sequences. The binding of regulatory proteins into CRE regions leads to chromatin high sensitive to DNase I digestion, which are termed as DNase I hypersensitive sites (DHSs). These DHSs can be efficiently detected through DNase I digestion followed by high-throughput DNA sequencing (DNase-seq). Thus, DNase-seq has become a powerful technique for DHSs mapping at whole-genome level in both plants and animals. Here we describe a DNase-seq procedure modified and developed for crop plants. These plants usually contain large amounts of repetitive sequences and complex organic constituents. With the main improvement in nuclei isolation, this method has been successfully used in mapping DHSs in cotton and sugarcane.

Unveiling the gene regulatory landscape in diseases through the identification of DNase I-hypersensitive sites

DNase I-hypersensitive sites (DHSs) serve key roles in the regulation of gene transcription as markers of cis-regulatory elements (CREs). Recent advances in next-generation sequencing have enabled the genome-wide location and annotation of DHSs in a variety of cells. Numerous studies have confirmed that DHSs are involved in several processes in cell fate decision and development. DHSs have also been indicated in cancer and inherited diseases as driver distal regulatory elements. Here, the definition of DHSs is reviewed, in addition to high-throughput methods of DHS identification. Furthermore, the function of DHSs in gene expression is probed. The roles of DHSs in disease occurrence are also reviewed and discussed. Concomitant advances in the identification of essential roles of DHSs will assist in disclosing the underlying molecular mechanisms, supplementing gene transcription and enlarging the molecular basis of DHS-related bioprocesses, phenotypes, distinct traits and diseases.

Epigenetic reprogramming of immune cells in injury, repair, and resolution

Immune cells are pivotal in the reaction to injury, whereupon, under ideal conditions, repair and resolution phases restore homeostasis following initial acute inflammation. Immune cell activation and reprogramming require transcriptional changes that can only be initiated if epigenetic alterations occur. Recently, accelerated deciphering of epigenetic mechanisms has extended knowledge of epigenetic regulation, including long-distance chromatin remodeling, DNA methylation, posttranslational histone modifications, and involvement of small and long noncoding RNAs. Epigenetic changes have been linked to aspects of immune cell development, activation, and differentiation. Furthermore, genome-wide epigenetic landscapes have been established for some immune cells, including tissue-resident macrophages, and blood-derived cells including T cells. The epigenetic mechanisms underlying developmental steps from hematopoietic stem cells to fully differentiated immune cells led to development of epigenetic technologies and insights into general rules of epigenetic regulation. Compared with more advanced research areas, epigenetic reprogramming of immune cells in injury remains in its infancy. While the early epigenetic mechanisms supporting activation of the immune response to injury have been studied, less is known about resolution and repair phases and cell type-specific changes. We review prominent recent findings concerning injury-mediated epigenetic reprogramming, particularly in stroke and myocardial infarction. Lastly, we illustrate how single-cell technologies will be crucial to understanding epigenetic reprogramming in the complex sequential processes following injury.

Navigating the non-coding genome in heart development and Congenital Heart Disease

Congenital Heart Disease (CHD) is characterised by a wide range of cardiac defects, from mild to life-threatening, which occur in babies worldwide. To date, there is no cure to CHD, however, progress in surgery has reduced its mortality allowing children affected by CHD to reach adulthood. In an effort to understand its genetic basis, several studies involving whole-genome sequencing (WGS) of patients with CHD have been undertaken and generated a great wealth of information. The majority of putative causative mutations identified in WGS studies fall into the non-coding part of the genome. Unfortunately, due to the lack of understanding of the function of these non-coding mutations, it is challenging to establish a causal link between the non-coding mutation and the disease. Thus, here we review the state-of-the-art approaches to interpret non-coding mutations in the context of CHD and address the following questions: What are the non-coding sequences important for cardiac function? Which technologies are used to identify them? Which resources are available to analyse them? What mutations are expected in these non-coding sequences? Learning from developmental process, what is their expected role in CHD?

Gene Therapy Leaves a Vicious Cycle

The human genetic code encrypted in thousands of genes holds the secret for synthesis of proteins that drive all biological processes necessary for normal life and death. Though the genetic ciphering remains unchanged through generations, some genes get disrupted, deleted and or mutated, manifesting diseases, and or disorders. Current treatment options—chemotherapy, protein therapy, radiotherapy, and surgery available for no more than 500 diseases—neither cure nor prevent genetic errors but often cause many side effects. However, gene therapy, colloquially called “living drug,” provides a one-time treatment option by rewriting or fixing errors in the natural genetic ciphering. Since gene therapy is predominantly a viral vector-based medicine, it has met with a fair bit of skepticism from both the science fraternity and patients. Now, thanks to advancements in gene editing and recombinant viral vector development, the interest of clinicians and pharmaceutical industries has been rekindled. With the advent of more than 12 different gene therapy drugs for curing cancer, blindness, immune, and neuronal disorders, this emerging experimental medicine has yet again come in the limelight. The present review article delves into the popular viral vectors used in gene therapy, advances, challenges, and perspectives.

Profiling of Pluripotency Factors in Single Cells and Early Embryos

Cell fate decisions are governed by sequence-specific transcription factors (TFs) that act in small populations of cells within developing embryos. To understand their functions in vivo, it is important to identify TF binding sites in these cells. However, current methods cannot profile TFs genome-wide at or near the single-cell level. Here we adapt the cleavage under targets and release using nuclease (CUT&RUN) method to profile TFs in low cell numbers, including single cells and individual pre-implantation embryos. Single-cell experiments suggest that only a fraction of TF binding sites are occupied in most cells, in a manner broadly consistent with measurements of peak intensity from multi-cell studies. We further show that chromatin binding by the pluripotency TF NANOG is highly dependent on the SWI/SNF chromatin remodeling complex in individual blastocysts but not in cultured cells. Ultra-low input CUT&RUN (uliCUT&RUN) therefore enables interrogation of TF binding from rare cell populations of particular importance in development or disease.

High-resolution interrogation of functional elements in the noncoding genome

The noncoding genome affects gene regulation and disease, yet we lack tools for rapid identification and manipulation of noncoding elements. We developed a CRISPR screen using ~18,000 single guide RNAs targeting >700 kilobases surrounding the genes NF1, NF2, and CUL3, which are involved in BRAF inhibitor resistance in melanoma. We find that noncoding locations that modulate drug resistance also harbor predictive hallmarks of noncoding function. With a subset of regions at the CUL3 locus, we demonstrate that engineered mutations alter transcription factor occupancy and long-range and local epigenetic environments, implicating these sites in gene regulation and chemotherapeutic resistance. Through our expansion of the potential of pooled CRISPR screens, we provide tools for genomic discovery and for elucidating biologically relevant mechanisms of gene regulation.

Human centromere repositioning within euchromatin after partial chromosome deletion

Centromeres are defined by a specialized chromatin organization that includes nucleosomes that contain the centromeric histone variant centromere protein A (CENP-A) instead of canonical histone H3. Studies in various organisms have shown that centromeric chromatin (i.e., CENP-A chromatin or centrochromatin) exhibits plasticity, in that it can assemble on different types of DNA sequences. However, once established on a chromosome, the centromere is maintained at the same position. In humans, this location is the highly homogeneous repetitive DNA alpha satellite. Mislocalization of centromeric chromatin to atypical locations can lead to genome instability, indicating that restriction of centromeres to a distinct genomic position is important for cell and organism viability. Here, we describe a rearrangement of Homo sapiens chromosome 17 (HSA17) that has placed alpha satellite DNA next to euchromatin. We show that on this mutant chromosome, CENP-A chromatin has spread from the alpha satellite into the short arm of HSA17, establishing a ∼700 kb hybrid centromeric domain that spans both repetitive and unique sequences and changes the expression of at least one gene over which it spreads. Our results illustrate the plasticity of human centromeric chromatin and suggest that heterochromatin normally constrains CENP-A chromatin onto alpha satellite DNA. This work highlights that chromosome rearrangements, particularly those that remove the pericentromere, create opportunities for centromeric nucleosomes to move into non-traditional genomic locations, potentially changing the surrounding chromatin environment and altering gene expression.

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Nuclease Footprints in Sperm Project Past and Future Chromatin Regulatory Events

Nuclear remodeling to a condensed state is a hallmark of spermatogenesis. This is achieved by replacement of histones with protamines. Regions retaining nucleosomes may be of functional significance. To determine their potential roles, sperm from wild type and transgenic mice harboring a single copy insert of the human protamine cluster were subjected to Micrococcal Nuclease-seq. CENTIPEDE, a hierarchical Bayesian model, was used to identify multiple spatial patterns, "footprints", of MNase-seq reads along the sperm genome. Regions predicted by CENTIPEDE analysis to be bound by a regulatory factor in sperm were correlated with genomic landmarks and higher order chromatin structure datasets to identify potential roles for these factors in regulating either prior or post spermatogenic, i.e., early embryonic events. This approach linked robust endogenous protamine transcription and transgene suppression to its chromatin environment within topologically associated domains. Of the candidate enhancer-bound regulatory proteins, Ctcf, was associated with chromatin domain boundaries in testes and embryonic stem cells. The continuity of Ctcf binding through the murine germline may permit rapid reconstitution of chromatin organization following fertilization. This likely reflects its preparation for early zygotic genome activation and comparatively accelerated preimplantation embryonic development program observed in mouse as compared to human and bull.

Integrative Tissue-Specific Functional Annotations in the Human Genome Provide Novel Insights on Many Complex Traits and Improve Signal Prioritization in Genome Wide Association Studies

Extensive efforts have been made to understand genomic function through both experimental and computational approaches, yet proper annotation still remains challenging, especially in non-coding regions. In this manuscript, we introduce GenoSkyline, an unsupervised learning framework to predict tissue-specific functional regions through integrating high-throughput epigenetic annotations. GenoSkyline successfully identified a variety of non-coding regulatory machinery including enhancers, regulatory miRNA, and hypomethylated transposable elements in extensive case studies. Integrative analysis of GenoSkyline annotations and results from genome-wide association studies (GWAS) led to novel biological insights on the etiologies of a number of human complex traits. We also explored using tissue-specific functional annotations to prioritize GWAS signals and predict relevant tissue types for each risk locus. Brain and blood-specific annotations led to better prioritization performance for schizophrenia than standard GWAS p-values and non-tissue-specific annotations. As for coronary artery disease, heart-specific functional regions was highly enriched of GWAS signals, but previously identified risk loci were found to be most functional in other tissues, suggesting a substantial proportion of still undetected heart-related loci. In summary, GenoSkyline annotations can guide genetic studies at multiple resolutions and provide valuable insights in understanding complex diseases. GenoSkyline is available at http://genocanyon.med.yale.edu/GenoSkyline.

Remodeling of retrotransposon elements during epigenetic induction of adult visual cortical plasticity by HDAC inhibitors

BACKGROUND

The capacity for plasticity in the adult brain is limited by the anatomical traces laid down during early postnatal life. Removing certain molecular brakes, such as histone deacetylases (HDACs), has proven to be effective in recapitulating juvenile plasticity in the mature visual cortex (V1). We investigated the chromatin structure and transcriptional control by genome-wide sequencing of DNase I hypersensitive sites (DHSS) and cap analysis of gene expression (CAGE) libraries after HDAC inhibition by valproic acid (VPA) in adult V1.

RESULTS

We found that VPA reliably reactivates the critical period plasticity and induces a dramatic change of chromatin organization in V1 yielding significantly greater accessibility distant from promoters, including at enhancer regions. VPA also induces nucleosome eviction specifically from retrotransposon (in particular SINE) elements. The transiently accessible SINE elements overlap with transcription factor-binding sites of the Fox family. Mapping of transcription start site activity using CAGE revealed transcription of epigenetic and neural plasticity-regulating genes following VPA treatment, which may help to re-program the genomic landscape and reactivate plasticity in the adult cortex.

CONCLUSIONS

Treatment with HDAC inhibitors increases accessibility to enhancers and repetitive elements underlying brain-specific gene expression and reactivation of visual cortical plasticity.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Genome-wide mapping of DNase I hypersensitive sites in plants

Genomic regions associated with regulatory proteins are known to be highly sensitive to DNase I digestion and are termed DNase I hypersensitive sites (DHSs). DHSs can be identified by DNase I digestion followed by high-throughput DNA sequencing (DNase-seq). DNase-seq has become a powerful technique for genome-wide mapping of chromatin accessibility in eukaryotes with a sequenced genome. We have developed a DNase-seq procedure in plants. This procedure was adapted from the protocol originally developed for mammalian cell lines. It includes plant nuclei isolation, digestion of purified nuclei with DNase I, recovery of DNase-trimmed DNA fragments, DNase-seq library development, Illumina sequencing and data analysis. We also introduce a barcoding system for library preparation. We have conducted DNase-seq in both Arabidopsis thaliana and rice, and developed genome-wide open chromatin maps in both species. These DHS datasets have been used to detect footprints from regulatory protein binding and to reveal genome-wide nucleosome positioning patterns.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Chromatin accessibility mapping identifies mediators of basal transcription and retinoid-induced repression of OTX2 in medulloblastoma

Despite an emerging understanding of the genetic alterations giving rise to various tumors, the mechanisms whereby most oncogenes are overexpressed remain unclear. Here we have utilized an integrated approach of genomewide regulatory element mapping via DNase-seq followed by conventional reporter assays and transcription factor binding site discovery to characterize the transcriptional regulation of the medulloblastoma oncogene Orthodenticle Homeobox 2 (OTX2). Through these studies we have revealed that OTX2 is differentially regulated in medulloblastoma at the level of chromatin accessibility, which is in part mediated by DNA methylation. In cell lines exhibiting chromatin accessibility of OTX2 regulatory regions, we found that autoregulation maintains OTX2 expression. Comparison of medulloblastoma regulatory elements with those of the developing brain reveals that these tumors engage a developmental regulatory program to drive OTX2 transcription. Finally, we have identified a transcriptional regulatory element mediating retinoid-induced OTX2 repression in these tumors. This work characterizes for the first time the mechanisms of OTX2 overexpression in medulloblastoma. Furthermore, this study establishes proof of principle for applying ENCODE datasets towards the characterization of upstream trans-acting factors mediating expression of individual genes.

A system-level model for the microbial regulatory genome

Microbes can tailor transcriptional responses to diverse environmental challenges despite having streamlined genomes and a limited number of regulators. Here, we present data-driven models that capture the dynamic interplay of the environment and genome-encoded regulatory programs of two types of prokaryotes: Escherichia coli (a bacterium) and Halobacterium salinarum (an archaeon). The models reveal how the genome-wide distributions of cis-acting gene regulatory elements and the conditional influences of transcription factors at each of those elements encode programs for eliciting a wide array of environment-specific responses. We demonstrate how these programs partition transcriptional regulation of genes within regulons and operons to re-organize gene-gene functional associations in each environment. The models capture fitness-relevant co-regulation by different transcriptional control mechanisms acting across the entire genome, to define a generalized, system-level organizing principle for prokaryotic gene regulatory networks that goes well beyond existing paradigms of gene regulation. An online resource (http://egrin2.systemsbiology.net) has been developed to facilitate multiscale exploration of conditional gene regulation in the two prokaryotes.

Genome-wide analysis of the chromatin composition of histone H2A and H3 variants in mouse embryonic stem cells

Genome-wide distribution of the majority of H2A and H3 variants (H2A, H2AX, H2AZ, macroH2A, H3.1, H3.2 and H3.3) was simultaneously investigated in mouse embryonic stem cells by chromatin immunoprecipitation sequencing. Around the transcription start site, histone variant distribution differed between genes possessing promoters of high and low CpG density, regardless of their expression levels. In the intergenic regions, regulatory elements were enriched in H2A.Z and H3.3, whereas repeat elements were abundant in H2A and macroH2A, and H3.1, respectively. Analysis of H2A and H3 variant combinations composing nucleosomes revealed that the H2A.Z and H3.3 combinations were present at a higher frequency throughout the genome than the other combinations, suggesting that H2A.Z and H3.3 associate preferentially with each other to comprise the nucleosomes independently of genome region. Finally, we found that chromatin was unstable only in regions where it was enriched in both H2A.Z and H3.3, but strongly quantified stable in regions in which only H3.3 was abundant. Therefore, histone variant composition is an important determinant of chromatin structure, which is associated with specific genomic functions.

FRA2 is a STAT5 target gene regulated by IL-2 in human CD4 T cells

Signal transducers and activators of transcription 5(STAT5) are cytokine induced signaling proteins, which regulate key immunological processes, such as tolerance induction, maintenance of homeostasis, and CD4 T-effector cell differentiation. In this study, transcriptional targets of STAT5 in CD4 T cells were studied by Chromatin Immunoprecipitation (ChIP). Genomic mapping of the sites cloned and identified in this study revealed the striking observation that the majority of STAT5-binding sites mapped to intergenic (>50 kb upstream) or intronic, rather than promoter proximal regions. Of the 105 STAT5 responsive binding sites identified, 94% contained the canonical (IFN-γ activation site) GAS motifs.

A number of putative target genes identified here are associated with tumor biology. Here, we identified Fos-related antigen 2 (FRA2) as a transcriptional target of IL-2 regulated STAT5. FRA2 is a basic -leucine zipper (bZIP) motif ‘Fos’ family transcription factor that is part of the AP-1 transcription factor complex and is also known to play a critical role in the progression of human tumours and more recently as a determinant of T cell plasticity. The binding site mapped to an internal intron within the FRA2 gene. The epigenetic architecture of FRA2, characterizes a transcriptionally active promoter as indicated by enrichment for histone methylation marks H3K4me1, H3K4me2, H3K4me3, and transcription/elongation associated marks H2BK5me1 and H4K20me1. FRA2 is regulated by IL-2 in activated CD4 T cells. Consistently, STAT5 bound to GAS sequence in the internal intron of FRA2 and reporter gene assays confirmed IL-2 induced STAT5 binding and transcriptional activation. Furthermore, addition of JAK3 inhibitor (R333) or Daclizumab inhibited the induction in TCR stimulated cells. Taken together, our data suggest that FRA2 is a novel STAT5 target gene, regulated by IL-2 in activated CD4 T cells.

Identification of active regulatory regions from DNA methylation data

We have recently shown that transcription factor binding leads to defined reduction in DNA methylation, allowing for the identification of active regulatory regions from high-resolution methylomes. Here, we present MethylSeekR, a computational tool to accurately identify such footprints from bisulfite-sequencing data. Applying our method to a large number of published human methylomes, we demonstrate its broad applicability and generalize our previous findings from a neuronal differentiation system to many cell types and tissues. MethylSeekR is available as an R package at www.bioconductor.org.

Measuring Arabidopsis chromatin accessibility using DNase I-polymerase chain reaction and DNase I-chip assays

DNA accessibility is an important layer of regulation of DNA-dependent processes. Methods that measure DNA accessibility at local and genome-wide scales have facilitated a rapid increase in the knowledge of chromatin architecture in animal and yeast systems. In contrast, much less is known about chromatin organization in plants. We developed a robust DNase I-polymerase chain reaction (PCR) protocol for the model plant Arabidopsis (Arabidopsis thaliana). DNA accessibility is probed by digesting nuclei with a gradient of DNase I followed by locus-specific PCR. The reduction in PCR product formation along the gradient of increasing DNase I concentrations is used to determine the accessibility of the chromatin DNA. We explain a strategy to calculate the decay constant of such signal reduction as a function of increasing DNase I concentration. This allows describing DNA accessibility using a single variable: the decay constant. We also used the protocol together with AGRONOMICS1 DNA tiling microarrays to establish genome-wide DNase I sensitivity landscapes.

Identification of a secondary promoter within the human B cell receptor component gene hCD79b

The human B cell-specific protein, CD79b (also known as Igβ and B29) constitutes an essential signal transduction component of the B cell receptor. Although its function is central to the triggering of B cell terminal differentiation in response to antigen stimulation, the transcriptional determinants that control CD79b gene expression remain poorly defined. In the present study, we explored these determinants using a series of hCD79b transgenic mouse models. Remarkably, we observed that the previously described hCD79b promoter along with its associated enhancer elements and first exon could be deleted without appreciable loss of hCD79b transcriptional activity or tissue specificity. In this deletion setting, a secondary promoter located within exon 2 maintained full levels and specificity of hCD79b transcription. Of note, this secondary promoter was also active, albeit at lower levels, in the wild-type hCD79b locus. The activity of the secondary promoter was dependent on the action(s) of a conserved sequence element mapping to a chromatin DNase I hypersensitive site located within intron 1. mRNA generated from this secondary promoter is predicted to encode an Igβ protein lacking a signal sequence and thus unable to serve normal B cell receptor function. Although the physiologic role of the hCD79b secondary promoter and its encoded protein remain unclear, the current data suggest that it has the capacity to play a role in normal as well as pathologic states in B cell proliferation and function.

Computational methodology for ChIP-seq analysis

Chromatin immunoprecipitation coupled with massive parallel sequencing (ChIP-seq) is a powerful technology to identify the genome-wide locations of DNA binding proteins such as transcription factors or modified histones. As more and more experimental laboratories are adopting ChIP-seq to unravel the transcriptional and epigenetic regulatory mechanisms, computational analyses of ChIP-seq also become increasingly comprehensive and sophisticated. In this article, we review current computational methodology for ChIP-seq analysis, recommend useful algorithms and workflows, and introduce quality control measures at different analytical steps. We also discuss how ChIP-seq could be integrated with other types of genomic assays, such as gene expression profiling and genome-wide association studies, to provide a more comprehensive view of gene regulatory mechanisms in important physiological and pathological processes.

Methylome of DNase I sensitive chromatin in Populus trichocarpa shoot apical meristematic cells: a simplified approach revealing characteristics of gene-body DNA methylation in open chromatin state

DNA methylation is involved in the control of plant development and adaptation to the environment through modifications of chromatin compaction and gene expression. In poplar (Populus trichocarpa), a perennial plant, variations in DNA methylation have been reported between genotypes and tissues or in response to drought. Nevertheless, the relationships between gene-body DNA methylation, gene expression and chromatin compaction still need clarification. Here, DNA methylation was mapped in the noncondensed chromatin fraction from P. trichocarpa shoot apical meristematic cells, the center of plant morphogenesis, where DNA methylation variations could influence the developmental trajectory. DNase I was used to isolate the noncondensed chromatin fraction. Methylated sequences were immunoprecipitated, sequenced using Illumina/Solexa technology and mapped on the v2.0 poplar genome. Bisulfite sequencing of candidate sequences was used to confirm mapping data and to assess cytosine contexts and methylation levels. While the methylated DNase I hypersensitive site fraction covered 1.9% of the poplar genome, it contained sequences corresponding to 74% of poplar gene models, mostly exons. The level and cytosine context of gene-body DNA methylation varied with the structural characteristics of the genes. Taken together, our data show that DNA methylation is widespread and variable among genes in open chromatin of meristematic cells, in agreement with a role in their developmental trajectory.

Correlation between DNase I hypersensitive site distribution and gene expression in HeLa S3 cells

Mapping DNase I hypersensitive sites (DHSs) within nuclear chromatin is a traditional and powerful method of identifying genetic regulatory elements. DHSs have been mapped by capturing the ends of long DNase I-cut fragments (>100,000 bp), or 100-1200 bp DNase I-double cleavage fragments (also called double-hit fragments). But next generation sequencing requires a DNA library containing DNA fragments of 100-500 bp. Therefore, we used short DNA fragments released by DNase I digestion to generate DNA libraries for next generation sequencing. The short segments are 100-300 bp and can be directly cloned and used for high-throughput sequencing. We identified 83,897 DHSs in 2,343,479 tags across the human genome. Our results indicate that the DHSs identified by this DHS assay are consistent with those identified by longer fragments in previous studies. We also found: (1) the distribution of DHSs in promoter and other gene regions of similarly expressed genes differs among different chromosomes; (2) silenced genes had a more open chromatin structure than previously thought; (3) DHSs in 3'untranslated regions (3'UTRs) are negatively correlated with level of gene expression.

Transcriptional and epigenetic control of T helper cell specification: molecular mechanisms underlying commitment and plasticity

T helper cell differentiation occurs in the context of the extracellular cytokine milieu evoked by diverse microbes and other pathogenic stimuli along with T cell receptor stimulation. The culmination of these signals results in specification of T helper lineages, which occurs through the combinatorial action of multiple transcription factors that establish distinctive transcriptomes. In this manner, inducible, but constitutively active, master regulators work in conjunction with factors such as the signal transducer and activator of transcriptions (STATs) that sense the extracellular environment. The acquisition of a distinctive transcriptome also depends on chromatin modifications that impact key cis elements as well as the changes in global genomic organization. Thus, signal transduction and epigenetics are linked in these processes of differentiation. In this review, recent advances in understanding T helper lineage specification and deciphering the action of transcription factors are summarized with emphasis on comprehensive views of the dynamic T cell epigenome.

Salt fractionation of nucleosomes for genome-wide profiling

Salt fractionation of nucleosomes, a classical method for defining "active" chromatin based on nucleosome solubility, has recently been adapted for genome-scale profiling. This method has several advantages for profiling chromatin dynamics, including general applicability to cell lines and tissues, quantitative recovery of chromatin, base-pair resolution of nucleosomes, and overall simplicity both in concept and execution. This chapter provides detailed protocols for nuclear isolation, chromatin fragmentation by micrococcal nuclease digestion, successive solubilization of chromatin fractions by addition of increasing concentrations of salt, and genome-wide analyses through microarray hybridization and next-generation sequencing.