相关论文
deepTools2: a next generation web server for deep-sequencing data analysis
Abstract:We present an update to our Galaxy-based web server for processing and visualizing deeply sequenced data. Its core tool set, deepTools, allows users to perform complete bioinformatic workflows ranging from quality controls and normalizations of aligned reads to integrative analyses, including clustering and visualization approaches. Since we first described our deepTools Galaxy server in 2014, we have implemented new solutions for many requests from the community and our users. Here, we introduce significant enhancements and new tools to further improve data visualization and interpretation. deepTools continue to be open to all users and freely available as a web service at deeptools.ie-freiburg.mpg.de. The new deepTools2 suite can be easily deployed within any Galaxy framework via the toolshed repository, and we also provide source code for command line usage under Linux and Mac OS X. A public and documented API for access to deepTools functionality is also available.
相关论文分享
暂无分享，去 创建一个
参考文献
[1] Kristian Vlahovicek,et al. Genomation: a Toolkit to Summarize, Annotate and Visualize Genomic Intervals , 2015, Bioinform..
[2] Gonçalo R. Abecasis,et al. The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..
[3] Eric Nestler,et al. ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases , 2014, BMC Genomics.
[4] Raymond K. Auerbach,et al. An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.
[5] ENCODEConsortium,et al. An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.
[6] Michael Q. Zhang,et al. Integrative analysis of 111 reference human epigenomes , 2015, Nature.
[7] Alexander S. Garruss,et al. Enhancer-associated H3K4 monomethylation by Trithorax-related, the Drosophila homolog of mammalian Mll3/Mll4. , 2012, Genes & development.
[8] Ho-Ryun Chung,et al. High-Affinity Sites Form an Interaction Network to Facilitate Spreading of the MSL Complex across the X Chromosome in Drosophila. , 2015, Molecular cell.
[9] Allon M. Klein,et al. Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells , 2015, Cell.
[10] P. Bickel,et al. Systematic evaluation of factors influencing ChIP-seq fidelity , 2012, Nature Methods.
[11] Steven J. M. Jones,et al. An Interactive Analysis and Exploration Tool for Epigenomic Data , 2013, Comput. Graph. Forum.
[12] Fidel Ramírez,et al. deepTools: a flexible platform for exploring deep-sequencing data , 2014, Nucleic Acids Res..
[13] Anton Nekrutenko,et al. Dissemination of scientific software with Galaxy ToolShed , 2014, Genome Biology.
[14] Anton Nekrutenko,et al. Using Galaxy to Perform Large‐Scale Interactive Data Analyses , 2007, Current protocols in bioinformatics.
[15] Data production leads,et al. An integrated encyclopedia of DNA elements in the human genome , 2012 .
[16] Yong Zhang,et al. Identifying ChIP-seq enrichment using MACS , 2012, Nature Protocols.
[17] D. Karolchik,et al. The UCSC Genome Browser database: 2016 update , 2015, bioRxiv.
[18] Daniel J. Blankenberg,et al. Using Galaxy to Perform Large‐Scale Interactive Data Analyses , 2012, Current protocols in bioinformatics.
[19] Olivier Cuvier,et al. Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing. , 2014, Molecular cell.
[20] C. Glass,et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.
[21] Marc D. Perry,et al. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia , 2012, Genome research.
Identification of piRNA Binding Sites Reveals the Argonaute Regulatory Landscape of the C. elegans Germline
Cell
2018
摘要
piRNAs (Piwi-interacting small RNAs) engage Piwi Argonautes to silence transposons and promote fertility in animal germlines. Genetic and computational studies have suggested that C. elegans piRNAs tolerate mismatched pairing and in principle could target every transcript. Here we employ in vivo cross-linking to identify transcriptome-wide interactions between piRNAs and target RNAs. We show that piRNAs engage all germline mRNAs and that piRNA binding follows microRNA-like pairing rules. Targeting correlates better with binding energy than with piRNA abundance, suggesting that piRNA concentration does not limit targeting. In mRNAs silenced by piRNAs, secondary small RNAs accumulate at the center and ends of piRNA binding sites. In germline-expressed mRNAs, however, targeting by the CSR-1 Argonaute correlates with reduced piRNA binding density and suppression of piRNA-associated secondary small RNAs. Our findings reveal physiologically important and nuanced regulation of individual piRNA targets and provide evidence for a comprehensive post-transcriptional regulatory step in germline gene expression.
Identification of the m6Am Methyltransferase PCIF1 Reveals the Location and Functions of m6Am in the Transcriptome.
Molecular cell
2019
摘要
mRNAs are regulated by nucleotide modifications that influence their cellular fate. Two of the most abundant modified nucleotides are N6-methyladenosine (m6A), found within mRNAs, and N6,2'-O-dimethyladenosine (m6Am), which is found at the first transcribed nucleotide. Distinguishing these modifications in mapping studies has been difficult. Here, we identify and biochemically characterize PCIF1, the methyltransferase that generates m6Am. We find that PCIF1 binds and is dependent on the m7G cap. By depleting PCIF1, we generated transcriptome-wide maps that distinguish m6Am and m6A. We find that m6A and m6Am misannotations arise from mRNA isoforms with alternative transcription start sites (TSSs). These isoforms contain m6Am that maps to "internal" sites, increasing the likelihood of misannotation. We find that depleting PCIF1 does not substantially affect mRNA translation but is associated with reduced stability of a subset of m6Am-annotated mRNAs. The discovery of PCIF1 and our accurate mapping technique will facilitate future studies to characterize m6Am's function.
The innate sensor ZBP1-IRF3 axis regulates cell proliferation in multiple myeloma
bioRxiv
2020
摘要
ZBP1 is an inducible, non-constitutively expressed cellular nucleic acid sensor that triggers type I interferon (IFN) responses via phosphorylation and activation of the transcription factor (TF) IRF3 by TBK1. However, the role of the ZBP1-IRF3 axis in cancer is not known. Here we show that ZBP1 is selectively and constitutively expressed in late B cell development and it is required for optimal T cell-dependent humoral immune responses. In the plasma cell (PC) cancer multiple myeloma, interaction of constitutively expressed ZBP1 with TBK1 and IRF3 results in IRF3 phosphorylation. Notably, rather than IFN type I response genes, IRF3 directly activates, in part through co-operation with the PC lineage-defining TF IRF4, cell cycle genes thus promoting myeloma cell proliferation. This generates a novel, potentially therapeutically targetable and relatively selective myeloma cell addiction to the ZBP1-IRF3 axis. These data expand our knowledge of the role of cellular immune sensors in cancer biology.
Predicting gene regulatory regions with a convolutional neural network for processing double-strand genome sequence information
PloS one
2020
摘要
With advances in sequencing technology, a vast amount of genomic sequence information has become available. However, annotating biological functions particularly of non-protein-coding regions in genome sequences without experiments is still a challenging task. Recently deep learning–based methods were shown to have the ability to predict gene regulatory regions from genome sequences, promising to aid the interpretation of genomic sequence data. Here, we report an improvement of the prediction accuracy for gene regulatory regions by using the design of convolution layers that efficiently process genomic sequence information, and developed a software, DeepGMAP, to train and compare different deep learning–based models (https://github.com/koonimaru/DeepGMAP). First, we demonstrate that our convolution layers, termed forward- and reverse-sequence scan (FRSS) layers, integrate both forward and reverse strand information, and enhance the power to predict gene regulatory regions. Second, we assessed previous studies and identified problems associated with data structures that caused overfitting. Finally, we introduce visualization methods to examine what the program learned. Together, our FRSS layers improve the prediction accuracy for gene regulatory regions.
Comprehensive benchmarking of software for mapping whole genome bisulfite data: from read alignment to DNA methylation analysis
bioRxiv
2020
摘要
Whole genome bisulfite sequencing is currently at the forefront of epigenetic analysis, facilitating the nucleotide-level resolution of 5-methylcytosine (5mC) on a genome-wide scale. Specialised software have been developed to accommodate the unique difficulties in aligning such sequencing reads to a given reference, building on the knowledge acquired from model organisms such as human, or Arabidopsis thaliana. As the field of epigenetics expands its purview to non-model plant species, new challenges arise which bring into question the suitability of previously established tools. Herein, nine short-read aligners are evaluated: Bismark, BS-Seeker2, BSMAP, BWA-meth, ERNE-BS5, GEM3, GSNAP, Last, and segemehl. Precision-recall of simulated alignments, in comparison to real sequencing data obtained from three natural accessions, reveals on-balance that BWA-meth and BSMAP are able to make the best use of the data during mapping. The influence of difficult-to-map regions, characterised by deviations in sequencing depth over repeat annotations, is evaluated in terms of the mean absolute deviation of the resulting methylation calls in comparison to a realistic methylome. Downstream methylation analysis is responsive to the handling of multi-mapping reads relative to mapping quality (MAPQ), and potentially susceptible to bias arising from the increased sequence complexity of densely-methylated reads.
Transcriptional network orchestrating regional patterning of cortical progenitors
bioRxiv
2020
摘要
We uncovered a transcription factor (TF) network that regulates cortical regional patterning. Screening the expression of hundreds of TFs in the developing mouse cortex identified 38 TFs that are expressed in gradients in the ventricular zone (VZ). We tested whether their cortical expression was altered in mutant mice with known patterning defects (Emx2, Nr2f1 and Pax6), which enabled us to define a cortical regionalization TF network (CRTFN). To identify genomic programming underlying this network, we performed TF ChIP-seq and chromatin-looping conformation to identify enhancer-gene interactions. To map enhancers involved in regional patterning of cortical progenitors, we performed assays for epigenomic marks and DNA accessibility in VZ cells purified from wild-type and patterning mutant mice. This integrated approach has identified a CRTFN and VZ enhancers involved in cortical regional patterning.
Efficient low-cost chromatin profiling with CUT&Tag
Nature Protocols
2020
摘要
We recently introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling strategy in which antibodies are bound to chromatin proteins in situ in permeabilized nuclei. These antibodies are then used to tether the cut-and-paste transposase Tn5. Activation of the transposase simultaneously cleaves DNA and adds adapters ('tagmentation') for paired-end DNA sequencing. Here, we introduce a streamlined CUT&Tag protocol that suppresses DNA accessibility artefacts to ensure high-fidelity mapping of the antibody-targeted protein and improves the signal-to-noise ratio over current chromatin profiling methods. Streamlined CUT&Tag can be performed in a single PCR tube, from cells to amplified libraries, providing low-cost genome-wide chromatin maps. By simplifying library preparation CUT&Tag requires less than a day at the bench, from live cells to sequencing-ready barcoded libraries. As a result of low background levels, barcoded and pooled CUT&Tag libraries can be sequenced for as little as $25 per sample. This enables routine genome-wide profiling of chromatin proteins and modifications and requires no special skills or equipment.
Peax: Interactive Visual Pattern Search in Sequential Data Using Unsupervised Deep Representation Learning
bioRxiv
2019
摘要
We present Peax, a novel feature-based technique for interactive visual pattern search in sequential data, like time series or data mapped to a genome sequence. Visually searching for patterns by similarity is often challenging because of the large search space, the visual complexity of patterns, and the user’s perception of similarity. For example, in genomics, researchers try to link patterns in multivariate sequential data to cellular or pathogenic processes, but a lack of ground truth and high variance makes automatic pattern detection unreliable. We have developed a convolutional autoencoder for unsupervised representation learning of regions in sequential data that can capture more visual details of complex patterns compared to existing similarity measures. Using this learned representation as features of the sequential data, our accompanying visual query system enables interactive feedback-driven adjustments of the pattern search to adapt to the users’ perceived similarity. Using an active learning sampling strategy, Peax collects user-generated binary relevance feedback. This feedback is used to train a model for binary classification, to ultimately find other regions that exhibit patterns similar to the search target. We demonstrate Peax’s features through a case study in genomics and report on a user study with eight domain experts to assess the usability and usefulness of Peax. Moreover, we evaluate the effectiveness of the learned feature representation for visual similarity search in two additional user studies. We find that our models retrieve significantly more similar patterns than other commonly used techniques.
Mutation detection in thousands of acute myeloid leukemia cells using single cell RNA-sequencing
bioRxiv
2018
摘要
Virtually all tumors are genetically heterogeneous, containing subclonal populations of cells that are defined by distinct mutations1. Subclones can have unique phenotypes that influence disease progression2, but these phenotypes are difficult to characterize: subclones usually cannot be physically purified, and bulk gene expression measurements obscure interclonal differences. Single-cell RNA-sequencing has revealed transcriptional heterogeneity within a variety of tumor types, but it is unclear how this expression heterogeneity relates to subclonal genetic events – for example, whether particular expression clusters correspond to mutationally defined subclones3,4,5,6-9. To address this question, we developed an approach that integrates enhanced whole genome sequencing (eWGS) with the 10x Genomics Chromium Single Cell 5’ Gene Expression workflow (scRNA-seq) to directly link expressed mutations with transcriptional profiles at single cell resolution. Using bone marrow samples from five cases of primary human Acute Myeloid Leukemia (AML), we generated WGS and scRNA-seq data for each case. Duplicate single cell libraries representing a median of 20,474 cells per case were generated from the bone marrow of each patient. Although the libraries were 5’ biased, we detected expressed mutations in cDNAs at distances up to 10 kbp from the 5’ ends of well-expressed genes, allowing us to identify hundreds to thousands of cells with AML-specific somatic mutations in every case. This data made it possible to distinguish AML cells (including normal-karyotype AML cells) from surrounding normal cells, to study tumor differentiation and intratumoral expression heterogeneity, to identify expression signatures associated with subclonal mutations, and to find cell surface markers that could be used to purify subclones for further study. The data also revealed transcriptional heterogeneity that occurred independently of subclonal mutations, suggesting that additional factors drive epigenetic heterogeneity. This integrative approach for connecting genotype to phenotype in AML cells is broadly applicable for analysis of any sample that is phenotypically and genetically heterogeneous.
Transcriptional and epigenomic landscapes of CNS and non-CNS vascular endothelial cells
eLife
2018
摘要
Vascular endothelial cell (EC) function depends on appropriate organ-specific molecular and cellular specializations. To explore genomic mechanisms that control this specialization, we have analyzed and compared the transcriptome, accessible chromatin, and DNA methylome landscapes from mouse brain, liver, lung, and kidney ECs. Analysis of transcription factor (TF) gene expression and TF motifs at candidate cis-regulatory elements reveals both shared and organ-specific EC regulatory networks. In the embryo, only those ECs that are adjacent to or within the central nervous system (CNS) exhibit canonical Wnt signaling, which correlates precisely with blood-brain barrier (BBB) differentiation and Zic3 expression. In the early postnatal brain, single-cell RNA-seq of purified ECs reveals (1) close relationships between veins and mitotic cells and between arteries and tip cells, (2) a division of capillary ECs into vein-like and artery-like classes, and (3) new endothelial subtype markers, including new validated tip cell markers.
TIMEOR: a web-based tool to uncover temporal regulatory mechanisms from multi-omics data
bioRxiv
2020
摘要
Uncovering how transcription factors (TFs) regulate their targets at the DNA, RNA and protein levels over time is critical to define gene regulatory networks (GRNs) in normal and diseased states. RNA-seq has become a standard method to measure gene regulation using an established set of analysis steps. However, none of the currently available pipeline methods for interpreting ordered genomic data (in time or space) use time series models to assign cause and effect relationships within GRNs, are adaptive to diverse experimental designs, or enable user interpretation through a web-based platform. Furthermore, methods which integrate ordered RNA-seq data with transcription factor binding data are urgently needed. Here, we present TIMEOR (Trajectory Inference and Mechanism Exploration with Omics data in R), the first web-based and adaptive time series multi-omics pipeline method which infers the relationship between gene regulatory events across time. TIMEOR addresses the critical need for methods to predict causal regulatory mechanism networks between TFs from time series multi-omics data. We used TIMEOR to identify a new link between insulin stimulation and the circadian rhythm cycle. TIMEOR is available at https://github.com/ashleymaeconard/TIMEOR.git.
Distinct roles of cohesin-SA1 and cohesin-SA2 in 3D chromosome organization
Nature Structural & Molecular Biology
2017
摘要
Two variant cohesin complexes containing SMC1, SMC3, RAD21 and either SA1 (also known as STAG1) or SA2 (also known as STAG2) are present in all cell types. We report here their genomic distribution and specific contributions to genome organization in human cells. Although both variants are found at CCCTC-binding factor (CTCF) sites, a distinct population of the SA2-containing cohesin complexes (hereafter referred to as cohesin-SA2) localize to enhancers lacking CTCF, are linked to tissue-specific transcription and cannot be replaced by the SA1-containing cohesin complex (cohesin-SA1) when SA2 is absent, a condition that has been observed in several tumors. Downregulation of each of these variants has different consequences for gene expression and genome architecture. Our results suggest that cohesin-SA1 preferentially contributes to the stabilization of topologically associating domain boundaries together with CTCF, whereas cohesin-SA2 promotes cell-type-specific contacts between enhancers and promoters independently of CTCF. Loss of cohesin-SA2 rewires local chromatin contacts and alters gene expression. These findings provide insights into how cohesin mediates chromosome folding and establish a novel framework to address the consequences of mutations in cohesin genes in cancer.Depletion, ChIP and Hi-C analyses of the genomic distribution of the two variant cohesin complexes in human cells reveal non-redundant functions and differential contributions to 3D genome organization.
DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome
Nature
2017
摘要
Transposable elements are viewed as ‘selfish genetic elements’, yet they contribute to gene regulation and genome evolution in diverse ways. More than half of the human genome consists of transposable elements. Alu elements belong to the short interspersed nuclear element (SINE) family of repetitive elements, and with over 1 million insertions they make up more than 10% of the human genome. Despite their abundance and the potential evolutionary advantages they confer, Alu elements can be mutagenic to the host as they can act as splice acceptors, inhibit translation of mRNAs and cause genomic instability. Alu elements are the main targets of the RNA-editing enzyme ADAR and the formation of Alu exons is suppressed by the nuclear ribonucleoprotein HNRNPC, but the broad effect of massive secondary structures formed by inverted-repeat Alu elements on RNA processing in the nucleus remains unknown. Here we show that DHX9, an abundant nuclear RNA helicase, binds specifically to inverted-repeat Alu elements that are transcribed as parts of genes. Loss of DHX9 leads to an increase in the number of circular-RNA-producing genes and amount of circular RNAs, translational repression of reporters containing inverted-repeat Alu elements, and transcriptional rewiring (the creation of mostly nonsensical novel connections between exons) of susceptible loci. Biochemical purifications of DHX9 identify the interferon-inducible isoform of ADAR (p150), but not the constitutively expressed ADAR isoform (p110), as an RNA-independent interaction partner. Co-depletion of ADAR and DHX9 augments the double-stranded RNA accumulation defects, leading to increased circular RNA production, revealing a functional link between these two enzymes. Our work uncovers an evolutionarily conserved function of DHX9. We propose that it acts as a nuclear RNA resolvase that neutralizes the immediate threat posed by transposon insertions and allows these elements to evolve as tools for the post-transcriptional regulation of gene expression.
Transposable elements drive reorganisation of 3D chromatin during early embryogenesis
bioRxiv
2019
摘要
Transposable elements are abundant genetic components of eukaryotic genomes with important regulatory features affecting transcription, splicing, and recombination, among others. Here we demonstrate that the Murine Endogenous Retroviral Element (MuERV-L/MERVL) family of transposable elements drives the 3D reorganisation of the genome in the early mouse embryo. By generating Hi-C data in 2-cell-like cells, we show that MERLV elements promote the formation of insulating domain boundaries through-out the genome in vivo and in vitro. The formation of these boundaries is coupled to the upregulation of directional transcription from MERVL, which results in the activation of a subset of the gene expression programme of the 2-cell stage embryo. Domain boundaries in the 2-cell stage embryo are transient and can be remodelled without undergoing cell division. Remarkably, we find extensive inter-strain MERVL variation, suggesting multiple non-overlapping rounds of recent genome invasion and a high regulatory plasticity of genome organisation. Our results demonstrate that MERVL drive chromatin organisation during early embryonic development shedding light into how nuclear organisation emerges during zygotic genome activation in mammals.
The dynamic, combinatorial cis-regulatory lexicon of epidermal differentiation
Nature Genetics
2020
摘要
Transcription factors (TFs) bind DNA sequence motif vocabularies in cis-regulatory elements (CREs) to modulate chromatin state and gene expression during cell state transitions. A quantitative understanding of how motif lexicons influence dynamic regulatory activity has been elusive due to the combinatorial nature of the cis-regulatory code. To address this, we undertook multi-omic data profiling of chromatin and expression dynamics across epidermal differentiation to identify 40,103 dynamic CREs associated with 3,609 dynamically expressed genes, then applied an interpretable deep learning framework to model the cis-regulatory logic of chromatin accessibility. This identified cooperative DNA sequence rules in dynamic CREs regulating synchronous gene modules with diverse roles in skin differentiation. Massively parallel reporter analysis validated temporal dynamics and cooperative cis-regulatory logic. Variants linked to human polygenic skin disease were enriched in these time-dependent combinatorial motif rules. This integrative approach reveals the combinatorial cis-regulatory lexicon of epidermal differentiation and represents a general framework for deciphering the organizational principles of the cis-regulatory code in dynamic gene regulation. HIGHLIGHTS An integrative multi-omic resource profiling chromatin and expression dynamics across keratinocyte differentiation Predictive deep learning models of chromatin dynamics reveal a high-resolution cis-regulatory DNA motif lexicon of epidermal differentiation Model interpretation enables discovery of combinatorial cis-regulatory logic of homotypic and heterotypic motif combinations Massively parallel reporter experiments validate temporal dynamics and cis-regulatory logic of the combinatorial motif lexicon
CTCF is a Barrier for Totipotent-like Reprogramming
bioRxiv
2020
摘要
Totipotent cells have the ability of generating embryonic and extra-embryonic tissues1,2. Interestingly, a rare population of cells with totipotent-like potential was identified within ESC cultures3. These cells, known as 2 cell (2C)-like cells, arise from ESC and display similar features to those found in the totipotent 2 cell embryo2-4. However, the molecular determinants of 2C-like conversion have not been completely elucidated. Here, we show that CTCF is a barrier for 2C-like reprogramming. Indeed, forced conversion to a 2C-like state by DUX expression was associated with DNA damage at a subset of CTCF binding sites. Endogenous or DUX-induced 2C-like ESC showed decreased CTCF enrichment at known binding sites, suggesting that acquisition of a totipotent-like state is associated with a highly dynamic chromatin architecture. Accordingly, depletion of CTCF in ESC efficiently promoted spontaneous and asynchronous conversion to a totipotent-like state. This phenotypic reprogramming was reversible upon restoration of CTCF levels. Furthermore, we showed that transcriptional activation of the ZSCAN4 cluster was necessary for successful 2C-like reprogramming. In summary, we revealed the intimate relation between CTCF and totipotent-like reprogramming.
CTCF Promotes Long-range Enhancer-promoter Interactions and Lineage-specific Gene Expression in Mammalian Cells
bioRxiv
2020
摘要
Topologically associating domains (TAD) and insulated neighborhoods (INs) have been proposed to constrain enhancer-promoter communications to enable cell-type specific transcription programs, but recent studies show that disruption of TADs and INs resulted in relatively mild changes in gene expression profiles. To better understand the role of chromatin architecture in dynamic enhancer-promoter contacts and lineage-specific gene expression, we have utilized the auxin-inducible degron system to acutely deplete CTCF, a key factor involved in TADs and IN formation, in mouse embryonic stem cells (mESCs) and examined chromatin architecture and gene regulation during neural differentiation. We find that while CTCF depletion leads to global weakening of TAD boundaries and loss of INs, only a minor fraction of enhancer-promoter contacts are lost, affecting a small subset of genes. The CTCF-dependent enhancer-promoter contacts tend to be long-range, spanning hundreds of kilobases, and are established directly by CTCF binding to promoters. Disruption of CTCF binding at the promoter reduces enhancer-promoter contacts and transcription, while artificial tethering of CTCF to the promoter restores the enhancer-promoter contacts and gene activation. Genome-wide analysis of CTCF binding and gene expression across multiple mouse tissues suggests that CTCF-dependent promoter-enhancer contacts may regulate expression of additional mouse genes, particularly those expressed in the brain. Our results uncover both CTCF-dependent and independent enhancer-promoter contacts, and highlight a distinct role for CTCF in promoting enhancer-promoter contacts and gene activation in addition to its insulator function.
Analysis of chromatin accessibility in human epidermis identifies putative barrier dysfunction-sensing enhancers
PloS one
2017
摘要
To identify putative gene regulatory regions that respond to epidermal injury, we mapped chromatin dynamics in a stratified human epidermis during barrier maturation and disruption. Engineered skin substitutes (ESS) cultured at the air-liquid interface were used as a model of developing human epidermis with incomplete barrier formation. The epidermal barrier stabilized following engraftment onto immunocompromised mice, and was compromised again upon injury. Modified formaldehyde-assisted isolation of regulatory elements (FAIRE) was used to identify accessible genomic regions characteristic of monolayer keratinocytes, ESS in vitro, grafted ESS, and tape-stripped ESS graft. We mapped differentiation- and maturation-associated changes in transcription factor binding sites enriched at each stage and observed overrepresentation of AP-1 gene family motifs in barrier-deficient samples. Transcription of TSLP, an important effector of immunological memory in response to allergen exposure, was dramatically elevated in our barrier-deficient samples. We identified dynamic DNA elements that correlated with TSLP induction and may contain enhancers that regulate TSLP. Two dynamic regions were located near the TSLP promoter and overlapped with allergy-associated SNPs rs17551370 and rs2289877, strongly implicating these loci in the regulation of TSLP expression in allergic disease. Additional dynamic chromatin regions ~250kb upstream of the TSLP promoter were found to be in high linkage disequilibrium with allergic disease SNPs. Taken together, these results define dynamic chromatin accessibility changes during epidermal development and dysfunction.
Spatial chromatin architecture alteration by structural variations in human genomes at the population scale
Genome Biology
2018
摘要
This genome-wide study is focused on the impact of structural variants identified in individuals from 26 human populations onto three-dimensional structures of their genomes. We assess the tendency of structural variants to accumulate in spatially interacting genomic segments and design a high-resolution computational algorithm to model the 3D conformational changes resulted by structural variations. We show that differential gene transcription is closely linked to variation in chromatin interaction networks mediated by RNA polymerase II. We also demonstrate that CTCF-mediated interactions are well conserved across population, but enriched with disease-associated SNPs. Altogether, this study assesses the critical impact of structural variants on the higher order organization of chromatin folding and provides unique insight into the mechanisms regulating gene transcription at the population scale, among which the local arrangement of chromatin loops seems to be the leading one. It is the first insight into the variability of the human 3D genome at the population scale.
Swr1 mediated H2A.ZPht1 incorporation designates centromere DNA for de novo CENP-ACnp1 assembly
bioRxiv
2017
摘要
The underlying hallmark of centromeres is the presence of specialized nucleosomes in which histone H3 is replaced by CENP-A. The events that mediate the installation of CENP-A in place of H3 remain poorly characterized. H2A.Z is linked to transcriptional competence and associates with mammalian centromeres. We find that H2A.ZPht1 and the Swr1 complex are enriched in fission yeast CENP-ACnp1 chromatin. Our analysis shows that Swr1, Msc1 and H2A.ZPht1 are required to maintain CENP-ACnp1 chromatin integrity. Cell cycle analyses demonstrate that H2A.ZPht1 is deposited in S phase, coincident with the deposition of placeholder H3, and prior to CENP-ACnp1 replenishment in G2. Establishment assays reveal that H2A.ZPht1 and Swr1 are required for de novo assembly of CENP-ACnp1 onto naïve centromere DNA. We propose that features akin to promoters within centromere DNA program the incorporation of H2A.ZPht1 via Swr1, and mediate the replacement of resident H3 nucleosomes with CENP-A nucleosomes thereby defining centromeres.