DNase I Hypersensitive Sites Research Articles

The esBAF and ISWI nucleosome remodeling complexes influence occupancy of overlapping dinucleosomes and fragile nucleosomes in murine embryonic stem cells

BackgroundNucleosome remodeling factors regulate the occupancy and positioning of nucleosomes genome-wide through ATP-driven DNA translocation. While many nucleosomes are consistently well-positioned, some nucleosomes and alternative nucleosome structures are more sensitive to nuclease digestion or are transitory. Fragile nucleosomes are nucleosome structures that are sensitive to nuclease digestion and may be composed of either six or eight histone proteins, making these either hexasomes or octasomes. Overlapping dinucleosomes are composed of two merged nucleosomes, lacking one H2A:H2B dimer, creating a 14-mer wrapped by ~ 250 bp of DNA. In vitro studies of nucleosome remodeling suggest that the collision of adjacent nucleosomes by sliding stimulates formation of overlapping dinucleosomes.ResultsTo better understand how nucleosome remodeling factors regulate alternative nucleosome structures, we depleted murine embryonic stem cells of the transcripts encoding remodeler ATPases BRG1 or SNF2H, then performed MNase-seq. We used high- and low-MNase digestion to assess the effects of nucleosome remodeling factors on nuclease-sensitive or “fragile” nucleosome occupancy. In parallel we gel-extracted MNase-digested fragments to enrich for overlapping dinucleosomes. We recapitulate prior identification of fragile nucleosomes and overlapping dinucleosomes near transcription start sites, and identify enrichment of these features around gene-distal DNaseI hypersensitive sites, CTCF binding sites, and pluripotency factor binding sites. We find that BRG1 stimulates occupancy of fragile nucleosomes but restricts occupancy of overlapping dinucleosomes.ConclusionsOverlapping dinucleosomes and fragile nucleosomes are prevalent within the ES cell genome, occurring at hotspots of gene regulation beyond their characterized existence at promoters. Although neither structure is fully dependent on either nucleosome remodeling factor, both fragile nucleosomes and overlapping dinucleosomes are affected by knockdown of BRG1, suggesting a role for the complex in creating or removing these structures.

NAGS, CPS1, and SLC25A13 (Citrin) at the Crossroads of Arginine and Pyrimidines Metabolism in Tumor Cells.

Urea cycle enzymes and transporters collectively convert ammonia into urea in the liver. Aberrant overexpression of carbamylphosphate synthetase 1 (CPS1) and SLC25A13 (citrin) genes has been associated with faster proliferation of tumor cells due to metabolic reprogramming that increases the activity of the CAD complex and pyrimidine biosynthesis. N-acetylglutamate (NAG), produced by NAG synthase (NAGS), is an essential activator of CPS1. Although NAGS is expressed in lung cancer derived cell lines, expression of the NAGS gene and its product was not evaluated in tumors with aberrant expression of CPS1 and citrin. We used data mining approaches to identify tumor types that exhibit aberrant overexpression of NAGS, CPS1, and citrin genes, and evaluated factors that may contribute to increased expression of the three genes and their products in tumors. Median expression of NAGS, CPS1, and citrin mRNA was higher in glioblastoma multiforme (GBM), glioma, and stomach adenocarcinoma (STAD) samples compared to the matched normal tissue. Median expression of CPS1 and citrin mRNA was higher in the lung adenocarcinoma (LUAD) sample while expression of NAGS mRNA did not differ. High NAGS expression was associated with an unfavorable outcome in patients with glioblastoma and GBM. Low NAGS expression was associated with an unfavorable outcome in patients with LUAD. Patterns of DNase hypersensitive sites and histone modifications in the upstream regulatory regions of NAGS, CPS1, and citrin genes were similar in liver tissue, lung tissue, and A549 lung adenocarcinoma cells despite different expression levels of the three genes in the liver and lung. Citrin gene copy numbers correlated with its mRNA expression in glioblastoma, GBM, LUAD, and STAD samples. There was little overlap between NAGS, CPS1, and citrin sequence variants found in patients with respective deficiencies, tumor samples, and individuals without known rare genetic diseases. The correlation between NAGS, CPS1, and citrin mRNA expression in the individual glioblastoma, GBM, LUAD, and STAD samples was very weak. These results suggest that the increased cytoplasmic supply of either carbamylphosphate, produced by CPS1, or aspartate may be sufficient to promote tumorigenesis, as well as the need for an alternative explanation of CPS1 activity in the absence of NAGS expression and NAG.

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Quantification and Mapping of Alkylation in the Human Genome Reveal Single Nucleotide Resolution Precursors of Mutational Signatures.

Chemical modifications to DNA bases, including DNA adducts arising from reactions with electrophilic chemicals, are well-known to impact cell growth, miscode during replication, and influence disease etiology. However, knowledge of how genomic sequences and structures influence the accumulation of alkylated DNA bases is not broadly characterized with high resolution, nor have these patterns been linked with overall quantities of modified bases in the genome. For benzo(a) pyrene (BaP), a ubiquitous environmental carcinogen, we developed a single-nucleotide resolution damage sequencing method to map in a human lung cell line the main mutagenic adduct arising from BaP. Furthermore, we combined this analysis with quantitative mass spectrometry to evaluate the dose-response profile of adduct formation. By comparing damage abundance with DNase hypersensitive sites, transcription levels, and other genome annotation data, we found that although overall adduct levels rose with increasing chemical exposure concentration, genomic distribution patterns consistently correlated with chromatin state and transcriptional status. Moreover, due to the single nucleotide resolution characteristics of this DNA damage map, we could determine preferred DNA triad sequence contexts for alkylation accumulation, revealing a characteristic DNA damage signature. This new BaP damage signature had a profile highly similar to mutational signatures identified previously in lung cancer genomes from smokers. Thus, these data provide insight on how genomic features shape the accumulation of alkylation products in the genome and predictive strategies for linking single-nucleotide resolution in vitro damage maps with human cancer mutations.

Analysis of endometrial carcinoma TCGA reveals differences in DNA methylation in tumors from Black and White women

ObjectiveRacial disparities exist in cancer patients both in incidence and death rates. In endometrial cancer, Black patients are reported to have higher incidence of aggressive endometrial cancer subtypes and higher death rates than White women. To date, diagnostic and prognostic biomarkers associated with race-specific methylation driven genes have yet to be identified.The objective of this study was to explore DNA methylation patterns in endometrial tumor samples from White and Black women. MethodsDifferentially methylated CpGs (DMCs) and differentially methylated regions (DMRs) were identified in White tumor samples compared to Black tumor samples using Endometrial Carcinoma (EC) methylation and clinical data from The Cancer Genome Atlas (TCGA). Survival analysis was performed using survival R package and results were visualized using Kaplan-Meier plots. To access the correlation between changes in methylation and gene expression, we downloaded raw RNA-sequencing by Expectation-Maximization (RSEM) counts data from The Cancer Genome Atlas (TCGA) using TCGABiolinks package (v2.18.0). ResultsOur analysis revealed 704 differentially methylated CpGs in tumors from Black and White women. These differentially methylated genes showed strong negative correlation with gene expression and statistically significant enrichment in regulatory regions such as DNase I hypersensitivity sites (DHSs) and transcription factor binding sites (TFBSs). Increased variability in methylation occurred in genes such as the insulin signaling pathway in Black tumor samples. ConclusionBy using two independent statistical method based on means (DMR, DMCs) and variances (DVCs) on the endometrial carcinoma TCGA data, we showed that endometrial tumors from Black women are hypomethylated and more hypervariable than tumors from White women. In-depth investigation of these methylation driven markers can aid in successful management of endometrial cancer disparities and improved overall survival in Black and White populations.

MYCN Regulates Cell Survival Via EIF4G1 in Inv(16) Acute Myeloid Leukemia

Transcription factors are frequently deregulated in leukemia and are pivotal to its pathogenesis. The core-binding factor (CBF) transcription factor is a heterodimer of CBFβ and RUNX proteins, is critical to hematopoiesis, and dysregulated in AML. The fusion oncogene CBFβ-SMMHC expressed in AML with the chromosome inversion inv(16)(p13q22) acts as a driver oncogene in hematopoietic stem cells by sequestering RUNX1 by protein-protein interaction to induce AML. The CBFβ-SMMHC inhibitor, AI-10-49, specifically disrupts its binding to RUNX1 resulting in genome-wide RUNX1 binding, apoptosis induction, and enhanced survival in inv(16) mouse model [Illendula* and Pulikkan* et al, Science 2015]. While previous studies have demonstrated that AI-10-49 induced cell death is partly mediated by repression of MYC expression in inv(16) AML cells [Pulikkan et al, Cell 2018], the function of other MYC family members, in inv(16) cell survival, and overall in AML, remains poorly understood. By conducting Real Time RT-PCR in CD34+ bone marrow cells from human AML samples, we identified MYCN upregulation in inv(16) AML. Analysis of RNA-seq data from inv(16) AML cells treated with AI-10-49 revealed that MYCN expression is downregulated (7-fold). MYCN silencing via CRISPR/Cas9 induced apoptosis in inv(16) AML cells. Furthermore, MYCN genetic targeting by CRISPR/Cas9-RNP approach in human primary inv(16) AML cells transplanted in humanized NSGS xenograft mouse model displayed significant delay in leukemia development, underlining the role of MYCN in inv(16) AML maintenance. ChIP-seq analysis conducted in inv(16) AML cells treated with AI-10-49 demonstrated increased RUNX1 occupancy at MYCN distal genomic element located at +26 Kb (named 'MYCN-e26') downstream from MYCN TSS without any RUNX1 binding at MYCN TSS. ChIP-seq analysis for histone marks such as H3K27ac and H3K4Me1, and chromatin accessibility by ATAC-seq suggested MYCN-e26 acts as a transcriptionally active enhancer element. Analysis of DNase I hypersensitive sites (DHSs) representing accessible chromatin regions in primary AML samples revealed that MYCN-e26 is a bona fide AML enhancer. Finally, deletion of a genomic element in MYCN-e26 surrounding the RUNX1 binding site (200 bp) by CRISPR/Cas9 genome editing downregulated MYCN transcript levels and inhibited inv(16) AML cell survival. Taken together, MYCN-e26 is instrumental in MYCN transcription and cell survival in inv(16) AML cells as a RUNX1 target. To identify the mechanism(s) by which MYCN regulates inv(16) AML cell survival, we performed CUT&Tag sequencing for MYCN and H3K27ac histone mark in inv(16) AML cells. Gene ontology analysis for MYCN binding peaks in transcriptionally active regions revealed mRNA metabolism and stability as the major pathways orchestrated by MYCN in inv(16) AML cells. By comparing the mRNA metabolism/ stability signature genes with transcripts significantly deregulated in human primary inv(16) AML samples, we identified eukaryotic translation initiation factor 4 gamma 1 (EIF4G1) as one of the major genes regulated by MYCN in inv(16) AML cell survival. Our CUT&Tag sequencing analysis demonstrated MYCN binds to EIF4G1 promoter and enhances its transcription. MYCN silencing by CRISPR/Cas9 approach in inv(16) AML cells showed significant downregulation of EIF4G1, suggesting MYCN is required for EIF4G1 expression in inv(16) AML cells. EIF4G1 genetic silencing via CRISPR/Cas9 approach induced apoptosis in inv(16) AML cells in vitro, and significantly delayed leukemia development in NSGS xenograft mouse model, suggesting EIF4G1 is required for the survival of inv(16) leukemia-initiating cells. Our studies demonstrate that CBFβ-SMMHC modulates MYCN levels by inhibition of RUNX1 binding to MYCN distal enhancer in inv(16) AML cells. In addition, our studies bring novel function of MYCN in inv(16) leukemic cell survival by regulating EIF4G1. Future therapeutic approaches targeting MYCN/EIF4G1 axis in combination with AI-10-49 warrant better clinical outcome for the treatment of inv(16) AML.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a "genome shock", leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton (Gossypium hirsutum and Gossypium barbadense, AADD) and its extant AA (Gossypium arboreum) and DD (Gossypium raimondii) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii. In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense, but absent from a resynthesized hybrid of G. arboreum and G. raimondii, suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis-regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

IDHS-FFLG: Identifying DNase I Hypersensitive Sites by Feature Fusion and Local-Global Feature Extraction Network.

The DNase I hypersensitive sites (DHSs) are active regions on chromatin that have been found to be highly sensitive to DNase I. These regions contain various cis-regulatory elements, including promoters, enhancers and silencers. Accurate identification of DHSs helps researchers better understand the transcriptional machinery of DNA and deepen the knowledge of functional DNA elements in non-coding sequences. Researchers have developed many methods based on traditional experiments and machine learning to identify DHSs. However, low prediction accuracy and robustness limit their application in genetics research. In this paper, a novel computational approach based on deep learning is proposed by feature fusion and local-global feature extraction network to identify DHSs in mouse, named iDHS-FFLG. First of all, multiple binary features of nucleotides are fused to better express sequence information. Then, a network consisting of the convolutional neural network (CNN), bidirectional long short-term memory (BiLSTM) and self-attention mechanism is designed to extract local features and global contextual associations. In the end, the prediction module is applied to distinguish between DHSs and non-DHSs. The results of several experiments demonstrate the superior performances of iDHS-FFLG compared to the latest methods.

COTTONOMICS: a comprehensive cotton multi-omics database.

The rapid advancement of sequencing technology, including next-generation sequencing (NGS), has greatly improved sequencing efficiency and decreased cost. Consequently, huge amounts of genomic, transcriptomic and epigenetic data concerning cotton species have been generated and released. These large-scale data provide immense opportunities for the study of cotton genomic structure and evolution, population genetic diversity and genome-wide mining of excellent genes for important traits. However, the complexity of NGS data also causes distress, as it cannot be utilized easily. Here, we presented the cotton omics data platform COTTONOMICS (http://cotton.zju.edu.cn/), an easily accessible web database that integrates 32.5 TB of omics data including seven assembled genomes, resequencing data from 1180 allotetraploid cotton accessions and RNA-sequencing (RNA-seq), small RNA-sequencing (smRNA-seq), Chromatin Immunoprecipitation sequencing (ChIP-seq), DNase hypersensitive sites sequencing (DNase-seq) and Bisulfite sequencing (BS-seq). COTTONOMICS allows users to employ various search scenarios and retrieve information concerning the cotton genomes, genomic variation (Single nucleotide polymorphisms (SNPs) and Insertion and Deletion (InDels)), gene expression, smRNA expression, epigenetic regulation and quantitative trait locus (QTLs). The user-friendly web interface offers a variety of modules for storing, retrieving, analyzing and visualizing cotton multi-omics data to diverse ends, thereby enabling users to decipher cotton population genetics and identify potential novel genes that influence agronomically beneficial traits. Database URL: http://cotton.zju.edu.cn

Fragmentomics of urinary cell-free DNA in nuclease knockout mouse models.

Urinary cell-free DNA (ucfDNA) is a potential biomarker for bladder cancer detection. However, the biological characteristics of ucfDNA are not well understood. We explored the roles of deoxyribonuclease 1 (DNASE1) and deoxyribonuclease 1-like 3 (DNASE1L3) in the fragmentation of ucfDNA using mouse models. The deletion of Dnase1 in mice (Dnase1-/-) caused aberrations in ucfDNA fragmentation, including a 24-fold increase in DNA concentration, and a 3-fold enrichment of long DNA molecules, with a relative decrease of fragments with thymine ends and reduction of jaggedness (i.e., the presence of single-stranded protruding ends). In contrast, such changes were not observed in mice with Dnase1l3 deletion (Dnase1l3-/-). These results suggested that DNASE1 was an important nuclease contributing to the ucfDNA fragmentation. Western blot analysis revealed that the concentration of DNASE1 protein was higher in urine than DNASE1L3. The native-polyacrylamide gel electrophoresis zymogram showed that DNASE1 activity in urine was higher than that in plasma. Furthermore, the proportion of ucfDNA fragment ends within DNase I hypersensitive sites (DHSs) was significantly increased in Dnase1-deficient mice. In humans, patients with bladder cancer had lower proportions of ucfDNA fragment ends within the DHSs when compared with participants without bladder cancer. The area under the curve (AUC) for differentiating patients with and without bladder cancer was 0.83, suggesting the analysis of ucfDNA fragmentation in the DHSs may have potential for bladder cancer detection. This work revealed the intrinsic links between the nucleases in urine and ucfDNA fragmentomics.

Comparison of chromatin accessibility landscapes during early development of prefrontal cortex between rhesus macaque and human

Epigenetic information regulates gene expression and development. However, our understanding of the evolution of epigenetic regulation on brain development in primates is limited. Here, we compared chromatin accessibility landscapes and transcriptomes during fetal prefrontal cortex (PFC) development between rhesus macaques and humans. A total of 304,761 divergent DNase I-hypersensitive sites (DHSs) are identified between rhesus macaques and humans, although many of these sites share conserved DNA sequences. Interestingly, most of the cis-elements linked to orthologous genes with dynamic expression are divergent DHSs. Orthologous genes expressed at earlier stages tend to have conserved cis-elements, whereas orthologous genes specifically expressed at later stages seldom have conserved cis-elements. These genes are enriched in synapse organization, learning and memory. Notably, DHSs in the PFC at early stages are linked to human educational attainment and cognitive performance. Collectively, the comparison of the chromatin epigenetic landscape between rhesus macaques and humans suggests a potential role for regulatory elements in the evolution of differences in cognitive ability between non-human primates and humans.

Chromatin accessibility shapes meiotic recombination in mouse primordial germ cells through assisting double-strand breaks and loop formation

Meiotic recombination is a driver of evolution, and aberrant recombination is a major contributor to aneuploidy in mammals. Mechanism of recombination remains elusive yet. Here, we present a computational analysis to explore recombination-related dynamics of chromatin accessibility in mouse primordial germ cells (PGCs). Our data reveals that: (1) recombination hotspots which get accessible at meiosis-specific DNase I-hypersensitive sites (DHSs) only when PGCs enter meiosis are located preferentially in intronic and distal intergenic regions; (2) stable DHSs maintained stably across PGC differentiation are enriched by CTCF motifs and CTCF binding and mediate chromatin loop formation; (3) compared with the specific DHSs aroused at meiotic stage, stable DHSs are largely encoded in DNA sequence and also enriched by epigenetic marks; (4) PRDM9 is likely to target nucleosome-occupied hotspot regions and remodels local chromatin structure to make them accessible for recombination machinery; and (5) cells undergoing meiotic recombination are deficient in TAD structure and chromatin loop arrays are organized regularly along the axis formed between homologous chromosomes. Taken together, by analyzing DHS-related DNA features, epigenetic marks and 3D genome structure, we revealed some specific roles of chromatin accessibility in recombination, which would expand our understanding of recombination mechanism.

Prediction of histone post-translational modification patterns based on nascent transcription data.

The role of histone modifications in transcription remains incompletely understood. Here we examine the relationship between histone modifications and transcription using experimental perturbations combined with sensitive machine-learning tools. Transcription predicted the variation in active histone marks and complex chromatin states, like bivalent promoters, down to single-nucleosome resolution and at an accuracy that rivaled the correspondence between independent ChIP-seq experiments. Blocking transcription rapidly removed two punctate marks, H3K4me3 and H3K27ac, from chromatin indicating that transcription is required for active histone modifications. Transcription was also required for maintenance of H3K27me3 consistent with a role for RNA in recruiting PRC2. A subset of DNase-I hypersensitive sites were refractory to prediction, precluding models where transcription initiates pervasively at any open chromatin. Our results, in combination with past literature, support a model in which active histone modifications serve a supportive, rather than an essential regulatory, role in transcription.

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.

Structural switching/polymorphism by sequential base substitution at quasi-palindromic SNP site (G → A) in LCR of human β-globin gene cluster

The human β-globin gene Locus Control Region (LCR), a dominant regulator of globin gene expression contains five tissue-specific DNase I-hypersensitive sites (HSs). A single nucleotide polymorphism (SNP) (A → G) present in HS4 region of locus control region (LCR), have shown a notable association between the G allele and the occurrence of β-thalassemia. This SNP site exhibiting a hairpin − duplex equilibrium manifested in A → B like DNA transition has previously been reported from this laboratory. Since, DNA is a dynamic and adaptable molecule, so any change of a single base within a primary DNA sequence can produce major biological consequences commonly manifested in genetic disorders such as sickle cell anemia and β-thalassemia. Herein, the differential behavior of sequential single base substitutions G → A on the quasi-palindromic sequence (d-TGGGGGCCCCA; HPG11) has been explored. A combination of native gel electrophoresis, circular dichroism (CD), and UV-thermal denaturation (Tm) techniques have been used to investigate the structural polymorphism associated with various variants of HPG11 i.e. HPG11A2 to HPG11A5. The CD spectra confirmed that all the HPG11 variants exhibit a hairpin − duplex equilibrium. Oligomer concentration dependence on CD spectra has been correlated with A → B DNA conformational transition. However, as revealed in gel electrophoresis, HPG11A2 → A5 exhibit the formation of a tetramolecular structure (four-way junction) at higher oligomer concentration. UV-melting studies also supported the melting of hairpin, duplex and four-way junction structure. This polymorphism pattern may possibly be significant for DNA-protein recognition, in the process of regulation of LCR in the β-globin gene.

Functional and HRE motifs count analysis of induction of selected hypoxia-responsive genes by HIF-1 and HIF-2 in human umbilical endothelial cells

We analyzed the effects of selective knockdown of either HIF-1α or HIF-2α on the transcriptional response to hypoxia of human umbilical endothelial cells at two time-points (2 h and 8 h) of hypoxia. We focused on 13 previously identified hypoxia-responsive genes, pre-selected to have different activation kinetics and different proportions of HRE motifs annotated to either HIF-1 or HIF-2 in open promoters – open chromatin DNase-hypersensitive sites (DHS) regions within ±1 kb of the gene start. We report that genes activated by both HIF-1 and 2 tend to be activated earlier than genes activated by HIF-1 only, which, in turn, tend to be activated earlier than genes activated by HIF-2 only. Moreover, for the 13 analyzed genes, we found that the effect of silencing HIF1A on the gene induction by hypoxia is greater for the genes with more HRE motifs annotated to HIF-1 in their promoter open chromatin DHS regions within ±1 kb and also within ±10 kb of the gene start. We corroborated and extended this finding by showing that among 232 genes previously identified as activated by hypoxia, the genes with ChIP-seq peak(s) for HIF-1α within a ±10 kb flank of the gene start contain more HRE motifs annotated to HIF-1 in the DHS regions within this flank than the genes with no ChIP-seq peaks. Also in the whole genome, the DHS regions intersecting ChIP-seq peaks for HIF-1α contain more HRE motifs annotated to HIF-1 than the DHS regions not intersecting the ChIP-seq peaks. This suggests a mechanism, by which higher promoter content of HRE motifs in DHS regions increases HIF-1 binding, which in turn increases gene induction by hypoxia.

Discrete regulatory modules instruct hematopoietic lineage commitment and differentiation

Lineage commitment and differentiation is driven by the concerted action of master transcriptional regulators at their target chromatin sites. Multiple efforts have characterized the key transcription factors (TFs) that determine the various hematopoietic lineages. However, the temporal interactions between individual TFs and their chromatin targets during differentiation and how these interactions dictate lineage commitment remains poorly understood. Here we perform dense, daily, temporal profiling of chromatin accessibility (DNase I-seq) and gene expression changes (total RNA-seq) along ex vivo human erythropoiesis to comprehensively define developmentally regulated DNase I hypersensitive sites (DHSs) and transcripts. We link both distal DHSs to their target gene promoters and individual TFs to their target DHSs, revealing that the regulatory landscape is organized in distinct sequential regulatory modules that regulate lineage restriction and maturation. Finally, direct comparison of transcriptional dynamics (bulk and single-cell) and lineage potential between erythropoiesis and megakaryopoiesis uncovers differential fate commitment dynamics between the two lineages as they exit the stem and progenitor stage. Collectively, these data provide insights into the temporally regulated synergy of the cis- and the trans-regulatory components underlying hematopoietic lineage commitment and differentiation.

IDHS-DT: Identifying DNase I hypersensitive sites by integrating DNA dinucleotide and trinucleotide information

DNase I hypersensitive sites (DHSs) is important for identifying the location of gene regulatory elements, such as promoters, enhancers, silencers, and so on. Thus, it is crucial for discriminating DHSs from non-DHSs. Although some traditional methods, such as Southern blots and DNase-seq technique, have the ability to identify DHSs, these approaches are time-consuming, laborious, and expensive. To address these issues, researchers paid their attention on computational approaches. Therefore, in this study, we developed a novel predictor called iDHS-DT to identify DHSs. In this predictor, the DNA sequences were firstly denoted by physicochemical properties (PC) of DNA dinucleotide and trinucleotide. Then, three different descriptors, including auto-covariance, cross-covariance, and discrete wavelet transform were used to collect related features from the PC matrix. Next, the least absolute shrinkage and selection operator (LASSO) algorithm was employed to remove these irrelevant and redundant features. Finally, these selected features were fed into support vector machine (SVM) for distinguishing DHSs from non-DHSs. The proposed method achieved 97.64% and 98.22% classification accuracy on dataset S1 and S2, respectively. Compared with the existing predictors, our proposed model has significantly improvement in classification performance. Experimental results demonstrated that the proposed method is powerful in identifying DHSs.

An epigenetic and transcriptomic signature of immune tolerance in human monocytes through multi-omics integration

BackgroundThe plasticity of monocytes enables them to exert multiple roles during an immune response, including promoting immune tolerance. How monocytes alter their functions to convey immune tolerance in the context of lower respiratory tract infections in humans is not well understood. Here, we sought to identify epigenetic and transcriptomic features of cytokine production capacity in circulating monocytes during community-acquired pneumonia (CAP).MethodsCirculating CD14+ monocytes were obtained from the blood of CAP patients included in a longitudinal, observational cohort study, on hospitalization (acute stage, n=75), and from the same patients after a 1-month follow-up (recovery stage, n=56). Age and sex-matched non-infectious participants were included as controls (n=41). Ex vivo cytokine production after lipopolysaccharide (LPS) exposure was assessed by multiplex assay. Transcriptomes of circulating monocytes were generated by RNA-sequencing, and DNA methylation levels in the same monocytes were measured by reduced representation bisulfite sequencing. Data were integrated by fitting projection-to-latent-structure models, and signatures derived by partial least squares discrimination.ResultsMonocytes captured during the acute stage exhibited impaired TNF, IL-1β, IL-6, and IL-10 production after ex vivo stimulation with LPS, relative to controls. IL-6 production was not resolved in recovery monocytes. Multivariate analysis of RNA-sequencing data identified 2938 significantly altered RNA transcripts in acute-stage monocytes (fold expression ≤−1.5 or ≥1.5; adjusted p ≤ 0.01), relative to controls. Comparing DNA methylation levels in circulating monocytes of CAP patients to controls revealed minimal differences, specifically in DNAse hypersensitive sites (HS) of acute-stage monocytes. Data integration identified a cholesterol biosynthesis gene signature and DNAse HS axis of IL-1β and IL-10 production (R2 =0.51).ConclusionsCirculating monocytes obtained from CAP patients during the acute stage exhibited impaired cytokine production capacities, indicative of reprogramming to a state of immune tolerance, which was not fully resolved after 1 month. Our split-sample study showed that 51% of the immune tolerance phenotype can be explained, at least in part, by coordinated shifts in cholesterol biosynthesis gene expression and DNAse HS methylation levels. A multi-scale model identified an epigenetic and transcriptomic signature of immune tolerance in monocytes, with implications for future interventions in immunosuppression.Trial registrationNCT number NCT02928367

Anatomy and evolution of a DNA replication origin.

DNA amplification occurs at the DNA puff II/9A locus in the fungus fly Sciara coprophila. As a foundation to study the molecular mechanism for the initiating events of II/9A DNA re-replication, we have sequenced 14kb spanning a DNase hypersensitive site (DHS) upstream of the 1kb amplification origin and through transcription units II/9-1 and II/9-2 downstream of the origin. These elements are annotated as well as the ORC binding site at the origin and the transition point (TP) between continuous and discontinuous DNA syntheses that marks the origin of bidirectional replication at the nucleotide level. A 9bp motif found at the TP is repeated near the other end of the 1 kb ORI and may identify a putative second TP. The steroid hormone ecdysone induces DNA amplification as well as transcription and puffing at locus II/9A. Within the 14 kb, several matches to the ecdysone response element (EcRE) consensus sequence were identified, including some in the amplification origin region. EcRE O-P is at a central axis of a remarkable symmetry, equidistant to the TPs that are themselves equidistant to EcRE O-1 and EcRE O-2. DNA sequence alterations have occurred throughout the II/9A region in a newly discovered polymorphism (#2). Polymorphism #2 is not specific to developmental stage, sex, or tissue, and it does not impair DNA amplification. The DHS, both 9bp TP sequences, and EcREs O-1, O-P, and O-2 are conserved between the polymorphism #1 and #2 sequences, suggesting their functional importance and retention during evolutionary selection. Moreover, a 72 bp sequence in the Sciara DHS at DNA puff II/9A is conserved in DNA puff C-3 of Rhynchosciara americana. Comparisons are discussed between the Sciara II/9A amplicon and the chorion locus amplicon on the third chromosome of Drosophila.